Include corrections by V.Uzhinskii JINR, Dubna, Russia from V1.37
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 13 Feb 2002 15:43:20 +0000 (15:43 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 13 Feb 2002 15:43:20 +0000 (15:43 +0000)
HIJING/hijing1_36/quench.F

index 675ac48..39b6263 100644 (file)
@@ -13,6 +13,12 @@ C
 #include "histrng.inc"
 C
        SAVE
+C
+        BB=HINT1(19)                                            ! Uzhi
+        PHI=HINT1(20)                                           ! Uzhi
+        BBX=BB*COS(PHI)                                         ! Uzhi
+        BBY=BB*SIN(PHI)                                         ! Uzhi
+c
        IF(NTP.EQ.2) GO TO 400
        IF(NTP.EQ.3) GO TO 2000 
 C*******************************************************
@@ -36,6 +42,7 @@ C******* find the wounded proj which can interact with jet***
              DY=YP(2,I2)-YP(2,JP)
              PHI=ULANGL_HIJING(DX,DY)
              DPHI=ABS(PHI-PHIP)
+              IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI      ! Uzhi
              IF(DPHI.GE.HIPR1(40)/2.0) GO TO 100
              RD0=SQRT(DX*DX+DY*DY)
              IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 100
@@ -62,6 +69,7 @@ C****** find wounded targ which can interact with jet********
                 DY=YT(2,I2)-YP(2,JP)
                 PHI=ULANGL_HIJING(DX,DY)
                 DPHI=ABS(PHI-PHIP)
+                IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi
                 IF(DPHI.GT.HIPR1(40)/2.0) GO TO 120
                 RD0=SQRT(DX*DX+DY*DY)
                 IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 120
@@ -202,6 +210,7 @@ C******* find the wounded proj which can interact with jet***
                    DY=YP(2,I2)-YT(2,JT)
                    PHI=ULANGL_HIJING(DX,DY)
                    DPHI=ABS(PHI-PHIT)
+                   IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi
                    IF(DPHI.GT.HIPR1(40)/2.0) GO TO 500
                    RD0=SQRT(DX*DX+DY*DY)
                    IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 500
@@ -228,6 +237,7 @@ C****** find wounded targ which can interact with jet********
                       DY=YT(2,I2)-YT(2,JT)
                       PHI=ULANGL_HIJING(DX,DY)
                       DPHI=ABS(PHI-PHIT)
+                      IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi
                       IF(DPHI.GT.HIPR1(40)/2.0) GO TO 520
                       RD0=SQRT(DX*DX+DY*DY)
                       IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 520
@@ -370,6 +380,7 @@ C
              DY=YP(2,I2)-YJ
              PHI=ULANGL_HIJING(DX,DY)
              DPHI=ABS(PHI-PHIQ)
+              IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI      ! Uzhi
              IF(DPHI.GT.HIPR1(40)/2.0) GO TO 2500
              RD0=SQRT(DX*DX+DY*DY)
              IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 2500
@@ -396,6 +407,7 @@ C****** find wounded targ which can interact with jet********
                 DY=YT(2,I2)-YJ
                 PHI=ULANGL_HIJING(DX,DY)
                 DPHI=ABS(PHI-PHIQ)
+                IF(DPHI.GE.HIPR1(40)) DPHI=2.*HIPR1(40)-DPHI ! Uzhi
                 IF(DPHI.GT.HIPR1(40)/2.0) GO TO 2520
                 RD0=SQRT(DX*DX+DY*DY)
                 IF(RD0*SIN(DPHI).GT.HIPR1(12)) GO TO 2520