Changes for #90817 Please commit PHOS trigger part in the AliAnalysisTaskESDfilter.cxx

[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / quench.F

diff --git a/HIJING/hijing1_36/quench.F b/HIJING/hijing1_36/quench.F
index 39b626376b485912facad921900f52e0ea302be8..11b344b0553fb2c9d7a4db86857e8e2933a46f94 100644 (file)
--- a/HIJING/hijing1_36/quench.F
+++ b/HIJING/hijing1_36/quench.F
@@ -1,10 +1,11 @@
-* $Id$
+C $Id$
 C
 C
 C
 C
        SUBROUTINE QUENCH(JPJT,NTP)
        DIMENSION RDP(300),LQP(300),RDT(300),LQT(300)
+#define BLANKET_SAVE
 #include "hijcrdn.inc"
 #include "hiparnt.inc"
 C
@@ -32,6 +33,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***
@@ -65,8 +71,8 @@ C****** find wounded targ which can interact with jet********
              KT=0
              DO 120 I2=1,IHNT2(3)
                 IF(NFT(I2,5).NE.3) GO TO 120
-                DX=YT(1,I2)-YP(1,JP)
-                DY=YT(2,I2)-YP(2,JP)
+                DX=YT(1,I2)-YP(1,JP)-BBX
+                DY=YT(2,I2)-YP(2,JP)-BBY
                 PHI=ULANGL_HIJING(DX,DY)
                 DPHI=ABS(PHI-PHIP)
                 IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi
@@ -107,7 +113,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 +151,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,13 +207,18 @@ 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
                 DO 500 I2=1,IHNT2(1)
                    IF(NFP(I2,5).NE.3) GO TO 500
-                   DX=YP(1,I2)-YT(1,JT)
-                   DY=YP(2,I2)-YT(2,JT)
+                   DX=YP(1,I2)+BBX-YT(1,JT)
+                   DY=YP(2,I2)+BBY-YT(2,JT)
                    PHI=ULANGL_HIJING(DX,DY)
                    DPHI=ABS(PHI-PHIT)
                    IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi
@@ -275,7 +286,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 +325,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
@@ -364,20 +375,25 @@ C********************************************************
 C
        JP=IASG(ISG,1)
        JT=IASG(ISG,2)
-       XJ=(YP(1,JP)+YT(1,JT))/2.0
-       YJ=(YP(2,JP)+YT(2,JT))/2.0
+       XJ=(YP(1,JP)+BBX+YT(1,JT))/2.0
+       YJ=(YP(2,JP)+BBY+YT(2,JT))/2.0
        DO 2690 I=1,NJSG(ISG)
           PTJET0=SQRT(PXSG(ISG,I)**2+PYSG(ISG,I)**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
           DO 2500 I2=1,IHNT2(1)
              IF(NFP(I2,5).NE.3.OR.I2.EQ.JP) GO TO 2500
-             DX=YP(1,I2)-XJ
-             DY=YP(2,I2)-YJ
+             DX=YP(1,I2)+BBX-XJ
+             DY=YP(2,I2)+BBY-YJ
              PHI=ULANGL_HIJING(DX,DY)
              DPHI=ABS(PHI-PHIQ)
               IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI      ! Uzhi
@@ -446,7 +462,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 +500,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 +544,10 @@ C          ********momentum and energy transfer from jet
  2690  CONTINUE
        RETURN
        END
+
+
+
+
+
+
+