6 SUBROUTINE QUENCH(JPJT,NTP)
7 DIMENSION RDP(300),LQP(300),RDT(300),LQT(300)
11 #include "hijjet1.inc"
12 #include "hijjet2.inc"
13 #include "histrng.inc"
16 IF(NTP.EQ.2) GO TO 400
17 IF(NTP.EQ.3) GO TO 2000
18 C*******************************************************
19 C Jet interaction for proj jet in the direction PHIP
20 C******************************************************
22 IF(NFP(JPJT,7).NE.1) RETURN
26 PTJET0=SQRT(PJPX(JP,I)**2+PJPY(JP,I)**2)
27 IF(PTJET0.LE.HIPR1(11)) GO TO 290
28 PTOT=SQRT(PTJET0*PTJET0+PJPZ(JP,I)**2)
29 IF(PTOT.LT.HIPR1(8)) GO TO 290
30 PHIP=ULANGL_HIJING(PJPX(JP,I),PJPY(JP,I))
31 C******* find the wounded proj which can interact with jet***
34 IF(NFP(I2,5).NE.3 .OR. I2.EQ.JP) GO TO 100
37 PHI=ULANGL_HIJING(DX,DY)
39 IF(DPHI.GE.HIPR1(40)/2.0) GO TO 100
41 IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 100
46 C******* rearrange according decending rd************
49 IF(RDP(I2).LT.RDP(J2)) GO TO 110
57 C****** find wounded targ which can interact with jet********
60 IF(NFT(I2,5).NE.3) GO TO 120
63 PHI=ULANGL_HIJING(DX,DY)
65 IF(DPHI.GT.HIPR1(40)/2.0) GO TO 120
67 IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 120
72 C******* rearrange according decending rd************
75 IF(RDT(I2).LT.RDT(J2)) GO TO 130
89 PTOT=SQRT(PJPX(JP,I)**2+PJPY(JP,I)**2+PJPZ(JP,I)**2)
95 210 IF(MT.GE.KT .AND. MP.GE.KP) GO TO 290
96 IF(MT.GE.KT) GO TO 220
97 IF(MP.GE.KP) GO TO 240
98 IF(RDP(MP+1).GT.RDT(MT+1)) GO TO 240
101 IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 210
103 IF(KFPJ(JP,I).NE.21) DP=0.5*DP
104 C ********string tension of quark jet is 0.5 of gluon's
105 IF(DP.LE.0.2) GO TO 210
106 IF(PTOT.LE.0.4) GO TO 290
107 IF(PTOT.LE.DP) DP=PTOT-0.2
110 IF(KFPJ(JP,I).NE.21) THEN
111 PRSHU=PP(LQP(MP),1)**2+PP(LQP(MP),2)**2
113 DE=SQRT(PJPM(JP,I)**2+PTOT**2)
114 & -SQRT(PJPM(JP,I)**2+(PTOT-DP)**2)
115 ERSHU=(PP(LQP(MP),4)+DE-DP)**2
117 IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 210
118 PP(LQP(MP),4)=SQRT(ERSHU)
119 PP(LQP(MP),5)=SQRT(AMSHU)
121 C ********reshuffle the energy when jet has mass
126 C ********momentum and energy transfer from jet
128 NPJ(LQP(MP))=NPJ(LQP(MP))+1
129 KFPJ(LQP(MP),NPJ(LQP(MP)))=21
130 PJPX(LQP(MP),NPJ(LQP(MP)))=DP1
131 PJPY(LQP(MP),NPJ(LQP(MP)))=DP2
132 PJPZ(LQP(MP),NPJ(LQP(MP)))=DP3
133 PJPE(LQP(MP),NPJ(LQP(MP)))=DP
134 PJPM(LQP(MP),NPJ(LQP(MP)))=0.0
139 IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 210
141 IF(DP.LE.0.2) GO TO 210
142 IF(PTOT.LE.0.4) GO TO 290
143 IF(PTOT.LE.DP) DP=PTOT-0.2
146 IF(KFPJ(JP,I).NE.21) THEN
147 PRSHU=PT(LQT(MT),1)**2+PT(LQT(MT),2)**2
149 DE=SQRT(PJPM(JP,I)**2+PTOT**2)
150 & -SQRT(PJPM(JP,I)**2+(PTOT-DP)**2)
151 ERSHU=(PT(LQT(MT),4)+DE-DP)**2
153 IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 210
154 PT(LQT(MT),4)=SQRT(ERSHU)
155 PT(LQT(MT),5)=SQRT(AMSHU)
157 C ********reshuffle the energy when jet has mass
163 C ********momentum and energy transfer from jet
164 NTJ(LQT(MT))=NTJ(LQT(MT))+1
165 KFTJ(LQT(MT),NTJ(LQT(MT)))=21
166 PJTX(LQT(MT),NTJ(LQT(MT)))=DP1
167 PJTY(LQT(MT),NTJ(LQT(MT)))=DP2
168 PJTZ(LQT(MT),NTJ(LQT(MT)))=DP3
169 PJTE(LQT(MT),NTJ(LQT(MT)))=DP
170 PJTM(LQT(MT),NTJ(LQT(MT)))=0.0
172 260 PJPX(JP,I)=(PTOT-DP)*V1
173 PJPY(JP,I)=(PTOT-DP)*V2
174 PJPZ(JP,I)=(PTOT-DP)*V3
175 PJPE(JP,I)=PJPE(JP,I)-DE
184 C*******************************************************
185 C Jet interaction for target jet in the direction PHIT
186 C******************************************************
188 C******* find the wounded proj which can interact with jet***
190 400 IF(NFT(JPJT,7).NE.1) RETURN
193 PTJET0=SQRT(PJTX(JT,I)**2+PJTY(JT,I)**2)
194 IF(PTJET0.LE.HIPR1(11)) GO TO 690
195 PTOT=SQRT(PTJET0*PTJET0+PJTZ(JT,I)**2)
196 IF(PTOT.LT.HIPR1(8)) GO TO 690
197 PHIT=ULANGL_HIJING(PJTX(JT,I),PJTY(JT,I))
200 IF(NFP(I2,5).NE.3) GO TO 500
203 PHI=ULANGL_HIJING(DX,DY)
205 IF(DPHI.GT.HIPR1(40)/2.0) GO TO 500
206 RD0=SQRT(DX*DX+DY*DY)
207 IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 500
210 RDP(KP)=COS(DPHI)*RD0
212 C******* rearrange according to decending rd************
215 IF(RDP(I2).LT.RDP(J2)) GO TO 510
223 C****** find wounded targ which can interact with jet********
226 IF(NFT(I2,5).NE.3 .OR. I2.EQ.JT) GO TO 520
229 PHI=ULANGL_HIJING(DX,DY)
231 IF(DPHI.GT.HIPR1(40)/2.0) GO TO 520
232 RD0=SQRT(DX*DX+DY*DY)
233 IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 520
236 RDT(KT)=COS(DPHI)*RD0
238 C******* rearrange according to decending rd************
241 IF(RDT(I2).LT.RDT(J2)) GO TO 530
255 PTOT=SQRT(PJTX(JT,I)**2+PJTY(JT,I)**2+PJTZ(JT,I)**2)
261 610 IF(MT.GE.KT .AND. MP.GE.KP) GO TO 690
262 IF(MT.GE.KT) GO TO 620
263 IF(MP.GE.KP) GO TO 640
264 IF(RDP(MP+1).GT.RDT(MT+1)) GO TO 640
267 IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 610
269 IF(KFTJ(JT,I).NE.21) DP=0.5*DP
270 C ********string tension of quark jet is 0.5 of gluon's
271 IF(DP.LE.0.2) GO TO 610
272 IF(PTOT.LE.0.4) GO TO 690
273 IF(PTOT.LE.DP) DP=PTOT-0.2
276 IF(KFTJ(JT,I).NE.21) THEN
277 PRSHU=PP(LQP(MP),1)**2+PP(LQP(MP),2)**2
279 DE=SQRT(PJTM(JT,I)**2+PTOT**2)
280 & -SQRT(PJTM(JT,I)**2+(PTOT-DP)**2)
281 ERSHU=(PP(LQP(MP),4)+DE-DP)**2
283 IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 610
284 PP(LQP(MP),4)=SQRT(ERSHU)
285 PP(LQP(MP),5)=SQRT(AMSHU)
287 C ********reshuffle the energy when jet has mass
293 C ********momentum and energy transfer from jet
294 NPJ(LQP(MP))=NPJ(LQP(MP))+1
295 KFPJ(LQP(MP),NPJ(LQP(MP)))=21
296 PJPX(LQP(MP),NPJ(LQP(MP)))=DP1
297 PJPY(LQP(MP),NPJ(LQP(MP)))=DP2
298 PJPZ(LQP(MP),NPJ(LQP(MP)))=DP3
299 PJPE(LQP(MP),NPJ(LQP(MP)))=DP
300 PJPM(LQP(MP),NPJ(LQP(MP)))=0.0
306 IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 610
308 IF(DP.LE.0.2) GO TO 610
309 IF(PTOT.LE.0.4) GO TO 690
310 IF(PTOT.LE.DP) DP=PTOT-0.2
313 IF(KFTJ(JT,I).NE.21) THEN
314 PRSHU=PT(LQT(MT),1)**2+PT(LQT(MT),2)**2
316 DE=SQRT(PJTM(JT,I)**2+PTOT**2)
317 & -SQRT(PJTM(JT,I)**2+(PTOT-DP)**2)
318 ERSHU=(PT(LQT(MT),4)+DE-DP)**2
320 IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 610
321 PT(LQT(MT),4)=SQRT(ERSHU)
322 PT(LQT(MT),5)=SQRT(AMSHU)
324 C ********reshuffle the energy when jet has mass
330 C ********momentum and energy transfer from jet
331 NTJ(LQT(MT))=NTJ(LQT(MT))+1
332 KFTJ(LQT(MT),NTJ(LQT(MT)))=21
333 PJTX(LQT(MT),NTJ(LQT(MT)))=DP1
334 PJTY(LQT(MT),NTJ(LQT(MT)))=DP2
335 PJTZ(LQT(MT),NTJ(LQT(MT)))=DP3
336 PJTE(LQT(MT),NTJ(LQT(MT)))=DP
337 PJTM(LQT(MT),NTJ(LQT(MT)))=0.0
339 660 PJTX(JT,I)=(PTOT-DP)*V1
340 PJTY(JT,I)=(PTOT-DP)*V2
341 PJTZ(JT,I)=(PTOT-DP)*V3
342 PJTE(JT,I)=PJTE(JT,I)-DE
349 C********************************************************
350 C Q-QBAR jet interaction
351 C********************************************************
353 IF(IASG(ISG,3).NE.1) RETURN
357 XJ=(YP(1,JP)+YT(1,JT))/2.0
358 YJ=(YP(2,JP)+YT(2,JT))/2.0
359 DO 2690 I=1,NJSG(ISG)
360 PTJET0=SQRT(PXSG(ISG,I)**2+PYSG(ISG,I)**2)
361 IF(PTJET0.LE.HIPR1(11).OR.PESG(ISG,I).LT.HIPR1(1))
363 PTOT=SQRT(PTJET0*PTJET0+PZSG(ISG,I)**2)
364 IF(PTOT.LT.MAX(HIPR1(1),HIPR1(8))) GO TO 2690
365 PHIQ=ULANGL_HIJING(PXSG(ISG,I),PYSG(ISG,I))
367 DO 2500 I2=1,IHNT2(1)
368 IF(NFP(I2,5).NE.3.OR.I2.EQ.JP) GO TO 2500
371 PHI=ULANGL_HIJING(DX,DY)
373 IF(DPHI.GT.HIPR1(40)/2.0) GO TO 2500
374 RD0=SQRT(DX*DX+DY*DY)
375 IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 2500
378 RDP(KP)=COS(DPHI)*RD0
380 C******* rearrange according to decending rd************
383 IF(RDP(I2).LT.RDP(J2)) GO TO 2510
391 C****** find wounded targ which can interact with jet********
393 DO 2520 I2=1,IHNT2(3)
394 IF(NFT(I2,5).NE.3 .OR. I2.EQ.JT) GO TO 2520
397 PHI=ULANGL_HIJING(DX,DY)
399 IF(DPHI.GT.HIPR1(40)/2.0) GO TO 2520
400 RD0=SQRT(DX*DX+DY*DY)
401 IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 2520
404 RDT(KT)=COS(DPHI)*RD0
406 C******* rearrange according to decending rd************
409 IF(RDT(I2).LT.RDT(J2)) GO TO 2530
423 PTOT=SQRT(PXSG(ISG,I)**2+PYSG(ISG,I)**2
429 2600 RN=RLU_HIJING(0)
430 2610 IF(MT.GE.KT .AND. MP.GE.KP) GO TO 2690
431 IF(MT.GE.KT) GO TO 2620
432 IF(MP.GE.KP) GO TO 2640
433 IF(RDP(MP+1).GT.RDT(MT+1)) GO TO 2640
436 IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 2610
438 IF(DP.LE.0.2) GO TO 2610
439 IF(PTOT.LE.0.4) GO TO 2690
440 IF(PTOT.LE.DP) DP=PTOT-0.2
443 IF(K2SG(ISG,I).NE.21) THEN
444 IF(PTOT.LT.DP+HIPR1(1)) GO TO 2690
445 PRSHU=PP(LQP(MP),1)**2+PP(LQP(MP),2)**2
447 DE=SQRT(PMSG(ISG,I)**2+PTOT**2)
448 & -SQRT(PMSG(ISG,I)**2+(PTOT-DP)**2)
449 ERSHU=(PP(LQP(MP),4)+DE-DP)**2
451 IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 2610
452 PP(LQP(MP),4)=SQRT(ERSHU)
453 PP(LQP(MP),5)=SQRT(AMSHU)
455 C ********reshuffle the energy when jet has mass
461 C ********momentum and energy transfer from jet
462 NPJ(LQP(MP))=NPJ(LQP(MP))+1
463 KFPJ(LQP(MP),NPJ(LQP(MP)))=21
464 PJPX(LQP(MP),NPJ(LQP(MP)))=DP1
465 PJPY(LQP(MP),NPJ(LQP(MP)))=DP2
466 PJPZ(LQP(MP),NPJ(LQP(MP)))=DP3
467 PJPE(LQP(MP),NPJ(LQP(MP)))=DP
468 PJPM(LQP(MP),NPJ(LQP(MP)))=0.0
474 IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 2610
476 IF(DP.LE.0.2) GO TO 2610
477 IF(PTOT.LE.0.4) GO TO 2690
478 IF(PTOT.LE.DP) DP=PTOT-0.2
481 IF(K2SG(ISG,I).NE.21) THEN
482 IF(PTOT.LT.DP+HIPR1(1)) GO TO 2690
483 PRSHU=PT(LQT(MT),1)**2+PT(LQT(MT),2)**2
485 DE=SQRT(PMSG(ISG,I)**2+PTOT**2)
486 & -SQRT(PMSG(ISG,I)**2+(PTOT-DP)**2)
487 ERSHU=(PT(LQT(MT),4)+DE-DP)**2
489 IF(AMSHU.LT.HIPR1(1)*HIPR1(1)) GO TO 2610
490 PT(LQT(MT),4)=SQRT(ERSHU)
491 PT(LQT(MT),5)=SQRT(AMSHU)
493 C ********reshuffle the energy when jet has mass
499 C ********momentum and energy transfer from jet
500 NTJ(LQT(MT))=NTJ(LQT(MT))+1
501 KFTJ(LQT(MT),NTJ(LQT(MT)))=21
502 PJTX(LQT(MT),NTJ(LQT(MT)))=DP1
503 PJTY(LQT(MT),NTJ(LQT(MT)))=DP2
504 PJTZ(LQT(MT),NTJ(LQT(MT)))=DP3
505 PJTE(LQT(MT),NTJ(LQT(MT)))=DP
506 PJTM(LQT(MT),NTJ(LQT(MT)))=0.0
508 2660 PXSG(ISG,I)=(PTOT-DP)*V1
509 PYSG(ISG,I)=(PTOT-DP)*V2
510 PZSG(ISG,I)=(PTOT-DP)*V3
511 PESG(ISG,I)=PESG(ISG,I)-DE