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 dimension tdi(3),tpr(3)
20 IF(
k(
i,1).LE.0.OR.
k(
i,1).GT.10) goto 100
21 IF(mstu(41).GE.2)
THEN
23 IF(kc.EQ.0.OR.kc.EQ.12.OR.kc.EQ.14.OR.kc.EQ.16.OR.
25 IF(mstu(41).GE.3.AND.kchg(kc,2).EQ.0.AND.
luchge(
k(
i,2)).EQ.0)
28 IF(
n+
np+mstu(44)+15.GE.mstu(4)-mstu(32)-5)
THEN
29 CALL
luerrm(11,
'(LUTHRU:) no more memory left in LUJETS')
41 IF(abs(paru(42)-1.).GT.0.001)
p(
n+
np,5)=
p(
n+
np,4)**(paru(42)-1.)
47 CALL
luerrm(8,
'(LUTHRU:) too few particles for analysis')
58 CALL ludbrb(
n+1,
n+
np+1,0.,-
phi,0d0,0d0,0d0)
60 CALL ludbrb(
n+1,
n+
np+1,-the,0.,0d0,0d0,0d0)
64 DO 110 ilf=
n+
np+4,
n+
np+mstu(44)+4
67 IF(ild.EQ.2)
p(
i,4)=sqrt(
p(
i,1)**2+
p(
i,2)**2)
68 DO 120 ilf=
n+
np+mstu(44)+3,
n+
np+4,-1
69 IF(
p(
i,4).LE.
p(ilf,4)) goto 130
71 120
p(ilf+1,
j)=
p(ilf,
j)
78 DO 160 ilg=
n+
np+mstu(44)+5,
n+
np+mstu(44)+15
80 nc=2**(min(mstu(44),
np)-1)
84 DO 180 ilf=1,min(mstu(44),
np)
86 IF(2**ilf*((ilc+2**(ilf-1)-1)/2**ilf).GE.ilc) sgn=-sgn
88 180 tdi(
j)=tdi(
j)+sgn*
p(
n+
np+ilf+3,
j)
89 tds=tdi(1)**2+tdi(2)**2+tdi(3)**2
90 DO 190 ilg=
n+
np+mstu(44)+min(ilc,10)+4,
n+
np+mstu(44)+5,-1
91 IF(tds.LE.
p(ilg,4)) goto 200
93 190
p(ilg+1,
j)=
p(ilg,
j)
107 IF(thp.LE.1
e-10) tdi(
j)=
p(
n+
np+mstu(44)+4+ilg,
j)
108 IF(thp.GT.1
e-10) tdi(
j)=tpr(
j)
111 sgn=
sign(
p(
i,5),tdi(1)*
p(
i,1)+tdi(2)*
p(
i,2)+tdi(3)*
p(
i,3))
113 260 tpr(
j)=tpr(
j)+sgn*
p(
i,
j)
114 thp=sqrt(tpr(1)**2+tpr(2)**2+tpr(3)**2)/
ps
115 IF(thp.GE.thps+paru(48)) goto 240
118 IF(thp.LT.
p(
n+
np+ild,4)-paru(48).AND.ilg.LT.min(10,nc)) goto 230
119 IF(thp.GT.
p(
n+
np+ild,4)+paru(48))
THEN
121 sgn=(-1.)**int(
rlu(0)+0.5)
123 270
p(
n+
np+ild,
j)=sgn*tpr(
j)/(
ps*thp)
128 280
IF(iagr.LT.mstu(45).AND.ilg.LT.min(10,nc)) goto 230
131 sgn=(-1.)**int(
rlu(0)+0.5)
151 CALL ludbrb(
n+1,
n+3,the,
phi,0d0,0d0,0d0)
156 IF(mstu(43).LE.1) mstu(3)=3
157 IF(mstu(43).GE.2)
n=
n+3
159 obl=
p(
n+2,4)-
p(
n+3,4)