12 IMPLICIT DOUBLE PRECISION(a-
h, o-
z)
16 parameter(ksusy1=1000000,ksusy2=2000000,ktechn=3000000,
17 &kexcit=4000000,kdimen=5000000)
20 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
21 common/
pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
22 common/
pydat3/mdcy(500,3),mdme(8000,2),brat(8000),kfdp(8000,5)
23 common/
pysubs/msel,mselpd,msub(500),kfin(2,-40:40),ckin(200)
26 common/
pyint2/iset(500),kfpr(500,2),coef(500,20),icol(40,4,2)
27 common/
pyint4/mwid(500),wids(500,5)
30 CHARACTER proc*28, chtmp*16
31 common/
pymssm/imss(0:99),rmss(0:99)
32 common/pymsrv/rvlam(3,3,3), rvlamp(3,3,3), rvlamb(3,3,3)
36 dimension wdtp(0:400),wdte(0:400,0:5),nmodes(0:20),pbrat(10)
37 CHARACTER proga(6)*28,chau*16,chkf*16,chd1*16,chd2*16,chd3*16,
38 &chin(2)*12,state(-1:5)*4,chkin(21)*18,disga(2)*28,
39 &progg9(13)*28,progg4(4)*28,progg2(2)*28,progp4(4)*28
40 CHARACTER*24 chd0, chdc(10)
43 &
'VMD/hadron * VMD ',
'VMD/hadron * direct ',
44 &
'VMD/hadron * anomalous ',
'direct * direct ',
45 &
'direct * anomalous ',
'anomalous * anomalous '/
46 DATA disga/
'e * VMD',
'e * anomalous'/
48 &
'direct * direct ',
'direct * VMD ',
49 &
'direct * anomalous ',
'VMD * direct ',
50 &
'VMD * VMD ',
'VMD * anomalous ',
51 &
'anomalous * direct ',
'anomalous * VMD ',
52 &
'anomalous * anomalous ',
'DIS * VMD ',
53 &
'DIS * anomalous ',
'VMD * DIS ',
56 &
'direct * direct ',
'direct * resolved ',
57 &
'resolved * direct ',
'resolved * resolved '/
59 &
'direct * hadron ',
'resolved * hadron '/
61 &
'VMD * hadron ',
'direct * hadron ',
62 &
'anomalous * hadron ',
'DIS * hadron '/
63 DATA state/
'----',
'off ',
'on ',
'on/+',
'on/-',
'on/1',
'on/2'/,
64 &chkin/
' m_hard (GeV/c^2) ',
' p_T_hard (GeV/c) ',
65 &
'm_finite (GeV/c^2)',
' y*_subsystem ',
' y*_large ',
66 &
' y*_small ',
' eta*_large ',
' eta*_small ',
67 &
'cos(theta*)_large ',
'cos(theta*)_small ',
' x_1 ',
68 &
' x_2 ',
' x_F ',
' cos(theta_hard) ',
69 &
'm''_hard (GeV/c^2) ',
' tau ',
' y* ',
70 &
'cos(theta_hard^-) ',
'cos(theta_hard^+) ',
' x_T^2 ',
72 DATA dname /
'q ',
'lepton',
'nu '/
79 WRITE(mstu(11),5200) 0,proc(0),
ngen(0,3),
ngen(0,1),
xsec(0,3)
81 IF(msub(
i).NE.1) goto 100
84 IF(
mint(121).GT.1)
THEN
86 DO 110 iga=1,
mint(121)
88 IF(
mint(121).EQ.2.AND.
mstp(14).EQ.10)
THEN
89 WRITE(mstu(11),5200) iga,disga(iga),
ngen(0,3),
ngen(0,1),
91 ELSEIF(
mint(121).EQ.9.OR.
mint(121).EQ.13)
THEN
92 WRITE(mstu(11),5200) iga,progg9(iga),
ngen(0,3),
ngen(0,1),
94 ELSEIF(
mint(121).EQ.4.AND.
mstp(14).EQ.30)
THEN
95 WRITE(mstu(11),5200) iga,progp4(iga),
ngen(0,3),
ngen(0,1),
97 ELSEIF(
mint(121).EQ.4)
THEN
98 WRITE(mstu(11),5200) iga,progg4(iga),
ngen(0,3),
ngen(0,1),
100 ELSEIF(
mint(121).EQ.2)
THEN
101 WRITE(mstu(11),5200) iga,progg2(iga),
ngen(0,3),
ngen(0,1),
104 WRITE(mstu(11),5200) iga,proga(iga),
ngen(0,3),
ngen(0,1),
110 WRITE(mstu(11),5400) mstu(23),mstu(30),mstu(27),
114 ELSEIF(mstat.EQ.2)
THEN
122 IF(kc.GT.2*
mstp(1).AND.kc.LE.10) goto 140
123 IF(kc.GT.10+2*
mstp(1).AND.kc.LE.20) goto 140
124 IF(kc.LE.5.OR.(kc.GE.11.AND.kc.LE.16)) ioff=1
125 IF(kc.EQ.18.AND.pmas(18,1).LT.1d0) ioff=1
126 IF(kc.EQ.21.OR.kc.EQ.22) ioff=1
128 IF(mwid(kc).LE.0) goto 140
129 IF(imss(1).LE.0.AND.(kf/ksusy1.EQ.1.OR.
130 & kf/ksusy1.EQ.2)) goto 140
135 IF(kc.LE.20) ngp=(mod(kc,10)+1)/2
136 IF(ngp.LE.
mstp(1))
WRITE(mstu(11),5700) kf,chkf(1:10),
137 & pmas(kc,1),0d0,0d0,state(mdcy(kc,1)),0d0
138 DO 120
j=1,mdcy(kc,3)
141 IF(iabs(kfdp(idc,1)).LE.20) ngp1=
142 & (mod(iabs(kfdp(idc,1)),10)+1)/2
144 IF(iabs(kfdp(idc,2)).LE.20) ngp2=
145 & (mod(iabs(kfdp(idc,2)),10)+1)/2
146 CALL
pyname(kfdp(idc,1),chd1)
147 CALL
pyname(kfdp(idc,2),chd2)
148 IF(kfdp(idc,3).EQ.0)
THEN
149 IF(mdme(idc,2).EQ.102.AND.ngp1.LE.
mstp(1).AND.
150 & ngp2.LE.
mstp(1))
WRITE(mstu(11),5800) idc,chd1(1:10),
151 & chd2(1:10),0d0,0d0,state(mdme(idc,1)),0d0
153 CALL
pyname(kfdp(idc,3),chd3)
154 IF(mdme(idc,2).EQ.102.AND.ngp1.LE.
mstp(1).AND.
155 & ngp2.LE.
mstp(1))
WRITE(mstu(11),5900) idc,chd1(1:10),
156 & chd2(1:10),chd3(1:10),0d0,0d0,state(mdme(idc,1)),0d0
161 CALL
pywidt(kf,pmas(kc,1)**2,wdtp,wdte)
163 IF(wdte(0,0).LE.0d0) brfin=0d0
164 WRITE(mstu(11),5700) kf,chkf(1:10),pmas(kc,1),wdtp(0),1d0,
165 & state(mdcy(kc,1)),brfin
166 DO 130
j=1,mdcy(kc,3)
169 IF(iabs(kfdp(idc,1)).LE.20) ngp1=
170 & (mod(iabs(kfdp(idc,1)),10)+1)/2
172 IF(iabs(kfdp(idc,2)).LE.20) ngp2=
173 & (mod(iabs(kfdp(idc,2)),10)+1)/2
175 IF(wdtp(0).GT.0d0) brpri=wdtp(
j)/wdtp(0)
177 IF(wdte(0,0).GT.0d0) brfin=wdte(
j,0)/wdte(0,0)
178 CALL
pyname(kfdp(idc,1),chd1)
179 CALL
pyname(kfdp(idc,2),chd2)
180 IF(kfdp(idc,3).EQ.0)
THEN
181 IF(ngp1.LE.
mstp(1).AND.ngp2.LE.
mstp(1))
182 &
WRITE(mstu(11),5800) idc,chd1(1:10),
183 & chd2(1:10),wdtp(
j),brpri,
184 & state(mdme(idc,1)),brfin
186 CALL
pyname(kfdp(idc,3),chd3)
187 IF(ngp1.LE.
mstp(1).AND.ngp2.LE.
mstp(1))
188 &
WRITE(mstu(11),5900) idc,chd1(1:10),
189 & chd2(1:10),chd3(1:10),wdtp(
j),brpri,
190 & state(mdme(idc,1)),brfin
198 ELSEIF(mstat.EQ.3)
THEN
204 WRITE(mstu(11),6200) chin(1),chin(2)
208 IF(ia.GT.
mstp(58).AND.ia.LE.10) goto 150
209 IF(ia.GT.10+2*
mstp(1).AND.ia.LE.20) goto 150
211 WRITE(mstu(11),6300) chau,state(kfin(1,
i)),chau,
217 ELSEIF(mstat.EQ.4)
THEN
221 IF(shrmax.LT.0d0) shrmax=
vint(1)
222 WRITE(mstu(11),6700) ckin(1),chkin(1),shrmax
223 pthmin=
max(ckin(3),ckin(5))
225 IF(pthmax.LT.0d0) pthmax=0.5d0*shrmax
226 WRITE(mstu(11),6800) ckin(3),pthmin,chkin(2),pthmax
227 WRITE(mstu(11),6900) chkin(3),ckin(6)
229 WRITE(mstu(11),6700) ckin(2*
i-1),chkin(
i),ckin(2*
i)
232 IF(sprmax.LT.0d0) sprmax=
vint(1)
233 WRITE(mstu(11),6700) ckin(31),chkin(15),sprmax
237 ELSEIF(mstat.EQ.5)
THEN
246 ELSEIF(mstat.EQ.6)
THEN
250 IF(iset(
i).LT.0) goto 180
251 WRITE(mstu(11),7600)
i,proc(
i),iset(
i),kfpr(
i,1),kfpr(
i,2)
255 ELSEIF(mstat.EQ.7)
THEN
256 WRITE (mstu(11),8000)
262 kfsusy=ilr*ksusy1+kfsm
265 IF (kfsm.EQ.1.OR.kfsm.EQ.3.OR.kfsm.EQ.5)
THEN
273 chdc(1)=dname(3) //
' + ' // dname(1)
274 chdc(2)=dname(2) //
' + ' // dname(1)
275 chdc(3)=dname(1) //
' + ' // dname(1)
277 DO 200
j=1,mdcy(kc,3)
279 id1=iabs(kfdp(idc,1))
280 id2=iabs(kfdp(idc,2))
281 IF (kfdp(idc,3).EQ.0)
THEN
282 IF ((id1.EQ.12.OR.id1.EQ.14.OR.id1.EQ.16).AND.(id2
283 & .EQ.1.OR.id2.EQ.3.OR.id2.EQ.5))
THEN
284 pbrat(1)=pbrat(1)+brat(idc)
285 nmodes(1)=nmodes(1)+1
286 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
287 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
288 ELSE IF ((id1.EQ.11.OR.id1.EQ.13.OR.id1.EQ.15).and
289 & .(id2.EQ.2.OR.id2.EQ.4.OR.id2.EQ.6))
THEN
290 pbrat(2)=pbrat(2)+brat(idc)
291 nmodes(2)=nmodes(2)+1
292 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
293 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
294 ELSE IF ((id1.EQ.2.OR.id1.EQ.4.OR.id1.EQ.6).and
295 & .(id2.EQ.1.OR.id2.EQ.3.OR.id2.EQ.5))
THEN
296 pbrat(3)=pbrat(3)+brat(idc)
297 nmodes(3)=nmodes(3)+1
298 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
299 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
305 IF (kfsm.EQ.2.OR.kfsm.EQ.4.OR.kfsm.EQ.6)
THEN
313 chdc(1)=dname(2) //
' + ' // dname(1)
314 chdc(2)=dname(1) //
' + ' // dname(1)
316 DO 220
j=1,mdcy(kc,3)
318 id1=iabs(kfdp(idc,1))
319 id2=iabs(kfdp(idc,2))
320 IF (kfdp(idc,3).EQ.0)
THEN
321 IF ((id1.EQ.11.OR.id1.EQ.13.OR.id1.EQ.15).AND.(id2
322 & .EQ.1.OR.id2.EQ.3.OR.id2.EQ.5))
THEN
323 pbrat(1)=pbrat(1)+brat(idc)
324 nmodes(1)=nmodes(1)+1
325 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
326 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
327 ELSE IF ((id1.EQ.1.OR.id1.EQ.3.OR.id1.EQ.5).AND.(id2
328 & .EQ.1.OR.id2.EQ.3.OR.id2.EQ.5))
THEN
329 pbrat(2)=pbrat(2)+brat(idc)
330 nmodes(2)=nmodes(2)+1
331 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
332 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
338 IF (kfsm.EQ.11.OR.kfsm.EQ.13.OR.kfsm.EQ.15)
THEN
346 chdc(1)=dname(3) //
' + ' // dname(2)
347 chdc(2)=dname(1) //
' + ' // dname(1)
349 DO 240
j=1,mdcy(kc,3)
351 id1=iabs(kfdp(idc,1))
352 id2=iabs(kfdp(idc,2))
353 IF (kfdp(idc,3).EQ.0)
THEN
354 IF ((id1.EQ.12.OR.id1.EQ.14.OR.id1.EQ.16).AND.(id2
355 & .EQ.11.OR.id2.EQ.13.OR.id2.EQ.15))
THEN
356 pbrat(1)=pbrat(1)+brat(idc)
357 nmodes(1)=nmodes(1)+1
358 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
359 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
361 IF ((id1.EQ.2.OR.id1.EQ.4.OR.id1.EQ.6).AND.(id2
362 & .EQ.1.OR.id2.EQ.3.OR.id2.EQ.5))
THEN
363 pbrat(2)=pbrat(2)+brat(idc)
364 nmodes(2)=nmodes(2)+1
365 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
366 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
372 IF ((kfsm.EQ.12.OR.kfsm.EQ.14.OR.kfsm.EQ.16).AND.ilr.EQ.1)
381 chdc(1)=dname(2) //
' + ' // dname(2)
382 chdc(2)=dname(1) //
' + ' // dname(1)
384 DO 260
j=1,mdcy(kc,3)
386 id1=iabs(kfdp(idc,1))
387 id2=iabs(kfdp(idc,2))
388 IF (kfdp(idc,3).EQ.0)
THEN
389 IF ((id1.EQ.11.OR.id1.EQ.13.OR.id1.EQ.15).AND.(id2
390 & .EQ.11.OR.id2.EQ.13.OR.id2.EQ.15))
THEN
391 pbrat(1)=pbrat(1)+brat(idc)
392 nmodes(1)=nmodes(1)+1
393 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
394 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
396 IF ((id1.EQ.1.OR.id1.EQ.3.OR.id1.EQ.5).AND.(id2
397 & .EQ.1.OR.id2.EQ.3.OR.id2.EQ.5))
THEN
398 nmodes(2)=nmodes(2)+1
399 pbrat(2)=pbrat(2)+brat(idc)
400 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
401 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
408 WRITE (mstu(11),8200) chd0, chdc(
i), pbrat(
i), nmodes(
i)
409 nmodes(0)=nmodes(0)+nmodes(
i)
418 IF (kfsm.EQ.22.OR.kfsm.EQ.23.OR.kfsm.EQ.25.OR.kfsm.EQ.35)
THEN
426 chdc(1)=dname(3) //
' + ' // dname(2) //
' + ' // dname(2)
427 chdc(2)=dname(3) //
' + ' // dname(1) //
' + ' // dname(1)
428 chdc(3)=dname(2) //
' + ' // dname(1) //
' + ' // dname(1)
429 chdc(4)=dname(1) //
' + ' // dname(1) //
' + ' // dname(1)
431 DO 310
j=1,mdcy(kc,3)
433 id1=iabs(kfdp(idc,1))
434 id2=iabs(kfdp(idc,2))
435 id3=iabs(kfdp(idc,3))
436 IF ((id1.EQ.12.OR.id1.EQ.14.OR.id1.EQ.16).AND.(id2
437 & .EQ.11.OR.id2.EQ.13.OR.id2.EQ.15).AND.(id3.EQ.11.or
438 & .id3.EQ.13.OR.id3.EQ.15))
THEN
439 pbrat(1)=pbrat(1)+brat(idc)
440 nmodes(1)=nmodes(1)+1
441 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
442 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
443 ELSE IF ((id1.EQ.12.OR.id1.EQ.14.OR.id1.EQ.16).and
444 & .(id2.EQ.1.OR.id2.EQ.3.OR.id2.EQ.5).AND.(id3.EQ.1
445 & .OR.id3.EQ.3.OR.id3.EQ.5))
THEN
446 pbrat(2)=pbrat(2)+brat(idc)
447 nmodes(2)=nmodes(2)+1
448 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
449 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
450 ELSE IF ((id1.EQ.11.OR.id1.EQ.13.OR.id1.EQ.15).and
451 & .(id2.EQ.2.OR.id2.EQ.4.OR.id2.EQ.6).AND.(id3.EQ.1
452 & .OR.id3.EQ.3.OR.id3.EQ.5))
THEN
453 pbrat(3)=pbrat(3)+brat(idc)
454 nmodes(3)=nmodes(3)+1
455 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
456 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
457 ELSE IF ((id1.EQ.2.OR.id1.EQ.4.OR.id1.EQ.6).and
458 & .(id2.EQ.1.OR.id2.EQ.3.OR.id2.EQ.5).AND.(id3.EQ.1
459 & .OR.id3.EQ.3.OR.id3.EQ.5))
THEN
460 pbrat(4)=pbrat(4)+brat(idc)
461 nmodes(4)=nmodes(4)+1
462 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
463 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
468 IF (kfsm.EQ.24.OR.kfsm.EQ.37)
THEN
476 chdc(1)=dname(3) //
' + ' // dname(3) //
' + ' // dname(2)
477 chdc(2)=dname(2) //
' + ' // dname(2) //
' + ' // dname(2)
478 chdc(3)=dname(3) //
' + ' // dname(1) //
' + ' // dname(1)
479 chdc(4)=dname(2) //
' + ' // dname(1) //
' + ' // dname(1)
480 chdc(5)=dname(1) //
' + ' // dname(1) //
' + ' // dname(1)
482 DO 330
j=1,mdcy(kc,3)
484 id1=iabs(kfdp(idc,1))
485 id2=iabs(kfdp(idc,2))
486 id3=iabs(kfdp(idc,3))
487 IF ((id1.EQ.12.OR.id1.EQ.14.OR.id1.EQ.16).AND.(id2
488 & .EQ.11.OR.id2.EQ.13.OR.id2.EQ.15).AND.(id3.EQ.12.or
489 & .id3.EQ.14.OR.id3.EQ.16))
THEN
490 pbrat(1)=pbrat(1)+brat(idc)
491 nmodes(1)=nmodes(1)+1
492 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
493 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
494 ELSE IF ((id1.EQ.12.OR.id1.EQ.14.OR.id1.EQ.16).and
495 & .(id2.EQ.12.OR.id2.EQ.14.OR.id2.EQ.16).AND.(id3.eq
496 & .11.OR.id3.EQ.13.OR.id3.EQ.15))
THEN
497 pbrat(1)=pbrat(1)+brat(idc)
498 nmodes(1)=nmodes(1)+1
499 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
500 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
501 ELSE IF ((id1.EQ.11.OR.id1.EQ.13.OR.id1.EQ.15).and
502 & .(id2.EQ.11.OR.id2.EQ.13.OR.id2.EQ.15).AND.(id3.eq
503 & .11.OR.id3.EQ.13.OR.id3.EQ.15))
THEN
504 pbrat(2)=pbrat(2)+brat(idc)
505 nmodes(2)=nmodes(2)+1
506 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
507 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
508 ELSE IF ((id1.EQ.12.OR.id1.EQ.14.OR.id1.EQ.16).and
509 & .(id2.EQ.1.OR.id2.EQ.3.OR.id2.EQ.5).AND.(id3.eq
510 & .2.OR.id3.EQ.4.OR.id3.EQ.6))
THEN
511 pbrat(3)=pbrat(3)+brat(idc)
512 nmodes(3)=nmodes(3)+1
513 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
514 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
515 ELSE IF ((id1.EQ.12.OR.id1.EQ.14.OR.id1.EQ.16).and
516 & .(id2.EQ.2.OR.id2.EQ.4.OR.id2.EQ.6).AND.(id3.eq
517 & .1.OR.id3.EQ.3.OR.id3.EQ.5))
THEN
518 pbrat(3)=pbrat(3)+brat(idc)
519 nmodes(3)=nmodes(3)+1
520 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
521 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
522 ELSE IF ((id1.EQ.11.OR.id1.EQ.13.OR.id1.EQ.15).and
523 & .(id2.EQ.2.OR.id2.EQ.4.OR.id2.EQ.6).AND.(id3.eq
524 & .2.OR.id3.EQ.4.OR.id3.EQ.6))
THEN
525 pbrat(4)=pbrat(4)+brat(idc)
526 nmodes(4)=nmodes(4)+1
527 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
528 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
529 ELSE IF ((id1.EQ.11.OR.id1.EQ.13.OR.id1.EQ.15).and
530 & .(id2.EQ.1.OR.id2.EQ.3.OR.id2.EQ.5).AND.(id3.eq
531 & .1.OR.id3.EQ.3.OR.id3.EQ.5))
THEN
532 pbrat(4)=pbrat(4)+brat(idc)
533 nmodes(4)=nmodes(4)+1
534 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
535 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
536 ELSE IF ((id1.EQ.2.OR.id1.EQ.4.OR.id1.EQ.6).and
537 & .(id2.EQ.2.OR.id2.EQ.4.OR.id2.EQ.6).AND.(id3.eq
538 & .1.OR.id3.EQ.3.OR.id3.EQ.5))
THEN
539 pbrat(5)=pbrat(5)+brat(idc)
540 nmodes(5)=nmodes(5)+1
541 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
542 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
543 ELSE IF ((id1.EQ.1.OR.id1.EQ.3.OR.id1.EQ.5).and
544 & .(id2.EQ.1.OR.id2.EQ.3.OR.id2.EQ.5).AND.(id3.eq
545 & .1.OR.id3.EQ.3.OR.id3.EQ.5))
THEN
546 pbrat(5)=pbrat(5)+brat(idc)
547 nmodes(5)=nmodes(5)+1
548 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
549 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
562 chdc(1)=dname(3) //
' + ' // dname(1) //
' + ' // dname(1)
563 chdc(2)=dname(2) //
' + ' // dname(1) //
' + ' // dname(1)
564 chdc(3)=dname(1) //
' + ' // dname(1) //
' + ' // dname(1)
566 DO 350
j=1,mdcy(kc,3)
568 id1=iabs(kfdp(idc,1))
569 id2=iabs(kfdp(idc,2))
570 id3=iabs(kfdp(idc,3))
571 IF ((id1.EQ.12.OR.id1.EQ.14.OR.id1.EQ.16).AND.(id2
572 & .EQ.1.OR.id2.EQ.3.OR.id2.EQ.5).AND.(id3.EQ.1.or
573 & .id3.EQ.3.OR.id3.EQ.5))
THEN
574 pbrat(1)=pbrat(1)+brat(idc)
575 nmodes(1)=nmodes(1)+1
576 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
577 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
578 ELSE IF ((id1.EQ.11.OR.id1.EQ.13.OR.id1.EQ.15).and
579 & .(id2.EQ.2.OR.id2.EQ.4.OR.id2.EQ.6).AND.(id3.EQ.1
580 & .OR.id3.EQ.3.OR.id3.EQ.5))
THEN
581 pbrat(2)=pbrat(2)+brat(idc)
582 nmodes(2)=nmodes(2)+1
583 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
584 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
585 ELSE IF ((id1.EQ.2.OR.id1.EQ.4.OR.id1.EQ.6).and
586 & .(id2.EQ.1.OR.id2.EQ.3.OR.id2.EQ.5).AND.(id3.EQ.1
587 & .OR.id3.EQ.3.OR.id3.EQ.5))
THEN
588 pbrat(3)=pbrat(3)+brat(idc)
589 nmodes(3)=nmodes(3)+1
590 IF (brat(idc).GT.0d0) nmodes(10)=nmodes(10)+1
591 IF (brat(idc).GT.
eps) nmodes(9)=nmodes(9)+1
598 WRITE (mstu(11),8200) chd0, chdc(
i), pbrat(
i), nmodes(
i)
599 nmodes(0)=nmodes(0)+nmodes(
i)
603 WRITE (mstu(11),8100) nmodes(0), nmodes(10), nmodes(9)
605 IF (imss(51).GE.1.OR.imss(52).GE.1.OR.imss(53).GE.1)
THEN
606 WRITE (mstu(11),8500)
610 WRITE (mstu(11),8700) irv,jrv,krv,rvlam(irv,jrv,krv)
611 & ,rvlamp(irv,jrv,krv),rvlamb(irv,jrv,krv)
615 WRITE (mstu(11),8600)
620 5000
FORMAT(
'1',9(
'*'),1
x,
'PYSTAT: Statistics on Number of ',
621 &
'Events and Cross-sections',1
x,9(
'*'))
622 5100
FORMAT(/1
x,78(
'=')/1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,
'I',12
x,
623 &
'Subprocess',12
x,
'I',6
x,
'Number of points',6
x,
'I',4
x,
'Sigma',3
x,
624 &
'I'/1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,
'I',34(
'-'),
'I',28(
'-'),
625 &
'I',4
x,
'(mb)',4
x,
'I'/1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,
'I',1
x,
626 &
'N:o',1
x,
'Type',25
x,
'I',4
x,
'Generated',9
x,
'Tried',1
x,
'I',12
x,
627 &
'I'/1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,78(
'=')/1
x,
'I',34
x,
'I',28
x,
629 5200
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,
631 5300
FORMAT(1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,78(
'=')/
632 &1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I')
633 5400
FORMAT(1
x,
'I',34
x,
'I',28
x,
'I',12
x,
'I'/1
x,78(
'=')//
634 &1
x,
'********* Total number of errors, excluding junctions =',
635 &1
x,i8,
' *************'/
636 &1
x,
'********* Total number of errors, including junctions =',
637 &1
x,i8,
' *************'/
638 &1
x,
'********* Total number of warnings = ',
639 &1
x,i8,
' *************'/
640 &1
x,
'********* Fraction of events that fail fragmentation ',
641 &
'cuts =',1
x,f8.5,
' *********'/)
642 5500
FORMAT(
'1',27(
'*'),1
x,
'PYSTAT: Decay Widths and Branching ',
643 &
'Ratios',1
x,27(
'*'))
644 5600
FORMAT(/1
x,98(
'=')/1
x,
'I',49
x,
'I',13
x,
'I',12
x,
'I',6
x,
'I',12
x,
'I'/
645 &1
x,
'I',5
x,
'Mother --> Branching/Decay Channel',8
x,
'I',1
x,
646 &
'Width (GeV)',1
x,
'I',7
x,
'B.R.',1
x,
'I',1
x,
'Stat',1
x,
'I',2
x,
647 &
'Eff. B.R.',1
x,
'I'/1
x,
'I',49
x,
'I',13
x,
'I',12
x,
'I',6
x,
'I',12
x,
'I'/
649 5700
FORMAT(1
x,
'I',49
x,
'I',13
x,
'I',12
x,
'I',6
x,
'I',12
x,
'I'/1
x,
'I',1
x,
650 &i8,2
x,a10,3
x,
'(m =',f10.3,
')',2
x,
'-->',5
x,
'I',2
x,1
p,d10.3,0
p,1
x,
651 &
'I',1
x,1
p,d10.3,0
p,1
x,
'I',1
x,a4,1
x,
'I',1
x,1
p,d10.3,0
p,1
x,
'I')
652 5800
FORMAT(1
x,
'I',1
x,i8,2
x,a10,1
x,
'+',1
x,a10,15
x,
'I',2
x,
653 &1
p,d10.3,0
p,1
x,
'I',1
x,1
p,d10.3,0
p,1
x,
'I',1
x,a4,1
x,
'I',1
x,
655 5900
FORMAT(1
x,
'I',1
x,i8,2
x,a10,1
x,
'+',1
x,a10,1
x,
'+',1
x,a10,2
x,
'I',2
x,
656 &1
p,d10.3,0
p,1
x,
'I',1
x,1
p,d10.3,0
p,1
x,
'I',1
x,a4,1
x,
'I',1
x,
658 6000
FORMAT(1
x,
'I',49
x,
'I',13
x,
'I',12
x,
'I',6
x,
'I',12
x,
'I'/1
x,98(
'='))
659 6100
FORMAT(
'1',7(
'*'),1
x,
'PYSTAT: Allowed Incoming Partons/',
660 &
'Particles at Hard Interaction',1
x,7(
'*'))
661 6200
FORMAT(/1
x,78(
'=')/1
x,
'I',38
x,
'I',37
x,
'I'/1
x,
'I',1
x,
662 &
'Beam particle:',1
x,a12,10
x,
'I',1
x,
'Target particle:',1
x,a12,7
x,
663 &
'I'/1
x,
'I',38
x,
'I',37
x,
'I'/1
x,
'I',1
x,
'Content',6
x,
'State',19
x,
664 &
'I',1
x,
'Content',6
x,
'State',18
x,
'I'/1
x,
'I',38
x,
'I',37
x,
'I'/1
x,
665 &78(
'=')/1
x,
'I',38
x,
'I',37
x,
'I')
666 6300
FORMAT(1
x,
'I',1
x,a9,5
x,a4,19
x,
'I',1
x,a9,5
x,a4,18
x,
'I')
667 6400
FORMAT(1
x,
'I',38
x,
'I',37
x,
'I'/1
x,78(
'='))
668 6500
FORMAT(
'1',12(
'*'),1
x,
'PYSTAT: User-Defined Limits on ',
669 &
'Kinematical Variables',1
x,12(
'*'))
670 6600
FORMAT(/1
x,78(
'=')/1
x,
'I',76
x,
'I')
671 6700
FORMAT(1
x,
'I',16
x,1
p,d10.3,0
p,1
x,
'<',1
x,a,1
x,
'<',1
x,1
p,d10.3,0
p,
673 6800
FORMAT(1
x,
'I',3
x,1
p,d10.3,0
p,1
x,
'(',1
p,d10.3,0
p,
')',1
x,
'<',1
x,a,
674 &1
x,
'<',1
x,1
p,d10.3,0
p,16
x,
'I')
675 6900
FORMAT(1
x,
'I',29
x,a,1
x,
'=',1
x,1
p,d10.3,0
p,16
x,
'I')
676 7000
FORMAT(1
x,
'I',76
x,
'I'/1
x,78(
'='))
677 7100
FORMAT(
'1',12(
'*'),1
x,
'PYSTAT: Summary of Status Codes and ',
678 &
'Parameter Values',1
x,12(
'*'))
679 7200
FORMAT(/3
x,
'I',4
x,
'MSTP(I)',9
x,
'PARP(I)',20
x,
'I',4
x,
'MSTP(I)',9
x,
681 7300
FORMAT(1
x,
i3,5
x,i6,6
x,1
p,d10.3,0
p,18
x,
i3,5
x,i6,6
x,1
p,d10.3)
682 7400
FORMAT(
'1',13(
'*'),1
x,
'PYSTAT: List of implemented processes',
684 7500
FORMAT(/1
x,65(
'=')/1
x,
'I',34
x,
'I',28
x,
'I'/1
x,
'I',12
x,
685 &
'Subprocess',12
x,
'I',1
x,
'ISET',2
x,
'KFPR(I,1)',2
x,
'KFPR(I,2)',1
x,
686 &
'I'/1
x,
'I',34
x,
'I',28
x,
'I'/1
x,65(
'=')/1
x,
'I',34
x,
'I',28
x,
'I')
687 7600
FORMAT(1
x,
'I',1
x,
i3,1
x,a28,1
x,
'I',1
x,
i4,1
x,i10,1
x,i10,1
x,
'I')
688 7700
FORMAT(1
x,
'I',34
x,
'I',28
x,
'I'/1
x,65(
'='))
690 & 17
x,
'Sums over R-Violating branching ratios',1
x/ 1
x
691 & /1
x,70(
'=')/1
x,
'I',50
x,
'I',11
x,
'I',5
x,
'I'/1
x,
'I',4
x
692 & ,
'Mother --> Sum over final state flavours',4
x,
'I',2
x
693 & ,
'BR(sum)',2
x,
'I',2
x,
'N',2
x,
'I'/1
x,
'I',50
x,
'I',11
x,
'I',5
x,
'I'
694 & /1
x,70(
'=')/1
x,
'I',50
x,
'I',11
x,
'I',5
x,
'I')
695 8100
FORMAT(1
x,
'I',50
x,
'I',11
x,
'I',5
x,
'I'/1
x,70(
'=')/1
x,
'I',1
x
696 & ,
'Total number of R-Violating modes :',3
x,
i5,24
x,
'I'/
697 & 1
x,
'I',1
x,
'Total number with non-vanishing BR :',2
x,
i5,24
x
698 & ,
'I'/1
x,
'I',1
x,
'Total number with BR > 0.001 :',8
x,
i5,24
x,
'I'
700 8200
FORMAT(1
x,
'I',1
x,a9,1
x,
'-->',1
x,a24,11
x,
701 &
'I',2
x,1
p,d8.2,0
p,1
x,
'I',2
x,
i2,1
x,
'I')
702 8300
FORMAT(1
x,
'I',50
x,
'I',11
x,
'I',5
x,
'I')
704 & 1
x,
'R-Violating couplings',1
x/ 1
x /
706 & 1
x,
'I',1
x,
'IJK',1
x,
'I',2
x,
'LAMBDA(IJK)',2
x,
'I',2
x
707 & ,
'LAMBDA''(IJK)',1
x,
'I',1
x,
"LAMBDA''(IJK)",1
x,
'I'/1
x,
'I',5
x
708 & ,
'I',15
x,
'I',15
x,
'I',15
x,
'I')
709 8600
FORMAT(1
x,55(
'='))
710 8700
FORMAT(1
x,
'I',1
x,
i1,
i1,
i1,1
x,
'I',1
x,1
p,d13.3,0
p,1
x,
'I',1
x,1
p
711 & ,d13.3,0
p,1
x,
'I',1
x,1
p,d13.3,0
p,1
x,
'I')