8 common/lujets/
n,
k(9000,5),
p(9000,5),
v(9000,5)
10 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
12 common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)
14 common/ludat3/mdcy(500,3),mdme(2000,2),brat(2000),kfdp(2000,5)
16 CHARACTER chap*16,chac*16,chan*16,chad(5)*16,chmo(12)*3,chdl(7)*4
18 DATA chmo/
'Jan',
'Feb',
'Mar',
'Apr',
'May',
'Jun',
'Jul',
'Aug',
'Sep',
19 &
'Oct',
'Nov',
'Dec'/,chdl/
'(())',
' ',
'()',
'!!',
'<>',
'==',
'(==)'/
30 IF(mlist.GE.1.AND.mlist.LE.3)
THEN
31 IF(mlist.EQ.1)
WRITE(mstu(11),1100)
32 IF(mlist.EQ.2)
WRITE(mstu(11),1200)
33 IF(mlist.EQ.3)
WRITE(mstu(11),1300)
38 IF(mstu(2).GT.0) imax=mstu(2)
40 IF((
i.GT.imax.AND.
i.LE.
n).OR.
k(
i,1).LT.0) goto 120
46 100
IF(chap(lem:lem).NE.
' ') len=lem
49 IF(mdl.EQ.2.OR.mdl.GE.8)
THEN
51 IF(len.GT.lmx) chac(lmx:lmx)=
'?'
54 IF(mdl.EQ.1.OR.mdl.EQ.7) ldl=2
56 chac=chdl(mdl)(1:2*ldl)//
' '
58 chac=chdl(mdl)(1:ldl)//chap(1:min(len,lmx-2*ldl))//
59 & chdl(mdl)(ldl+1:2*ldl)//
' '
60 IF(len+2*ldl.GT.lmx) chac(lmx:lmx)=
'?'
65 IF(
k(
i,1).EQ.1.OR.
k(
i,1).EQ.2.OR.
k(
i,1).EQ.11.OR.
k(
i,1).EQ.12)
69 IF(kc.NE.0) kcc=kchg(kc,2)
70 IF(kcc.NE.0.AND.istr.EQ.0)
THEN
72 IF(len+2*ldl+3.LE.lmx) chac(lmx-1:lmx-1)=
'A'
73 ELSEIF(kcc.NE.0.AND.(
k(
i,1).EQ.2.OR.
k(
i,1).EQ.12))
THEN
74 IF(len+2*ldl+3.LE.lmx) chac(lmx-1:lmx-1)=
'I'
77 IF(len+2*ldl+3.LE.lmx) chac(lmx-1:lmx-1)=
'V'
82 IF(mlist.EQ.1.AND.abs(
p(
i,4)).LT.9999.)
THEN
83 WRITE(mstu(11),1400)
i,chac(1:12),(
k(
i,j1),j1=1,3),
85 ELSEIF(mlist.EQ.1.AND.abs(
p(
i,4)).LT.99999.)
THEN
86 WRITE(mstu(11),1500)
i,chac(1:12),(
k(
i,j1),j1=1,3),
88 ELSEIF(mlist.EQ.1)
THEN
89 WRITE(mstu(11),1600)
i,chac(1:12),(
k(
i,j1),j1=1,3),
91 ELSEIF(mstu(5).EQ.10000.AND.(
k(
i,1).EQ.3.OR.
k(
i,1).EQ.13.OR.
93 WRITE(mstu(11),1700)
i,chac,(
k(
i,j1),j1=1,3),
94 &
k(
i,4)/100000000,mod(
k(
i,4)/10000,10000),mod(
k(
i,4),10000),
95 &
k(
i,5)/100000000,mod(
k(
i,5)/10000,10000),mod(
k(
i,5),10000),
98 WRITE(mstu(11),1800)
i,chac,(
k(
i,j1),j1=1,5),(
p(
i,j2),j2=1,5)
100 IF(mlist.EQ.3)
WRITE(mstu(11),1900) (
v(
i,
j),
j=1,5)
103 IF(mstu(70).GE.1)
THEN
105 DO 110
j=1,min(10,mstu(70))
106 110
IF(
i.EQ.mstu(70+
j)) isep=1
107 IF(isep.EQ.1.AND.mlist.EQ.1)
WRITE(mstu(11),2000)
108 IF(isep.EQ.1.AND.mlist.GE.2)
WRITE(mstu(11),2100)
115 IF(mlist.EQ.1.AND.abs(
ps(4)).LT.9999.)
THEN
116 WRITE(mstu(11),2200)
ps(6),(
ps(
j),
j=1,5)
117 ELSEIF(mlist.EQ.1.AND.abs(
ps(4)).LT.99999.)
THEN
118 WRITE(mstu(11),2300)
ps(6),(
ps(
j),
j=1,5)
119 ELSEIF(mlist.EQ.1)
THEN
120 WRITE(mstu(11),2400)
ps(6),(
ps(
j),
j=1,5)
122 WRITE(mstu(11),2500)
ps(6),(
ps(
j),
j=1,5)
126 ELSEIF(mlist.EQ.11)
THEN
131 IF(chap.NE.
' '.AND.chan.EQ.
' ')
WRITE(mstu(11),2700) kf,chap
132 140
IF(chan.NE.
' ')
WRITE(mstu(11),2700) kf,chap,-kf,chan
135 DO 150 kflb=1,kfla-(3-kfls)/2
136 kf=1000*kfla+100*kflb+kfls
139 150
WRITE(mstu(11),2700) kf,chap,-kf,chan
142 IF(kmul.EQ.0.OR.kmul.EQ.3) kfls=1
145 IF(kmul.EQ.2.OR.kmul.EQ.3) kflr=1
149 kf=10000*kflr+100*kflb+10*kflc+kfls
152 160
WRITE(mstu(11),2700) kf,chap,-kf,chan
153 kf=10000*kflr+110*kflb+kfls
155 170
WRITE(mstu(11),2700) kf,chap
158 WRITE(mstu(11),2700) kf,chap
161 WRITE(mstu(11),2700) kf,chap
167 IF(kflsp.EQ.1.AND.(kfla.EQ.kflb.OR.kflb.EQ.kflc)) goto 180
168 IF(kflsp.EQ.2.AND.kfla.EQ.kflc) goto 180
169 IF(kflsp.EQ.1) kf=1000*kfla+100*kflc+10*kflb+kfls
170 IF(kflsp.GE.2) kf=1000*kfla+100*kflb+10*kflc+kfls
173 WRITE(mstu(11),2700) kf,chap,-kf,chan
178 ELSEIF(mlist.EQ.12)
THEN
183 IF(mstu(2).NE.0) kfmax=mstu(2)
184 DO 220 kf=
max(1,mstu(1)),kfmax
187 IF(mstu(14).EQ.0.AND.kf.GT.100.AND.kc.LE.100) goto 220
188 IF(mstu(14).GT.0.AND.kf.GT.100.AND.
max(mod(kf/1000,10),
189 & mod(kf/100,10)).GT.mstu(14)) goto 220
193 IF(kf.LE.100.AND.chap.EQ.
' '.AND.mdcy(kc,2).EQ.0) goto 220
196 WRITE(mstu(11),2900) kf,kc,chap,chan,kchg(kc,1),kchg(kc,2),
197 & kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
201 IF(kf.GT.100.AND.kc.LE.100) goto 220
202 DO 210 idc=mdcy(kc,2),mdcy(kc,2)+mdcy(kc,3)-1
204 200 CALL
luname(kfdp(idc,
j),chad(
j))
205 210
WRITE(mstu(11),3000) idc,mdme(idc,1),mdme(idc,2),brat(idc),
211 ELSEIF(mlist.EQ.13)
THEN
214 230
WRITE(mstu(11),3200)
i,mstu(
i),paru(
i),mstj(
i),parj(
i),parf(
i)
218 1000
FORMAT(///20
x,
'The Lund Monte Carlo - JETSET version ',
i1,
'.',
i1/
219 &20
x,
'** Last date of change: ',
i2,1
x,a3,1
x,
i4,
' **'/)
220 1100
FORMAT(///28
x,
'Event listing (summary)'//4
x,
'I particle/jet KS',
221 &5
x,
'KF orig p_x p_y p_z E m'/)
222 1200
FORMAT(///28
x,
'Event listing (standard)'//4
x,
'I particle/jet',
223 &
' K(I,1) K(I,2) K(I,3) K(I,4) K(I,5) P(I,1)',
224 &
' P(I,2) P(I,3) P(I,4) P(I,5)'/)
225 1300
FORMAT(///28
x,
'Event listing (with vertices)'//4
x,
'I particle/j',
226 &
'et K(I,1) K(I,2) K(I,3) K(I,4) K(I,5) P(I,1)',
227 &
' P(I,2) P(I,3) P(I,4) P(I,5)'/73
x,
228 &
'V(I,1) V(I,2) V(I,3) V(I,4) V(I,5)'/)
229 1400
FORMAT(1
x,
i4,2
x,a12,1
x,
i2,1
x,i6,1
x,
i4,5f9.3)
230 1500
FORMAT(1
x,
i4,2
x,a12,1
x,
i2,1
x,i6,1
x,
i4,5f9.2)
231 1600
FORMAT(1
x,
i4,2
x,a12,1
x,
i2,1
x,i6,1
x,
i4,5f9.1)
232 1700
FORMAT(1
x,
i4,2
x,a16,1
x,
i3,1
x,i8,2
x,
i4,2(3
x,
i1,2
i4),5f13.5)
233 1800
FORMAT(1
x,
i4,2
x,a16,1
x,
i3,1
x,i8,2
x,
i4,2(3
x,i9),5f13.5)
234 1900
FORMAT(66
x,5(1
x,f12.3))
235 2000
FORMAT(1
x,78(
'='))
236 2100
FORMAT(1
x,130(
'='))
237 2200
FORMAT(19
x,
'sum:',f6.2,5
x,5f9.3)
238 2300
FORMAT(19
x,
'sum:',f6.2,5
x,5f9.2)
239 2400
FORMAT(19
x,
'sum:',f6.2,5
x,5f9.1)
240 2500
FORMAT(19
x,
'sum charge:',f6.2,3
x,
'sum momentum and inv. mass:',
242 2600
FORMAT(///20
x,
'List of KF codes in program'/)
243 2700
FORMAT(4
x,i6,4
x,a16,6
x,i6,4
x,a16)
244 2800
FORMAT(///30
x,
'Particle/parton data table'//5
x,
'KF',5
x,
'KC',4
x,
245 &
'particle',8
x,
'antiparticle',6
x,
'chg col anti',8
x,
'mass',7
x,
246 &
'width',7
x,
'w-cut',5
x,
'lifetime',1
x,
'decay'/11
x,
'IDC',1
x,
'on/off',
247 &1
x,
'ME',3
x,
'Br.rat.',4
x,
'decay products')
248 2900
FORMAT(/1
x,i6,3
x,
i4,4
x,a16,a16,3
i5,1
x,f12.5,2(1
x,f11.5),
251 3100
FORMAT(///20
x,
'Parameter value table'//4
x,
'I',3
x,
'MSTU(I)',
252 &8
x,
'PARU(I)',3
x,
'MSTJ(I)',8
x,
'PARJ(I)',8
x,
'PARF(I)')
253 3200
FORMAT(1
x,
i4,1
x,i9,1
x,f14.5,1
x,i9,1
x,f14.5,1
x,f14.5)