Main_HIJING.f

Go to the documentation of this file.

      program test
      CHARACTER frame*8,proj*8,targ*8
      common/hiparnt/hipr1(100),ihpr2(50),hint1(100),ihnt2(50)
      SAVE  /hiparnt/
      common/himain1/natt,eatt,jatt,nt,np,n0,n01,n10,n11
      SAVE  /himain1/
      common/himain2/katt(130000,4),patt(130000,4)
      SAVE  /himain2/
      common/ranseed/nseed
      SAVE  /ranseed/
C
C       
C....initialize HIJING for Au+Au collisions at c.m. energy of 200 GeV:
c        FRAME="CMS"
c        PROJ="A"
c   TARG="A"
        iap=1     
    izp=1       
    iat=197      
    izt=79    
        efrm=200.0   !GeV 
        n_event=100
        ihpr2(10)=0
c       ---------------------------------------------------------
c        initialize HIJING for geant4
    ihnt2(1)=iap
    ihnt2(2)=izp
    ihnt2(3)=iat
    ihnt2(4)=izt
    ihnt2(5)=211   !pi+
    ihnt2(6)=0
cc     --------------Ap>1----------------------
         
        hint1(8)=max(ulmass(2112),ulmass(2212))
        hint1(9)=hint1(8)
        write(*,*)'mass at the beginning',hint1(8), hint1(9)
cc---------------------------------------------
        if(ihnt2(5).ne.0) then
        hint1(8)=ulmass(ihnt2(5))
        endif
    
c        write(*,*)'modified mass',HINT1(8)
c
c
c         ranseed=1097569630
c         print *,"seed",ranseed
c         call sseed(ranseed)

c         do i=1, 10
c         write(*,*)'ran',NSEED, RAN(NSEED), RLU(NSEED)
c         end do
C
C
C*** initialize HIJING
c        CALL HIJSET(EFRM,FRAME,PROJ,TARG,IAP,IZP,IAT,IZT)
         CALL hijset(efrm,iap,izp,iat,izt)
c         CALL HIJSET(EFRM)
        
          write(*,*)'Sjet=', hint1(11),'mb','Stot=',hint1(13),'mb'
C
C....set BMIN=0 and BMAX=0.0 for central interaction
      bmin=0.0
      bmax=0.0    !HIPR1(34)+HIPR1(35)
C....generating N_EVENT events of central AA interaction:
      DO 200 ie=1,n_event
c         CALL HIJING(FRAME,BMIN,BMAX)
         CALL hijing(bmin,bmax)
         WRITE(*,*) ie,natt,eatt
 200  continue
      stop
      END 

c   FUNCTION RAN(NSEED)
c   RAN=RLU(NSEED)
c   RETURN
c   END