Corrected number of collisions

Chiara

diff --git a/DPMJET/dpmjet3.0-5F.f b/DPMJET/dpmjet3.0-5F.f
index 388b0c79b24132678074b2fa31de7df4d041b2ce..2c4897c43c6581324556bba96bdf6de9af9cecec 100644 (file)
--- a/DPMJET/dpmjet3.0-5F.f
+++ b/DPMJET/dpmjet3.0-5F.f
@@ -5532,7 +5532,8 @@ C           ENDIF
 
 * Glauber formalism: collision properties
       COMMON /DTGLCP/ RPROJ,RTARG,BIMPAC,
-     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC
+     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC,
+     &                NCP,NCT
 
 * central particle production, impact parameter biasing
       COMMON /DTIMPA/ BIMIN,BIMAX,XSFRAC,ICENTR
@@ -5547,6 +5548,8 @@ C           ENDIF
       IREJ   = 0
       ICREQU = ICREQU+1
       NC     = 0
+      NCP    = 0
+      NCT    = 0
 
     1 CONTINUE
       ICSAMP = ICSAMP+1
@@ -5586,6 +5589,14 @@ C           ENDIF
          ITOLD  = IT
          JJPOLD = JJPROJ
          EPROLD = EPROJ
+     DO 8 I=1, IP
+        NCP = NCP+JSSH(I)
+*        WRITE(6,*)' PROJ.NUCL. ',I,' NCOLL = ',NCP
+    8 CONTINUE
+     DO 9 I=1, IT
+        NCT = NCT+JTSH(I)
+*        WRITE(6,*)' TAR.NUCL. ',I,' NCOLL = ',NCT
+    9 CONTINUE
       ENDIF
 
 * force diffractive particle production in h-K interactions