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 dimension psum(5),pini(6),pfin(6)
15 IF(mtest.GE.1) CALL
lutabu(20)
34 IF(iev.EQ.301.OR.iev.EQ.351.OR.iev.EQ.401) mstj(116)=3
40 IF(ity.EQ.3.OR.ity.EQ.4) mstj(11)=2
41 IF(ity.EQ.1) CALL
lu1ent(1,1,15.,0.,0.)
42 IF(ity.EQ.2) CALL
lu1ent(1,3101,15.,0.,0.)
43 IF(ity.EQ.3) CALL
lu1ent(1,-2203,15.,0.,0.)
44 IF(ity.EQ.4) CALL
lu1ent(1,-4,30.,0.,0.)
45 IF(ity.EQ.5) CALL
lu1ent(1,21,15.,0.,0.)
48 ELSEIF(iev.LE.130)
THEN
50 IF(ity.EQ.1) CALL
lu2ent(1,1,-1,40.)
51 IF(ity.EQ.2) CALL
lu2ent(1,4,-4,30.)
52 IF(ity.EQ.3) CALL
lu2ent(1,2,2103,100.)
53 IF(ity.EQ.4) CALL
lu2ent(1,21,21,40.)
54 IF(ity.EQ.5) CALL
lu3ent(1,2101,21,-3203,30.,0.6,0.8)
55 IF(ity.EQ.6) CALL
lu3ent(1,5,21,-5,40.,0.9,0.8)
56 IF(ity.EQ.7) CALL
lu3ent(1,21,21,21,60.,0.7,0.5)
57 IF(ity.EQ.8) CALL
lu4ent(1,2,21,21,-2,40.,0.4,0.64,0.6,0.12,0.2)
60 ELSEIF(iev.LE.200)
THEN
62 mstj(2)=1+mod(iev-131,4)
63 mstj(3)=1+mod((iev-131)/4,4)
64 IF(ity.EQ.1) CALL
lu2ent(1,4,-5,40.)
65 IF(ity.EQ.2) CALL
lu3ent(1,3,21,-3,40.,0.9,0.4)
66 IF(ity.EQ.3) CALL
lu4ent(1,2,21,21,-2,40.,0.4,0.64,0.6,0.12,0.2)
67 IF(ity.GE.4) CALL
lu4ent(1,2,-3,3,-2,40.,0.4,0.64,0.6,0.12,0.2)
70 ELSEIF(iev.LE.300)
THEN
76 IF(
i.EQ.1) kfl=int(1.+4.*
rlu(0))
77 IF(
i.EQ.njet) kfl=-int(1.+4.*
rlu(0))
79 theta=acos(2.*
rlu(0)-1.)
83 IF(
i.EQ.1.OR.
i.EQ.njet) psum(5)=psum(5)+
ulmass(kfl)
85 120 psum(
j)=psum(
j)+
p(
i,
j)
86 IF(psum(4)**2-psum(1)**2-psum(2)**2-psum(3)**2.LT.
87 & (psum(5)+parj(32))**2) goto 100
90 ELSEIF(iev.LE.350)
THEN
95 ELSEIF(iev.LE.400)
THEN
97 mstj(43)=1+mod(iev/2,4)
102 ELSEIF(iev.LE.450)
THEN
107 ELSEIF(iev.LE.500)
THEN
111 ELSEIF(iev.LE.560)
THEN
116 CALL
lu1ent(1,100*kflb+10*kflc+kfls,0.,0.,0.)
119 ELSEIF(iev.LE.600)
THEN
125 CALL
lu1ent(1,1000*kfla+100*kflb+10*kflc+kfls,0.,0.,0.)
141 IF((pfin(1)-pini(1))**2+(pfin(2)-pini(2))**2.GE.4.) merr=merr+1
142 epzrem=pini(4)+pini(3)-pfin(4)-pfin(3)
143 IF(epzrem.LT.0..OR.epzrem.GT.2.*parj(31)) merr=merr+1
144 IF(abs(pfin(6)-pini(6)).GT.2.1) merr=merr+1
147 150
IF(abs(pfin(
j)-pini(
j)).GT.0001*pini(4)) merr=merr+1
148 IF(abs(pfin(6)-pini(6)).GT.0.1) merr=merr+1
150 IF(merr.NE.0)
WRITE(mstu(11),1000) (pini(
j),
j=1,4),pini(6),
151 &(pfin(
j),
j=1,4),pfin(6)
156 IF(
k(
i,1).GT.20) goto 160
158 WRITE(mstu(11),1100)
i
161 pd=
p(
i,4)**2-
p(
i,1)**2-
p(
i,2)**2-
p(
i,3)**2-
p(
i,5)**2
162 IF(abs(pd).GT.
max(0.1,0.001*
p(
i,4)**2).OR.
p(
i,4).LT.0.)
THEN
163 WRITE(mstu(11),1200)
i
167 IF(mtest.GE.1) CALL
lutabu(21)
170 IF(merr.NE.0.OR.mstu(24).NE.0.OR.mstu(28).NE.0)
THEN
172 ELSEIF(mtest.GE.1.AND.mod(iev-5,100).EQ.0)
THEN
177 IF(merr.NE.0) nerr=nerr+1
179 WRITE(mstu(11),1300) iev
183 IF(mtest.GE.1) CALL
lutabu(22)
184 WRITE(mstu(11),1400) nerr
196 1000
FORMAT(/
' Momentum, energy and/or charge were not conserved ',
197 &
'in following event'/
' sum of',9
x,
'px',11
x,
'py',11
x,
'pz',11
x,
198 &
'E',8
x,
'charge'/
' before',2
x,4(1
x,f12.5),1
x,f8.2/
' after',3
x,
199 &4(1
x,f12.5),1
x,f8.2)
200 1100
FORMAT(/5
x,
'Entry no.',
i4,
' in following event not known code')
201 1200
FORMAT(/5
x,
'Entry no.',
i4,
' in following event has faulty ',
203 1300
FORMAT(/5
x,
'Ten errors experienced by event ',
i3/
204 &5
x,
'Something is seriously wrong! Execution stopped now!')
205 1400
FORMAT(/5
x,
'Number of erroneous or suspect events in run:',
i3/
206 &5
x,
'(0 fine, 1 acceptable if a single jet, ',
207 &
'>=2 something is wrong)')