cosmetics

[u/mrichter/AliRoot.git] / DPMJET / dpmjet3.0-5F.f

diff --git a/DPMJET/dpmjet3.0-5F.f b/DPMJET/dpmjet3.0-5F.f
index 9061bf8c9372eeca35c24e07ef0f6f32ffaaeb3d..863d66b1cb5727cc72f4110cf021a7a5f1bab754 100644 (file)
--- a/DPMJET/dpmjet3.0-5F.f
+++ b/DPMJET/dpmjet3.0-5F.f
@@ -3750,8 +3750,8 @@ C #include "dtu_dtevtp.inc"
 
 * Glauber formalism: collision properties
       COMMON /DTGLCP/ RPROJ,RTARG,BIMPAC,
-     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC
-
+     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC,
+     &                NCP,NCT
 * flags for diffractive interactions (DTUNUC 1.x)
       COMMON /DTFLG3/ ISINGD,IDOUBD,IFLAGD,IDIFF
 
@@ -5873,7 +5873,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
 
 * flavors of partons (DTUNUC 1.x)
       COMMON /DTDPMF/ IPVQ(MAXVQU),IPPV1(MAXVQU),IPPV2(MAXVQU),
@@ -6125,7 +6126,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
 **temporary
 
 * statistics: Glauber-formalism
@@ -9203,8 +9205,8 @@ C     PA2(3) = SIGN(SQRT(PZ2NSQ),PA2(3))
 
 * Glauber formalism: collision properties
       COMMON /DTGLCP/ RPROJ,RTARG,BIMPAC,
-     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC
-
+     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC,
+     &                NCP,NCT
       DIMENSION R(3),PIN(4),POUT(4),DEV(4)
 
       DO 1 K=1,4
@@ -9314,7 +9316,8 @@ C              PHKK(K,I) = PHKK(K,I)-DEV(K)/DBLE(NVAL)
 
 * Glauber formalism: collision properties
       COMMON /DTGLCP/ RPROJ,RTARG,BIMPAC,
-     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC
+     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC,
+     &                NCP,NCT
 
 * various options for treatment of partons (DTUNUC 1.x)
 * (chain recombination, Cronin,..)
@@ -11255,8 +11258,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
 * Glauber formalism: flags and parameters for statistics
       LOGICAL LPROD
       CHARACTER*8 CGLB
@@ -16911,7 +16914,8 @@ C           QARJ(I) = PARJ(I)
 
 * Glauber formalism: collision properties
       COMMON /DTGLCP/ RPROJ,RTARG,BIMPAC,
-     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC
+     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC,
+     &                NCP,NCT
 
 * flags for input different options
       LOGICAL LEMCCK,LHADRO,LSEADI,LEVAPO
@@ -17110,8 +17114,8 @@ C9990 CONTINUE
 
 * Glauber formalism: collision properties
       COMMON /DTGLCP/ RPROJ,RTARG,BIMPAC,
-     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC
-
+     &                NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC,
+     &                NCP,NCT
 * nuclear potential
       LOGICAL LFERMI
       COMMON /DTNPOT/ PFERMP(2),PFERMN(2),FERMOD,
@@ -18854,37 +18858,37 @@ C     EPOT(2,8) = PFERMN(2)**2/(2.0D0*AAM(8)) + EBINDN(2)
 
 * correction of projectile 4-momentum for effective target pot.
 * and Coulomb-energy (in case of hadron-nucleus interaction only)
-      IF ((IP.EQ.1).AND.(IT.GT.1).AND.LFERMI) THEN
-         EPNI = EPN
+*      IF ((IP.EQ.1).AND.(IT.GT.1).AND.LFERMI) THEN
+*         EPNI = EPN
 *   Coulomb-energy:
 *     positively charged hadron - check energy for Coloumb pot.
-         IF (IICH(IJPROJ).EQ.1) THEN
-            THRESH = ETACOU(2)+AAM(IJPROJ)
-            IF (EPNI.LE.THRESH) THEN
-               WRITE(LOUT,1000)
- 1000          FORMAT(/,1X,'KKINC:  WARNING!  projectile energy',
-     &                ' below Coulomb threshold - event rejected',/)
-               ISTHKK(1) = 1
-               RETURN
-            ENDIF
+*         IF (IICH(IJPROJ).EQ.1) THEN
+*            THRESH = ETACOU(2)+AAM(IJPROJ)
+*            IF (EPNI.LE.THRESH) THEN
+*               WRITE(LOUT,1000)
+* 1000          FORMAT(/,1X,'KKINC:  WARNING!  projectile energy',
+*     &                ' below Coulomb threshold - event rejected',/)
+*               ISTHKK(1) = 1
+*               RETURN
+*            ENDIF
 *     negatively charged hadron - increase energy by Coulomb energy
-         ELSEIF (IICH(IJPROJ).EQ.-1) THEN
-            EPNI = EPNI+ETACOU(2)
-         ENDIF
-         IF ((IJPROJ.EQ.1).OR.(IJPROJ.EQ.8)) THEN
+*         ELSEIF (IICH(IJPROJ).EQ.-1) THEN
+*            EPNI = EPNI+ETACOU(2)
+*         ENDIF
+*         IF ((IJPROJ.EQ.1).OR.(IJPROJ.EQ.8)) THEN
 *   Effective target potential
 *sr 6.6. binding energy only (to avoid negative exc. energies)
 C           EPNI = EPNI+EPOT(2,IJPROJ)
-            EBIPOT = EBINDP(2)
-            IF ((IJPROJ.NE.1).AND.(ABS(EPOT(2,IJPROJ)).GT.5.0D-3))
-     &         EBIPOT = EBINDN(2)
-            EPNI = EPNI+ABS(EBIPOT)
+*            EBIPOT = EBINDP(2)
+*            IF ((IJPROJ.NE.1).AND.(ABS(EPOT(2,IJPROJ)).GT.5.0D-3))
+*     &         EBIPOT = EBINDN(2)
+*            EPNI = EPNI+ABS(EBIPOT)
 * re-initialization of DTLTRA
-            DUM1 = ZERO
-            DUM2 = ZERO
-            CALL DT_LTINI(IJPROJ,IJTARG,EPNI,DUM1,DUM2,0)
-         ENDIF
-      ENDIF
+*            DUM1 = ZERO
+*            DUM2 = ZERO
+*            CALL DT_LTINI(IJPROJ,IJTARG,EPNI,DUM1,DUM2,0)
+*         ENDIF
+*      ENDIF
 
 * projectile in n-n cms
       IF ((IP.LE.1).AND.(IT.GT.1)) THEN