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 sm(3,3),sax(3),
ps(3,5)
27 IF(
k(
i,1).LE.0.OR.
k(
i,1).GT.10) goto 150
28 IF(mstu(41).GE.2)
THEN
30 IF(kc.EQ.0.OR.kc.EQ.12.OR.kc.EQ.14.OR.kc.EQ.16.OR.
32 IF(mstu(41).GE.3.AND.kchg(kc,2).EQ.0.AND.
luchge(
k(
i,2)).EQ.0)
35 IF(
n+
np+1.GE.mstu(4)-mstu(32)-5)
THEN
36 CALL
luerrm(11,
'(LUJMAS:) no more memory left in LUJETS')
44 IF(mstu(42).EQ.0)
p(
n+
np,5)=0.
45 IF(mstu(42).EQ.1.AND.
k(
i,2).NE.22)
p(
n+
np,5)=pmas(101,1)
51 130 sm(j1,j2)=sm(j1,j2)+
p(
i,j1)*
p(
i,j2)
52 pss=pss+(
p(
i,1)**2+
p(
i,2)**2+
p(
i,3)**2)
59 CALL
luerrm(8,
'(LUJMAS:) too few particles for analysis')
64 paru(61)=sqrt(
max(0.,
ps(3,4)**2-
ps(3,1)**2-
ps(3,2)**2-
ps(3,3)**2))
69 160 sm(j1,j2)=sm(j1,j2)/pss
70 sq=(sm(1,1)*sm(2,2)+sm(1,1)*sm(3,3)+sm(2,2)*sm(3,3)-sm(1,2)**2-
71 &sm(1,3)**2-sm(2,3)**2)/3.-1./9.
72 sr=-0.5*(sq+1./9.+sm(1,1)*sm(2,3)**2+sm(2,2)*sm(1,3)**2+sm(3,3)*
73 &sm(1,2)**2-sm(1,1)*sm(2,2)*sm(3,3))+sm(1,2)*sm(1,3)*sm(2,3)+1./27.
74 sp=cos(acos(
max(min(sr/sqrt(-sq**3),1.),-1.))/3.)
75 sma=1./3.+sqrt(-sq)*
max(2.*
sp,sqrt(3.*(1.-
sp**2))-
sp)
79 sm(j1,j1)=sm(j1,j1)-sma
81 170 sm(j2,j1)=sm(j1,j2)
85 IF(abs(sm(j1,j2)).LE.smax) goto 180
93 rl=sm(j1,jb)/sm(ja,jb)
95 sm(j1,j2)=sm(j1,j2)-rl*sm(ja,j2)
96 IF(abs(sm(j1,j2)).LE.smax) goto 190
101 jb2=jb+2-3*((jb+1)/3)
104 sax(jb)=-(sm(ja,jb1)*sax(jb1)+sm(ja,jb2)*sax(jb2))/sm(ja,jb)
108 psax=
p(
i,1)*sax(1)+
p(
i,2)*sax(2)+
p(
i,3)*sax(3)
114 pms=(
ps(1,4)**2-
ps(1,1)**2-
ps(1,2)**2-
ps(1,3)**2)+
115 &(
ps(2,4)**2-
ps(2,1)**2-
ps(2,2)**2-
ps(2,3)**2)
124 IF(
k(
i,3).EQ.1) pmdi=2.*(
p(
i,5)**2-pps)
125 IF(
k(
i,3).EQ.2) pmdi=2.*(
p(
i,5)**2+pps)
133 IF(pmd.LT.-paru(48)*pms)
THEN
146 ps(1,5)=sqrt(
max(0.,
ps(1,4)**2-
ps(1,1)**2-
ps(1,2)**2-
ps(1,3)**2))
147 ps(2,5)=sqrt(
max(0.,
ps(2,4)**2-
ps(2,1)**2-
ps(2,2)**2-
ps(2,3)**2))
149 pml=min(
ps(1,5),
ps(2,5))