hijing/examples/example_4.f

0001 C****************************************************************************
0002 C          Program # 3 from Comp. Phys. Commun. 83 (1994) 307
0003 C                  by M. Gyulassy and X-.N. Wang
0004 C             Modified by V.Uzhinsky, CERN, Oct. 2003
0005 C***************************************************************************
0006
0007       CHARACTER FRAME*8,PROJ*8,TARG*8
0008
0009       COMMON/HIPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)
0010       SAVE  /HIPARNT/
0011
0012       COMMON/RANSEED/NSEED
0013       SAVE  /RANSEED/
0014
0015       CHARACTER *1 KEY
0016
0017       NSEED=0
0018
0019 C....switch on triggered jet production:
0020       IHPR2(3)=1
0021
0022       FRAME='CMS'
0023
0024       write(6,*)'===================================================='
0025       write(6,*)' Simulation of events with triggered hard processes '
0026       write(6,*)'          in hh-, hA- and AA-collisions             '
0027       write(6,*)'    Calculation will be performed in CM System      '
0028       write(6,*)'===================================================='
0029
0030       write(6,*)
0031       write(6,*)'Enter the energy per NN-collision (GeV)'
0032       read(5,*)  EFRM
0033
0034       write(6,*)
0035       write(6,*)'Enter a type of the "projectile" particle'
0036       write(6,*)
0037       write(6,*)' P proton,            PBAR anti-proton,'
0038       write(6,*)' N neutron,           NBAR anti-neutron,'
0039       write(6,*)' PI+ - positive pion, PI- negative pion,'
0040       write(6,*)' K+ positive kaon,    K- negative kaon'
0041       write(6,*)
0042       write(6,*)' A - nucleus --------------------------'
0043
0044       read(5,1) PROJ
0045 1     format(A8)
0046
0047       if(PROJ.ne.'A') then
0048         IAP=1
0049         if(PROJ.eq.'P'   ) IZP= 1
0050         if(PROJ.eq.'PBAR') IZP=-1
0051         if(PROJ.eq.'N'   ) IZP= 0
0052         if(PROJ.eq.'NBAR') IZP= 0
0053         if(PROJ.eq.'PI+' ) IZP= 1
0054         if(PROJ.eq.'PI-' ) IZP=-1
0055         if(PROJ.eq.'K+'  ) IZP= 1
0056         if(PROJ.eq.'K-'  ) IZP=-1
0057       else
0058         write(6,*)
0059         write(6,*)'Enter mass number and charge of the proj. nucleus'
0060         read(5,*)  IAP, IZP
0061       endif
0062
0063       write(6,*)
0064       write(6,*)'Enter a type of the "target" particle (same notations)'
0065       read(5,1) TARG
0066
0067       if(TARG.ne.'A') then
0068         IAT=1
0069         if(TARG.eq.'P'   ) IZT= 1
0070         if(TARG.eq.'PBAR') IZT=-1
0071         if(TARG.eq.'N'   ) IZT= 0
0072         if(TARG.eq.'NBAR') IZT= 0
0073         if(TARG.eq.'PI+' ) IZT= 1
0074         if(TARG.eq.'PI-' ) IZT=-1
0075         if(TARG.eq.'K+'  ) IZT= 1
0076         if(TARG.eq.'K-'  ) IZT=-1
0077       else
0078         write(6,*)
0079         write(6,*)'Enter mass number and charge of the target nucleus'
0080         read(5,*)  IAT, IZT
0081       endif
0082
0083       if(PROJ.eq.'A'.or.TARG.eq.'A') then
0084         write(6,*)'Enter Min. and Max. values of impact parameter (fm)'
0085         read(5,*)  BMIN, BMAX
0086       else
0087         BMIN=0.0
0088         BMAX=0.0
0089       endif
0090
0091       if(PROJ.eq.'A'.or.TARG.eq.'A') then
0092         write(6,*)
0093         write(6,*)' Would you like to take jet quenching into account?'
0094         write(6,*)'                 Y - Yes, N - No?   '
0095
0096         read(5,2) KEY
0097 2       format(A1)
0098
0099         if(KEY.eq.'Y'.or.KEY.eq.'y') then
0100           IHPR2(4)=1
0101         else
0102           IHPR2(4)=0
0103         endif
0104       endif
0105
0106 C....set the pt range of the triggered jets:
0107       write(6,*)
0108       write(6,*)' Set the Pt of triggered jets (in GeV/c)'
0109       read(5,*)   Pt_trigger
0110
0111       write(6,*)
0112       write(6,*)'Enter number of events'
0113       read(5,*)  N_events
0114
0115 C....initialize HIJING for requested interactions
0116       CALL HIJSET(EFRM,FRAME,PROJ,TARG,IAP,IZP,IAT,IZT)
0117
0118       write(6,*)' Simulation of interactions with'
0119       write(6,*)
0120       write(6,*)' Proj = ',PROJ,' and  Targ = ',TARG
0121       write(6,*)' IAP  =',IAP  ,'            IAT  =',IAT
0122       write(6,*)' IZP  =',IZP  ,'             IZT  =',IZT
0123       write(6,*)
0124       write(6,*)' Reference frame -   ',FRAME
0125       write(6,*)' ENERGY            ',EFRM,' GeV'
0126       write(6,*)' Number of generated events',N_events
0127       write(6,*)' Triggered jets Pt =>',Pt_trigger,' (GeV/c)'
0128       write(6,*)
0129
0130       DO 100 I_event=1,N_events
0131
0132         write(6,*)'Event # ',I_event,' -------------------------'
0133         write(6,*)
0134
0135         HIPR1(10)=-Pt_trigger
0136
0137         CALL HIJING(FRAME,BMIN,BMAX)
0138
0139 C....print out flavor code of the first jet:
0140         write(6,*)'Flavor code of the first jet: ',IHNT2(9)
0141
0142 C....and its four momentum:
0143         write(6,*)'     Px           Py           Pz',
0144      ,'            E     (GeV/c, GeV)'
0145         write(6,3)HINT1(21),HINT1(22),HINT1(23),HINT1(24)
0146 3       format(4(e11.4,2x))
0147         write(6,*)
0148
0149 C....print out flavor code of the second jet:
0150         write(6,*)'Flavor code of the second jet:',IHNT2(10)
0151
0152 C....and its four momentum:
0153         write(6,*)'     Px           Py           Pz',
0154      ,'            E     (GeV/c, GeV)'
0155         write(6,3)HINT1(31),HINT1(32),HINT1(33),HINT1(34)
0156         write(6,*)
0157
0158         PAUSE
0159 100   CONTINUE
0160
0161       STOP
0162       END
0163
0164       FUNCTION RAN(NSEED)
0165       RAN=RLU(NSEED)
0166       RETURN
0167       END