Use log(E) dependence for energy loss if IHPR2(50) == 1.

diff --git a/HIJING/hijing1_36/quench.F b/HIJING/hijing1_36/quench.F
index 39b6263..d993d01 100644 (file)
--- a/HIJING/hijing1_36/quench.F
+++ b/HIJING/hijing1_36/quench.F
@@ -1,4 +1,4 @@
-* $Id$
+C $Id$
 C
 C
 C
@@ -32,6 +32,11 @@ 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***
@@ -107,7 +112,7 @@ C*******    rearrange according decending rd************
  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
@@ -145,7 +150,7 @@ C           ********momentum and energy transfer from jet
  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
@@ -201,6 +206,11 @@ C******* find the wounded proj which can interact with jet***
                 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
@@ -275,7 +285,7 @@ C*******    rearrange according to decending rd************
 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
@@ -314,7 +324,7 @@ C           ********momentum and energy transfer from jet
 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
@@ -371,6 +381,11 @@ C
           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
@@ -446,7 +461,7 @@ C*******    rearrange according to decending rd************
  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
@@ -484,7 +499,7 @@ C           ********momentum and energy transfer from jet
  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
@@ -528,3 +543,10 @@ C          ********momentum and energy transfer from jet
  2690  CONTINUE
        RETURN
        END
+
+
+
+
+
+
+