7 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
9 common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)
22 IF((mstj(93).EQ.1.OR.mstj(93).EQ.2).AND.kfa.LE.10)
THEN
27 ELSEIF(kfa.LE.100.OR.kc.LE.80.OR.kc.GT.100)
THEN
36 kflr=mod(kfa/10000,10)
42 IF(kfla.EQ.0.AND.kflr.EQ.0.AND.kfls.LE.3)
THEN
43 IF(kfls.EQ.1) pmspl=-3./(pmb*pmc)
44 IF(kfls.GE.3) pmspl=1./(pmb*pmc)
45 ulmass=parf(111)+pmb+pmc+parf(113)*parf(101)**2*pmspl
46 ELSEIF(kfla.EQ.0)
THEN
51 ulmass=parf(113+kmul)+pmb+pmc
52 ELSEIF(kflc.EQ.0)
THEN
53 IF(kfls.EQ.1) pmspl=-3./(pma*pmb)
54 IF(kfls.EQ.3) pmspl=1./(pma*pmb)
55 ulmass=2.*parf(112)/3.+pma+pmb+parf(114)*parf(101)**2*pmspl
56 IF(mstj(93).EQ.1)
ulmass=pma+pmb
60 IF(kfls.EQ.2.AND.kfla.EQ.kflb)
THEN
61 pmspl=1./(pma*pmb)-2./(pma*pmc)-2./(pmb*pmc)
62 ELSEIF(kfls.EQ.2.AND.kflb.GE.kflc)
THEN
63 pmspl=-2./(pma*pmb)-2./(pma*pmc)+1./(pmb*pmc)
64 ELSEIF(kfls.EQ.2)
THEN
67 pmspl=1./(pma*pmb)+1./(pma*pmc)+1./(pmb*pmc)
69 ulmass=parf(112)+pma+pmb+pmc+parf(114)*parf(101)**2*pmspl
75 IF(mstj(24).GE.1.AND.pmas(kc,2).GT.1
e-4)
THEN
76 IF(mstj(24).EQ.1.OR.(mstj(24).EQ.2.AND.kfa.GT.100))
THEN
78 & atan(2.*pmas(kc,3)/pmas(kc,2)))
81 pmlow=atan((
max(0.,pm0-pmas(kc,3))**2-pm0**2)/
83 pmupp=atan((pm0+pmas(kc,3))**2-pm0**2)/(pm0*pmas(kc,2))
84 ulmass=sqrt(
max(0.,pm0**2+pm0*pmas(kc,2)*tan(pmlow+
85 & (pmupp-pmlow)*
rlu(0))))