Corrected number of collisions
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 12 Feb 2013 07:36:52 +0000 (07:36 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 12 Feb 2013 07:36:52 +0000 (07:36 +0000)
Chiara

DPMJET/dpmjet3.0-5F.f

index 388b0c7..2c4897c 100644 (file)
@@ -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