9 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
11 common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)
13 common/ludat3/mdcy(500,3),mdme(2000,2),brat(2000),kfdp(2000,5)
15 common/pyhisubs/msel,msub(200),kfin(2,-40:40),ckin(200)
21 common/pyhiint2/iset(200),kfpr(200,2),coef(200,20),icol(40,4,2)
23 common/pyhiint4/widp(21:40,0:40),wide(21:40,0:40),wids(21:40,3)
25 common/pyhiint6/proc(0:200)
28 dimension wdtp(0:40),wdte(0:40,0:5)
40 wfac=aem/(24.*xw)*wmas
42 wids(24,1)=((wdte(0,1)+wdte(0,2))*(wdte(0,1)+wdte(0,3))+
43 &(wdte(0,1)+wdte(0,2)+wdte(0,1)+wdte(0,3))*(wdte(0,4)+wdte(0,5))+
44 &2.*wdte(0,4)*wdte(0,5))/wdtp(0)**2
45 wids(24,2)=(wdte(0,1)+wdte(0,2)+wdte(0,4))/wdtp(0)
46 wids(24,3)=(wdte(0,1)+wdte(0,3)+wdte(0,4))/wdtp(0)
48 widp(24,
i)=wfac*wdtp(
i)
49 110 wide(24,
i)=wfac*wdte(
i,0)
53 hcfac=aem/(8.*xw)*(hcmas/wmas)**2*hcmas
55 wids(37,1)=((wdte(0,1)+wdte(0,2))*(wdte(0,1)+wdte(0,3))+
56 &(wdte(0,1)+wdte(0,2)+wdte(0,1)+wdte(0,3))*(wdte(0,4)+wdte(0,5))+
57 &2.*wdte(0,4)*wdte(0,5))/wdtp(0)**2
58 wids(37,2)=(wdte(0,1)+wdte(0,2)+wdte(0,4))/wdtp(0)
59 wids(37,3)=(wdte(0,1)+wdte(0,3)+wdte(0,4))/wdtp(0)
61 widp(37,
i)=hcfac*wdtp(
i)
62 120 wide(37,
i)=hcfac*wdte(
i,0)
66 zfac=aem/(48.*xw*(1.-xw))*zmas
68 wids(23,1)=((wdte(0,1)+wdte(0,2))**2+
69 &2.*(wdte(0,1)+wdte(0,2))*(wdte(0,4)+wdte(0,5))+
70 &2.*wdte(0,4)*wdte(0,5))/wdtp(0)**2
71 wids(23,2)=(wdte(0,1)+wdte(0,2)+wdte(0,4))/wdtp(0)
74 widp(23,
i)=zfac*wdtp(
i)
75 130 wide(23,
i)=zfac*wdte(
i,0)
79 hfac=aem/(8.*xw)*(hmas/wmas)**2*hmas
81 wids(25,1)=((wdte(0,1)+wdte(0,2))**2+
82 &2.*(wdte(0,1)+wdte(0,2))*(wdte(0,4)+wdte(0,5))+
83 &2.*wdte(0,4)*wdte(0,5))/wdtp(0)**2
84 wids(25,2)=(wdte(0,1)+wdte(0,2)+wdte(0,4))/wdtp(0)
87 widp(25,
i)=hfac*wdtp(
i)
88 140 wide(25,
i)=hfac*wdte(
i,0)
92 zpfac=aem/(48.*xw*(1.-xw))*zpmas
94 wids(32,1)=((wdte(0,1)+wdte(0,2)+wdte(0,3))**2+
95 &2.*(wdte(0,1)+wdte(0,2))*(wdte(0,4)+wdte(0,5))+
96 &2.*wdte(0,4)*wdte(0,5))/wdtp(0)**2
97 wids(32,2)=(wdte(0,1)+wdte(0,2)+wdte(0,4))/wdtp(0)
100 widp(32,
i)=zpfac*wdtp(
i)
101 150 wide(32,
i)=zpfac*wdte(
i,0)
105 rfac=0.08*rmas/((
mstp(1)-1)*(1.+6.*(1.+
ulalps(rmas**2)/paru(1))))
107 wids(40,1)=((wdte(0,1)+wdte(0,2))*(wdte(0,1)+wdte(0,3))+
108 &(wdte(0,1)+wdte(0,2)+wdte(0,1)+wdte(0,3))*(wdte(0,4)+wdte(0,5))+
109 &2.*wdte(0,4)*wdte(0,5))/wdtp(0)**2
110 wids(40,2)=(wdte(0,1)+wdte(0,2)+wdte(0,4))/wdtp(0)
111 wids(40,3)=(wdte(0,1)+wdte(0,3)+wdte(0,4))/wdtp(0)
113 widp(40,
i)=wfac*wdtp(
i)
114 160 wide(40,
i)=wfac*wdte(
i,0)
118 DO 170
i=1,min(8,mdcy(21,3))
120 IF(mdme(idc,1).LE.0) goto 170
135 pmas(kc,2)=widp(kc,0)
136 pmas(kc,3)=min(0.9*pmas(kc,1),10.*pmas(kc,2))
137 DO 180
j=1,mdcy(kc,3)
139 brat(idc)=wide(kc,
j)/wide(kc,0)
143 IF(
mstp(43).EQ.1)
THEN
144 proc(1)=
'f + fb -> gamma*'
145 ELSEIF(
mstp(43).EQ.2)
THEN
146 proc(1)=
'f + fb -> Z0'
147 ELSEIF(
mstp(43).EQ.3)
THEN
148 proc(1)=
'f + fb -> gamma*/Z0'
152 IF(
mstp(44).EQ.1)
THEN
153 proc(141)=
'f + fb -> gamma*'
154 ELSEIF(
mstp(44).EQ.2)
THEN
155 proc(141)=
'f + fb -> Z0'
156 ELSEIF(
mstp(44).EQ.3)
THEN
157 proc(141)=
'f + fb -> Z''0'
158 ELSEIF(
mstp(44).EQ.4)
THEN
159 proc(141)=
'f + fb -> gamma*/Z0'
160 ELSEIF(
mstp(44).EQ.5)
THEN
161 proc(141)=
'f + fb -> gamma*/Z''0'
162 ELSEIF(
mstp(44).EQ.6)
THEN
163 proc(141)=
'f + fb -> Z0/Z''0'
164 ELSEIF(
mstp(44).EQ.7)
THEN
165 proc(141)=
'f + fb -> gamma*/Z0/Z''0'