]> git.uio.no Git - u/mrichter/AliRoot.git/blame - HIJING/hijing1_36/jetini.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
changed the order of call of endofcycle so that images are produced
[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / jetini.F
CommitLineData
e74335a4 1* $Id$
2C
3C
4C
5C
6C
7 SUBROUTINE JETINI(JP,JT,I_TRIG)
8C*******Initialize PYTHIA for jet production**********************
9C I_TRIG=0: for normal processes
10C I_TRIG=1: for triggered processes
11C JP: sequence number of the projectile
12C JT: sequence number of the target
13C For A+A collisions, one has to initilize pythia
14C separately for each type of collisions, pp, pn,np and nn,
15C or hp and hn for hA collisions. In this subroutine we use the following
16C catalogue for different type of collisions:
17C h+h: h+h (I_TYPE=1)
18C h+A: h+p (I_TYPE=1), h+n (I_TYPE=2)
19C A+h: p+h (I_TYPE=1), n+h (I_TYPE=2)
20C A+A: p+p (I_TYPE=1), p+n (I_TYPE=2), n+p (I_TYPE=3), n+n (I_TYPE=4)
21C*****************************************************************
22 CHARACTER BEAM*16,TARG*16
23 DIMENSION XSEC0(8,0:200),COEF0(8,200,20),INI(8),
24 & MINT44(8),MINT45(8)
bc676b8e 25#define BLANKET_SAVE
e74335a4 26#include "hijcrdn.inc"
27#include "hiparnt.inc"
28#include "histrng.inc"
29#include "hipyint.inc"
30C
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/
39 SAVE
40C
41 IHNT2(11)=JP
42 IHNT2(12)=JT
43 IF(IHNT2(5).NE.0 .AND. IHNT2(6).NE.0) THEN
44 I_TYPE=1
45 ELSE IF(IHNT2(5).NE.0 .AND. IHNT2(6).EQ.0) THEN
46 I_TYPE=1
47 IF(NFT(JT,4).EQ.2112) I_TYPE=2
48 ELSE IF(IHNT2(5).EQ.0 .AND. IHNT2(6).NE.0) THEN
49 I_TYPE=1
50 IF(NFP(JP,4).EQ.2112) I_TYPE=2
51 ELSE
52 IF(NFP(JP,4).EQ.2212 .AND. NFT(JT,4).EQ.2212) THEN
53 I_TYPE=1
54 ELSE IF(NFP(JP,4).EQ.2212 .AND. NFT(JT,4).EQ.2112) THEN
55 I_TYPE=2
56 ELSE IF(NFP(JP,4).EQ.2112 .AND. NFT(JT,4).EQ.2212) THEN
57 I_TYPE=3
58 ELSE
59 I_TYPE=4
60 ENDIF
61 ENDIF
62c
385f8e32 63
e74335a4 64 IF(I_TRIG.NE.0) GO TO 160
65 IF(I_TRIG.EQ.I_LAST) GO TO 150
66 MSTP(2)=2
67c ********second order running alpha_strong
68 MSTP(33)=1
69 PARP(31)=HIPR1(17)
70C ********inclusion of K factor
71 MSTP(51)=3
72C ********Duke-Owens set 1 structure functions
73 MSTP(61)=1
74C ********INITIAL STATE RADIATION
75 MSTP(71)=1
76C ********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
25d2dbca 79 MSTP(91) = 1
e74335a4 80c
81 MSTP(81)=0
82C ******** NO MULTIPLE INTERACTION
83 MSTP(82)=1
84C *******STRUCTURE OF MUTLIPLE INTERACTION
85 MSTP(111)=0
86C ********frag off(have to be done by local call)
87 IF(IHPR2(10).EQ.0) MSTP(122)=0
88C ********No printout of initialization information
89 PARP(81)=HIPR1(8)
90 CKIN(5)=HIPR1(8)
91 CKIN(3)=HIPR1(8)
92 CKIN(4)=HIPR1(9)
93 IF(HIPR1(9).LE.HIPR1(8)) CKIN(4)=-1.0
94 CKIN(9)=-10.0
95 CKIN(10)=10.0
96 MSEL=0
97 DO 100 ISUB=1,200
98 MSUB(ISUB)=0
99 100 CONTINUE
100 MSUB(11)=1
101 MSUB(12)=1
102 MSUB(13)=1
103 MSUB(28)=1
104 MSUB(53)=1
105 MSUB(68)=1
385f8e32 106 MSUB(81)=1
107 MSUB(82)=1
e74335a4 108 DO 110 J=1,MIN(8,MDCY(21,3))
109 110 MDME(MDCY(21,2)+J-1,1)=0
110 ISEL=4
111 IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
112 MDME(MDCY(21,2)+ISEL-1,1)=1
113C ********QCD subprocesses
114 MSUB(14)=1
115 MSUB(18)=1
116 MSUB(29)=1
117C ******* direct photon production
118 150 IF(INI(I_TYPE).NE.0) GO TO 800
119 GO TO 400
120C
121C *****triggered subprocesses, jet, photon, heavy quark and DY
122C
123 160 I_TYPE=4+I_TYPE
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
25d2dbca 130C WRITE(6,*) "TRIGGERED", CKIN(3), CKIN(4), CKIN(5), IHPR2(3)
131C
132C ALICE
133C
134 IF (HIPR1(50) .EQ. 1) THEN
135 MSTP(61) = 0
136 MSTP(71) = 0
137 MSTP(91) = 0
138 ENDIF
139
e74335a4 140c
141 MSEL=0
142 DO 101 ISUB=1,200
143 MSUB(ISUB)=0
144 101 CONTINUE
145 IF(IHPR2(3).EQ.1) THEN
146 MSUB(11)=1
147 MSUB(12)=1
148 MSUB(13)=1
149 MSUB(28)=1
150 MSUB(53)=1
151 MSUB(68)=1
152 MSUB(81)=1
153 MSUB(82)=1
154 MSUB(14)=1
155 MSUB(18)=1
156 MSUB(29)=1
157 DO 102 J=1,MIN(8,MDCY(21,3))
158 102 MDME(MDCY(21,2)+J-1,1)=0
159 ISEL=4
160 IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
161 MDME(MDCY(21,2)+ISEL-1,1)=1
162C ********QCD subprocesses
163 ELSE IF(IHPR2(3).EQ.2) THEN
25d2dbca 164C WRITE(6,*) "DIRECT PHOTON PRODUCTION"
e74335a4 165 MSUB(14)=1
166 MSUB(18)=1
167 MSUB(29)=1
168C ********Direct photon production
169c 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))
172 CKIN(5)=HIPR1(8)
173 PARP(81)=HIPR1(8)
174 MSUB(81)=1
175 MSUB(82)=1
176 DO 105 J=1,MIN(8,MDCY(21,3))
177 105 MDME(MDCY(21,2)+J-1,1)=0
178 ISEL=4
179 IF(HINT1(1).GE.20.0 .and. IHPR2(18).EQ.1) ISEL=5
180 MDME(MDCY(21,2)+ISEL-1,1)=1
181C **********Heavy quark production
182 ENDIF
183260 IF(INI(I_TYPE).NE.0) GO TO 800
184C
185C
186400 INI(I_TYPE)=1
187 IF(IHPR2(10).EQ.0) MSTP(122)=0
188 IF(NFP(JP,4).EQ.2212) THEN
189 BEAM='P'
190 ELSE IF(NFP(JP,4).EQ.-2212) THEN
191 BEAM='P~'
192 ELSE IF(NFP(JP,4).EQ.2112) THEN
193 BEAM='N'
194 ELSE IF(NFP(JP,4).EQ.-2112) THEN
195 BEAM='N~'
196 ELSE IF(NFP(JP,4).EQ.211) THEN
197 BEAM='PI+'
198 ELSE IF(NFP(JP,4).EQ.-211) THEN
199 BEAM='PI-'
200 ELSE IF(NFP(JP,4).EQ.321) THEN
201 BEAM='PI+'
202 ELSE IF(NFP(JP,4).EQ.-321) THEN
203 BEAM='PI-'
204 ELSE
205 WRITE(6,*) 'unavailable beam type', NFP(JP,4)
206 ENDIF
207 IF(NFT(JT,4).EQ.2212) THEN
208 TARG='P'
209 ELSE IF(NFT(JT,4).EQ.-2212) THEN
210 TARG='P~'
211 ELSE IF(NFT(JT,4).EQ.2112) THEN
212 TARG='N'
213 ELSE IF(NFT(JT,4).EQ.-2112) THEN
214 TARG='N~'
215 ELSE IF(NFT(JT,4).EQ.211) THEN
216 TARG='PI+'
217 ELSE IF(NFT(JT,4).EQ.-211) THEN
218 TARG='PI-'
219 ELSE IF(NFT(JT,4).EQ.321) THEN
220 TARG='PI+'
221 ELSE IF(NFT(JT,4).EQ.-321) THEN
222 TARG='PI-'
223 ELSE
224 WRITE(6,*) 'unavailable target type', NFT(JT,4)
225 ENDIF
226C
227 IHNT2(16)=1
228C ******************indicate for initialization use when
229C structure functions are called in PYTHIA
230C
25d2dbca 231C WRITE(6,*) "CALL PYINIT WITH", CKIN(3), CKIN(4), MSTP(61), MSTP(71), MSTP(91)
385f8e32 232C ALICE ->
233 IF (IHPR2(49) .EQ. 1) THEN
234C SWITCH OFF MASSIVE FINAL QUARKS
235 MSEL = 0
236C QCD on
237 MSUB(11) = 1
238 MSUB(12) = 1
239 MSUB(13) = 1
240 MSUB(28) = 1
241 MSUB(53) = 1
242 MSUB(68) = 1
243C HF OFF
244 MSUB(81) = 0
245 MSUB(82) = 0
246 MDME(127, 1) = 1
247 MDME(128, 1) = 1
248 MDME(129, 1) = 1
249C g->QQbar off
250 MDME(130, 1) = 0
251 MDME(131, 1) = 0
252C highest flavor in final state shower
253 MSTJ(45) = 3
254C highest flavor in initial state shower and pdf
255 MSTP(54) = 3
256 ENDIF
257C <- ALICE
e74335a4 258 CALL PYINIT_HIJING('CMS',BEAM,TARG,HINT1(1))
259 MINT4=MINT(44)
260 MINT5=MINT(45)
261 MINT44(I_TYPE)=MINT(44)
262 MINT45(I_TYPE)=MINT(45)
263 ATXS(0)=XSEC(0,1)
264 XSEC0(I_TYPE,0)=XSEC(0,1)
265 DO 500 I=1,200
266 ATXS(I)=XSEC(I,1)
267 XSEC0(I_TYPE,I)=XSEC(I,1)
268 DO 500 J=1,20
269 ATCO(I,J)=COEF(I,J)
270 COEF0(I_TYPE,I,J)=COEF(I,J)
271500 CONTINUE
272C
273 IHNT2(16)=0
274C
57367a3d 275 I_LAST=I_TRIG
e74335a4 276 RETURN
277C ********Store the initialization information for
278C late use
279C
280C
281800 MINT(44)=MINT44(I_TYPE)
282 MINT(45)=MINT45(I_TYPE)
283 MINT4=MINT(44)
284 MINT5=MINT(45)
285 XSEC(0,1)=XSEC0(I_TYPE,0)
286 ATXS(0)=XSEC(0,1)
287 DO 900 I=1,200
288 XSEC(I,1)=XSEC0(I_TYPE,I)
289 ATXS(I)=XSEC(I,1)
290 DO 900 J=1,20
291 COEF(I,J)=COEF0(I_TYPE,I,J)
292 ATCO(I,J)=COEF(I,J)
293900 CONTINUE
294 I_LAST=I_TRIG
295 MINT(11)=NFP(JP,4)
296 MINT(12)=NFT(JT,4)
297 RETURN
298 END
25d2dbca 299
300
301
302
303
304
305
306
307