Hadr02/hijing1_383/Main_HIJING.f

0001       program test
0002       CHARACTER FRAME*8,PROJ*8,TARG*8
0003       COMMON/HIPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)
0004       SAVE  /HIPARNT/
0005       COMMON/HIMAIN1/NATT,EATT,JATT,NT,NP,N0,N01,N10,N11
0006       SAVE  /HIMAIN1/
0007       COMMON/HIMAIN2/KATT(130000,4),PATT(130000,4)
0008       SAVE  /HIMAIN2/
0009       COMMON/RANSEED/NSEED
0010       SAVE  /RANSEED/
0011 C
0012 C
0013 C....initialize HIJING for Au+Au collisions at c.m. energy of 200 GeV:
0014 c        FRAME="CMS"
0015 c        PROJ="A"
0016 c       TARG="A"
0017         IAP=1
0018         IZP=1
0019         IAT=197
0020         IZT=79
0021         EFRM=200.0   !GeV
0022         N_EVENT=100
0023         IHPR2(10)=0
0024 c       ---------------------------------------------------------
0025 c        initialize HIJING for geant4
0026         IHNT2(1)=IAP
0027         IHNT2(2)=IZP
0028         IHNT2(3)=IAT
0029         IHNT2(4)=IZT
0030         IHNT2(5)=211   !pi+
0031         IHNT2(6)=0
0032 cc     --------------Ap>1----------------------
0033
0034         HINT1(8)=MAX(ULMASS(2112),ULMASS(2212))
0035         HINT1(9)=HINT1(8)
0036         write(*,*)'mass at the beginning',HINT1(8), HINT1(9)
0037 cc---------------------------------------------
0038         if(IHNT2(5).ne.0) then
0039         HINT1(8)=ULMASS(IHNT2(5))
0040         endif
0041
0042 c        write(*,*)'modified mass',HINT1(8)
0043 c
0044 c
0045 c         ranseed=1097569630
0046 c         print *,"seed",ranseed
0047 c         call sseed(ranseed)
0048
0049 c         do i=1, 10
0050 c         write(*,*)'ran',NSEED, RAN(NSEED), RLU(NSEED)
0051 c         end do
0052 C
0053 C
0054 C*** initialize HIJING
0055 c        CALL HIJSET(EFRM,FRAME,PROJ,TARG,IAP,IZP,IAT,IZT)
0056          CALL HIJSET(EFRM,IAP,IZP,IAT,IZT)
0057 c         CALL HIJSET(EFRM)
0058
0059           write(*,*)'Sjet=', HINT1(11),'mb','Stot=',HINT1(13),'mb'
0060 C
0061 C....set BMIN=0 and BMAX=0.0 for central interaction
0062       BMIN=0.0
0063       BMAX=0.0    !HIPR1(34)+HIPR1(35)
0064 C....generating N_EVENT events of central AA interaction:
0065       DO 200 IE=1,N_EVENT
0066 c         CALL HIJING(FRAME,BMIN,BMAX)
0067          CALL HIJING(BMIN,BMAX)
0068          WRITE(*,*) IE,NATT,EATT
0069  200  continue
0070       STOP
0071       END
0072
0073 c       FUNCTION RAN(NSEED)
0074 c       RAN=RLU(NSEED)
0075 c       RETURN
0076 c       END