7 SUBROUTINE JETINI(JP,JT,I_TRIG)
8 C*******Initialize PYTHIA for jet production**********************
9 C I_TRIG=0: for normal processes
10 C I_TRIG=1: for triggered processes
11 C JP: sequence number of the projectile
12 C JT: sequence number of the target
13 C For A+A collisions, one has to initilize pythia
14 C separately for each type of collisions, pp, pn,np and nn,
15 C or hp and hn for hA collisions. In this subroutine we use the following
16 C catalogue for different type of collisions:
18 C h+A: h+p (I_TYPE=1), h+n (I_TYPE=2)
19 C A+h: p+h (I_TYPE=1), n+h (I_TYPE=2)
20 C A+A: p+p (I_TYPE=1), p+n (I_TYPE=2), n+p (I_TYPE=3), n+n (I_TYPE=4)
21 C*****************************************************************
22 CHARACTER BEAM*16,TARG*16
23 DIMENSION XSEC0(8,0:200),COEF0(8,200,20),INI(8),
26 #include "hijcrdn.inc"
27 #include "hiparnt.inc"
28 #include "histrng.inc"
29 #include "hipyint.inc"
31 #include "ludat1_hijing.inc"
32 #include "ludat3_hijing.inc"
33 #include "pysubs_hijing.inc"
34 #include "pypars_hijing.inc"
35 #include "pyint1_hijing.inc"
36 #include "pyint2_hijing.inc"
37 #include "pyint5_hijing.inc"
38 DATA INI/8*0/I_LAST/-1/
43 IF(IHNT2(5).NE.0 .AND. IHNT2(6).NE.0) THEN
45 ELSE IF(IHNT2(5).NE.0 .AND. IHNT2(6).EQ.0) THEN
47 IF(NFT(JT,4).EQ.2112) I_TYPE=2
48 ELSE IF(IHNT2(5).EQ.0 .AND. IHNT2(6).NE.0) THEN
50 IF(NFP(JP,4).EQ.2112) I_TYPE=2
52 IF(NFP(JP,4).EQ.2212 .AND. NFT(JT,4).EQ.2212) THEN
54 ELSE IF(NFP(JP,4).EQ.2212 .AND. NFT(JT,4).EQ.2112) THEN
56 ELSE IF(NFP(JP,4).EQ.2112 .AND. NFT(JT,4).EQ.2212) THEN
64 IF(I_TRIG.NE.0) GO TO 160
65 IF(I_TRIG.EQ.I_LAST) GO TO 150
67 c ********second order running alpha_strong
70 C ********inclusion of K factor
72 C ********Duke-Owens set 1 structure functions
74 C ********INITIAL STATE RADIATION
76 C ********FINAL STATE RADIATION
77 IF(IHPR2(2).EQ.0.OR.IHPR2(2).EQ.2) MSTP(61)=0
78 IF(IHPR2(2).EQ.0.OR.IHPR2(2).EQ.1) MSTP(71)=0
82 C ******** NO MULTIPLE INTERACTION
84 C *******STRUCTURE OF MUTLIPLE INTERACTION
86 C ********frag off(have to be done by local call)
87 IF(IHPR2(10).EQ.0) MSTP(122)=0
88 C ********No printout of initialization information
93 IF(HIPR1(9).LE.HIPR1(8)) CKIN(4)=-1.0
108 DO 110 J=1,MIN(8,MDCY(21,3))
109 110 MDME(MDCY(21,2)+J-1,1)=0
111 IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
112 MDME(MDCY(21,2)+ISEL-1,1)=1
113 C ********QCD subprocesses
117 C ******* direct photon production
118 150 IF(INI(I_TYPE).NE.0) GO TO 800
121 C *****triggered subprocesses, jet, photon, heavy quark and DY
124 IF(I_TRIG.EQ.I_LAST) GO TO 260
125 PARP(81)=ABS(HIPR1(10))-0.25
126 CKIN(5)=ABS(HIPR1(10))-0.25
127 CKIN(3)=ABS(HIPR1(10))-0.25
128 CKIN(4)=ABS(HIPR1(10))+0.25
129 IF(HIPR1(10).LT.HIPR1(8)) CKIN(4)=-1.0
130 C WRITE(6,*) "TRIGGERED", CKIN(3), CKIN(4), CKIN(5), IHPR2(3)
134 IF (HIPR1(50) .EQ. 1) THEN
145 IF(IHPR2(3).EQ.1) THEN
157 DO 102 J=1,MIN(8,MDCY(21,3))
158 102 MDME(MDCY(21,2)+J-1,1)=0
160 IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
161 MDME(MDCY(21,2)+ISEL-1,1)=1
162 C ********QCD subprocesses
163 ELSE IF(IHPR2(3).EQ.2) THEN
164 C WRITE(6,*) "DIRECT PHOTON PRODUCTION"
168 C ********Direct photon production
169 c q+qbar->g+gamma,q+qbar->gamma+gamma, q+g->q+gamma
170 ELSE IF(IHPR2(3).EQ.3) THEN
171 CKIN(3)=MAX(0.0,HIPR1(10))
176 DO 105 J=1,MIN(8,MDCY(21,3))
177 105 MDME(MDCY(21,2)+J-1,1)=0
179 IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
180 MDME(MDCY(21,2)+ISEL-1,1)=1
181 C **********Heavy quark production
183 260 IF(INI(I_TYPE).NE.0) GO TO 800
187 IF(IHPR2(10).EQ.0) MSTP(122)=0
188 IF(NFP(JP,4).EQ.2212) THEN
190 ELSE IF(NFP(JP,4).EQ.-2212) THEN
192 ELSE IF(NFP(JP,4).EQ.2112) THEN
194 ELSE IF(NFP(JP,4).EQ.-2112) THEN
196 ELSE IF(NFP(JP,4).EQ.211) THEN
198 ELSE IF(NFP(JP,4).EQ.-211) THEN
200 ELSE IF(NFP(JP,4).EQ.321) THEN
202 ELSE IF(NFP(JP,4).EQ.-321) THEN
205 WRITE(6,*) 'unavailable beam type', NFP(JP,4)
207 IF(NFT(JT,4).EQ.2212) THEN
209 ELSE IF(NFT(JT,4).EQ.-2212) THEN
211 ELSE IF(NFT(JT,4).EQ.2112) THEN
213 ELSE IF(NFT(JT,4).EQ.-2112) THEN
215 ELSE IF(NFT(JT,4).EQ.211) THEN
217 ELSE IF(NFT(JT,4).EQ.-211) THEN
219 ELSE IF(NFT(JT,4).EQ.321) THEN
221 ELSE IF(NFT(JT,4).EQ.-321) THEN
224 WRITE(6,*) 'unavailable target type', NFT(JT,4)
228 C ******************indicate for initialization use when
229 C structure functions are called in PYTHIA
231 C WRITE(6,*) "CALL PYINIT WITH", CKIN(3), CKIN(4), MSTP(61), MSTP(71), MSTP(91)
233 IF (IHPR2(49) .EQ. 1) THEN
234 C SWITCH OFF MASSIVE FINAL QUARKS
252 C highest flavor in final state shower
254 C highest flavor in initial state shower and pdf
258 CALL PYINIT_HIJING('CMS',BEAM,TARG,HINT1(1))
261 MINT44(I_TYPE)=MINT(44)
262 MINT45(I_TYPE)=MINT(45)
264 XSEC0(I_TYPE,0)=XSEC(0,1)
267 XSEC0(I_TYPE,I)=XSEC(I,1)
270 COEF0(I_TYPE,I,J)=COEF(I,J)
277 C ********Store the initialization information for
281 800 MINT(44)=MINT44(I_TYPE)
282 MINT(45)=MINT45(I_TYPE)
285 XSEC(0,1)=XSEC0(I_TYPE,0)
288 XSEC(I,1)=XSEC0(I_TYPE,I)
291 COEF(I,J)=COEF0(I_TYPE,I,J)