Substracting the x-projection of impact parameter
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 20 Feb 2004 16:29:28 +0000 (16:29 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 20 Feb 2004 16:29:28 +0000 (16:29 +0000)
HIJING/hijing1_36/hijsft.F

index a9c0e5d..54fa8eb 100644 (file)
@@ -113,8 +113,8 @@ C
              IF(ABS(COS(PHI)).LT.1.0E-5) THEN
 C AM            DD1=R1
 C AM            DD2=R1
-                 DD1 = ABS(SQRT(R1**2-XT0**2) - YT0)
-                 DD2 = ABS(SQRT(R1**2-XT0**2) + YT0)
+                 DD1 = ABS(SQRT(R1**2-(XT0-BX)**2) - YT0)
+                 DD2 = ABS(SQRT(R1**2-(XT0-BX)**2) + YT0)
 C PH            DD3=ABS(BY+SQRT(R2**2-(XP0-BX)**2)-YP0)
 C PH            DD4=ABS(BY-SQRT(R2**2-(XP0-BX)**2)-YP0)
                 DD3=ABS(BY+SQRT(R2**2-(XT0-BX)**2)-YP0)