Do not update impact parameter in case of event repetition
authorAndreas Morsch <morsch@pcalice102b.(none)>
Thu, 12 Dec 2013 14:05:24 +0000 (15:05 +0100)
committerAndreas Morsch <morsch@pcalice102b.(none)>
Thu, 12 Dec 2013 14:05:24 +0000 (15:05 +0100)
except for Ncoll==0

HIJING/hijing1_36/hijing.F

index 07db9ce..23c319a 100644 (file)
@@ -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