3 C*********************************************************************
5 SUBROUTINE PYOVLY_HIJING(MOVLY)
7 C...Initializes multiplicity distribution and selects mutliplicity
8 C...of overlayed events, i.e. several events occuring at the same
10 #include "ludat1_hijing.inc"
11 #include "pypars_hijing.inc"
12 #include "pyint1_hijing.inc"
16 C...Sum of allowed cross-sections for overlayed events.
19 IF(MSTP(132).GE.2) VINT(131)=VINT(131)+VINT(104)
20 IF(MSTP(132).GE.3) VINT(131)=VINT(131)+VINT(103)
21 IF(MSTP(132).GE.4) VINT(131)=VINT(131)+VINT(102)
23 C...Initialize multiplicity distribution for unbiased events.
24 IF(MSTP(133).EQ.1) THEN
25 XNAVE=VINT(131)*PARP(131)
26 IF(XNAVE.GT.40.) WRITE(MSTU(11),1000) XNAVE
27 WTI(0)=EXP(-MIN(50.,XNAVE))
31 WTI(I)=WTI(I-1)*XNAVE/I
32 IF(I-2.5.GT.XNAVE.AND.WTI(I).LT.1E-6) GOTO 110
40 C...Initialize mutiplicity distribution for biased events.
41 ELSEIF(MSTP(133).EQ.2) THEN
42 XNAVE=VINT(131)*PARP(131)
43 IF(XNAVE.GT.40.) WRITE(MSTU(11),1000) XNAVE
44 WTI(1)=EXP(-MIN(50.,XNAVE))*XNAVE
48 WTI(I)=WTI(I-1)*XNAVE/(I-1)
49 IF(I-2.5.GT.XNAVE.AND.WTI(I).LT.1E-6) GOTO 130
58 C...Pick multiplicity of overlayed events.
60 IF(MSTP(133).EQ.0) THEN
61 MINT(81)=MAX(1,MSTP(134))
67 IF(WTR.LE.0.) GOTO 150
73 C...Format statement for error message.
74 1000 FORMAT(1X,'Warning: requested average number of events per bunch',
75 &'crossing too large, ',1P,E12.4)