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 IF(mstu(12).GE.1) CALL
lulist(0)
17 IF(kfl.LT.0.OR.kfl.GT.8)
THEN
18 CALL
luerrm(16,
'(LUONIA:) called with unknown flavour code')
19 IF(mstu(21).GE.1)
RETURN
21 IF(ecm.LT.parj(127)+2.02*parf(101))
THEN
22 CALL
luerrm(16,
'(LUONIA:) called with too small CM energy')
23 IF(mstu(21).GE.1)
RETURN
28 IF(mstj(115).GE.2)
THEN
30 CALL
lu1ent(nc-1,11,0.5*ecm,0.,0.)
32 CALL
lu1ent(nc,-11,0.5*ecm,paru(1),0.)
36 IF(mstj(115).GE.3.AND.kflc.GE.5)
THEN
42 CALL
lu1ent(nc,kf,ecm,0.,0.)
53 IF(x3.GE.1..OR.((1.-x1)/(x2*x3))**2+((1.-x2)/(x1*x3))**2+
54 &((1.-x3)/(x1*x2))**2.LE.2.*
rlu(0)) goto 100
57 IF(mstj(101).LE.4) CALL
lu3ent(nc+1,21,21,21,ecm,x1,x3)
58 IF(mstj(101).GE.5) CALL
lu3ent(-(nc+1),21,21,21,ecm,x1,x3)
62 IF(mstj(108).EQ.2.AND.(mstj(101).EQ.0.OR.mstj(101).EQ.1))
65 IF(mstu(111).EQ.2) paru(112)=parj(122)
67 IF(kflc.NE.0) qf=kchg(kflc,1)/3.
68 rgam=7.2*qf**2*paru(101)/
ulalps(ecm**2)
71 IF(
rlu(0).GT.rgam/(1.+rgam))
THEN
72 IF(1.-
max(x1,x2,x3).LE.
max((parj(126)/ecm)**2,parj(125)))
74 IF(njet.EQ.2.AND.mstj(101).LE.4) CALL
lu2ent(nc+1,21,21,ecm)
75 IF(njet.EQ.2.AND.mstj(101).GE.5) CALL
lu2ent(-(nc+1),21,21,ecm)
79 IF(ecmc.LT.2.*parj(127)) goto 100
84 IF(mstj(101).GE.5)
k(nc+2,4)=mstu(5)*(nc+3)
85 IF(mstj(101).GE.5)
k(nc+2,5)=mstu(5)*(nc+3)
86 IF(mstj(101).GE.5)
k(nc+3,4)=mstu(5)*(nc+2)
87 IF(mstj(101).GE.5)
k(nc+3,5)=mstu(5)*(nc+2)
89 IF(ecmc.LT.4.*parj(127))
THEN
93 CALL
lu1ent(nc+2,83,0.5*(x2+x3)*ecm,paru(1),0.)
99 110
k(ip,3)=
k(ip,3)+(mstj(115)/2)+(kflc/5)*(mstj(115)/3)*(nc-1)
102 IF(mstj(106).EQ.1)
THEN
104 hf1=1.-parj(131)*parj(132)
106 ct13=(x1*x3-2.*x1-2.*x3+2.)/(x1*x3)
107 st13=sqrt(1.-ct13**2)
108 sigl=0.5*x3**2*((1.-x2)**2+(1.-x3)**2)*st13**2
109 sigu=(x1*(1.-x1))**2+(x2*(1.-x2))**2+(x3*(1.-x3))**2-sigl
111 sigi=(sigl*ct13/st13+0.5*x1*x3*(1.-x2)**2*st13)/sq2
112 sigmax=(2.*hf1+hf3)*abs(sigu)+2.*(hf1+hf3)*abs(sigl)+2.*(hf1+
113 & 2.*hf3)*abs(sigt)+2.*sq2*(hf1+2.*hf3)*abs(sigi)
116 120 chi=paru(2)*
rlu(0)
125 c2phi=cos(2.*(
phi-parj(134)))
126 s2phi=sin(2.*(
phi-parj(134)))
127 sig=((1.+cthe**2)*hf1+sthe**2*c2phi*hf3)*sigu+2.*(sthe**2*hf1-
128 & sthe**2*c2phi*hf3)*sigl+2.*(sthe**2*c2chi*hf1+((1.+cthe**2)*
129 & c2chi*c2phi-2.*cthe*s2chi*s2phi)*hf3)*sigt-2.*sq2*(2.*sthe*cthe*
130 & cchi*hf1-2.*sthe*(cthe*cchi*c2phi-schi*s2phi)*hf3)*sigi
131 IF(sig.LT.sigmax*
rlu(0)) goto 120
132 CALL ludbrb(nc+1,
n,0.,chi,0d0,0d0,0d0)
133 CALL ludbrb(nc+1,
n,the,
phi,0d0,0d0,0d0)
137 IF(mstj(101).GE.5.AND.njet.GE.2)
THEN
138 CALL
lushow(nc+mk+1,-njet,ecmc)
140 IF(mstj(105).EQ.-1) mstj(14)=0
141 IF(mstj(105).GE.0) mstu(28)=0
144 IF(mstj(105).GE.0.AND.mstu(28).NE.0) goto 100
148 IF(mstj(105).EQ.1) CALL
luexec