Do not update impact parameter in case of event repetition

except for Ncoll==0

diff --git a/HIJING/hijing1_36/hijing.F b/HIJING/hijing1_36/hijing.F
index 07db9ceb025dd87e8454cdf3b87c414cd853e03f..23c319a89d98f90b150d9481093cfb88d30db372 100644 (file)
--- a/HIJING/hijing1_36/hijing.F
+++ b/HIJING/hijing1_36/hijing.F
@@ -296,6 +296,12 @@ c*********************************
 
 C********************
 24     MISS=-1
+C****  BB IS THE ABSOLUTE VALUE OF IMPACT PARAMETER,BB**2 IS 
+C       RANDOMLY GENERATED AND ITS ORIENTATION IS RANDOMLY SET 
+C       BY THE ANGLE PHI  FOR EACH COLLISION.******************
+C
+       BB=SQRT(BMIN**2+RLU_HIJING(0)*(BMAX**2-BMIN**2))
+       PHI=2.0*HIPR1(40)*RLU_HIJING(0)
 
 50     MISS=MISS+1
        IF(MISS.GT.50) THEN
@@ -326,12 +332,6 @@ C                  ********Initialize for a new event
         NPSPECT=0
         NNSPECT=0 
 
-C****  BB IS THE ABSOLUTE VALUE OF IMPACT PARAMETER,BB**2 IS 
-C       RANDOMLY GENERATED AND ITS ORIENTATION IS RANDOMLY SET 
-C       BY THE ANGLE PHI  FOR EACH COLLISION.******************
-C
-       BB=SQRT(BMIN**2+RLU_HIJING(0)*(BMAX**2-BMIN**2))
-       PHI=2.0*HIPR1(40)*RLU_HIJING(0)
        BBX=BB*COS(PHI)
        BBY=BB*SIN(PHI)
        HINT1(19)=BB
@@ -467,8 +467,10 @@ C                          suffered
        IF(NCOLT.EQ.0) THEN
           NLOP=NLOP+1
            IF(NLOP.LE.20.OR.
-     &           (IHNT2(1).EQ.1.AND.IHNT2(3).EQ.1)) GO TO 60
-
+     &          (IHNT2(1).EQ.1.AND.IHNT2(3).EQ.1)) THEN 
+              BB=SQRT(BMIN**2+RLU_HIJING(0)*(BMAX**2-BMIN**2))
+              GO TO 60
+           ENDIF
 
            RETURN
        ENDIF