Correction V.Uzhinskii (v1.37) introduced.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 8 Oct 2003 11:59:21 +0000 (11:59 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 8 Oct 2003 11:59:21 +0000 (11:59 +0000)
HIJING/hijing1_36/quench.F

index d993d01..7d2fdc0 100644 (file)
@@ -70,8 +70,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
@@ -216,8 +216,8 @@ C******* find the wounded proj which can interact with jet***
                 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
@@ -374,8 +374,8 @@ 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))
@@ -391,8 +391,8 @@ C
           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