]>
Commit | Line | Data |
---|---|---|
79a689ca | 1 | C $Id$ |
e74335a4 | 2 | C |
3 | C | |
4 | C | |
5 | C | |
6 | SUBROUTINE QUENCH(JPJT,NTP) | |
7 | DIMENSION RDP(300),LQP(300),RDT(300),LQT(300) | |
8 | #include "hijcrdn.inc" | |
9 | #include "hiparnt.inc" | |
10 | C | |
11 | #include "hijjet1.inc" | |
12 | #include "hijjet2.inc" | |
13 | #include "histrng.inc" | |
14 | C | |
15 | SAVE | |
9804fb1d | 16 | C |
17 | BB=HINT1(19) ! Uzhi | |
18 | PHI=HINT1(20) ! Uzhi | |
19 | BBX=BB*COS(PHI) ! Uzhi | |
20 | BBY=BB*SIN(PHI) ! Uzhi | |
21 | c | |
e74335a4 | 22 | IF(NTP.EQ.2) GO TO 400 |
23 | IF(NTP.EQ.3) GO TO 2000 | |
24 | C******************************************************* | |
25 | C Jet interaction for proj jet in the direction PHIP | |
26 | C****************************************************** | |
27 | C | |
28 | IF(NFP(JPJT,7).NE.1) RETURN | |
29 | ||
30 | JP=JPJT | |
31 | DO 290 I=1,NPJ(JP) | |
32 | PTJET0=SQRT(PJPX(JP,I)**2+PJPY(JP,I)**2) | |
33 | IF(PTJET0.LE.HIPR1(11)) GO TO 290 | |
34 | PTOT=SQRT(PTJET0*PTJET0+PJPZ(JP,I)**2) | |
79a689ca | 35 | IF (IHPR2(50) .EQ. 1) THEN |
36 | DEDX0 = HIPR1(14)*LOG10(PTJET0)/LOG10(5.) | |
37 | ELSE | |
38 | DEDX0 = HIPR1(14) | |
39 | ENDIF | |
e74335a4 | 40 | IF(PTOT.LT.HIPR1(8)) GO TO 290 |
41 | PHIP=ULANGL_HIJING(PJPX(JP,I),PJPY(JP,I)) | |
42 | C******* find the wounded proj which can interact with jet*** | |
43 | KP=0 | |
44 | DO 100 I2=1,IHNT2(1) | |
45 | IF(NFP(I2,5).NE.3 .OR. I2.EQ.JP) GO TO 100 | |
46 | DX=YP(1,I2)-YP(1,JP) | |
47 | DY=YP(2,I2)-YP(2,JP) | |
48 | PHI=ULANGL_HIJING(DX,DY) | |
49 | DPHI=ABS(PHI-PHIP) | |
9804fb1d | 50 | IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi |
e74335a4 | 51 | IF(DPHI.GE.HIPR1(40)/2.0) GO TO 100 |
52 | RD0=SQRT(DX*DX+DY*DY) | |
53 | IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 100 | |
54 | KP=KP+1 | |
55 | LQP(KP)=I2 | |
56 | RDP(KP)=COS(DPHI)*RD0 | |
57 | 100 CONTINUE | |
58 | C******* rearrange according decending rd************ | |
59 | DO 110 I2=1,KP-1 | |
60 | DO 110 J2=I2+1,KP | |
61 | IF(RDP(I2).LT.RDP(J2)) GO TO 110 | |
62 | RD=RDP(I2) | |
63 | LQ=LQP(I2) | |
64 | RDP(I2)=RDP(J2) | |
65 | LQP(I2)=LQP(J2) | |
66 | RDP(J2)=RD | |
67 | LQP(J2)=LQ | |
68 | 110 CONTINUE | |
69 | C****** find wounded targ which can interact with jet******** | |
70 | KT=0 | |
71 | DO 120 I2=1,IHNT2(3) | |
72 | IF(NFT(I2,5).NE.3) GO TO 120 | |
73 | DX=YT(1,I2)-YP(1,JP) | |
74 | DY=YT(2,I2)-YP(2,JP) | |
75 | PHI=ULANGL_HIJING(DX,DY) | |
76 | DPHI=ABS(PHI-PHIP) | |
9804fb1d | 77 | IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi |
e74335a4 | 78 | IF(DPHI.GT.HIPR1(40)/2.0) GO TO 120 |
79 | RD0=SQRT(DX*DX+DY*DY) | |
80 | IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 120 | |
81 | KT=KT+1 | |
82 | LQT(KT)=I2 | |
83 | RDT(KT)=COS(DPHI)*RD0 | |
84 | 120 CONTINUE | |
85 | C******* rearrange according decending rd************ | |
86 | DO 130 I2=1,KT-1 | |
87 | DO 130 J2=I2+1,KT | |
88 | IF(RDT(I2).LT.RDT(J2)) GO TO 130 | |
89 | RD=RDT(I2) | |
90 | LQ=LQT(I2) | |
91 | RDT(I2)=RDT(J2) | |
92 | LQT(I2)=LQT(J2) | |
93 | RDT(J2)=RD | |
94 | LQT(J2)=LQ | |
95 | 130 CONTINUE | |
96 | ||
97 | MP=0 | |
98 | MT=0 | |
99 | R0=0.0 | |
100 | NQ=0 | |
101 | DP=0.0 | |
102 | PTOT=SQRT(PJPX(JP,I)**2+PJPY(JP,I)**2+PJPZ(JP,I)**2) | |
103 | V1=PJPX(JP,I)/PTOT | |
104 | V2=PJPY(JP,I)/PTOT | |
105 | V3=PJPZ(JP,I)/PTOT | |
106 | ||
107 | 200 RN=RLU_HIJING(0) | |
108 | 210 IF(MT.GE.KT .AND. MP.GE.KP) GO TO 290 | |
109 | IF(MT.GE.KT) GO TO 220 | |
110 | IF(MP.GE.KP) GO TO 240 | |
111 | IF(RDP(MP+1).GT.RDT(MT+1)) GO TO 240 | |
112 | 220 MP=MP+1 | |
113 | DRR=RDP(MP)-R0 | |
114 | IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 210 | |
79a689ca | 115 | DP=DRR*DEDX0 |
e74335a4 | 116 | IF(KFPJ(JP,I).NE.21) DP=0.5*DP |
117 | C ********string tension of quark jet is 0.5 of gluon's | |
118 | IF(DP.LE.0.2) GO TO 210 | |
119 | IF(PTOT.LE.0.4) GO TO 290 | |
120 | IF(PTOT.LE.DP) DP=PTOT-0.2 | |
121 | DE=DP | |
122 | ||
123 | IF(KFPJ(JP,I).NE.21) THEN | |
124 | PRSHU=PP(LQP(MP),1)**2+PP(LQP(MP),2)**2 | |
125 | & +PP(LQP(MP),3)**2 | |
126 | DE=SQRT(PJPM(JP,I)**2+PTOT**2) | |
127 | & -SQRT(PJPM(JP,I)**2+(PTOT-DP)**2) | |
128 | ERSHU=(PP(LQP(MP),4)+DE-DP)**2 | |
129 | AMSHU=ERSHU-PRSHU | |
130 | IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 210 | |
131 | PP(LQP(MP),4)=SQRT(ERSHU) | |
132 | PP(LQP(MP),5)=SQRT(AMSHU) | |
133 | ENDIF | |
134 | C ********reshuffle the energy when jet has mass | |
135 | R0=RDP(MP) | |
136 | DP1=DP*V1 | |
137 | DP2=DP*V2 | |
138 | DP3=DP*V3 | |
139 | C ********momentum and energy transfer from jet | |
140 | ||
141 | NPJ(LQP(MP))=NPJ(LQP(MP))+1 | |
142 | KFPJ(LQP(MP),NPJ(LQP(MP)))=21 | |
143 | PJPX(LQP(MP),NPJ(LQP(MP)))=DP1 | |
144 | PJPY(LQP(MP),NPJ(LQP(MP)))=DP2 | |
145 | PJPZ(LQP(MP),NPJ(LQP(MP)))=DP3 | |
146 | PJPE(LQP(MP),NPJ(LQP(MP)))=DP | |
147 | PJPM(LQP(MP),NPJ(LQP(MP)))=0.0 | |
148 | GO TO 260 | |
149 | ||
150 | 240 MT=MT+1 | |
151 | DRR=RDT(MT)-R0 | |
152 | IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 210 | |
79a689ca | 153 | DP=DRR*DEDX0 |
e74335a4 | 154 | IF(DP.LE.0.2) GO TO 210 |
155 | IF(PTOT.LE.0.4) GO TO 290 | |
156 | IF(PTOT.LE.DP) DP=PTOT-0.2 | |
157 | DE=DP | |
158 | ||
159 | IF(KFPJ(JP,I).NE.21) THEN | |
160 | PRSHU=PT(LQT(MT),1)**2+PT(LQT(MT),2)**2 | |
161 | & +PT(LQT(MT),3)**2 | |
162 | DE=SQRT(PJPM(JP,I)**2+PTOT**2) | |
163 | & -SQRT(PJPM(JP,I)**2+(PTOT-DP)**2) | |
164 | ERSHU=(PT(LQT(MT),4)+DE-DP)**2 | |
165 | AMSHU=ERSHU-PRSHU | |
166 | IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 210 | |
167 | PT(LQT(MT),4)=SQRT(ERSHU) | |
168 | PT(LQT(MT),5)=SQRT(AMSHU) | |
169 | ENDIF | |
170 | C ********reshuffle the energy when jet has mass | |
171 | ||
172 | R0=RDT(MT) | |
173 | DP1=DP*V1 | |
174 | DP2=DP*V2 | |
175 | DP3=DP*V3 | |
176 | C ********momentum and energy transfer from jet | |
177 | NTJ(LQT(MT))=NTJ(LQT(MT))+1 | |
178 | KFTJ(LQT(MT),NTJ(LQT(MT)))=21 | |
179 | PJTX(LQT(MT),NTJ(LQT(MT)))=DP1 | |
180 | PJTY(LQT(MT),NTJ(LQT(MT)))=DP2 | |
181 | PJTZ(LQT(MT),NTJ(LQT(MT)))=DP3 | |
182 | PJTE(LQT(MT),NTJ(LQT(MT)))=DP | |
183 | PJTM(LQT(MT),NTJ(LQT(MT)))=0.0 | |
184 | ||
185 | 260 PJPX(JP,I)=(PTOT-DP)*V1 | |
186 | PJPY(JP,I)=(PTOT-DP)*V2 | |
187 | PJPZ(JP,I)=(PTOT-DP)*V3 | |
188 | PJPE(JP,I)=PJPE(JP,I)-DE | |
189 | ||
190 | PTOT=PTOT-DP | |
191 | NQ=NQ+1 | |
192 | GO TO 200 | |
193 | 290 CONTINUE | |
194 | ||
195 | RETURN | |
196 | ||
197 | C******************************************************* | |
198 | C Jet interaction for target jet in the direction PHIT | |
199 | C****************************************************** | |
200 | C | |
201 | C******* find the wounded proj which can interact with jet*** | |
202 | ||
203 | 400 IF(NFT(JPJT,7).NE.1) RETURN | |
204 | JT=JPJT | |
205 | DO 690 I=1,NTJ(JT) | |
206 | PTJET0=SQRT(PJTX(JT,I)**2+PJTY(JT,I)**2) | |
207 | IF(PTJET0.LE.HIPR1(11)) GO TO 690 | |
208 | PTOT=SQRT(PTJET0*PTJET0+PJTZ(JT,I)**2) | |
79a689ca | 209 | IF (IHPR2(50) .EQ. 1) THEN |
210 | DEDX0 = HIPR1(14)*LOG10(PTJET0)/LOG10(5.) | |
211 | ELSE | |
212 | DEDX0 = HIPR1(14) | |
213 | ENDIF | |
e74335a4 | 214 | IF(PTOT.LT.HIPR1(8)) GO TO 690 |
215 | PHIT=ULANGL_HIJING(PJTX(JT,I),PJTY(JT,I)) | |
216 | KP=0 | |
217 | DO 500 I2=1,IHNT2(1) | |
218 | IF(NFP(I2,5).NE.3) GO TO 500 | |
219 | DX=YP(1,I2)-YT(1,JT) | |
220 | DY=YP(2,I2)-YT(2,JT) | |
221 | PHI=ULANGL_HIJING(DX,DY) | |
222 | DPHI=ABS(PHI-PHIT) | |
9804fb1d | 223 | IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi |
e74335a4 | 224 | IF(DPHI.GT.HIPR1(40)/2.0) GO TO 500 |
225 | RD0=SQRT(DX*DX+DY*DY) | |
226 | IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 500 | |
227 | KP=KP+1 | |
228 | LQP(KP)=I2 | |
229 | RDP(KP)=COS(DPHI)*RD0 | |
230 | 500 CONTINUE | |
231 | C******* rearrange according to decending rd************ | |
232 | DO 510 I2=1,KP-1 | |
233 | DO 510 J2=I2+1,KP | |
234 | IF(RDP(I2).LT.RDP(J2)) GO TO 510 | |
235 | RD=RDP(I2) | |
236 | LQ=LQP(I2) | |
237 | RDP(I2)=RDP(J2) | |
238 | LQP(I2)=LQP(J2) | |
239 | RDP(J2)=RD | |
240 | LQP(J2)=LQ | |
241 | 510 CONTINUE | |
242 | C****** find wounded targ which can interact with jet******** | |
243 | KT=0 | |
244 | DO 520 I2=1,IHNT2(3) | |
245 | IF(NFT(I2,5).NE.3 .OR. I2.EQ.JT) GO TO 520 | |
246 | DX=YT(1,I2)-YT(1,JT) | |
247 | DY=YT(2,I2)-YT(2,JT) | |
248 | PHI=ULANGL_HIJING(DX,DY) | |
249 | DPHI=ABS(PHI-PHIT) | |
9804fb1d | 250 | IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi |
e74335a4 | 251 | IF(DPHI.GT.HIPR1(40)/2.0) GO TO 520 |
252 | RD0=SQRT(DX*DX+DY*DY) | |
253 | IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 520 | |
254 | KT=KT+1 | |
255 | LQT(KT)=I2 | |
256 | RDT(KT)=COS(DPHI)*RD0 | |
257 | 520 CONTINUE | |
258 | C******* rearrange according to decending rd************ | |
259 | DO 530 I2=1,KT-1 | |
260 | DO 530 J2=I2+1,KT | |
261 | IF(RDT(I2).LT.RDT(J2)) GO TO 530 | |
262 | RD=RDT(I2) | |
263 | LQ=LQT(I2) | |
264 | RDT(I2)=RDT(J2) | |
265 | LQT(I2)=LQT(J2) | |
266 | RDT(J2)=RD | |
267 | LQT(J2)=LQ | |
268 | 530 CONTINUE | |
269 | ||
270 | MP=0 | |
271 | MT=0 | |
272 | NQ=0 | |
273 | DP=0.0 | |
274 | R0=0.0 | |
275 | PTOT=SQRT(PJTX(JT,I)**2+PJTY(JT,I)**2+PJTZ(JT,I)**2) | |
276 | V1=PJTX(JT,I)/PTOT | |
277 | V2=PJTY(JT,I)/PTOT | |
278 | V3=PJTZ(JT,I)/PTOT | |
279 | ||
280 | 600 RN=RLU_HIJING(0) | |
281 | 610 IF(MT.GE.KT .AND. MP.GE.KP) GO TO 690 | |
282 | IF(MT.GE.KT) GO TO 620 | |
283 | IF(MP.GE.KP) GO TO 640 | |
284 | IF(RDP(MP+1).GT.RDT(MT+1)) GO TO 640 | |
285 | 620 MP=MP+1 | |
286 | DRR=RDP(MP)-R0 | |
287 | IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 610 | |
79a689ca | 288 | DP=DRR*DEDX0 |
e74335a4 | 289 | IF(KFTJ(JT,I).NE.21) DP=0.5*DP |
290 | C ********string tension of quark jet is 0.5 of gluon's | |
291 | IF(DP.LE.0.2) GO TO 610 | |
292 | IF(PTOT.LE.0.4) GO TO 690 | |
293 | IF(PTOT.LE.DP) DP=PTOT-0.2 | |
294 | DE=DP | |
295 | C | |
296 | IF(KFTJ(JT,I).NE.21) THEN | |
297 | PRSHU=PP(LQP(MP),1)**2+PP(LQP(MP),2)**2 | |
298 | & +PP(LQP(MP),3)**2 | |
299 | DE=SQRT(PJTM(JT,I)**2+PTOT**2) | |
300 | & -SQRT(PJTM(JT,I)**2+(PTOT-DP)**2) | |
301 | ERSHU=(PP(LQP(MP),4)+DE-DP)**2 | |
302 | AMSHU=ERSHU-PRSHU | |
303 | IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 610 | |
304 | PP(LQP(MP),4)=SQRT(ERSHU) | |
305 | PP(LQP(MP),5)=SQRT(AMSHU) | |
306 | ENDIF | |
307 | C ********reshuffle the energy when jet has mass | |
308 | C | |
309 | R0=RDP(MP) | |
310 | DP1=DP*V1 | |
311 | DP2=DP*V2 | |
312 | DP3=DP*V3 | |
313 | C ********momentum and energy transfer from jet | |
314 | NPJ(LQP(MP))=NPJ(LQP(MP))+1 | |
315 | KFPJ(LQP(MP),NPJ(LQP(MP)))=21 | |
316 | PJPX(LQP(MP),NPJ(LQP(MP)))=DP1 | |
317 | PJPY(LQP(MP),NPJ(LQP(MP)))=DP2 | |
318 | PJPZ(LQP(MP),NPJ(LQP(MP)))=DP3 | |
319 | PJPE(LQP(MP),NPJ(LQP(MP)))=DP | |
320 | PJPM(LQP(MP),NPJ(LQP(MP)))=0.0 | |
321 | ||
322 | GO TO 660 | |
323 | ||
324 | 640 MT=MT+1 | |
325 | DRR=RDT(MT)-R0 | |
326 | IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 610 | |
79a689ca | 327 | DP=DRR*DEDX0 |
e74335a4 | 328 | IF(DP.LE.0.2) GO TO 610 |
329 | IF(PTOT.LE.0.4) GO TO 690 | |
330 | IF(PTOT.LE.DP) DP=PTOT-0.2 | |
331 | DE=DP | |
332 | ||
333 | IF(KFTJ(JT,I).NE.21) THEN | |
334 | PRSHU=PT(LQT(MT),1)**2+PT(LQT(MT),2)**2 | |
335 | & +PT(LQT(MT),3)**2 | |
336 | DE=SQRT(PJTM(JT,I)**2+PTOT**2) | |
337 | & -SQRT(PJTM(JT,I)**2+(PTOT-DP)**2) | |
338 | ERSHU=(PT(LQT(MT),4)+DE-DP)**2 | |
339 | AMSHU=ERSHU-PRSHU | |
340 | IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 610 | |
341 | PT(LQT(MT),4)=SQRT(ERSHU) | |
342 | PT(LQT(MT),5)=SQRT(AMSHU) | |
343 | ENDIF | |
344 | C ********reshuffle the energy when jet has mass | |
345 | ||
346 | R0=RDT(MT) | |
347 | DP1=DP*V1 | |
348 | DP2=DP*V2 | |
349 | DP3=DP*V3 | |
350 | C ********momentum and energy transfer from jet | |
351 | NTJ(LQT(MT))=NTJ(LQT(MT))+1 | |
352 | KFTJ(LQT(MT),NTJ(LQT(MT)))=21 | |
353 | PJTX(LQT(MT),NTJ(LQT(MT)))=DP1 | |
354 | PJTY(LQT(MT),NTJ(LQT(MT)))=DP2 | |
355 | PJTZ(LQT(MT),NTJ(LQT(MT)))=DP3 | |
356 | PJTE(LQT(MT),NTJ(LQT(MT)))=DP | |
357 | PJTM(LQT(MT),NTJ(LQT(MT)))=0.0 | |
358 | ||
359 | 660 PJTX(JT,I)=(PTOT-DP)*V1 | |
360 | PJTY(JT,I)=(PTOT-DP)*V2 | |
361 | PJTZ(JT,I)=(PTOT-DP)*V3 | |
362 | PJTE(JT,I)=PJTE(JT,I)-DE | |
363 | ||
364 | PTOT=PTOT-DP | |
365 | NQ=NQ+1 | |
366 | GO TO 600 | |
367 | 690 CONTINUE | |
368 | RETURN | |
369 | C******************************************************** | |
370 | C Q-QBAR jet interaction | |
371 | C******************************************************** | |
372 | 2000 ISG=JPJT | |
373 | IF(IASG(ISG,3).NE.1) RETURN | |
374 | C | |
375 | JP=IASG(ISG,1) | |
376 | JT=IASG(ISG,2) | |
377 | XJ=(YP(1,JP)+YT(1,JT))/2.0 | |
378 | YJ=(YP(2,JP)+YT(2,JT))/2.0 | |
379 | DO 2690 I=1,NJSG(ISG) | |
380 | PTJET0=SQRT(PXSG(ISG,I)**2+PYSG(ISG,I)**2) | |
381 | IF(PTJET0.LE.HIPR1(11).OR.PESG(ISG,I).LT.HIPR1(1)) | |
382 | & GO TO 2690 | |
383 | PTOT=SQRT(PTJET0*PTJET0+PZSG(ISG,I)**2) | |
79a689ca | 384 | IF (IHPR2(50) .EQ. 1) THEN |
385 | DEDX0 = HIPR1(14)*LOG10(PTJET0)/LOG10(5.) | |
386 | ELSE | |
387 | DEDX0 = HIPR1(14) | |
388 | ENDIF | |
e74335a4 | 389 | IF(PTOT.LT.MAX(HIPR1(1),HIPR1(8))) GO TO 2690 |
390 | PHIQ=ULANGL_HIJING(PXSG(ISG,I),PYSG(ISG,I)) | |
391 | KP=0 | |
392 | DO 2500 I2=1,IHNT2(1) | |
393 | IF(NFP(I2,5).NE.3.OR.I2.EQ.JP) GO TO 2500 | |
394 | DX=YP(1,I2)-XJ | |
395 | DY=YP(2,I2)-YJ | |
396 | PHI=ULANGL_HIJING(DX,DY) | |
397 | DPHI=ABS(PHI-PHIQ) | |
9804fb1d | 398 | IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi |
e74335a4 | 399 | IF(DPHI.GT.HIPR1(40)/2.0) GO TO 2500 |
400 | RD0=SQRT(DX*DX+DY*DY) | |
401 | IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 2500 | |
402 | KP=KP+1 | |
403 | LQP(KP)=I2 | |
404 | RDP(KP)=COS(DPHI)*RD0 | |
405 | 2500 CONTINUE | |
406 | C******* rearrange according to decending rd************ | |
407 | DO 2510 I2=1,KP-1 | |
408 | DO 2510 J2=I2+1,KP | |
409 | IF(RDP(I2).LT.RDP(J2)) GO TO 2510 | |
410 | RD=RDP(I2) | |
411 | LQ=LQP(I2) | |
412 | RDP(I2)=RDP(J2) | |
413 | LQP(I2)=LQP(J2) | |
414 | RDP(J2)=RD | |
415 | LQP(J2)=LQ | |
416 | 2510 CONTINUE | |
417 | C****** find wounded targ which can interact with jet******** | |
418 | KT=0 | |
419 | DO 2520 I2=1,IHNT2(3) | |
420 | IF(NFT(I2,5).NE.3 .OR. I2.EQ.JT) GO TO 2520 | |
421 | DX=YT(1,I2)-XJ | |
422 | DY=YT(2,I2)-YJ | |
423 | PHI=ULANGL_HIJING(DX,DY) | |
424 | DPHI=ABS(PHI-PHIQ) | |
9804fb1d | 425 | IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi |
e74335a4 | 426 | IF(DPHI.GT.HIPR1(40)/2.0) GO TO 2520 |
427 | RD0=SQRT(DX*DX+DY*DY) | |
428 | IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 2520 | |
429 | KT=KT+1 | |
430 | LQT(KT)=I2 | |
431 | RDT(KT)=COS(DPHI)*RD0 | |
432 | 2520 CONTINUE | |
433 | C******* rearrange according to decending rd************ | |
434 | DO 2530 I2=1,KT-1 | |
435 | DO 2530 J2=I2+1,KT | |
436 | IF(RDT(I2).LT.RDT(J2)) GO TO 2530 | |
437 | RD=RDT(I2) | |
438 | LQ=LQT(I2) | |
439 | RDT(I2)=RDT(J2) | |
440 | LQT(I2)=LQT(J2) | |
441 | RDT(J2)=RD | |
442 | LQT(J2)=LQ | |
443 | 2530 CONTINUE | |
444 | ||
445 | MP=0 | |
446 | MT=0 | |
447 | NQ=0 | |
448 | DP=0.0 | |
449 | R0=0.0 | |
450 | PTOT=SQRT(PXSG(ISG,I)**2+PYSG(ISG,I)**2 | |
451 | & +PZSG(ISG,I)**2) | |
452 | V1=PXSG(ISG,I)/PTOT | |
453 | V2=PYSG(ISG,I)/PTOT | |
454 | V3=PZSG(ISG,I)/PTOT | |
455 | ||
456 | 2600 RN=RLU_HIJING(0) | |
457 | 2610 IF(MT.GE.KT .AND. MP.GE.KP) GO TO 2690 | |
458 | IF(MT.GE.KT) GO TO 2620 | |
459 | IF(MP.GE.KP) GO TO 2640 | |
460 | IF(RDP(MP+1).GT.RDT(MT+1)) GO TO 2640 | |
461 | 2620 MP=MP+1 | |
462 | DRR=RDP(MP)-R0 | |
463 | IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 2610 | |
79a689ca | 464 | DP=DRR*DEDX0/2.0 |
e74335a4 | 465 | IF(DP.LE.0.2) GO TO 2610 |
466 | IF(PTOT.LE.0.4) GO TO 2690 | |
467 | IF(PTOT.LE.DP) DP=PTOT-0.2 | |
468 | DE=DP | |
469 | C | |
470 | IF(K2SG(ISG,I).NE.21) THEN | |
471 | IF(PTOT.LT.DP+HIPR1(1)) GO TO 2690 | |
472 | PRSHU=PP(LQP(MP),1)**2+PP(LQP(MP),2)**2 | |
473 | & +PP(LQP(MP),3)**2 | |
474 | DE=SQRT(PMSG(ISG,I)**2+PTOT**2) | |
475 | & -SQRT(PMSG(ISG,I)**2+(PTOT-DP)**2) | |
476 | ERSHU=(PP(LQP(MP),4)+DE-DP)**2 | |
477 | AMSHU=ERSHU-PRSHU | |
478 | IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 2610 | |
479 | PP(LQP(MP),4)=SQRT(ERSHU) | |
480 | PP(LQP(MP),5)=SQRT(AMSHU) | |
481 | ENDIF | |
482 | C ********reshuffle the energy when jet has mass | |
483 | C | |
484 | R0=RDP(MP) | |
485 | DP1=DP*V1 | |
486 | DP2=DP*V2 | |
487 | DP3=DP*V3 | |
488 | C ********momentum and energy transfer from jet | |
489 | NPJ(LQP(MP))=NPJ(LQP(MP))+1 | |
490 | KFPJ(LQP(MP),NPJ(LQP(MP)))=21 | |
491 | PJPX(LQP(MP),NPJ(LQP(MP)))=DP1 | |
492 | PJPY(LQP(MP),NPJ(LQP(MP)))=DP2 | |
493 | PJPZ(LQP(MP),NPJ(LQP(MP)))=DP3 | |
494 | PJPE(LQP(MP),NPJ(LQP(MP)))=DP | |
495 | PJPM(LQP(MP),NPJ(LQP(MP)))=0.0 | |
496 | ||
497 | GO TO 2660 | |
498 | ||
499 | 2640 MT=MT+1 | |
500 | DRR=RDT(MT)-R0 | |
501 | IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 2610 | |
79a689ca | 502 | DP=DRR*DEDX0 |
e74335a4 | 503 | IF(DP.LE.0.2) GO TO 2610 |
504 | IF(PTOT.LE.0.4) GO TO 2690 | |
505 | IF(PTOT.LE.DP) DP=PTOT-0.2 | |
506 | DE=DP | |
507 | ||
508 | IF(K2SG(ISG,I).NE.21) THEN | |
509 | IF(PTOT.LT.DP+HIPR1(1)) GO TO 2690 | |
510 | PRSHU=PT(LQT(MT),1)**2+PT(LQT(MT),2)**2 | |
511 | & +PT(LQT(MT),3)**2 | |
512 | DE=SQRT(PMSG(ISG,I)**2+PTOT**2) | |
513 | & -SQRT(PMSG(ISG,I)**2+(PTOT-DP)**2) | |
514 | ERSHU=(PT(LQT(MT),4)+DE-DP)**2 | |
515 | AMSHU=ERSHU-PRSHU | |
516 | IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 2610 | |
517 | PT(LQT(MT),4)=SQRT(ERSHU) | |
518 | PT(LQT(MT),5)=SQRT(AMSHU) | |
519 | ENDIF | |
520 | C ********reshuffle the energy when jet has mass | |
521 | ||
522 | R0=RDT(MT) | |
523 | DP1=DP*V1 | |
524 | DP2=DP*V2 | |
525 | DP3=DP*V3 | |
526 | C ********momentum and energy transfer from jet | |
527 | NTJ(LQT(MT))=NTJ(LQT(MT))+1 | |
528 | KFTJ(LQT(MT),NTJ(LQT(MT)))=21 | |
529 | PJTX(LQT(MT),NTJ(LQT(MT)))=DP1 | |
530 | PJTY(LQT(MT),NTJ(LQT(MT)))=DP2 | |
531 | PJTZ(LQT(MT),NTJ(LQT(MT)))=DP3 | |
532 | PJTE(LQT(MT),NTJ(LQT(MT)))=DP | |
533 | PJTM(LQT(MT),NTJ(LQT(MT)))=0.0 | |
534 | ||
535 | 2660 PXSG(ISG,I)=(PTOT-DP)*V1 | |
536 | PYSG(ISG,I)=(PTOT-DP)*V2 | |
537 | PZSG(ISG,I)=(PTOT-DP)*V3 | |
538 | PESG(ISG,I)=PESG(ISG,I)-DE | |
539 | ||
540 | PTOT=PTOT-DP | |
541 | NQ=NQ+1 | |
542 | GO TO 2600 | |
543 | 2690 CONTINUE | |
544 | RETURN | |
545 | END | |
79a689ca | 546 | |
547 | ||
548 | ||
549 | ||
550 | ||
551 | ||
552 |