7 common/lujets/
n,
k(9000,5),
p(9000,5),
v(9000,5)
9 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
11 common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)
13 common/ludat3/mdcy(500,3),mdme(2000,2),brat(2000),kfdp(2000,5)
19 IF(mstu(12).GE.1) CALL
lulist(0)
32 IF(
k(
i,1).LE.0.OR.
k(
i,1).GT.10) goto 120
48 IF(
k(ip,1).GT.0.AND.
k(ip,1).LE.10) kc=
lucomp(
k(ip,2))
53 ELSEIF(kchg(kc,2).EQ.0)
THEN
54 IF(mstj(21).GE.1.AND.mdcy(kc,1).GE.1.AND.(mstj(51).LE.0.OR.mbe.
55 & eq.2.OR.pmas(kc,2).GE.parj(91).OR.iabs(
k(ip,2)).EQ.311))
59 IF(mstj(92).GT.0)
THEN
61 qmax=sqrt(
max(0.,(
p(ip1,4)+
p(ip1+1,4))**2-(
p(ip1,1)+
p(ip1+1,
62 & 1))**2-(
p(ip1,2)+
p(ip1+1,2))**2-(
p(ip1,3)+
p(ip1+1,3))**2))
63 CALL
lushow(ip1,ip1+1,qmax)
66 ELSEIF(mstj(92).LT.0)
THEN
74 ELSEIF(
k(ip,1).EQ.1.OR.
k(ip,1).EQ.2)
THEN
76 IF(mfrag.GE.1.AND.
k(ip,1).EQ.1) mfrag=2
77 IF(mstj(21).GE.2.AND.
k(ip,1).EQ.2.AND.
n.GT.ip)
THEN
78 IF(
k(ip+1,1).EQ.1.AND.
k(ip+1,3).EQ.
k(ip,3).AND.
79 &
k(ip,3).GT.0.AND.
k(ip,3).LT.ip)
THEN
80 IF(kchg(
lucomp(
k(
k(ip,3),2)),2).EQ.0) mfrag=min(1,mfrag)
83 IF(mfrag.EQ.1) CALL
lustrf(ip)
84 IF(mfrag.EQ.2) CALL
luindf(ip)
85 IF(mfrag.EQ.2.AND.
k(ip,1).EQ.1) mcons=0
86 IF(mfrag.EQ.2.AND.(mstj(3).LE.0.OR.mod(mstj(3),5).EQ.0)) mcons=0
90 IF(mstu(24).NE.0.AND.mstu(21).GE.2)
THEN
91 ELSEIF(ip.LT.
n.AND.
n.LT.mstu(4)-20-mstu(32))
THEN
94 CALL
luerrm(11,
'(LUEXEC:) no more memory left in LUJETS')
98 IF(mbe.EQ.1.AND.mstj(51).GE.1)
THEN
105 IF(
k(
i,1).LE.0.OR.
k(
i,1).GT.10) goto 160
110 pdev=(abs(
ps(2,1)-
ps(1,1))+abs(
ps(2,2)-
ps(1,2))+abs(
ps(2,3)-
111 &
ps(1,3))+abs(
ps(2,4)-
ps(1,4)))/(1.+abs(
ps(2,4))+abs(
ps(1,4)))
112 IF(mcons.EQ.1.AND.pdev.GT.paru(11)) CALL
luerrm(15,
113 &
'(LUEXEC:) four-momentum was not conserved')
114 IF(mcons.EQ.1.AND.abs(
ps(2,6)-
ps(1,6)).GT.0.1) CALL
luerrm(15,
115 &
'(LUEXEC:) charge was not conserved')