Warning, /coresoftware/generators/hijing/doc/hijing1.3.txt is written in an unsupported language. File is not indexed.
0001 ********************************************************************
0002 * HIJING Monte Carlo Program *
0003 * For Relativisitic Heavy Ion Collision *
0004 * Version 1.30 *
0005 * Miklos Gyulassy and Xin-Nian Wang *
0006 * Nuclear Science Division, MS 70A-3307 *
0007 * Lawrence Berkeley Laboratory and University of California *
0008 * Berkeley, CA 94720 USA *
0009 * email: XNWANG@nsdssd.lbl.gov *
0010 ********************************************************************
0011
0012 INSTRUCTION TO USERS
0013 *****************************************
0014 1)General information
0015 HIJING basically is designed to simulate multiple jets
0016 and particle production in pp, pA or AA collisions. PYTHIA
0017 subroutines and Lund jet fragmentation scheme are used. Therefore,
0018 this program uses the same particle identification code as JETSET.72.
0019 Users can also use subroutines in JETSET.72 and change the values
0020 of parameters in JETSET.72 and PYTHIA therein. But there are some
0021 special parameters in HIJING which users have to specify or change.
0022 This is the first version. Things may change in the future.
0023 If you find any bugs or problems please report them to the authors.
0024 This file contains a example program and all the HIJING subroutines.
0025 To save compiling time, JETSET.7.2 is not included in HIJING. To
0026 meet some specific needs of HIJING, JETSET.72 has been modified and
0027 renamed as HIPYSET. One should link HIJING with HIPYSET and the main
0028 program:
0029 ************************************************************************
0030 * $link yourmain,hijing,hipyset *
0031 ************************************************************************
0032
0033 2)Main subroutines
0034 There are two subroutines in HIJING which users have to call:
0035
0036 SUBROUTINE HIJSET(EFRM,FRAME,PROJ,TARG,IAP,IZP,IAT,IZT)
0037 Purpose: This subroutine is to initialize HIJING which should
0038 be called first before any subroutine.
0039 EFRM(GeV): is the colliding energy in the frame specified by
0040 FRAME.
0041 FRAME: (must be 'LAB' or 'CMS') to specify the frame of the
0042 collision
0043 PROJ,TARG:Character variables of proj and targ particles.
0044 ='P' proton
0045 ='PBAR' anti-proton
0046 ='N' neutron
0047 ='NBAR' anti-neutron
0048 ='PI+' positive pion
0049 ='PI-' negative pion
0050 ='K+' positive kaon
0051 ='K-' negative kaon
0052 ='A' nucleus
0053 IAP,IAT: mass number of the projectile and target nuclei, set
0054 to 1 for hadrons
0055 IZP,IZT: charge number of the projectile and target nuclei,for
0056 hadrons it is the charge of that hadron(=1,0,-1)
0057
0058 SUBROUTINE HIJING(FRAME,BMIN,BMAX)
0059 Purpose: This is the main program of HIJING which can be
0060 called(many times) only after HIJSET is called once.
0061 FRAME: (must be 'LAB' or 'CMS') to specify the frame of the
0062 collision
0063 BMIN,BMAX: are the low and up limits bewteen which the impact
0064 parameters are randomly but evenly selected for pA and
0065 AA collisions. For pp collisions, the events are
0066 averaged over all impact parameters.
0067
0068
0069
0070 5)Main Common blocks
0071 There are mainly three common blocks which provide users with
0072 important information of the generated events.
0073
0074 COMMON /HIMAIN1/ NATT,EATT,JATT,NT,NP,N0,N01,N10,N11
0075 COMMON /HIMAIN2/KATT(130000,4),PATT(130000,4)
0076 These give information of interaction and produced particles.
0077 NATT: total number of produced stable particles of the current
0078 event.
0079 EATT: the total c.m.s energy of the produced particles to check
0080 energy conservation.
0081 JATT: the total number of jet-pairs in the current event.
0082 NT,NP: number of wounded target and projetile nucleons.
0083 N0,N01,N10,N11: number of N-N,N-N_wounded,N_wounded-N, and
0084 N_wounded-N_wounded collisions in the current event.
0085 KATT(I,1)(I=1,NATT): partilce ID of the I_th produced particle.
0086 Users have to refer to JETSET7.2 for identifying
0087 particles with their ID's.(See the end of this
0088 documentation)
0089 KATT(I,2)(I=1,NATT):This is a code to identify the sources
0090 from which the particle comes.
0091 =0 projectile which has not interacted at all
0092 =1 projectile which only suffers an elastic collision
0093 =2 from the diffractive projectile in a single
0094 diffractive interaction
0095 =3 from the fragmentation of projectile strings
0096 =10 target which has not interacted at all
0097 =11 target which only suffers an elastic collision
0098 =12 from the diffractive target in a single diffractive
0099 interaction
0100 =13 from the fragmentation of target strings
0101 =20 from hard scattered partons which form string
0102 systems themselves.
0103 =30 from soft sea-quarks strings
0104 =40 from direct production in the hard processes(
0105 currently,only direct photons are included).
0106 KATT(I,3)(I=1,NATT): The line number of the parent particle.
0107 For finally produced or directly produced (not from
0108 the decay of another particle) particles, it is set
0109 to 0 (The option to retain the information of all
0110 particles including the decayed ones is IHPR2(21)=1).
0111 KATT(I,4)(I=1,NATT): Status number of the particle:
0112 =1 finally or directly produced particles.
0113 =11 particles which has already decayed.
0114 PATT(I,1:4)(I=1,NATT): four-momentum of I_th produced particle.
0115
0116
0117 COMMON/HIJJET1/NPJ(300),KFPJ(300,500),PJPX(300,500),PJPY(300,500)
0118 & ,PJPZ(300,500),PJPE(300,500),PJPM(300,500)
0119 & ,NTJ(300),KFTJ(300,500),PJTX(300,500),PJTY(300,500)
0120 & ,PJTZ(300,500),PJTE(300,500),PJTM(300,500)
0121 This common block contains information about produced
0122 partons which are connected with the leading quarks
0123 and diquarks of projectile or target to form a string
0124 system for fragmentation.
0125 NPJ(I)(I=1,IAP): number of partons associated with I_th
0126 projectile
0127 KFPJ(I,J)(J=1,NPJ(I)): parton ID of the J_th parton associated
0128 with I_th projectile.
0129 PJPX(I,J),PJPY(I,J),PJPZ(I,J),PJPE(I,J),PJPM(I,J): the four
0130 momentum and rest mass of J_th parton associated with
0131 the I_th projectile
0132 NTJ(I)(I=1,IAT): number of partons associated with I_th
0133 target
0134 KFTJ(I,J)(J=1,NTJ(I)): parton ID of the J_th parton associated
0135 with I_th target.
0136 PJTX(I,J),PJTY(I,J),PJTZ(I,J),PJTE(I,J),PJTM(I,J): the four
0137 momentum and rest mass of J_th parton associated with
0138 the I_th target.
0139
0140 COMMON/HIJJET2/NSG,NJSG(900),IASG(900,3),K1SG(900,100),K2SG(900,100)
0141 & ,PXSG(900,100),PYSG(900,100),PZSG(900,100),PESG(900,100)
0142 & ,PMSG(900,100)
0143
0144 This common block contains information about the partons which
0145 form string systems themselves without be connected to leading
0146 quarks and diquarks.
0147
0148 NSG: the total number of such string systems.
0149 NJSG(I)(I=1,NSG): number of partons in the I_th string system.
0150 IASG(I,1),IASG(I,2): specify which projetile and target produce
0151 the I_th string system
0152 IASG(I,3):=0 the jets will be quenched
0153 1 the jets will not be quenched
0154 K1SG(I,J)(J=1,NJSG(I)): color flow information of J_th parton
0155 in the I_th string system.
0156 K2SG(I,J)(J=1,NJSG(I)): ID code of J_th parton in I_th string.
0157 PXSG(I,J),PYSG(I,J),PZSG(I,J),PESG(I,J),PMSG(I,J): four
0158 momentum and rest mass of the J_th parton in the I_th
0159 string system.
0160
0161
0162 COMMON/HISTRNG/NFP(300,15),PP(300,15),NFT(300,15),PT(300,15)
0163
0164 This common block contains information about the projectile
0165 and target and the corresponding constituent quarks, diquarks.
0166
0167 NFP(K,1):flavor of quark in projectile.
0168 NFP(K,2):flavor of diquark in projectile.
0169 NFP(K,3):present particle ID of the projectile.
0170 NFP(K,4):original particle ID of projectile.
0171 NFP(K,5):collision status of projectile
0172 =0:no collision;
0173 =1:elastic;
0174 =2:the diffractive one in a single-diffractive
0175 =3:beccomes excited string after an inelastic
0176 NFP(K,6):is the total number of jet production associated
0177 with ths projectile. if NFP(K,6)<0 it can not
0178 produce jet anymore due to energy conservation.
0179 NFP(K,10):to indicate whether the valence quarks or
0180 diquarks has had a hard scattering,
0181 =0: have no hard scattering;
0182 =1,-1: has or had a hard scattering.
0183 NFP(K,11):total number of interactions this projectile
0184 has suffered so far.
0185 PP(K,1),PP(K,2),PP(K,3),PP(K,4),PP(K,5): four momentum and
0186 the invariant mass (Px,Py,Pz,E,M) of the projectile.
0187 PP(K,6),PP(K,7):transverse momentum (Px,Py) of the quark.
0188 PP(K,8),PP(K,9):transverse momentum (Px,Py) of the diquark.
0189 PP(K,10),PP(K,11),PP(K,12):three momentum (Px,Py,Pz) the
0190 quark or diquark has from the last hard scattering
0191 it suffered.
0192 PP(K,14):the mass of the projectile quark.
0193 PP(K,15):the mass of the projectile diquark.
0194 ......
0195 NFT(K,1-15),PT(K,1-15), are about the same information of
0196 target and target quarks, diquarks.
0197
0198 4)Event Options and parameters
0199
0200 COMMON/HIPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)
0201
0202 This common block contains most of HIJING input parameters,
0203 switchs(HIPR1,IHPR2) and event information(HINT1,IHNT2). The
0204 defualt values are given by D.
0205
0206 HIPR1(1):(D=1.5 GeV)Minimum value for the invariant mass of
0207 the excited string system in a hadron-hadron
0208 interaction.
0209 HIPR1(2):(D=0.35 GeV)width of the Gaussian pt distribution of
0210 hadron in string fragmenation.
0211 HIPR1(3),HIPR1(4):(D=0.5,0.9 GeV**-2) give the a and b
0212 parameters of the symmetric LUND fragmentatin fuction.
0213 HIPR1(5):(D=2.0) cut-off for the QCD radiation, string system
0214 with invariant mass below it will not radiate soft
0215 gluons.
0216 HIPR1(6):(D=0.1) The shadowing of structure fucntion at x=0
0217 is parametrized as HIPR1(6)*(A**(1/3)-1.0)
0218 HIPR1(7):not used
0219 HIPR1(8):(D=2.0 GeV) Minimum pt of hard or semihard scatterings.
0220 HIPR1(9):(D=-1.0 GeV) Maximum pt of hard or semihard
0221 scatterings, if negative there is no up-limit.
0222 HIPR1(10):(D=-2.25 GeV) To specify the value of PT for one hard
0223 scattering generated per event(we call them sample
0224 jets), which could be used to study jet in AA
0225 collision. If HIPR1(10) is negtive, its absolute value
0226 gives the low limit of the PT of the sample jets.
0227 HIPR1(11):(D=2.0 GeV) Minimum pt of a jet which will interact
0228 with excited nuclear matter, or when the pt of a jet
0229 is smaller than HIPR1(11) it will stop interacting
0230 further more.
0231 HIPR1(12):(D=1.0 fm)The transverse distance between a
0232 traversing jet and an excited nucleon(or string system)
0233 below which they will interact and the jet will lose
0234 energy and momentum to the nucleon.
0235 HIPR1(13):(D=1.0 fm)The mean free path of a jet to interact
0236 when it goes through an excited nuclear matter.
0237 HIPR1(14):(D=2.0 GeV/fm) The string tension of a gluon jet.
0238 It is the energy loss dE/dx=HIPR1(14) of the gluon
0239 jet when the string is streched. For quark jet
0240 dE/dx=0.5*HIPR1(14).
0241 HIPR1(15):(D=0.2 GeV) the momentum scale lamda in the
0242 calculation of alpha_strong.
0243 HIPR1(16):(D=2.0 GeV) The initial mometum scale Q0 for the
0244 evolution of the structure fuction.
0245 HIPR1(17):(D=2.0) K factor for the differential jet cross
0246 sections.
0247 HIPR1(18):(D=0.2) Fractional probability for color flip
0248 occuring to sea quarks among inelatic events in
0249 the DPM-like option.
0250 HIPR1(19),HIPR1(20):(D=0.1,1.6 GeV) The pt kick from the
0251 initial and final state interaction will follow
0252 1/(HIPR1(19)**2+pt**2)/(HIPR1(20)**2+pt**2)
0253 distribution. HIPR1(19) regularizes the singularity
0254 at pt=0 and HIPR1(20) gives the scale beyond which
0255 the distribution will be like 1/pt**4.
0256 HIPR1(21):(D=1.6 GeV) The maximum pt for the pt kick per
0257 initial or final state interaction. Beyond this value
0258 the pt_kick-HIPR1(21) will follow the same Gaussion
0259 distribution given by HIPR1(2).
0260 HIPR1(22):(D=2.0 GeV) The form factor parameter to supress
0261 the diquarks from having a pt transfer relative to the
0262 quarks,1/(1+pt**2/HIPR1(22)**2).
0263 HIPR1(23)-HIPR1(28):not used.
0264 HIPR1(29):(D=0.4 fm) the minimum distance between two nucleons
0265 inside a nucleus when the coordinates of all the
0266 nucleons in a nucleus are initialized.
0267 HIPR1(30):(D=2*HIPR1(31)=57.0 bm) the inclusive cross
0268 section sigma_s for soft interaction. The default
0269 value sigma_s=sigma_0 is used to insure the geometrical
0270 scaling of pp interaction cross section at low
0271 energies.
0272 HIPR1(31):(D=28.5 bm) the cross section sigam_0 which
0273 characterizes the geometrical size of a nucleon
0274 (HIPR1(31)=pi*b0**2). The default value is only
0275 for high-energy limite(sqrt{s}>200 GeV). At lower
0276 energies, a slight decrease which depends on energy
0277 is parametrized in the program. The default value
0278 of the above two parameters HIPR1(30),HIPR1(31) are
0279 only for NN type interactions. For other kinds
0280 of projectile or target hadrons users should
0281 change these values so that a correct inelastic
0282 and total cross section HINT1(12),HINT1(13) are
0283 obtained by the program.
0284 HIPR1(32):(D=3.90) The parameter in the scaled eikonal function
0285 of nucleon Xi(R)=A**2*(A*R)**3*K3(A*R),A=HIPR1(32)
0286 which is a Fourier transform of a dipole formtor
0287 1/(1+Q**2/mu**2)**2 and A=mu*b0 and pi*b0**2=HIPR1(31).
0288 HIPR1(33):fractional cross section of single-diffractive
0289 interaction.(parametrized)
0290 HIPR1(34):Maximum radial coordinator for projectile
0291 nucleons to be given by the initialzation program
0292 HIJSET.
0293 HIPR1(35):Maximum radial coordinator for target
0294 nucleons to be given by the initialzation program
0295 HIJSET.
0296 HIPR1(36)-HIPR1(39):not used
0297 HIPR1(40):(D=3.141592654) value of pi.
0298 HIPR1(41):not used
0299 HIPR1(42):(D=0.4)parameter in the expotnetial pt distribution
0300 of sea quarks production in DPM option
0301 HIPR1(43):(D=0.01) The fractional energy error relative to the
0302 colliding energy permitted per nucleon-nucleon
0303 collisions.
0304 HIPR1(44),HIPR1(45),HIPR1(46):(D=1.5,0.1,0.25) These
0305 parameters give the x distribution of valence quarks
0306 (1-x)**HIPR1(44)/(x**2+HIPR1(45)**2/s)**HIPR1(46)
0307 for baryons, and 1/(x**2+HIPR1(45)**2/s)**HIPR1(46)/
0308 ((1-x)**2+HIPR1(45)**2/s)**HIPR(46) for mesons, in
0309 soft interactions.
0310 HIPR1(47),HIPR1(48):(D=0.0,0.5) These two parameters give
0311 x distribution of the sea quarks
0312 (1-x)**HIPR1(47)/(x**2+HIPR1(45)**2/s)**HIPR1(48),
0313 in a soft interaction. The same is used for the
0314 the distribution of the invariant mass of the
0315 disassociated nucleon in a single diffractive
0316 collision.
0317 HIPR1(49)-HIPR1(100):not used.
0318
0319
0320 IHPR2(1):(D=1)switch for dipole-approximated QCD radiation
0321 of the string system in soft interaction.
0322 IHPR2(2):(D=3)switch for initial and final state radiation in
0323 the hard scattering.
0324 =0: both initial and final radiation off;
0325 =1: initial radiation on and final radiation off;
0326 =2: initial radiation off and final radiation on;
0327 =3: both initial and final radiation on.
0328 IHPR2(3):(D=0)switch for one hard scattering with specified PT
0329 PT=HIPR1(10) per event.
0330 =1:ordinary hard processes;
0331 =2:only direct photon production.
0332 IHPR2(4):(D=1)switch for jet quenching in the the exicited
0333 nuclear matter,i.e. jets only interact with those
0334 nuclear matter that has suffered collisions.
0335 IHPR2(5):(D=1)switch for the pt kick due to soft initial state
0336 interaction both for valence quarks or diquarks and
0337 gluons.
0338 IHPR2(6):(D=1)switch for the nuclear effect on the parton
0339 distribution function.
0340 IHPR2(7):(D=1)(1 or 2) selection of Duke-Owns parametrization
0341 set of nucleon structure function.
0342 IHPR2(8):(D=10)The maximum number of jet production per
0343 nucleon-nucleon interaction. When IHPR2(8)=0, jet
0344 production will be turned off. When IHPR2(8)<0, the
0345 number of jets production will be fixed for each
0346 NN collision at a value of abs(IHPR2(8)).
0347 IHPR2(9):(D=0)The switch to guarrantee at leat one minijet
0348 production per event(pp,pA or AA).
0349 IHPR2(10):(D=0) Option to print error warning messages that
0350 might happen. When a fatal error happens the current
0351 event will be abandoned and a new one is generated.
0352 IHPR2(11):(D=1) Choice of baryon production model.
0353 =0 no baryon-antibaryon pair production,
0354 initial diquark treated as a unit.
0355 =1 diquark-antidiquark pair production allowed, initial
0356 diquark treated as a unit.
0357 =2 diquark-antidiquark pair production allowed, with the
0358 possibility for diquark to be split according to the
0359 "popcorn" scheme.(see the menu of jetset72)
0360 IHPR2(12):(D=1)Option to turn off the decay of partiles
0361 such as pi0,K0S,D+-,Lampda,Sigma+-,X-0,Omega-.
0362 IHPR2(13):(D=1)Option to turn on single diffractive reactions.
0363 IHPR2(14):(D=1)Option to turn on elastic scattering
0364 IHPR2(15 - 17):not used
0365 IHPR2(18):(D=0)option to switch on B-quark production. Charm
0366 production is the default. When B-quark production is
0367 on, charm quark production is automatically off.
0368 IHPR2(19):(D=1)Option to turn on initial state interaction
0369 IHPR2(20):(D=1)Option to turn off(IHPR2(20)=0) the
0370 fragmentation.
0371 IHPR2(21):(D=0)Option to retain the information of all
0372 particles including those which have decayed
0373 and the decay chain. The line number of the
0374 parent particle is KATT(I,3). The status
0375 of a partcile, whether it is a finally produced
0376 particle (KATT(I,4)=1) or a decayed particle
0377 (KATT(I,4)=11).
0378 IHPR2(22)-IHPR2(50): not used.
0379
0380
0381
0382 HINT1(1):(GeV) The center-of-mass energy per nucleon-nucleon
0383 collision.
0384 HINT1(2): The Lorenze transform variable beta from LAB to CM
0385 frame of nucelon nucleon collisions.
0386 HINT1(3):The rapidity of the CM frame beta=tang(ycm).
0387 HINT1(4):The rapidity of projectile y_proj.
0388 HINT1(5):The rapidity of target y_targ.
0389 HINT1(6):(GeV)The energy of projectile in the chosen frame.
0390 HINT1(7):(GeV)The energy of target in the chosen frame.
0391 HINT1(8):(GeV)The rest mass of projectile particles
0392 HINT1(9):(GeV)The rest mass of target particles.
0393 HINT1(10):(mb) The averaged cross section for jet production
0394 per nucleon-nucleon collisions of the nuclear
0395 interaction(When very small compared to the total
0396 nucleon-nucleon cross section this is very close
0397 to the jet inclusive cross section HINT1(11)).
0398 HINT1(11):(mb) The averaged inclusive cross section for jet
0399 production per nucleon-nucleon collisions of the
0400 nuclear interaction.
0401 HINT1(12):(mb) The averaged inelastic cross section of
0402 nucleon-nucleon collisions.
0403 HINT1(13):(mb) The averaged total cross section of nucleon
0404 nucleon collisions.
0405 HINT1(14):(mb) The jet production cross section without
0406 shadowing effect.
0407 HINT1(15):(mb)The cross section to account for the projectile
0408 shadowing correction term in the jet cross section.
0409 HINT1(16):(mb)The cross section to account for the target
0410 shadowing correction term in the jet cross section.
0411 HINT1(17):(mb)The cross section to account for the cross
0412 term of shadowing correction in the jet cross section.
0413 HINT1(18):(mb)The cross section for jet production for
0414 the latest nucleon-nucleon collisions which
0415 depends on the impact parameters(i.e. HINT1(18)=
0416 HINT1(14)-a1*HINT1(15)-a2*HINT1(16)+a1*a2*HINT1(17)).
0417 HINT1(19):(fm) The imapct parameter of the latest event.
0418 HINT1(20):(rad)The orientation of the latest collision.
0419 HINT1(21)-HINT1(25):the three momenta, energy and mass of the
0420 first scattered parton in the specified hard collision
0421 per event with specified PT (This is before the final
0422 state radiation but after the initial state radiation).
0423 HINT1(26)-HINT1(29): The total momentum and energy of the
0424 partons in the specified hard collision from the
0425 final state radiation of the first
0426 scattered parton which are a little different from
0427 HINT1(21)-HINT1(24) due to the rearrangement of
0428 the parton showers with energy and momentum
0429 conservation.
0430 HINT1(30):not used.
0431 HINT1(31)-HINT1(35):the three momenta, energy and mass of the
0432 second scattered parton in the specified hard collision
0433 per event with specified PT (This is before the final
0434 state radiation but after the initial state radiation).
0435 HINT1(36)-HINT1(39):The total momentum and energy of the
0436 partons in the specified hard collision from the
0437 final state radiation of the second
0438 scattered parton which are a little different from
0439 HINT1(31)-HINT1(34) due to the rearrangement of
0440 the parton showers with energy and momentum
0441 conservation.
0442 HINT1(40):not used.
0443 HINT1(41)-HINT1(45)the three momenta, energy and mass of the
0444 first scattered parton in the latest hard collision
0445 of the latest event.
0446 HINT1(46): PT of the first scattered parton in the latest hard
0447 collision of the latest event.
0448 HINT1(47)-HINT1(50):The total momentum and energy of the
0449 partons in the latest hard collision from the
0450 final state radiation of the second
0451 scattered parton which are a little different from
0452 HINT1(41)-HINT1(44) due to the rearrangement of
0453 the parton showers with energy and momentum
0454 conservation.
0455 HINT1(51)-HINT1(55)the three momenta, energy and mass of the
0456 second scattered parton in the latest hard collision
0457 of the latest event.
0458 HINT1(56): PT of the second scattered parton in the latest
0459 hard collision of the latest event.
0460 HINT1(57):Correction factor for cross section of parton_nucleon
0461 interactins due to processes below pt<pt_cut(HIPR1(8)).
0462 HINT1(58):(mb) The averaged cross section of a parton
0463 scattering off a nucleon which is used to simulate
0464 initial state interaction.
0465 HINT1(59):(mb) The averaged total cross section of the
0466 triggered jet production with specified P_T by
0467 HIPR1(10)(and switch by IHPR2(3))(usually, it is
0468 very close to the inclusive cross section HINT1(60)
0469 when the later is very small).
0470 HINT1(60):(mb)The averaged inclusive cross section of the
0471 triggered jet production with specified P_T by
0472 HIPR1(10)(and switch by IHPR2(3)).
0473 HINT1(61):(mb) The trigger jet production cross section without
0474 shadowing effect.
0475 HINT1(62):(mb)The cross section to account for the projectile
0476 shadowing correction term in the trigger jet cross
0477 section.
0478 HINT1(63):(mb)The cross section to account for the target
0479 shadowing correction term in the trigger jet cross
0480 section.
0481 HINT1(64):(mb)The cross section to account for the cross
0482 term of shadowing correction in the trigger jet
0483 cross section.
0484 HINT1(65):(mb)The inclusive cross section for latest trigger
0485 jet production which depends on the impact parameters
0486 (i.e. HINT1(65)=HINT1(61)-a1*HINT1(62)-a2*HINT1(63)
0487 +a1*a2*HINT1(64)).
0488 HINT1(67)-HINT1(70):The total momentum and energy of the
0489 partons in the latest hard collision from the
0490 final state radiation of the second
0491 scattered parton which are a little different from
0492 HINT1(51)-HINT1(54) due to the rearrangement of
0493 the parton showers with energy and momentum
0494 conservation.
0495 HINT1(71):not used.
0496 HINT1(72)-HINT1(75): Three parameters for the wood-saxon
0497 projectile nuclear distribution and the normalization,
0498 rho(x)=FNORM*(1.+W*(X/R)**2)/(1+EXP((X-R)/D)
0499 R=HINT1(72), D=HINT1(73), W=HINT1(74), FNORM=HINT1(75).
0500 HINT1(76)-HINT1(79): Three parameters for the wood-saxon
0501 target nuclear distribution and the normalization,
0502 rho(x)=FNORM*(1.+W*(X/R)**2)/(1+EXP((X-R)/D)
0503 R=HINT1(76), D=HINT1(77), W=HINT1(78), FNORM=HINT1(79).
0504 HINT1(80)-HINT1(100): The probability of 0-20 jets production
0505 per pp interaction.
0506
0507
0508
0509 IHNT2(1):The atomic number of the projectile nucleus
0510 (=1 for hadron).
0511 IHNT2(2):The charge number of the projectile nucleus, if the
0512 projectile is a hadron it gives the charge of that
0513 hadron.
0514 IHNT2(3):The atomic number of the target nucleus(=1 for hadron).
0515 IHNT2(4):The charge number of the terget nucleus,if the
0516 target is a hadron it gives the charge of
0517 that hadron.
0518 IHNT2(5):The ID number of projectile hadrons(=0 for nucleus).
0519 IHNT2(6):The ID number of target hadrons(=0 for nucleus).
0520 IHNT2(7)-IHNT2(8):not used
0521 IHNT2(9):The parton id of the first scattered parton in the
0522 specified hard collision per event with fixed PT.
0523 IHNT2(10):The parton id of the second scattered parton in the
0524 specified hard collision per event with fixed PT.
0525 IHNT2(11):The sequence number of projectile nucleon of the
0526 latest nucleon-nucleon interaction in the latest event.
0527 IHNT2(12):The sequence number of target nucleon of the latest
0528 nucleon-nucleon interaction in the present event.
0529 IHNT2(13):status of the latest soft LUND type exicitation:
0530 =1 double diffractive;
0531 =2 single diffractive
0532 =3 non-single diffractive.
0533 IHNT2(14):The parton id of the first scattered parton in the
0534 latest hard collision of the latest event.
0535 IHNT2(15):The parton id of the second scattered parton in the
0536 latest hard collision of the latest event.
0537 IHNT2(16):parameter to indicate the initialization of the
0538 the cross sections(=1). When given value 1
0539 the initilization of pythia will used un-shadowed
0540 structure functions to evaluate the maximum
0541 cross sections.
0542 IHNT2(17)-IHNT2(50):not used.
0543
0544
0545 (5)Example program
0546 C****************************************************************************
0547 C
0548 C
0549 C
0550 C The following is an example program for calling HIJING. one should
0551 C include all the common blocks and the data values which are listed
0552 C below in his own program.
0553 C***************************************************************************
0554 CHARACTER FRAME*8,PROJ*8,TARG*8
0555 COMMON/HIMAIN1/ NATT,EATT,JATT,NT,NP,N0,N01,N10,N11
0556 COMMON/HIMAIN2/KATT(130000,4),PATT(130000,4)
0557 C ********information of produced particles
0558 C
0559 COMMON/HIJJET1/NPJ(300),KFPJ(300,500),PJPX(300,500),PJPY(300,500)
0560 & ,PJPZ(300,500),PJPE(300,500),PJPM(300,500)
0561 & ,NTJ(300),KFTJ(300,500),PJTX(300,500),PJTY(300,500)
0562 & ,PJTZ(300,500),PJTE(300,500),PJTM(300,500)
0563 COMMON/HIJJET2/NSG,NJSG(900),IASG(900,3),K1SG(900,100),K2SG(900,100)
0564 & ,PXSG(900,100),PYSG(900,100),PZSG(900,100),PESG(900,100),PMSG(900,100)
0565 C ********information of produced partons
0566 EFRM=200.0
0567 FRAME='CMS'
0568 PROJ='A'
0569 TARG='A'
0570 IAP=197
0571 IZP=79
0572 IAT=197
0573 IZT=79
0574 CALL HIJSET(EFRM,FRAME,PROJ,TARG,IAP,IZP,IAT,IZT)
0575 C ********Initialize HIJING
0576
0577 BMIN=0.0
0578 BMAX=0.0
0579 DO 2000 I_event=1,100
0580
0581 CALL HIJING(FRAME,BMIN,BMAX)
0582 C
0583 DO 1000 I=1,NATT
0584
0585 IF(LUCHGE(KATT(I,1)).NE.0) THEN
0586 C this select charged particles only
0587 C !information of produced particles
0588 C !is stored in common blocks HIMAIN1 and
0589 C HIMAIN2
0590 ENDIF
0591 1000 CONTINUE
0592
0593 2000 CONTINUE
0594 STOP
0595 END
0596 c****************************************************************
0597
0598 If one want to trig on jets with specified high P_T, the
0599 following should be included in the program
0600
0601 IHPR2(3)=1(2,for direct photon production)
0602 HIPR1(10)=P_T(-P_T for jets above P_T)
0603
0604 Because large P_T jet has very small cross sections, this will
0605 improve the simulation efficiency. The absolute cross section
0606 can be obtained by multiplying the results with the total
0607 cross section of such jet production in pp which is stored in
0608 HINT1(59).
0609
0610
0611
0612
0613 (6)Appendix: For user's reference a selection of particle ID codes from
0614 JETSET7.2 are listed below. For full list please check JETSET users
0615 instruction. The codes for anti-particles are just the negative
0616 values of the corresponding particles.
0617
0618 Quarks and leptons
0619
0620 1 d 11 e-
0621 2 u 12 nu_e
0622 3 s 13 mu-
0623 4 c 14 nu_mu
0624 5 b 15 tau-
0625 6 t 16 nu_tau
0626
0627 Gauge bosons
0628
0629 21 gluon
0630 22 gamma
0631
0632 Diquarks
0633 1103 dd_1
0634 2101 ud_0 2103 ud_1
0635 2203 uu_1
0636
0637 3101 sd_0 3103 sd_1
0638 3201 su_0 3203 su_1
0639 3303 ss_1
0640
0641 Mesons
0642 211 pi+ 213 rho+
0643 311 K0 313 K*0
0644 321 K+ 323 K*+
0645 411 D+ 413 D*+
0646 421 D0 423 D*0
0647 431 D_s+ 433 D*_s+
0648 511 B0 513 B*0
0649 521 B+ 523 B*+
0650 531 B_s0 533 B*_s0
0651
0652 111 pi0 113 rho0
0653 221 eta 223 omega
0654 331 eta' 333 phi
0655 441 eta_c 443 J/psi
0656 551 eta_b 553 Upsilon
0657 661 eta_t 663 Theta
0658
0659 130 K_L0
0660 310 K_S0
0661
0662 Baryons
0663 1114 Delta-
0664 2112 n 2114 Delta0
0665 2212 p 2214 Delta+
0666 2224 Delta++
0667 3112 Sigma- 3114 Sigma*-
0668 3122 Lambda0
0669 3212 Sigma0 3214 Sigma*0
0670 3222 Sigma+ 3224 Sigma*+
0671 3312 Xi- 3314 Xi*-
0672 3322 Xi0 3324 Xi*0
0673 3334 Omega-
0674 4112 Sigma_c0 4114 Sigma*_c0
0675 4122 Lambda_c+
0676 4212 Sigma_c+ 4214 Sigma*_c+
0677 4222 Sigma_c++ 4224 Sigma*_c++
0678 4132 Xi_c0
0679 4312 Xi'_c0 4314 Xi*_c0
0680 4232 Xi_c+
0681 4322 Xi'_c+ 4324 Xi*_c+
0682 4332 Omega_c0 4334 Omega*_c0
0683 5112 Sigma_b- 5114 Sigma*_b-
0684 5122 Lambda_b0
0685 5212 Sigma_b0 5214 Sigma*_b0
0686 5222 Sigma_b+ 5224 Sigma*_b+
0687
