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)
16 nce2=2*mstu(51)*mstu(52)
17 ptlrat=1./sinh(paru(51))**2
21 IF(
k(
i,1).LE.0.OR.
k(
i,1).GT.10) goto 110
22 IF(
p(
i,1)**2+
p(
i,2)**2.LE.ptlrat*
p(
i,3)**2) goto 110
23 IF(mstu(41).GE.2)
THEN
25 IF(kc.EQ.0.OR.kc.EQ.12.OR.kc.EQ.14.OR.kc.EQ.16.OR.
27 IF(mstu(41).GE.3.AND.kchg(kc,2).EQ.0.AND.
luchge(
k(
i,2)).EQ.0)
31 pt=sqrt(
p(
i,1)**2+
p(
i,2)**2)
33 ieta=
max(1,min(mstu(51),1+int(mstu(51)*0.5*(
eta/paru(51)+1.))))
35 iphi=
max(1,min(mstu(52),1+int(mstu(52)*0.5*(
phi/paru(1)+1.))))
36 ietph=mstu(52)*ieta+iphi
40 IF(ietph.EQ.
k(ic,3))
THEN
46 IF(nc.GE.mstu(4)-mstu(32)-5)
THEN
47 CALL
luerrm(11,
'(LUCELL:) no more memory left in LUJETS')
55 p(nc,1)=(paru(51)/mstu(51))*(2*ieta-1-mstu(51))
56 p(nc,2)=(paru(1)/mstu(52))*(2*iphi-1-mstu(52))
61 IF(mstu(53).GE.1)
THEN
64 IF(mstu(53).EQ.2) pei=
p(ic,5)/cosh(
p(ic,1))
65 120 pef=pei+paru(55)*sqrt(-2.*
log(
max(1
e-10,
rlu(0)))*pei)*
67 IF(pef.LT.0..OR.pef.GT.paru(56)*pei) goto 120
69 130
IF(mstu(53).EQ.2)
p(ic,5)=pef*cosh(
p(ic,1))
76 IF(
k(ic,5).NE.2) goto 150
77 IF(
p(ic,5).LE.etmax) goto 150
83 IF(etmax.LT.paru(52)) goto 210
84 IF(nj.GE.mstu(4)-mstu(32)-5)
THEN
85 CALL
luerrm(11,
'(LUCELL:) no more memory left in LUJETS')
101 IF(
k(ic,5).EQ.0) goto 160
102 IF(abs(
p(ic,1)-
eta).GT.paru(54)) goto 160
103 dphia=abs(
p(ic,2)-
phi)
104 IF(dphia.GT.paru(54).AND.dphia.LT.paru(2)-paru(54)) goto 160
106 IF(dphia.GT.paru(1)) phic=phic+
sign(paru(2),
phi)
107 IF((
p(ic,1)-
eta)**2+(phic-
phi)**2.GT.paru(54)**2) goto 160
109 k(nj,4)=
k(nj,4)+
k(ic,4)
110 p(nj,3)=
p(nj,3)+
p(ic,5)*
p(ic,1)
111 p(nj,4)=
p(nj,4)+
p(ic,5)*phic
112 p(nj,5)=
p(nj,5)+
p(ic,5)
116 IF(
p(nj,5).LT.paru(53))
THEN
119 170
IF(
k(ic,5).LT.0)
k(ic,5)=-
k(ic,5)
120 ELSEIF(mstu(54).LE.2)
THEN
121 p(nj,3)=
p(nj,3)/
p(nj,5)
122 p(nj,4)=
p(nj,4)/
p(nj,5)
123 IF(abs(
p(nj,4)).GT.paru(1))
p(nj,4)=
p(nj,4)-
sign(paru(2),
126 180
IF(
k(ic,1).LT.0)
k(ic,1)=0
131 IF(
k(ic,5).GE.0) goto 200
132 p(nj,1)=
p(nj,1)+
p(ic,5)*cos(
p(ic,2))
133 p(nj,2)=
p(nj,2)+
p(ic,5)*sin(
p(ic,2))
134 p(nj,3)=
p(nj,3)+
p(ic,5)*sinh(
p(ic,1))
135 p(nj,4)=
p(nj,4)+
p(ic,5)*cosh(
p(ic,1))
145 IF(
k(ij,5).EQ.0) goto 220
146 IF(
p(ij,5).LT.etmax) goto 220
162 IF(mstu(54).EQ.2)
THEN
170 ELSEIF(mstu(54).GE.3)
THEN
172 250
p(
i,5)=sqrt(
max(0.,
p(
i,4)**2-
p(
i,1)**2-
p(
i,2)**2-
p(
i,3)**2))
179 IF(mstu(43).LE.1) mstu(3)=njet
180 IF(mstu(43).GE.2)
n=
n+njet