I discovered a nasty bug in CORRFW/AliCFGridSparse.cxx:

[u/mrichter/AliRoot.git] / DPMJET / phojet1.12-35c3.f

diff --git a/DPMJET/phojet1.12-35c3.f b/DPMJET/phojet1.12-35c3.f
index c3bd70f62d7fdbe3dc8c2aaa69ec82eb53c5768f..8e793debb6f5c67991f5e2501701059436051443 100644 (file)
--- a/DPMJET/phojet1.12-35c3.f
+++ b/DPMJET/phojet1.12-35c3.f
@@ -494,7 +494,7 @@ C  hard cross sections and MC selection weights
       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
       INTEGER MDCY,MDME,KFDP
       DOUBLE PRECISION  BRAT
-      COMMON/PYDAT3/MDCY(500,3),MDME(4000,2),BRAT(4000),KFDP(4000,5)
+      COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
 
       INTEGER PYCOMP
 
@@ -2496,7 +2496,7 @@ C  standard particle data interface
       DOUBLE PRECISION PHEP,VHEP
       COMMON /POEVT1/ NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
      &                JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),
-     &                VHEP(4,NMXHEP)
+     &                VHEP(4,NMXHEP), NSD1, NSD2, NDD
 C  extension to standard particle data interface (PHOJET specific)
       INTEGER IMPART,IPHIST,ICOLOR
       COMMON /POEVT2/ IMPART(NMXHEP),IPHIST(2,NMXHEP),ICOLOR(2,NMXHEP)
@@ -11836,6 +11836,11 @@ C  event debugging information
      &        KHLOO,KSDPO,KHDPO,KEVENT,KSOFT,KHARD
       COMMON /PODEBG/ IDEB(NMAXD),KSPOM,KHPOM,KSREG,KHDIR,KACCEP,KSTRG,
      &                KHTRG,KSLOO,KHLOO,KSDPO,KHDPO,KEVENT,KSOFT,KHARD
+      PARAMETER (NMXHEP=4000)
+      COMMON /POEVT1/ NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
+     &                JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),
+     &                VHEP(4,NMXHEP), NSD1, NSD2, NDD
+
 C  model switches and parameters
       CHARACTER*8 MDLNA
       INTEGER ISWMDL,IPAMDL
@@ -12204,10 +12209,18 @@ C  sample strings to prepare fragmentation
           CALL PHO_PREVNT(-1)
           RETURN
         ENDIF
-        IF(IPROC.EQ.5) ID1A = ID1A+1
-        IF(IPROC.EQ.6) ID2A = ID2A+1
-        IF(IPROC.EQ.7) ID3A = ID3A+1
-
+        IF(IPROC.EQ.5) THEN 
+           ID1A = ID1A+1
+           NSD1 = NSD1 +1
+           ENDIF
+        IF(IPROC.EQ.6) THEN
+           ID2A = ID2A+1
+           NSD2 = NSD2 + 1
+        ENDIF
+        IF(IPROC.EQ.7) THEN
+           ID3A = ID3A+1
+           NDD = NDD + 1
+        ENDIF
 C-----------------------------------------------------------------------
 C  single / double direct processes
 
@@ -12260,9 +12273,18 @@ C  sample strings to prepare fragmentation
           CALL PHO_PREVNT(-1)
           RETURN
         ENDIF
-        IF(IPROC.EQ.5) ID1A = ID1A+1
-        IF(IPROC.EQ.6) ID2A = ID2A+1
-        IF(IPROC.EQ.7) ID3A = ID3A+1
+        IF(IPROC.EQ.5) THEN 
+           ID1A = ID1A+1
+           NSD1 = NSD1 +1
+        ENDIF
+        IF(IPROC.EQ.6) THEN
+           ID2A = ID2A+1
+           NSD2 = NSD2 + 1
+        ENDIF
+        IF(IPROC.EQ.7) THEN
+           ID3A = ID3A+1
+           NDD = NDD + 1
+        ENDIF
         IDIA = IDIA+MHDIR
 
 C-----------------------------------------------------------------------
@@ -33870,7 +33892,7 @@ C  input/output channels
       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
       INTEGER MDCY,MDME,KFDP
       DOUBLE PRECISION  BRAT
-      COMMON/PYDAT3/MDCY(500,3),MDME(4000,2),BRAT(4000),KFDP(4000,5)
+      COMMON/PYDAT3/MDCY(500,3),MDME(8000,2),BRAT(8000),KFDP(8000,5)
 
       INTEGER PYCOMP
 
@@ -33889,7 +33911,7 @@ C  defaults
       DEF21 = PARJ(21)
 
 C  declare stable particles
-      IF(IDEFAB.GE.2) MSTJ(22) = 2
+c     IF(IDEFAB.GE.2) MSTJ(22) = 2
 
 C  load optimized parameters
       IF(IDEFAB.GE.3) THEN