12 IMPLICIT DOUBLE PRECISION(a-
h, o-
z)
16 parameter(ksusy1=1000000,ksusy2=2000000,ktechn=3000000,
17 &kexcit=4000000,kdimen=5000000)
19 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
20 common/
pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
21 common/
pydat3/mdcy(500,3),mdme(8000,2),brat(8000),kfdp(8000,5)
24 common/
pyint2/iset(500),kfpr(500,2),coef(500,20),icol(40,4,2)
25 common/
pyint3/xsfx(2,-40:40),isig(1000,3),sigh(1000)
26 common/
pyint4/mwid(500),wids(500,5)
27 common/
pyint7/sigt(0:6,0:6,0:5)
28 common/pysgcm/isub,isubsv,mmin1,mmax1,mmin2,mmax2,mmina,mmaxa,
29 &kfac(2,-40:40),comfac,fack,faca,sh,th,uh,sh2,th2,uh2,sqm3,sqm4,
30 &shr,sqpth,taup,be34,cth,
x(2),sqmz,sqmw,gmmz,gmmw,
31 &aem,as,xw,xw1,xwc,xwv,poll,polr,polll,polrr
35 dimension wdtp(0:400),wdte(0:400,0:5)
42 facggf=comfac*aem**2*2d0*(sh2+uh2)/th2
43 facgzf=comfac*aem**2*xwc*4d0*sh2/(th*(th-sqmz))
44 faczzf=comfac*(aem*xwc)**2*2d0*sh2/(th-sqmz)**2
45 facwwf=comfac*(0.5d0*aem/xw)**2*sh2/(th-sqmw)**2
47 IF(
i.EQ.0.OR.kfac(1,
i).EQ.0) goto 110
50 IF(
j.EQ.0.OR.kfac(2,
j).EQ.0) goto 100
53 ei=kchg(ia,1)*isign(1,
i)/3d0
54 ai=
sign(1d0,kchg(ia,1)+0.5d0)*isign(1,
i)
56 ej=kchg(ja,1)*isign(1,
j)/3d0
57 aj=
sign(1d0,kchg(ja,1)+0.5d0)*isign(1,
j)
61 IF(
mstp(21).GE.1.AND.
mstp(21).LE.4)
THEN
62 IF(
mstp(21).EQ.1.OR.
mstp(21).EQ.4)
THEN
63 facncf=facggf*
ei**2*ej**2+facgzf*
ei*ej*
64 & (vi*vj*(1d0+uh2/sh2)+ai*aj*epsij*(1d0-uh2/sh2))+
65 & faczzf*((vi**2+ai**2)*(vj**2+aj**2)*(1d0+uh2/sh2)+
66 & 4d0*vi*vj*ai*aj*epsij*(1d0-uh2/sh2))
67 ELSEIF(
mstp(21).EQ.2)
THEN
68 facncf=facggf*
ei**2*ej**2
70 facncf=faczzf*((vi**2+ai**2)*(vj**2+aj**2)*
71 & (1d0+uh2/sh2)+4d0*vi*vj*ai*aj*epsij*(1d0-uh2/sh2))
74 IF(ia.GT.10.AND.mod(ia,2).EQ.0) facncf=2d0*facncf
75 IF(ja.GT.10.AND.mod(ja,2).EQ.0) facncf=2d0*facncf
83 IF((
mstp(21).EQ.1.OR.
mstp(21).EQ.5).AND.ai*aj.LT.0d0)
THEN
85 IF(epsij.LT.0d0) facccf=facccf*uh2/sh2
86 IF(ia.GT.10.AND.mod(ia,2).EQ.0) facccf=2d0*facccf
87 IF(ja.GT.10.AND.mod(ja,2).EQ.0) facccf=2d0*facccf
97 ELSEIF(isub.EQ.11)
THEN
99 facqq1=comfac*as**2*4d0/9d0*(sh2+uh2)/th2
100 facqqb=comfac*as**2*4d0/9d0*((sh2+uh2)/th2*faca-
101 &
mstp(34)*2d0/3d0*uh2/(sh*th))
102 facqq2=comfac*as**2*4d0/9d0*((sh2+th2)/uh2-
103 &
mstp(34)*2d0/3d0*sh2/(th*uh))
106 IF(
i.EQ.0.OR.ia.GT.
mstp(58).OR.kfac(1,
i).EQ.0) goto 130
109 IF(
j.EQ.0.OR.ja.GT.
mstp(58).OR.kfac(2,
j).EQ.0) goto 120
115 IF(
i.EQ.-
j) sigh(nchn)=facqqb
117 sigh(nchn)=0.5d0*sigh(nchn)
122 sigh(nchn)=0.5d0*facqq2
127 ELSEIF(isub.EQ.12)
THEN
129 CALL
pywidt(21,sh,wdtp,wdte)
130 facqqb=comfac*as**2*4d0/9d0*(th2+uh2)/sh2*
131 & (wdte(0,1)+wdte(0,2)+wdte(0,4))
133 IF(
i.EQ.0.OR.iabs(
i).GT.
mstp(58).OR.
134 & kfac(1,
i)*kfac(2,-
i).EQ.0) goto 140
142 ELSEIF(isub.EQ.13)
THEN
144 facgg1=comfac*as**2*32d0/27d0*(uh/th-(2d0+
mstp(34)*1d0/4d0)*
146 facgg2=comfac*as**2*32d0/27d0*(th/uh-(2d0+
mstp(34)*1d0/4d0)*
149 IF(
i.EQ.0.OR.iabs(
i).GT.
mstp(58).OR.
150 & kfac(1,
i)*kfac(2,-
i).EQ.0) goto 150
155 sigh(nchn)=0.5d0*facgg1
160 sigh(nchn)=0.5d0*facgg2
163 ELSEIF(isub.EQ.14)
THEN
165 facgg=comfac*as*aem*8d0/9d0*(th2+uh2)/(th*uh)
167 IF(
i.EQ.0.OR.iabs(
i).GT.
mstp(58).OR.
168 & kfac(1,
i)*kfac(2,-
i).EQ.0) goto 160
169 ei=kchg(iabs(
i),1)/3d0
174 sigh(nchn)=facgg*
ei**2
177 ELSEIF(isub.EQ.18)
THEN
179 facgg=comfac*aem**2*2d0*(th2+uh2)/(th*uh)
181 IF(
i.EQ.0.OR.kfac(1,
i)*kfac(2,-
i).EQ.0) goto 170
182 ei=kchg(iabs(
i),1)/3d0
184 IF(iabs(
i).LE.10) fcoi=faca/3d0
189 sigh(nchn)=0.5d0*facgg*fcoi*
ei**4
193 ELSEIF(isub.LE.40)
THEN
196 facqg1=comfac*as**2*4d0/9d0*((2d0+
mstp(34)*1d0/4d0)*uh2/th2-
198 facqg2=comfac*as**2*4d0/9d0*((2d0+
mstp(34)*1d0/4d0)*sh2/th2-
201 IF(
i.EQ.0.OR.iabs(
i).GT.10) goto 190
203 IF(isde.EQ.1.AND.kfac(1,
i)*kfac(2,21).EQ.0) goto 180
204 IF(isde.EQ.2.AND.kfac(1,21)*kfac(2,
i).EQ.0) goto 180
218 ELSEIF(isub.EQ.29)
THEN
220 fgq=comfac*faca*as*aem*1d0/3d0*(sh2+uh2)/(-sh*uh)
222 IF(
i.EQ.0.OR.iabs(
i).GT.
mstp(58)) goto 210
223 ei=kchg(iabs(
i),1)/3d0
226 IF(isde.EQ.1.AND.kfac(1,
i)*kfac(2,21).EQ.0) goto 200
227 IF(isde.EQ.2.AND.kfac(1,21)*kfac(2,
i).EQ.0) goto 200
236 ELSEIF(isub.EQ.33)
THEN
238 fgq=comfac*as*aem*8d0/3d0*(sh2+uh2)/(-sh*uh)
240 IF(
i.EQ.0.OR.iabs(
i).GT.
mstp(58)) goto 230
241 ei=kchg(iabs(
i),1)/3d0
244 IF(isde.EQ.1.AND.kfac(1,
i)*kfac(2,22).EQ.0) goto 220
245 IF(isde.EQ.2.AND.kfac(1,22)*kfac(2,
i).EQ.0) goto 220
254 ELSEIF(isub.EQ.34)
THEN
256 fgq=comfac*aem**2*2d0*(sh2+uh2)/(-sh*uh)
259 ei=kchg(iabs(
i),1)/3d0
262 IF(isde.EQ.1.AND.kfac(1,
i)*kfac(2,22).EQ.0) goto 240
263 IF(isde.EQ.2.AND.kfac(1,22)*kfac(2,
i).EQ.0) goto 240
273 ELSEIF(isub.LE.80)
THEN
276 IF(kfac(1,21)*kfac(2,21).EQ.0) goto 270
280 IF(mdme(idc0+1,1).GE.1) flavwt=flavwt+
281 & sqrt(
max(0d0,1d0-4d0*pmas(1,1)**2/sh))
282 IF(mdme(idc0+2,1).GE.1) flavwt=flavwt+
283 & sqrt(
max(0d0,1d0-4d0*pmas(2,1)**2/sh))
284 IF(mdme(idc0+3,1).GE.1) flavwt=flavwt+
285 & sqrt(
max(0d0,1d0-4d0*pmas(3,1)**2/sh))
286 facqq1=comfac*as**2*1d0/6d0*(uh/th-(2d0+
mstp(34)*1d0/4d0)*
287 & uh2/sh2)*flavwt*faca
288 facqq2=comfac*as**2*1d0/6d0*(th/uh-(2d0+
mstp(34)*1d0/4d0)*
289 & th2/sh2)*flavwt*faca
303 sqmavg=pmas(ifl,1)**2
304 IF(mdme(idc0+ifl,1).GE.1.AND.sh.GT.4.04d0*sqmavg)
THEN
305 be34=sqrt(1d0-4d0*sqmavg/sh)
306 thq=-0.5d0*sh*(1d0-be34*cth)
307 uhq=-0.5d0*sh*(1d0+be34*cth)
308 thuhq=thq*uhq-sqmavg*sh
309 IF(
mstp(34).EQ.0)
THEN
310 facqq1=uhq/thq-2d0*uhq**2/sh2+4d0*(sqmavg/sh)*thuhq/thq**2
311 facqq2=thq/uhq-2d0*thq**2/sh2+4d0*(sqmavg/sh)*thuhq/uhq**2
313 facqq1=uhq/thq-2.25d0*uhq**2/sh2+4.5d0*(sqmavg/sh)*thuhq/
314 & thq**2+0.5d0*sqmavg*(thq+sqmavg)/thq**2-sqmavg**2/(sh*thq)
315 facqq2=thq/uhq-2.25d0*thq**2/sh2+4.5d0*(sqmavg/sh)*thuhq/
316 & uhq**2+0.5d0*sqmavg*(uhq+sqmavg)/uhq**2-sqmavg**2/(sh*uhq)
318 facqq1=comfac*faca*as**2*(1d0/6d0)*facqq1*be34
319 facqq2=comfac*faca*as**2*(1d0/6d0)*facqq2*be34
323 isig(nchn,3)=1+2*(ifl-3)
328 isig(nchn,3)=2+2*(ifl-3)
334 ELSEIF(isub.EQ.54)
THEN
336 CALL
pywidt(21,sh,wdtp,wdte)
338 DO 280
i=1,min(8,mdcy(21,3))
340 wdtesu=wdtesu+ef**2*(wdte(
i,1)+wdte(
i,2)+wdte(
i,3)+
343 facqq=comfac*aem*as*wdtesu*(th2+uh2)/(th*uh)
344 IF(kfac(1,21)*kfac(2,22).NE.0)
THEN
351 IF(kfac(1,22)*kfac(2,21).NE.0)
THEN
359 ELSEIF(isub.EQ.58)
THEN
361 CALL
pywidt(22,sh,wdtp,wdte)
363 DO 290
i=1,min(12,mdcy(22,3))
364 IF(
i.LE.8) ef= kchg(
i,1)/3d0
365 IF(
i.GE.9) ef= kchg(9+2*(
i-8),1)/3d0
366 wdtesu=wdtesu+ef**2*(wdte(
i,1)+wdte(
i,2)+wdte(
i,3)+
369 facff=comfac*aem**2*wdtesu*2d0*(th2+uh2)/(th*uh)
370 IF(kfac(1,22)*kfac(2,22).NE.0)
THEN
378 ELSEIF(isub.EQ.68)
THEN
380 IF(kfac(1,21)*kfac(2,21).EQ.0) goto 300
381 facgg1=comfac*as**2*9d0/4d0*(sh2/th2+2d0*sh/th+3d0+2d0*th/sh+
383 facgg2=comfac*as**2*9d0/4d0*(uh2/sh2+2d0*uh/sh+3d0+2d0*sh/uh+
385 facgg3=comfac*as**2*9d0/4d0*(th2/uh2+2d0*th/uh+3d0+2d0*uh/th+
391 sigh(nchn)=0.5d0*facgg1
396 sigh(nchn)=0.5d0*facgg2
401 sigh(nchn)=0.5d0*facgg3
404 ELSEIF(isub.EQ.80)
THEN
406 fqpi=comfac*(2d0*aem/9d0)*(-sh/th)*(1d0/sh2+1d0/th2)
408 q2fpsh=0.55d0/
log(
max(2d0,2d0*sh))
409 delsh=uh*sqrt(assh*q2fpsh)
411 q2fpuh=0.55d0/
log(
max(2d0,-2d0*uh))
412 deluh=sh*sqrt(asuh*q2fpuh)
413 DO 320
i=
max(-2,mmina),min(2,mmaxa)
415 ei=kchg(iabs(
i),1)/3d0
418 IF(isde.EQ.1.AND.kfac(1,
i)*kfac(2,22).EQ.0) goto 310
419 IF(isde.EQ.2.AND.kfac(1,22)*kfac(2,
i).EQ.0) goto 310
424 sigh(nchn)=fqpi*(
ei*delsh+ej*deluh)**2
429 ELSEIF(isub.LE.100)
THEN
432 sigs=
vint(315)*
vint(316)*sigt(0,0,1)
434 ELSEIF(isub.EQ.92)
THEN
436 sigs=
vint(315)*
vint(316)*sigt(0,0,2)
438 ELSEIF(isub.EQ.93)
THEN
440 sigs=
vint(315)*
vint(316)*sigt(0,0,3)
442 ELSEIF(isub.EQ.94)
THEN
444 sigs=
vint(315)*
vint(316)*sigt(0,0,4)
446 ELSEIF(isub.EQ.95)
THEN
448 sigs=
vint(315)*
vint(316)*sigt(0,0,5)
450 ELSEIF(isub.EQ.96)
THEN
452 CALL
pywidt(21,sh,wdtp,wdte)
455 facqq1=comfac*as**2*4d0/9d0*(sh2+uh2)/th2
456 facqqb=comfac*as**2*4d0/9d0*((sh2+uh2)/th2*faca-
457 &
mstp(34)*2d0/3d0*uh2/(sh*th))
458 facqq2=comfac*as**2*4d0/9d0*(sh2+th2)/uh2
459 facqqi=-comfac*as**2*4d0/9d0*
mstp(34)*2d0/3d0*sh2/(th*uh)
460 ratqqi=(facqq1+facqq2+facqqi)/(facqq1+facqq2)
470 IF(
i.EQ.-
j) sigh(nchn)=facqqb
472 sigh(nchn)=0.5d0*facqq1*ratqqi
477 sigh(nchn)=0.5d0*facqq2*ratqqi
483 facqqb=comfac*as**2*4d0/9d0*(th2+uh2)/sh2*
484 & (wdte(0,1)+wdte(0,2)+wdte(0,3)+wdte(0,4))
485 facgg1=comfac*as**2*32d0/27d0*(uh/th-(2d0+
mstp(34)*1d0/4d0)*
487 facgg2=comfac*as**2*32d0/27d0*(th/uh-(2d0+
mstp(34)*1d0/4d0)*
500 sigh(nchn)=0.5d0*facgg1
505 sigh(nchn)=0.5d0*facgg2
509 facqg1=comfac*as**2*4d0/9d0*((2d0+
mstp(34)*1d0/4d0)*uh2/th2-
511 facqg2=comfac*as**2*4d0/9d0*((2d0+
mstp(34)*1d0/4d0)*sh2/th2-
532 IF(mdme(idc0+1,1).GE.1) flavwt=flavwt+
533 & sqrt(
max(0d0,1d0-4d0*pmas(1,1)**2/sh))
534 IF(mdme(idc0+2,1).GE.1) flavwt=flavwt+
535 & sqrt(
max(0d0,1d0-4d0*pmas(2,1)**2/sh))
536 IF(mdme(idc0+3,1).GE.1) flavwt=flavwt+
537 & sqrt(
max(0d0,1d0-4d0*pmas(3,1)**2/sh))
538 facqq1=comfac*as**2*1d0/6d0*(uh/th-(2d0+
mstp(34)*1d0/4d0)*
539 & uh2/sh2)*flavwt*faca
540 facqq2=comfac*as**2*1d0/6d0*(th/uh-(2d0+
mstp(34)*1d0/4d0)*
541 & th2/sh2)*flavwt*faca
556 sqmavg=pmas(ifl,1)**2
557 IF(mdme(idc0+ifl,1).GE.1.AND.sh.GT.4.04d0*sqmavg)
THEN
558 be34=sqrt(1d0-4d0*sqmavg/sh)
559 thq=-0.5d0*sh*(1d0-be34*cth)
560 uhq=-0.5d0*sh*(1d0+be34*cth)
561 thuhq=thq*uhq-sqmavg*sh
562 IF(
mstp(34).EQ.0)
THEN
563 facqq1=uhq/thq-2d0*uhq**2/sh2+4d0*(sqmavg/sh)*thuhq/thq**2
564 facqq2=thq/uhq-2d0*thq**2/sh2+4d0*(sqmavg/sh)*thuhq/uhq**2
566 facqq1=uhq/thq-2.25d0*uhq**2/sh2+4.5d0*(sqmavg/sh)*thuhq/
567 & thq**2+0.5d0*sqmavg*(thq+sqmavg)/thq**2-sqmavg**2/(sh*thq)
568 facqq2=thq/uhq-2.25d0*thq**2/sh2+4.5d0*(sqmavg/sh)*thuhq/
569 & uhq**2+0.5d0*sqmavg*(uhq+sqmavg)/uhq**2-sqmavg**2/(sh*uhq)
571 facqq1=comfac*faca*as**2*(1d0/6d0)*facqq1*be34
572 facqq2=comfac*faca*as**2*(1d0/6d0)*facqq2*be34
576 isig(nchn,3)=531+2*(ifl-3)
581 isig(nchn,3)=532+2*(ifl-3)
587 facgg1=comfac*as**2*9d0/4d0*(sh2/th2+2d0*sh/th+3d0+
588 & 2d0*th/sh+th2/sh2)*faca
589 facgg2=comfac*as**2*9d0/4d0*(uh2/sh2+2d0*uh/sh+3d0+
590 & 2d0*sh/uh+sh2/uh2)*faca
591 facgg3=comfac*as**2*9d0/4d0*(th2/uh2+2d0*th/uh+3+
597 sigh(nchn)=0.5d0*facgg1
602 sigh(nchn)=0.5d0*facgg2
607 sigh(nchn)=0.5d0*facgg3
609 ELSEIF(isub.EQ.99)
THEN
611 IF(
mint(107).EQ.4)
THEN
620 comfac=paru(5)*4d0*paru(1)**2*paru(101)*
vint(315)*
vint(316)
621 pm2rho=pmas(
pycomp(113),1)**2
622 IF(
mstp(19).EQ.0)
THEN
625 ELSEIF(
mstp(19).EQ.1)
THEN
626 comfac=comfac/(q2ga+pm2rho)
628 IF(
mint(11).EQ.22.AND.
mint(12).EQ.22)
THEN
631 xga=q2ga/(w2ga+q2ga-pmas(
pycomp(
mint(10+isde)),1)**2)
633 comfac=comfac/
max(1d-2,1d0-xga)
635 ELSEIF(
mstp(19).EQ.2)
THEN
636 comfac=comfac*q2ga/(q2ga+pm2rho)**2
638 comfac=comfac*q2ga/(q2ga+pm2rho)**2
640 IF(
mint(11).EQ.22.AND.
mint(12).EQ.22)
THEN
641 rdrds=4.1d-3*w2ga**2.167d0/((q2ga+0.15d0*w2ga)**2*
642 & q2ga**0.75d0)*(1d0+0.11d0*q2ga*p2ga/(1d0+0.02d0*p2ga**2))
645 rdrds=1.5d-4*w2ga**2.167d0/((q2ga+0.041d0*w2ga)**2*
647 xga=q2ga/(w2ga+q2ga-pmas(
pycomp(
mint(10+isde)),1)**2)
649 comfac=comfac*exp(-
max(1d-10,rdrds))
650 IF(
mstp(19).EQ.4) comfac=comfac/
max(1d-2,1d0-xga)
653 IF(
i.EQ.0.OR.kfac(isde,
i).EQ.0) goto 390
654 IF(iabs(
i).LT.10.AND.iabs(
i).GT.
mstp(58)) goto 390
655 ei=kchg(iabs(
i),1)/3d0
660 sigh(nchn)=comfac*
ei**2
665 IF(isub.EQ.114.OR.isub.EQ.115)
THEN
683 ei=kchg(iabs(
i),1)/3d0
685 IF(isub.EQ.115) eiwt=
ei
690 IF((
mstp(38).GE.1.AND.
mstp(38).LE.8).OR.epss.LT.1d-4)
THEN
691 b0stur=1d0+(th-uh)/sh*altu+0.5d0*(th2+uh2)/sh2*(altu**2+
694 b0tsur=1d0+(sh-uh)/th*alsu+0.5d0*(sh2+uh2)/th2*alsu**2
695 b0tsui=-paru(1)*((sh-uh)/th+(sh2+uh2)/th2*alsu)
696 b0utsr=1d0+(sh-th)/uh*alst+0.5d0*(sh2+th2)/uh2*alst**2
697 b0utsi=-paru(1)*((sh-th)/uh+(sh2+th2)/uh2*alst)
703 CALL
pywaux(1,epss,w1sr,w1si)
704 CALL
pywaux(1,epst,w1tr,w1ti)
705 CALL
pywaux(1,epsu,w1ur,w1ui)
706 CALL
pywaux(2,epss,w2sr,w2si)
707 CALL
pywaux(2,epst,w2tr,w2ti)
708 CALL
pywaux(2,epsu,w2ur,w2ui)
709 CALL
pyi3au(epss,th/uh,y3stur,y3stui)
710 CALL
pyi3au(epss,uh/th,y3sutr,y3suti)
711 CALL
pyi3au(epst,sh/uh,y3tsur,y3tsui)
712 CALL
pyi3au(epst,uh/sh,y3tusr,y3tusi)
713 CALL
pyi3au(epsu,sh/th,y3ustr,y3usti)
714 CALL
pyi3au(epsu,th/sh,y3utsr,y3utsi)
715 b0stur=1d0+(1d0+2d0*th/sh)*w1tr+(1d0+2d0*uh/sh)*w1ur+
716 & 0.5d0*((th2+uh2)/sh2-epss)*(w2tr+w2ur)-
717 & 0.25d0*epst*(1d0-0.5d0*epss)*(y3sutr+y3tusr)-
718 & 0.25d0*epsu*(1d0-0.5d0*epss)*(y3stur+y3utsr)+
719 & 0.25d0*(-2d0*(th2+uh2)/sh2+4d0*epss+epst+epsu+
720 & 0.5d0*epst*epsu)*(y3tsur+y3ustr)
721 b0stui=(1d0+2d0*th/sh)*w1ti+(1d0+2d0*uh/sh)*w1ui+
722 & 0.5d0*((th2+uh2)/sh2-epss)*(w2ti+w2ui)-
723 & 0.25d0*epst*(1d0-0.5d0*epss)*(y3suti+y3tusi)-
724 & 0.25d0*epsu*(1d0-0.5d0*epss)*(y3stui+y3utsi)+
725 & 0.25d0*(-2d0*(th2+uh2)/sh2+4d0*epss+epst+epsu+
726 & 0.5d0*epst*epsu)*(y3tsui+y3usti)
727 b0tsur=1d0+(1d0+2d0*sh/th)*w1sr+(1d0+2d0*uh/th)*w1ur+
728 & 0.5d0*((sh2+uh2)/th2-epst)*(w2sr+w2ur)-
729 & 0.25d0*epss*(1d0-0.5d0*epst)*(y3tusr+y3sutr)-
730 & 0.25d0*epsu*(1d0-0.5d0*epst)*(y3tsur+y3ustr)+
731 & 0.25d0*(-2d0*(sh2+uh2)/th2+4d0*epst+epss+epsu+
732 & 0.5d0*epss*epsu)*(y3stur+y3utsr)
733 b0tsui=(1d0+2d0*sh/th)*w1si+(1d0+2d0*uh/th)*w1ui+
734 & 0.5d0*((sh2+uh2)/th2-epst)*(w2si+w2ui)-
735 & 0.25d0*epss*(1d0-0.5d0*epst)*(y3tusi+y3suti)-
736 & 0.25d0*epsu*(1d0-0.5d0*epst)*(y3tsui+y3usti)+
737 & 0.25d0*(-2d0*(sh2+uh2)/th2+4d0*epst+epss+epsu+
738 & 0.5d0*epss*epsu)*(y3stui+y3utsi)
739 b0utsr=1d0+(1d0+2d0*th/uh)*w1tr+(1d0+2d0*sh/uh)*w1sr+
740 & 0.5d0*((th2+sh2)/uh2-epsu)*(w2tr+w2sr)-
741 & 0.25d0*epst*(1d0-0.5d0*epsu)*(y3ustr+y3tsur)-
742 & 0.25d0*epss*(1d0-0.5d0*epsu)*(y3utsr+y3stur)+
743 & 0.25d0*(-2d0*(th2+sh2)/uh2+4d0*epsu+epst+epss+
744 & 0.5d0*epst*epss)*(y3tusr+y3sutr)
745 b0utsi=(1d0+2d0*th/uh)*w1ti+(1d0+2d0*sh/uh)*w1si+
746 & 0.5d0*((th2+sh2)/uh2-epsu)*(w2ti+w2si)-
747 & 0.25d0*epst*(1d0-0.5d0*epsu)*(y3usti+y3tsui)-
748 & 0.25d0*epss*(1d0-0.5d0*epsu)*(y3utsi+y3stui)+
749 & 0.25d0*(-2d0*(th2+sh2)/uh2+4d0*epsu+epst+epss+
750 & 0.5d0*epst*epss)*(y3tusi+y3suti)
751 b1stur=-1d0-0.25d0*(epss+epst+epsu)*(w2sr+w2tr+w2ur)+
752 & 0.25d0*(epsu+0.5d0*epss*epst)*(y3sutr+y3tusr)+
753 & 0.25d0*(epst+0.5d0*epss*epsu)*(y3stur+y3utsr)+
754 & 0.25d0*(epss+0.5d0*epst*epsu)*(y3tsur+y3ustr)
755 b1stui=-0.25d0*(epss+epst+epsu)*(w2si+w2ti+w2ui)+
756 & 0.25d0*(epsu+0.5d0*epss*epst)*(y3suti+y3tusi)+
757 & 0.25d0*(epst+0.5d0*epss*epsu)*(y3stui+y3utsi)+
758 & 0.25d0*(epss+0.5d0*epst*epsu)*(y3tsui+y3usti)
759 b2stur=-1d0+0.125d0*epss*epst*(y3sutr+y3tusr)+
760 & 0.125d0*epss*epsu*(y3stur+y3utsr)+
761 & 0.125d0*epst*epsu*(y3tsur+y3ustr)
762 b2stui=0.125d0*epss*epst*(y3suti+y3tusi)+
763 & 0.125d0*epss*epsu*(y3stui+y3utsi)+
764 & 0.125d0*epst*epsu*(y3tsui+y3usti)
766 a0stur=a0stur+eiwt*b0stur
767 a0stui=a0stui+eiwt*b0stui
768 a0tsur=a0tsur+eiwt*b0tsur
769 a0tsui=a0tsui+eiwt*b0tsui
770 a0utsr=a0utsr+eiwt*b0utsr
771 a0utsi=a0utsi+eiwt*b0utsi
772 a1stur=a1stur+eiwt*b1stur
773 a1stui=a1stui+eiwt*b1stui
774 a2stur=a2stur+eiwt*b2stur
775 a2stui=a2stui+eiwt*b2stui
777 asqsum=a0stur**2+a0stui**2+a0tsur**2+a0tsui**2+a0utsr**2+
778 & a0utsi**2+4d0*a1stur**2+4d0*a1stui**2+a2stur**2+a2stui**2
779 facgg=comfac*faca/(16d0*paru(1)**2)*as**2*aem**2*asqsum
780 facgp=comfac*faca*5d0/(192d0*paru(1)**2)*as**3*aem*asqsum
781 IF(kfac(1,21)*kfac(2,21).EQ.0) goto 410
786 IF(isub.EQ.114) sigh(nchn)=0.5d0*facgg
787 IF(isub.EQ.115) sigh(nchn)=facgp
790 ELSEIF(isub.EQ.131.OR.isub.EQ.132)
THEN
793 IF(
mint(15).EQ.22.AND.
mint(107).EQ.0.AND.
vint(3).LT.0d0)
795 IF(
mint(16).EQ.22.AND.
mint(108).EQ.0.AND.
vint(4).LT.0d0)
798 fgq=comfac*as*aem*8d0/3d0*sh**2/(sh+ph)**2*
799 & ((sh2+uh2-2d0*ph*th)/(-sh*uh)-2d0*ph*th/(sh+ph)**2)
801 fgq=comfac*as*aem*8d0/3d0*sh**2/(sh+ph)**4*(-4d0*ph*th)
804 IF(
i.EQ.0.OR.iabs(
i).GT.
mstp(58)) goto 430
805 ei=kchg(iabs(
i),1)/3d0
808 IF(isde.EQ.1.AND.kfac(1,
i)*kfac(2,22).EQ.0) goto 420
809 IF(isde.EQ.2.AND.kfac(1,22)*kfac(2,
i).EQ.0) goto 420
818 ELSEIF(isub.EQ.133.OR.isub.EQ.134)
THEN
821 IF(
mint(15).EQ.22.AND.
mint(107).EQ.0.AND.
vint(3).LT.0d0)
823 IF(
mint(16).EQ.22.AND.
mint(108).EQ.0.AND.
vint(4).LT.0d0)
826 fgq=comfac*aem**2*2d0*sh**2/(sh+ph)**2*
827 & ((sh2+uh2-2d0*ph*th)/(-sh*uh)-2d0*ph*th/(sh+ph)**2)
829 fgq=comfac*aem**2*2d0*sh**2/(sh+ph)**4*(-4d0*ph*th)
833 ei=kchg(iabs(
i),1)/3d0
836 IF(isde.EQ.1.AND.kfac(1,
i)*kfac(2,22).EQ.0) goto 440
837 IF(isde.EQ.2.AND.kfac(1,22)*kfac(2,
i).EQ.0) goto 440
846 ELSEIF(isub.EQ.135.OR.isub.EQ.136)
THEN
849 IF(
mint(15).EQ.22.AND.
mint(107).EQ.0.AND.
vint(3).LT.0d0)
851 IF(
mint(16).EQ.22.AND.
mint(108).EQ.0.AND.
vint(4).LT.0d0)
853 CALL
pywidt(21,sh,wdtp,wdte)
855 DO 460
i=1,min(8,mdcy(21,3))
857 wdtesu=wdtesu+ef**2*(wdte(
i,1)+wdte(
i,2)+wdte(
i,3)+
861 facqq=comfac*aem*as*wdtesu*sh**2/(sh+ph)**2*
862 & ((th2+uh2-2d0*ph*sh)/(th*uh)+4d0*ph*sh/(sh+ph)**2)
864 facqq=comfac*aem*as*wdtesu*sh**2/(sh+ph)**4*8d0*ph*sh
866 IF(kfac(1,21)*kfac(2,22).NE.0)
THEN
873 IF(kfac(1,22)*kfac(2,21).NE.0)
THEN
881 ELSEIF(isub.GE.137.AND.isub.LE.140)
THEN
887 CALL
pywidt(22,sh,wdtp,wdte)
889 DO 470
i=1,min(12,mdcy(22,3))
890 IF(
i.LE.8) ef= kchg(
i,1)/3d0
891 IF(
i.GE.9) ef= kchg(9+2*(
i-8),1)/3d0
892 wdtesu=wdtesu+ef**2*(wdte(
i,1)+wdte(
i,2)+wdte(
i,3)+
895 dlamb2=(th+uh)**2-4d0*ph1*ph2
897 fparam=-sh*(th+uh)/dlamb2
898 facff=comfac*aem**2*wdtesu*2d0*sh2/(dlamb2*th2*uh2)*
899 & (th*uh-ph1*ph2)*((th2+uh2)*(1d0-2d0*fparam*(1d0-fparam))-
900 & 2d0*ph1*ph2*fparam**2)
901 ELSEIF(isub.EQ.138)
THEN
902 facff=comfac*aem**2*wdtesu*4d0*sh2*sh/(dlamb2**2*th2*uh2)*
903 & ph2*(4d0*(th*uh-ph1*ph2)*(th*uh+ph1*sh*(th-uh)**2/dlamb2)+
904 & 2d0*ph1**2*(th-uh)**2)
905 ELSEIF(isub.EQ.139)
THEN
906 facff=comfac*aem**2*wdtesu*4d0*sh2*sh/(dlamb2**2*th2*uh2)*
907 & ph1*(4d0*(th*uh-ph1*ph2)*(th*uh+ph2*sh*(th-uh)**2/dlamb2)+
908 & 2d0*ph2**2*(th-uh)**2)
910 facff=comfac*aem**2*wdtesu*32d0*sh2**2/(dlamb2**3*th2*uh2)*
911 & ph1*ph2*(th*uh-ph1*ph2)*(th-uh)**2
913 IF(kfac(1,22)*kfac(2,22).NE.0)
THEN