8 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
10 common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)
12 common/ludat3/mdcy(500,3),mdme(2000,2),brat(2000),kfdp(2000,5)
14 common/pyhisubs/msel,msub(200),kfin(2,-40:40),ckin(200)
20 common/pyhiint4/widp(21:40,0:40),wide(21:40,0:40),wids(21:40,3)
22 common/pyhiint5/
ngen(0:200,3),
xsec(0:200,3)
24 common/pyhiint6/proc(0:200)
27 CHARACTER chau*16,chpa(-40:40)*12,chin(2)*12,
28 &state(-1:5)*4,chkin(21)*18
29 DATA state/
'----',
'off ',
'on ',
'on/+',
'on/-',
'on/1',
'on/2'/,
30 &chkin/
' m_hard (GeV/c^2) ',
' p_T_hard (GeV/c) ',
31 &
'm_finite (GeV/c^2)',
' y*_subsystem ',
' y*_large ',
32 &
' y*_small ',
' eta*_large ',
' eta*_small ',
33 &
'cos(theta*)_large ',
'cos(theta*)_small ',
' x_1 ',
34 &
' x_2 ',
' x_F ',
' cos(theta_hard) ',
35 &
'm''_hard (GeV/c^2) ',
' tau ',
' y* ',
36 &
'cos(theta_hard^-) ',
'cos(theta_hard^+) ',
' x_T^2 ',
43 WRITE(mstu(11),1200) 0,proc(0),
ngen(0,3),
ngen(0,1),
xsec(0,3)
45 IF(msub(
i).NE.1) goto 100
48 WRITE(mstu(11),1300) 1.-float(
ngen(0,3))/
52 ELSEIF(mstat.EQ.2)
THEN
55 110 chpa(kf)=chau(1:12)
63 WRITE(mstu(11),1600) chpa(kc),0.,0.,state(mdcy(kc,1)),0.
66 120
IF(mdme(idc,2).EQ.102)
WRITE(mstu(11),1700) chpa(kfdp(idc,1)),
67 & chpa(kfdp(idc,2)),0.,0.,state(mdme(idc,1)),0.
75 IF(wide(kc,0).GT.0.)
THEN
76 WRITE(mstu(11),1600) chpa(kc),widp(kc,0),1.,
77 & state(mdcy(kc,1)),1.
80 140
WRITE(mstu(11),1700) chpa(kfdp(idc,1)),chpa(kfdp(idc,2)),
81 & widp(kc,
j),widp(kc,
j)/widp(kc,0),state(mdme(idc,1)),
82 & wide(kc,
j)/wide(kc,0)
84 WRITE(mstu(11),1600) chpa(kc),widp(kc,0),1.,
85 & state(mdcy(kc,1)),0.
91 ELSEIF(mstat.EQ.3)
THEN
97 WRITE(mstu(11),2000) chin(1),chin(2)
100 160 chpa(kf)=chau(1:12)
101 IF(
mint(43).EQ.1)
THEN
102 WRITE(mstu(11),2100) chpa(
mint(11)),state(kfin(1,
mint(11))),
103 & chpa(
mint(12)),state(kfin(2,
mint(12)))
104 ELSEIF(
mint(43).EQ.2)
THEN
105 WRITE(mstu(11),2100) chpa(
mint(11)),state(kfin(1,
mint(11))),
106 & chpa(-
mstp(54)),state(kfin(2,-
mstp(54)))
107 DO 170
i=-
mstp(54)+1,-1
108 170
WRITE(mstu(11),2200) chpa(
i),state(kfin(2,
i))
110 180
WRITE(mstu(11),2200) chpa(
i),state(kfin(2,
i))
111 WRITE(mstu(11),2200) chpa(21),state(kfin(2,21))
112 ELSEIF(
mint(43).EQ.3)
THEN
113 WRITE(mstu(11),2100) chpa(-
mstp(54)),state(kfin(1,-
mstp(54))),
114 & chpa(
mint(12)),state(kfin(2,
mint(12)))
115 DO 190
i=-
mstp(54)+1,-1
116 190
WRITE(mstu(11),2300) chpa(
i),state(kfin(1,
i))
118 200
WRITE(mstu(11),2300) chpa(
i),state(kfin(1,
i))
119 WRITE(mstu(11),2300) chpa(21),state(kfin(1,21))
120 ELSEIF(
mint(43).EQ.4)
THEN
121 DO 210
i=-
mstp(54),-1
122 210
WRITE(mstu(11),2100) chpa(
i),state(kfin(1,
i)),chpa(
i),
125 220
WRITE(mstu(11),2100) chpa(
i),state(kfin(1,
i)),chpa(
i),
127 WRITE(mstu(11),2100) chpa(21),state(kfin(1,21)),chpa(21),
133 ELSEIF(mstat.EQ.4)
THEN
137 IF(shrmax.LT.0.) shrmax=
vint(1)
138 WRITE(mstu(11),2700) ckin(1),chkin(1),shrmax
139 pthmin=
max(ckin(3),ckin(5))
141 IF(pthmax.LT.0.) pthmax=0.5*shrmax
142 WRITE(mstu(11),2800) ckin(3),pthmin,chkin(2),pthmax
143 WRITE(mstu(11),2900) chkin(3),ckin(6)
145 230
WRITE(mstu(11),2700) ckin(2*
i-1),chkin(
i),ckin(2*
i)
147 IF(sprmax.LT.0.) sprmax=
vint(1)
148 WRITE(mstu(11),2700) ckin(31),chkin(13),sprmax
153 240
WRITE(mstu(11),2700)
vint(
i-5),chkin(
i),
vint(
i+15)
157 ELSEIF(mstat.EQ.5)
THEN
166 1000
FORMAT(
'1',9(
'*'),1
x,
'PYHISTAT: Statistics on Number of ',
167 &
'Events and Cross-sections',1
x,9(
'*'))
168 1100
FORMAT(/1
x,78(
'=')/1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,
'I',12
x,
169 &
'Subprocess',12
x,
'I',6
x,
'Number of points',6
x,
'I',4
x,
'Sigma',3
x,
170 &
'I'/1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,
'I',34(
'-'),
'I',28(
'-'),
171 &
'I',4
x,
'(mb)',4
x,
'I'/1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,
'I',1
x,
172 &
'N:o',1
x,
'Type',25
x,
'I',4
x,
'Generated',9
x,
'Tried',1
x,
'I',12
x,
173 &
'I'/1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,78(
'=')/1
x,
'I',34
x,
'I',28
x,
175 1200
FORMAT(1
x,
'I',1
x,
i3,1
x,a28,1
x,
'I',1
x,i12,1
x,i13,1
x,
'I',1
x,1
p,
177 1300
FORMAT(1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,78(
'=')//
178 &1
x,
'********* Fraction of events that fail fragmentation ',
179 &
'cuts =',1
x,f8.5,
' *********'/)
180 1400
FORMAT(
'1',17(
'*'),1
x,
'PYHISTAT: Decay Widths and Branching ',
181 &
'Ratios',1
x,17(
'*'))
182 1500
FORMAT(/1
x,78(
'=')/1
x,
'I',29
x,
'I',13
x,
'I',12
x,
'I',6
x,
'I',12
x,
'I'/
183 &1
x,
'I',1
x,
'Branching/Decay Channel',5
x,
'I',1
x,
'Width (GeV)',1
x,
184 &
'I',7
x,
'B.R.',1
x,
'I',1
x,
'Stat',1
x,
'I',2
x,
'Eff. B.R.',1
x,
'I'/1
x,
185 &
'I',29
x,
'I',13
x,
'I',12
x,
'I',6
x,
'I',12
x,
'I'/1
x,78(
'='))
186 1600
FORMAT(1
x,
'I',29
x,
'I',13
x,
'I',12
x,
'I',6
x,
'I',12
x,
'I'/1
x,
'I',1
x,
187 &a12,1
x,
'->',13
x,
'I',2
x,1
p,e10.3,0
p,1
x,
'I',1
x,1
p,e10.3,0
p,1
x,
'I',
188 &1
x,a4,1
x,
'I',1
x,1
p,e10.3,0
p,1
x,
'I')
189 1700
FORMAT(1
x,
'I',1
x,a12,1
x,
'+',1
x,a12,1
x,
'I',2
x,1
p,e10.3,0
p,1
x,
'I',
190 &1
x,1
p,e10.3,0
p,1
x,
'I',1
x,a4,1
x,
'I',1
x,1
p,e10.3,0
p,1
x,
'I')
191 1800
FORMAT(1
x,
'I',29
x,
'I',13
x,
'I',12
x,
'I',6
x,
'I',12
x,
'I'/1
x,78(
'='))
192 1900
FORMAT(
'1',7(
'*'),1
x,
'PYHISTAT: Allowed Incoming Partons/',
193 &
'Particles at Hard Interaction',1
x,7(
'*'))
194 2000
FORMAT(/1
x,78(
'=')/1
x,
'I',38
x,
'I',37
x,
'I'/1
x,
'I',1
x,
195 &
'Beam particle:',1
x,a,10
x,
'I',1
x,
'Target particle:',1
x,a,7
x,
196 &
'I'/1
x,
'I',38
x,
'I',37
x,
'I'/1
x,
'I',1
x,
'Content',9
x,
'State',16
x,
197 &
'I',1
x,
'Content',9
x,
'State',15
x,
'I'/1
x,
'I',38
x,
'I',37
x,
'I'/1
x,
198 &78(
'=')/1
x,
'I',38
x,
'I',37
x,
'I')
199 2100
FORMAT(1
x,
'I',1
x,a,5
x,a,16
x,
'I',1
x,a,5
x,a,15
x,
'I')
200 2200
FORMAT(1
x,
'I',38
x,
'I',1
x,a,5
x,a,15
x,
'I')
201 2300
FORMAT(1
x,
'I',1
x,a,5
x,a,16
x,
'I',37
x,
'I')
202 2400
FORMAT(1
x,
'I',38
x,
'I',37
x,
'I'/1
x,78(
'='))
203 2500
FORMAT(
'1',12(
'*'),1
x,
'PYHISTAT: User-Defined Limits on ',
204 &
'Kinematical Variables',1
x,12(
'*'))
205 2600
FORMAT(/1
x,78(
'=')/1
x,
'I',76
x,
'I')
206 2700
FORMAT(1
x,
'I',16
x,1
p,e10.3,0
p,1
x,
'<',1
x,a,1
x,
'<',1
x,1
p,e10.3,0
p,
208 2800
FORMAT(1
x,
'I',3
x,1
p,e10.3,0
p,1
x,
'(',1
p,e10.3,0
p,
')',1
x,
'<',1
x,a,
209 &1
x,
'<',1
x,1
p,e10.3,0
p,16
x,
'I')
210 2900
FORMAT(1
x,
'I',29
x,a,1
x,
'=',1
x,1
p,e10.3,0
p,16
x,
'I')
211 3000
FORMAT(1
x,
'I',76
x,
'I'/1
x,78(
'='))
212 3100
FORMAT(////1
x,5(
'*'),1
x,
213 &
'PYHISTAT: Derived Limits on Kinematical ',
214 &
'Variables Used in Generation',1
x,5(
'*'))
215 3200
FORMAT(
'1',12(
'*'),1
x,
'PYHISTAT: Summary of Status Codes and ',
216 &
'Parameter Values',1
x,12(
'*'))
217 3300
FORMAT(/3
x,
'I',4
x,
'MSTP(I)',9
x,
'PARP(I)',20
x,
'I',4
x,
'MSTP(I)',9
x,
219 3400
FORMAT(1
x,
i3,5
x,i6,6
x,1
p,e10.3,0
p,18
x,
i3,5
x,i6,6
x,1
p,e10.3)