Logarithmic dependence of energy loss added as option.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 19 Jul 2002 14:29:43 +0000 (14:29 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 19 Jul 2002 14:29:43 +0000 (14:29 +0000)
HIJING/hijing1_36/hijsft.F

index 779e93b..120a992 100644 (file)
@@ -138,7 +138,12 @@ C
              IF(DD1.LT.HIPR1(13)) DD1=0.0
              IF(DD2.LT.HIPR1(13)) DD2=0.0
              IF(NFP(JP,10).EQ.1.AND.PPJET.GT.HIPR1(11)) THEN
-                DP1=DD1*HIPR1(14)/2.0
+                IF (IHPR2(50) .EQ. 1) THEN 
+                   DEDX0 = HIPR1(14)*LOG10(PPJET)/LOG10(5.)
+                ELSE
+                   DEDX0 = HIPR1(14)
+                ENDIF
+                DP1=DD1*DEDX0/2.0
                 DP1=MIN(DP1,PPJET-HIPR1(11))
                 PKC1=PPJET-DP1
                 DPX1=COS(PHI1)*DP1
@@ -168,7 +173,12 @@ C
                 ENDIF
              ENDIF
  15          IF(NFT(JT,10).EQ.1.AND.PTJET.GT.HIPR1(11)) THEN
-                DP2=DD2*HIPR1(14)/2.0
+                IF (IHPR2(50) .EQ. 1) THEN 
+                   DEDX0 = HIPR1(14)*LOG10(PTJET)/LOG10(5.)
+                ELSE
+                   DEDX0 = HIPR1(14)
+                ENDIF
+                DP2=DD2*DEDX0/2.0
                 DP2=MIN(DP2,PTJET-HIPR1(11))
                 PKC2=PTJET-DP2
                 DPX2=COS(PHI2)*DP2