COMMON /DTHIS3/ IHISPP(50),IHISXS(50),IXSTBL
COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
- 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)
* LEPTO
**LUND single / double precision
& IDBAM(NMXHKK),IDCH(NMXHKK),NPOINT(10),
& IHIST(2,NMXHKK)
- 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)
COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
PARAMETER (MAXLND=4000)
COMMON/PYJETS/N,NPAD,K(MAXLND,5),P(MAXLND,5),V(MAXLND,5)
COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
- 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)
* flags for particle decays
COMMON /DTFRPA/ MSTUX(20),PARUX(20),MSTJX(20),PARJX(20),
WRITE(LOUT,1000)
1000 FORMAT(/,/,1X,'STATIS:',20X,'statistics of the run',/,
& 28X,'---------------------')
+ IF (ICREQU.GT.0) THEN
WRITE(LOUT,1001) ICREQU,ICSAMP,DBLE(ICSAMP)/DBLE(ICREQU)
1001 FORMAT(/,1X,'number of events requested / sampled',13X,
& I8,' / ',I8,/,1X,'number of samp. evts per requested ',
& 'event',11X,F9.1)
+ ENDIF
IF (ICDIFF(1).NE.0) THEN
WRITE(LOUT,1009) ICDIFF
1009 FORMAT(/,1X,'diffractive events: total ',I8,/,49X,
& 'low mass high mass',/,24X,'single diffraction',
& 7X,I8,4X,I8,/,24X,'double diffraction',7X,I8,4X,I8)
ENDIF
- IF (ICENTR.GT.0) THEN
+ IF (ICENTR.GT.0.AND.ICSAMP.GT.0.AND.ICCPRO.GT.0) THEN
WRITE(LOUT,1002) DBLE(ICCPRO)/DBLE(ICSAMP),
& DBLE(ICSAMP)/DBLE(ICCPRO)
1002 FORMAT(/,1X,'central production:',/,2X,'mean number',
& ' of sampled Glauber-events per event',9X,F9.1,/,
& 2X,'fraction of production cross section',21X,F10.6)
ENDIF
+ IF (ICSAMP.GT.0) THEN
WRITE(LOUT,1003) DBLE(ICDPR)/DBLE(ICSAMP),
& DBLE(ICDTA)/DBLE(ICSAMP)
1003 FORMAT(/,54X,'proj. targ.',/,1X,'average number of wounded',
& ' nucleons after x-sampling',2(4X,F6.2))
+ ENDIF
IF (MCGENE.EQ.1) THEN
+ IF (ICSAMP.GT.0) THEN
WRITE(LOUT,1004) DBLE(ICRJSS)/DBLE(ICSAMP)
1004 FORMAT(/,1X,'mean number of sea-sea chain rejections per',
& ' event',3X,F9.1)
1005 FORMAT(/,1X,'Reggeon contribution:',/,1X,'mean number ',
& 'of single chains per event',13X,F9.1)
ENDIF
+ ENDIF
+ IF (ICSAMP.GT.0.AND.ICREQU.GT.0) THEN
WRITE(LOUT,1006)
1006 FORMAT(/,1X,'chain system statistics: (per event)',/,
& 23X,'mean number of chains mean number of chains',/,
& 'Treatment of final nucleon conf.',10X,'IREXCI(1)/NEVT = '
& ,F7.2,/,43X,'IREXCI(2)/NEVT = ',F7.2,/,48X,
& 'IREXCI(3) = ',I5,/)
+ ENDIF
ELSEIF (MCGENE.EQ.2) THEN
WRITE(LOUT,1010) ELOJET
1010 FORMAT(/,/,1X,'PHOJET-treatment of chain systems above ',
1013 FORMAT(/,1X,'2. chain system statistics -',
& ' mean numbers per evt:',/,30X,'---------------------',
& /,/,16X,'s-s',7X,'d-s',7X,'s-d')
+ IF (ICSAMP.GT.0) THEN
WRITE(LOUT,1014)
& ((DBLE(ICEVTG(I,J))/DBLE(ICSAMP),I=1,3),J=0,1),
& (DBLE(ICCHAI(2,I))/(2.0D0*DBLE(ICSAMP)),I=1,3),
& ' dbl-dir ',3E10.2,/,' s-Pom. ',3E10.2,/,
& ' h-Pom. ',3E10.2,/,' s-Reg. ',3E10.2,/,
& ' enh-trg ',3E10.2,/,' enh-log ',3E10.2)
+ ENDIF
WRITE(LOUT,1015)
1015 FORMAT(/,16X,'s-v',7X,'d-v',7X,'v-s',7X,'v-d',7X,'v-v')
+ IF (ICSAMP.GT.0) THEN
WRITE(LOUT,1016)
& ((DBLE(ICEVTG(I,J))/DBLE(ICSAMP),I=4,8),J=0,1),
& (DBLE(ICCHAI(2,I))/(2.0D0*DBLE(ICSAMP)),I=4,8),
& ' dbl-dir ',5E10.2,/,' s-Pom. ',5E10.2,/,
& ' h-Pom. ',5E10.2,/,' s-Reg. ',5E10.2,/,
& ' enh-trg ',5E10.2,/,' enh-log ',5E10.2)
+ ENDIF
ENDIF
CALL DT_CHASTA(1)
*
DIMENSION XX(61),YY(61),LL(101,10)
DIMENSION X(N),Y(N),L(40),LOV(40,10)
+ INTEGER*4 LL, L, LOV
DATA L/
11H*,1H2,1H3,1H4,1H5,1H6,1H7,1H8,1H9,1HZ,
21H+,1HA,1HO,1HB,1HC,1HD,1HE,1HF,1HG,1HH,