-* $Id$
+C $Id$
C
C
C
PTJET0=SQRT(PJPX(JP,I)**2+PJPY(JP,I)**2)
IF(PTJET0.LE.HIPR1(11)) GO TO 290
PTOT=SQRT(PTJET0*PTJET0+PJPZ(JP,I)**2)
+ IF (IHPR2(50) .EQ. 1) THEN
+ DEDX0 = HIPR1(14)*LOG10(PTJET0)/LOG10(5.)
+ ELSE
+ DEDX0 = HIPR1(14)
+ ENDIF
IF(PTOT.LT.HIPR1(8)) GO TO 290
PHIP=ULANGL_HIJING(PJPX(JP,I),PJPY(JP,I))
C******* find the wounded proj which can interact with jet***
220 MP=MP+1
DRR=RDP(MP)-R0
IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 210
- DP=DRR*HIPR1(14)
+ DP=DRR*DEDX0
IF(KFPJ(JP,I).NE.21) DP=0.5*DP
C ********string tension of quark jet is 0.5 of gluon's
IF(DP.LE.0.2) GO TO 210
240 MT=MT+1
DRR=RDT(MT)-R0
IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 210
- DP=DRR*HIPR1(14)
+ DP=DRR*DEDX0
IF(DP.LE.0.2) GO TO 210
IF(PTOT.LE.0.4) GO TO 290
IF(PTOT.LE.DP) DP=PTOT-0.2
PTJET0=SQRT(PJTX(JT,I)**2+PJTY(JT,I)**2)
IF(PTJET0.LE.HIPR1(11)) GO TO 690
PTOT=SQRT(PTJET0*PTJET0+PJTZ(JT,I)**2)
+ IF (IHPR2(50) .EQ. 1) THEN
+ DEDX0 = HIPR1(14)*LOG10(PTJET0)/LOG10(5.)
+ ELSE
+ DEDX0 = HIPR1(14)
+ ENDIF
IF(PTOT.LT.HIPR1(8)) GO TO 690
PHIT=ULANGL_HIJING(PJTX(JT,I),PJTY(JT,I))
KP=0
620 MP=MP+1
DRR=RDP(MP)-R0
IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 610
- DP=DRR*HIPR1(14)
+ DP=DRR*DEDX0
IF(KFTJ(JT,I).NE.21) DP=0.5*DP
C ********string tension of quark jet is 0.5 of gluon's
IF(DP.LE.0.2) GO TO 610
640 MT=MT+1
DRR=RDT(MT)-R0
IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 610
- DP=DRR*HIPR1(14)
+ DP=DRR*DEDX0
IF(DP.LE.0.2) GO TO 610
IF(PTOT.LE.0.4) GO TO 690
IF(PTOT.LE.DP) DP=PTOT-0.2
IF(PTJET0.LE.HIPR1(11).OR.PESG(ISG,I).LT.HIPR1(1))
& GO TO 2690
PTOT=SQRT(PTJET0*PTJET0+PZSG(ISG,I)**2)
+ IF (IHPR2(50) .EQ. 1) THEN
+ DEDX0 = HIPR1(14)*LOG10(PTJET0)/LOG10(5.)
+ ELSE
+ DEDX0 = HIPR1(14)
+ ENDIF
IF(PTOT.LT.MAX(HIPR1(1),HIPR1(8))) GO TO 2690
PHIQ=ULANGL_HIJING(PXSG(ISG,I),PYSG(ISG,I))
KP=0
2620 MP=MP+1
DRR=RDP(MP)-R0
IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 2610
- DP=DRR*HIPR1(14)/2.0
+ DP=DRR*DEDX0/2.0
IF(DP.LE.0.2) GO TO 2610
IF(PTOT.LE.0.4) GO TO 2690
IF(PTOT.LE.DP) DP=PTOT-0.2
2640 MT=MT+1
DRR=RDT(MT)-R0
IF(RN.GE.1.0-EXP(-DRR/HIPR1(13))) GO TO 2610
- DP=DRR*HIPR1(14)
+ DP=DRR*DEDX0
IF(DP.LE.0.2) GO TO 2610
IF(PTOT.LE.0.4) GO TO 2690
IF(PTOT.LE.DP) DP=PTOT-0.2
2690 CONTINUE
RETURN
END
+
+
+
+
+
+
+