7 SUBROUTINE pywidt(KFLR,SH,WDTP,WDTE)
10 IMPLICIT DOUBLE PRECISION(a-
h, o-
z)
14 parameter(ksusy1=1000000,ksusy2=2000000,ktechn=3000000,
15 &kexcit=4000000,kdimen=5000000)
17 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
18 common/
pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
19 common/
pydat3/mdcy(500,3),mdme(8000,2),brat(8000),kfdp(8000,5)
20 common/
pysubs/msel,mselpd,msub(500),kfin(2,-40:40),ckin(200)
23 common/
pyint4/mwid(500),wids(500,5)
24 common/
pymssm/imss(0:99),rmss(0:99)
25 common/
pyssmt/zmix(4,4),umix(2,2),vmix(2,2),smz(4),smw(2),
26 &sfmix(16,4),zmixi(4,4),umixi(2,2),vmixi(2,2)
27 common/pytcsm/itcm(0:99),rtcm(0:99)
31 COMPLEX*16 zmixc(4,4),al,bl,ar,br,fl,fr
32 dimension wdtp(0:400),wdte(0:400,0:5),mofsv(3,2),widwsv(3,2),
33 &wid2sv(3,2),wdtpp(0:400),wdtep(0:400,0:5)
34 SAVE mofsv,widwsv,wid2sv
35 DATA mofsv/6*0/,widwsv/6*0d0/,wid2sv/6*0d0/
54 IF(
mstp(110).NE.0.AND.(mwid(kc).EQ.1.OR.mwid(kc).EQ.2.OR.
55 &(mwid(kc).EQ.3.AND.
mint(63).EQ.1)))
THEN
56 IF(
mstp(110).EQ.kfla)
THEN
58 ELSEIF(
mstp(110).EQ.-1)
THEN
59 IF(kfla.NE.6.AND.kfla.NE.23.AND.kfla.NE.24) fudge=
parp(110)
60 ELSEIF(
mstp(110).EQ.-2)
THEN
66 IF((mwid(kc).LE.0.OR.mwid(kc).GE.4).AND.kfla.NE.21.AND.
76 ELSEIF(mwid(kc).EQ.2.OR.(mwid(kc).EQ.3.AND.
mint(63).EQ.0))
THEN
80 IF(mdme(idc,1).LT.0) goto 120
85 IF(kchg(kfc1,3).EQ.1) kfd1=kfls*kfd1
89 IF(kchg(kfc2,3).EQ.1) kfd2=kfls*kfd2
95 IF(kchg(kfc3,3).EQ.1) kfd3=kfls*kfd3
100 wdtp(
i)=pmas(kc,2)*brat(idc)*(shr/pmr)
101 IF(mdme(idc,2).GE.51.AND.mdme(idc,2).LE.53.AND.
102 & pm1+pm2+pm3.GE.shr)
THEN
104 ELSEIF(mdme(idc,2).EQ.52.AND.kfd3.EQ.0)
THEN
105 wdtp(
i)=wdtp(
i)*sqrt(
max(0d0,(sh-pm1**2-pm2**2)**2-
106 & 4d0*pm1**2*pm2**2))/sh
107 ELSEIF(mdme(idc,2).EQ.52)
THEN
110 pmb=pm1+pm2+pm3-pma-pmc
111 pmbc=pmb+pmc+0.5d0*(shr-pma-pmc-pmc)
116 wdtp(
i)=wdtp(
i)*sqrt(
max(0d0,
117 & ((1d0-pman-pmbcn)**2-4d0*pman*pmbcn)*
118 & ((pmbcn-pmbn-pmcn)**2-4d0*pmbn*pmcn)))*
119 & ((shr-pma)**2-(pmb+pmc)**2)*
120 & (1d0+0.25d0*(pma+pmb+pmc)/shr)/
121 & ((1d0-pmbcn)*pmbcn*sh)
122 ELSEIF(mdme(idc,2).EQ.53.AND.kfd3.EQ.0)
THEN
123 wdtp(
i)=wdtp(
i)*sqrt(
124 &
max(0d0,(sh-pm1**2-pm2**2)**2-4d0*pm1**2*pm2**2)/
125 &
max(1d-4,(pmr**2-pm1**2-pm2**2)**2-4d0*pm1**2*pm2**2))
126 ELSEIF(mdme(idc,2).EQ.53)
THEN
129 pmb=pm1+pm2+pm3-pma-pmc
130 pmbc=pmb+pmc+0.5d0*(shr-pma-pmb-pmc)
136 & ((1d0-pman-pmbcn)**2-4d0*pman*pmbcn)*
137 & ((pmbcn-pmbn-pmcn)**2-4d0*pmbn*pmcn)))*
138 & ((shr-pma)**2-(pmb+pmc)**2)*
139 & (1d0+0.25d0*(pma+pmb+pmc)/shr)/
140 & ((1d0-pmbcn)*pmbcn*sh)
141 pmbc=pmb+pmc+0.5d0*(pmr-pma-pmb-pmc)
147 & ((1d0-pman-pmbcn)**2-4d0*pman*pmbcn)*
148 & ((pmbcn-pmbn-pmcn)**2-4d0*pmbn*pmcn)))*
149 & ((pmr-pma)**2-(pmb+pmc)**2)*
150 & (1d0+0.25d0*(pma+pmb+pmc)/pmr)/
151 & ((1d0-pmbcn)*pmbcn*pmr**2)
152 wdtp(
i)=wdtp(
i)*facact/
max(1d-6,facnom)
154 wdtp(
i)=fudge*wdtp(
i)
155 wdtp(0)=wdtp(0)+wdtp(
i)
159 IF(mdme(idc,1).GT.0)
THEN
160 IF(kfd2.EQ.kfd1)
THEN
161 IF(kchg(kfc1,3).EQ.0)
THEN
163 ELSEIF(kfd1.GT.0)
THEN
169 wid2=wid2*wids(kfc3,2)
170 ELSEIF(kfd3.LT.0)
THEN
171 wid2=wid2*wids(kfc3,3)
173 ELSEIF(kfd2.EQ.-kfd1)
THEN
176 wid2=wid2*wids(kfc3,2)
177 ELSEIF(kfd3.LT.0)
THEN
178 wid2=wid2*wids(kfc3,3)
180 ELSEIF(kfd3.EQ.kfd1)
THEN
181 IF(kchg(kfc1,3).EQ.0)
THEN
183 ELSEIF(kfd1.GT.0)
THEN
189 wid2=wid2*wids(kfc2,2)
190 ELSEIF(kfd2.LT.0)
THEN
191 wid2=wid2*wids(kfc2,3)
193 ELSEIF(kfd3.EQ.-kfd1)
THEN
196 wid2=wid2*wids(kfc2,2)
197 ELSEIF(kfd2.LT.0)
THEN
198 wid2=wid2*wids(kfc2,3)
200 ELSEIF(kfd3.EQ.kfd2)
THEN
201 IF(kchg(kfc2,3).EQ.0)
THEN
203 ELSEIF(kfd2.GT.0)
THEN
209 wid2=wid2*wids(kfc1,2)
210 ELSEIF(kfd1.LT.0)
THEN
211 wid2=wid2*wids(kfc1,3)
213 ELSEIF(kfd3.EQ.-kfd2)
THEN
216 wid2=wid2*wids(kfc1,2)
217 ELSEIF(kfd1.LT.0)
THEN
218 wid2=wid2*wids(kfc1,3)
227 wid2=wid2*wids(kfc2,2)
229 wid2=wid2*wids(kfc2,3)
232 wid2=wid2*wids(kfc3,2)
233 ELSEIF(kfd3.LT.0)
THEN
234 wid2=wid2*wids(kfc3,3)
239 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
240 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
241 wdte(
i,0)=wdte(
i,mdme(idc,1))
242 wdte(0,0)=wdte(0,0)+wdte(
i,0)
256 IF(kfla.EQ.35.OR.kfla.EQ.36)
THEN
264 IF(
mstp(8).GE.2) xw=1d0-(pmas(24,1)/pmas(23,1))**2
267 IF(
mstp(8).GE.1) aem=sqrt(2d0)*paru(105)*pmas(24,1)**2*xw/paru(1)
273 fac=(aem/(16d0*xw))*(sh/pmas(24,1)**2)*shr
274 radct=1d0-2.5d0*as/paru(1)
275 DO 140
i=1,mdcy(kc,3)
277 IF(mdme(idc,1).LT.0) goto 140
278 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
279 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
280 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 140
282 IF(
i.GE.4.AND.
i.LE.7)
THEN
284 wdtp(
i)=fac*vckm(3,
i-3)*radct*
285 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
286 & ((1d0-rm2)**2+(1d0+rm2)*rm1-2d0*rm1**2)
289 IF(
i.EQ.7) wid2=wid2*wids(7,2)
292 IF(
i.EQ.7) wid2=wid2*wids(7,3)
296 rm2r=
pymrun(kfdp(idc,2),sh)**2/sh
297 wdtp(
i)=fac*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
298 & ((1d0+rm2-rm1)*(rm2r*paru(141)**2+1d0/paru(141)**2)+
299 & 4d0*sqrt(rm2r*rm2))
301 IF(kflr.LT.0) wid2=wids(37,3)
303 ELSEIF(
i.GE.10.AND.
i.LE.13.AND.imss(1).NE.0)
THEN
307 tanw=sqrt(paru(102)/(1d0-paru(102)))
314 IF(shr.GT.pmnchi+pmstop)
THEN
317 zmixc(iz,ik)=dcmplx(zmix(iz,ik),zmixi(iz,ik))
319 al=shr*dconjg(zmixc(iz,4))/(2.0d0*pmas(24,1)*sinb)
320 ar=-et*zmixc(iz,1)*tanw
321 bl=t3l*(zmixc(iz,2)-zmixc(iz,1)*tanw)-ar
323 fl=sfmix(6,1)*al+sfmix(6,2)*ar
324 fr=sfmix(6,1)*bl+sfmix(6,2)*br
325 pcm=sqrt((sh-(pmnchi+pmstop)**2)*
326 & (sh-(pmnchi-pmstop)**2))/(2d0*shr)
327 wdtp(
i)=(0.5d0*
pyalem(sh)/paru(102))*pcm*
328 & ((abs(fl)**2+abs(fr)**2)*(sh+pmnchi**2-pmstop**2)+
329 & smz(iz)*4d0*shr*dble(fl*dconjg(fr)))/sh
331 wid2=wids(kfc1,2)*wids(kfc2,2)
333 wid2=wids(kfc1,2)*wids(kfc2,3)
336 ELSEIF(
i.EQ.14.AND.imss(1).NE.0)
THEN
342 IF(shr.GT.pmnchi+pmstop)
THEN
345 pcm=sqrt((sh-(pmnchi+pmstop)**2)*
346 & (sh-(pmnchi-pmstop)**2))/(2d0*shr)
347 wdtp(
i)=4d0/3d0*0.5d0*
pyalps(sh)*pcm*((rl**2+rr**2)*
348 & (sh+pmnchi**2-pmstop**2)+pmnchi*4d0*shr*rl*rr)/sh
350 wid2=wids(kfc1,2)*wids(kfc2,2)
352 wid2=wids(kfc1,2)*wids(kfc2,3)
355 ELSEIF(
i.EQ.15.AND.imss(1).NE.0)
THEN
359 xmgr2=pmas(kfc1,1)**2
360 wdtp(
i)=sh**2*shr/(96d0*paru(1)*xmp2*xmgr2)*(1d0-rm2)**4
363 IF(kflr.LT.0) wid2=wids(kfc2,3)
366 wdtp(
i)=fudge*wdtp(
i)
367 wdtp(0)=wdtp(0)+wdtp(
i)
368 IF(mdme(idc,1).GT.0)
THEN
369 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
370 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
371 wdte(
i,0)=wdte(
i,mdme(idc,1))
372 wdte(0,0)=wdte(0,0)+wdte(
i,0)
376 ELSEIF(kfla.EQ.7)
THEN
378 fac=(aem/(16d0*xw))*(sh/pmas(24,1)**2)*shr
379 DO 150
i=1,mdcy(kc,3)
381 IF(mdme(idc,1).LT.0) goto 150
382 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
383 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
384 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 150
386 IF(
i.GE.4.AND.
i.LE.7)
THEN
388 wdtp(
i)=fac*vckm(
i-3,4)*
389 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
390 & ((1d0-rm2)**2+(1d0+rm2)*rm1-2d0*rm1**2)
393 IF(
i.EQ.6) wid2=wid2*wids(6,2)
394 IF(
i.EQ.7) wid2=wid2*wids(8,2)
397 IF(
i.EQ.6) wid2=wid2*wids(6,3)
398 IF(
i.EQ.7) wid2=wid2*wids(8,3)
401 IF(kflr.LT.0) wid2=wids(24,2)
402 ELSEIF(
i.EQ.9.OR.
i.EQ.10)
THEN
404 wdtp(
i)=fac*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
405 & ((1d0+rm2-rm1)*(paru(141)**2+rm2/paru(141)**2)+4d0*rm2)
408 IF(
i.EQ.10) wid2=wid2*wids(6,2)
411 IF(
i.EQ.10) wid2=wid2*wids(6,3)
414 wdtp(
i)=fudge*wdtp(
i)
415 wdtp(0)=wdtp(0)+wdtp(
i)
416 IF(mdme(idc,1).GT.0)
THEN
417 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
418 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
419 wdte(
i,0)=wdte(
i,mdme(idc,1))
420 wdte(0,0)=wdte(0,0)+wdte(
i,0)
424 ELSEIF(kfla.EQ.8)
THEN
426 fac=(aem/(16d0*xw))*(sh/pmas(24,1)**2)*shr
427 DO 160
i=1,mdcy(kc,3)
429 IF(mdme(idc,1).LT.0) goto 160
430 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
431 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
432 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 160
434 IF(
i.GE.4.AND.
i.LE.7)
THEN
436 wdtp(
i)=fac*vckm(4,
i-3)*
437 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
438 & ((1d0-rm2)**2+(1d0+rm2)*rm1-2d0*rm1**2)
441 IF(
i.EQ.7) wid2=wid2*wids(7,2)
444 IF(
i.EQ.7) wid2=wid2*wids(7,3)
446 ELSEIF(
i.EQ.9.OR.
i.EQ.10)
THEN
448 wdtp(
i)=fac*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
449 & ((1d0+rm2-rm1)*(rm2*paru(141)**2+1d0/paru(141)**2)+4d0*rm2)
452 IF(
i.EQ.10) wid2=wid2*wids(7,2)
455 IF(
i.EQ.10) wid2=wid2*wids(7,3)
458 wdtp(
i)=fudge*wdtp(
i)
459 wdtp(0)=wdtp(0)+wdtp(
i)
460 IF(mdme(idc,1).GT.0)
THEN
461 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
462 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
463 wdte(
i,0)=wdte(
i,mdme(idc,1))
464 wdte(0,0)=wdte(0,0)+wdte(
i,0)
468 ELSEIF(kfla.EQ.17)
THEN
470 fac=(aem/(16d0*xw))*(sh/pmas(24,1)**2)*shr
471 DO 170
i=1,mdcy(kc,3)
473 IF(mdme(idc,1).LT.0) goto 170
474 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
475 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
476 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 170
480 wdtp(
i)=fac*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
481 & ((1d0-rm2)**2+(1d0+rm2)*rm1-2d0*rm1**2)
491 wdtp(
i)=fac*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
492 & ((1d0+rm2-rm1)*(paru(141)**2+rm2/paru(141)**2)+4d0*rm2)
501 wdtp(
i)=fudge*wdtp(
i)
502 wdtp(0)=wdtp(0)+wdtp(
i)
503 IF(mdme(idc,1).GT.0)
THEN
504 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
505 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
506 wdte(
i,0)=wdte(
i,mdme(idc,1))
507 wdte(0,0)=wdte(0,0)+wdte(
i,0)
511 ELSEIF(kfla.EQ.18)
THEN
513 fac=(aem/(16d0*xw))*(sh/pmas(24,1)**2)*shr
514 DO 180
i=1,mdcy(kc,3)
516 IF(mdme(idc,1).LT.0) goto 180
517 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
518 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
519 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 180
523 wdtp(
i)=fac*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
524 & ((1d0-rm2)**2+(1d0+rm2)*rm1-2d0*rm1**2)
534 wdtp(
i)=fac*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
535 & ((1d0+rm2-rm1)*(rm2*paru(141)**2+1d0/paru(141)**2)+4d0*rm2)
544 wdtp(
i)=fudge*wdtp(
i)
545 wdtp(0)=wdtp(0)+wdtp(
i)
546 IF(mdme(idc,1).GT.0)
THEN
547 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
548 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
549 wdte(
i,0)=wdte(
i,mdme(idc,1))
550 wdte(0,0)=wdte(0,0)+wdte(
i,0)
554 ELSEIF(kfla.EQ.21)
THEN
557 DO 190
i=1,mdcy(kc,3)
559 IF(mdme(idc,1).LT.0) goto 190
560 rm1=pmas(iabs(kfdp(idc,1)),1)**2/sh
561 rm2=pmas(iabs(kfdp(idc,2)),1)**2/sh
562 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 190
566 wdtp(
i)=(1d0+2d0*rm1)*sqrt(
max(0d0,1d0-4d0*rm1))
567 IF(
i.EQ.6) wid2=wids(6,1)
568 IF((
i.EQ.7.OR.
i.EQ.8)) wid2=wids(
i,1)
570 wdtp(
i)=fudge*wdtp(
i)
571 wdtp(0)=wdtp(0)+wdtp(
i)
572 IF(mdme(idc,1).GT.0)
THEN
573 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
574 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
575 wdte(
i,0)=wdte(
i,mdme(idc,1))
576 wdte(0,0)=wdte(0,0)+wdte(
i,0)
580 ELSEIF(kfla.EQ.22)
THEN
583 DO 200
i=1,mdcy(kc,3)
585 IF(mdme(idc,1).LT.0) goto 200
586 rm1=pmas(iabs(kfdp(idc,1)),1)**2/sh
587 rm2=pmas(iabs(kfdp(idc,2)),1)**2/sh
588 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 200
594 IF(
i.GE.6.AND.
mstp(35).GE.1) fcof=fcof*
pyhfth(sh,sh*rm1,1d0)
595 wdtp(
i)=fcof*ef**2*(1d0+2d0*rm1)*sqrt(
max(0d0,1d0-4d0*rm1))
596 IF(
i.EQ.6) wid2=wids(6,1)
597 IF((
i.EQ.7.OR.
i.EQ.8)) wid2=wids(
i,1)
600 ef=kchg(9+2*(
i-8),1)/3d0
601 wdtp(
i)=ef**2*(1d0+2d0*rm1)*sqrt(
max(0d0,1d0-4d0*rm1))
602 IF(
i.EQ.12) wid2=wids(17,1)
604 wdtp(
i)=fudge*wdtp(
i)
605 wdtp(0)=wdtp(0)+wdtp(
i)
606 IF(mdme(idc,1).GT.0)
THEN
607 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
608 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
609 wdte(
i,0)=wdte(
i,mdme(idc,1))
610 wdte(0,0)=wdte(0,0)+wdte(
i,0)
614 ELSEIF(kfla.EQ.23)
THEN
617 xwc=1d0/(16d0*xw*xw1)
618 fac=(aem*xwc/3d0)*shr
620 IF(
mint(61).GE.1.AND.icase.EQ.2)
THEN
625 IF(
mint(61).EQ.1.AND.icase.EQ.2)
THEN
627 IF(kfi.GT.20) kfi=iabs(
mint(16))
635 & 2d0*xwc*sh*(sh-sqmz)/((sh-sqmz)**2+hz**2)
637 & xwc**2*sh**2/((sh-sqmz)**2+hz**2)
639 DO 220
i=1,mdcy(kc,3)
641 IF(mdme(idc,1).LT.0) goto 220
642 rm1=pmas(iabs(kfdp(idc,1)),1)**2/sh
643 rm2=pmas(iabs(kfdp(idc,2)),1)**2/sh
644 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 220
649 af=
sign(1d0,ef+0.1d0)
652 IF(
i.GE.6.AND.
mstp(35).GE.1) fcof=fcof*
pyhfth(sh,sh*rm1,1d0)
653 IF(
i.EQ.6) wid2=wids(6,1)
654 IF((
i.EQ.7.OR.
i.EQ.8)) wid2=wids(
i,1)
658 af=
sign(1d0,ef+0.1d0)
661 IF((
i.EQ.15.OR.
i.EQ.16)) wid2=wids(2+
i,1)
663 be34=sqrt(
max(0d0,1d0-4d0*rm1))
665 wdtp(
i)=fac*fcof*(vf**2*(1d0+2d0*rm1)+af**2*(1d0-4d0*rm1))*
667 ELSEIF(
mint(61).EQ.1.AND.icase.EQ.2)
THEN
669 & ef*vf+(vi**2+ai**2)*
vint(114)*vf**2)*(1d0+2d0*rm1)+
670 & (vi**2+ai**2)*
vint(114)*af**2*(1d0-4d0*rm1))*be34
671 ELSEIF(
mint(61).EQ.2.AND.icase.EQ.2)
THEN
672 fggf=fcof*ef**2*(1d0+2d0*rm1)*be34
673 fgzf=fcof*ef*vf*(1d0+2d0*rm1)*be34
674 fzzf=fcof*(vf**2*(1d0+2d0*rm1)+af**2*(1d0-4d0*rm1))*be34
676 IF(icase.EQ.1) wdtp(
i)=fudge*wdtp(
i)
677 IF(icase.EQ.1) wdtp(0)=wdtp(0)+wdtp(
i)
678 IF(mdme(idc,1).GT.0)
THEN
679 IF((icase.EQ.1.AND.
mint(61).NE.1).OR.
680 & (icase.EQ.2.AND.
mint(61).EQ.1))
THEN
681 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
682 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+
683 & wdte(
i,mdme(idc,1))
684 wdte(
i,0)=wdte(
i,mdme(idc,1))
685 wdte(0,0)=wdte(0,0)+wdte(
i,0)
687 IF(
mint(61).EQ.2.AND.icase.EQ.2)
THEN
689 &
vint(111)+fggf*wid2
692 &
vint(114)+fzzf*wid2
696 IF(
mint(61).GE.1) icase=3-icase
697 IF(icase.EQ.2) goto 210
699 ELSEIF(kfla.EQ.24)
THEN
701 fac=(aem/(24d0*xw))*shr
702 DO 230
i=1,mdcy(kc,3)
704 IF(mdme(idc,1).LT.0) goto 230
705 rm1=pmas(iabs(kfdp(idc,1)),1)**2/sh
706 rm2=pmas(iabs(kfdp(idc,2)),1)**2/sh
707 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 230
711 fcof=3d0*radc*vckm((
i-1)/4+1,mod(
i-1,4)+1)
713 IF(mod(
i,4).EQ.3) wid2=wids(6,2)
714 IF(mod(
i,4).EQ.0) wid2=wids(8,2)
715 IF(
i.GE.13) wid2=wid2*wids(7,3)
717 IF(mod(
i,4).EQ.3) wid2=wids(6,3)
718 IF(mod(
i,4).EQ.0) wid2=wids(8,3)
719 IF(
i.GE.13) wid2=wid2*wids(7,2)
725 IF(
i.EQ.20) wid2=wids(17,3)*wids(18,2)
727 IF(
i.EQ.20) wid2=wids(17,2)*wids(18,3)
730 wdtp(
i)=fac*fcof*(2d0-rm1-rm2-(rm1-rm2)**2)*
731 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
732 wdtp(
i)=fudge*wdtp(
i)
733 wdtp(0)=wdtp(0)+wdtp(
i)
734 IF(mdme(idc,1).GT.0)
THEN
735 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
736 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
737 wdte(
i,0)=wdte(
i,mdme(idc,1))
738 wdte(0,0)=wdte(0,0)+wdte(
i,0)
742 ELSEIF(kfla.EQ.25.OR.kfla.EQ.35.OR.kfla.EQ.36)
THEN
745 fac=(aem/(8d0*xw))*(shfs/pmas(24,1)**2)*shr
746 DO 270
i=1,mdcy(kfhigg,3)
747 idc=
i+mdcy(kfhigg,2)-1
748 IF(mdme(idc,1).LT.0) goto 270
751 rm1=pmas(kfc1,1)**2/sh
752 rm2=pmas(kfc2,1)**2/sh
753 IF(
i.NE.16.AND.
i.NE.17.AND.sqrt(rm1)+sqrt(rm2).GT.1d0)
759 wdtp(
i)=fac*3d0*(
pymrun(kfdp(idc,1),sh)**2/shfs)*
760 & sqrt(
max(0d0,1d0-4d0*rm1))*radc
762 IF(ihigg.NE.3) wdtp(
i)=wdtp(
i)*(1d0-4d0*rm1)
763 IF(
mstp(4).GE.1.OR.ihigg.GE.2)
THEN
764 IF(mod(
i,2).EQ.1) wdtp(
i)=wdtp(
i)*paru(151+10*ihigg)**2
765 IF(mod(
i,2).EQ.0) wdtp(
i)=wdtp(
i)*paru(152+10*ihigg)**2
766 IF(imss(1).NE.0.AND.kfc1.EQ.5)
THEN
767 wdtp(
i)=wdtp(
i)/(1d0+rmss(41))**2
769 wdtp(
i)=wdtp(
i)*(1d0+rmss(41)*paru(152+10*ihigg)/
770 & paru(151+10*ihigg))**2
774 IF(
i.EQ.6) wid2=wids(6,1)
775 IF((
i.EQ.7.OR.
i.EQ.8)) wid2=wids(
i,1)
778 wdtp(
i)=fac*rm1*sqrt(
max(0d0,1d0-4d0*rm1))*(sh/shfs)
780 IF(ihigg.NE.3) wdtp(
i)=wdtp(
i)*(1d0-4d0*rm1)
781 IF(
mstp(4).GE.1.OR.ihigg.GE.2) wdtp(
i)=wdtp(
i)*
782 & paru(153+10*ihigg)**2
783 IF(
i.EQ.12) wid2=wids(17,1)
790 eps=(2d0*pmas(
j,1))**2/sh
795 rln=
log((1d0+root)/(1d0-root))
799 phire=-0.25d0*(rln**2-paru(1)**2)
800 phiim=0.5d0*paru(1)*rln
802 phire=(asin(1d0/sqrt(
eps)))**2
806 etarej=-0.5d0*
eps*(1d0+(1d0-
eps)*phire)
807 etaimj=-0.5d0*
eps*(1d0-
eps)*phiim
809 etarej=-0.5d0*
eps*phire
810 etaimj=-0.5d0*
eps*phiim
813 IF(
mstp(4).GE.1.OR.ihigg.GE.2)
THEN
814 IF(mod(
j,2).EQ.1)
THEN
815 etarej=etarej*paru(151+10*ihigg)
816 etaimj=etaimj*paru(151+10*ihigg)
818 etarej=etarej*paru(152+10*ihigg)
819 etaimj=etaimj*paru(152+10*ihigg)
825 eta2=etare**2+etaim**2
826 wdtp(
i)=fac*(as/paru(1))**2*eta2
833 IF(
mstp(4).GE.1.OR.ihigg.GE.2) jmax=jmax+1
835 IF(
j.LE.2*
mstp(1))
THEN
837 eps=(2d0*pmas(
j,1))**2/sh
838 ELSEIF(
j.LE.3*
mstp(1))
THEN
841 eps=(2d0*pmas(10+jl,1))**2/sh
842 ELSEIF(
j.EQ.3*
mstp(1)+1)
THEN
843 eps=(2d0*pmas(24,1))**2/sh
845 eps=(2d0*pmas(37,1))**2/sh
851 rln=
log((1d0+root)/(1d0-root))
855 phire=-0.25d0*(rln**2-paru(1)**2)
856 phiim=0.5d0*paru(1)*rln
858 phire=(asin(1d0/sqrt(
eps)))**2
861 IF(
j.LE.3*
mstp(1))
THEN
864 phipre=-0.5d0*
eps*(1d0+(1d0-
eps)*phire)
865 phipim=-0.5d0*
eps*(1d0-
eps)*phiim
867 phipre=-0.5d0*
eps*phire
868 phipim=-0.5d0*
eps*phiim
870 IF(
j.LE.2*
mstp(1).AND.mod(
j,2).EQ.1)
THEN
872 ejh=paru(151+10*ihigg)
873 ELSEIF(
j.LE.2*
mstp(1))
THEN
875 ejh=paru(152+10*ihigg)
878 ejh=paru(153+10*ihigg)
880 IF(
mstp(4).EQ.0.AND.ihigg.EQ.1) ejh=1d0
881 etarej=ejc*ejh*phipre
882 etaimj=ejc*ejh*phipim
883 ELSEIF(
j.EQ.3*
mstp(1)+1)
THEN
885 etarej=0.5d0+0.75d0*
eps*(1d0+(2d0-
eps)*phire)
886 etaimj=0.75d0*
eps*(2d0-
eps)*phiim
887 IF(
mstp(4).GE.1.OR.ihigg.GE.2)
THEN
888 etarej=etarej*paru(155+10*ihigg)
889 etaimj=etaimj*paru(155+10*ihigg)
893 fachhh=(pmas(24,1)/pmas(37,1))**2*
894 & paru(158+10*ihigg+2*(ihigg/3))
895 etarej=
eps*(1d0-
eps*phire)*fachhh
896 etaimj=-
eps**2*phiim*fachhh
901 eta2=etare**2+etaim**2
902 wdtp(
i)=fac*(aem/paru(1))**2*0.5d0*eta2
909 IF(
mstp(4).GE.1.OR.ihigg.GE.2) jmax=jmax+1
911 IF(
j.LE.2*
mstp(1))
THEN
913 aj=
sign(1d0,ej+0.1d0)
915 eps=(2d0*pmas(
j,1))**2/sh
916 epsp=(2d0*pmas(
j,1)/pmas(23,1))**2
917 ELSEIF(
j.LE.3*
mstp(1))
THEN
920 aj=
sign(1d0,ej+0.1d0)
922 eps=(2d0*pmas(10+jl,1))**2/sh
923 epsp=(2d0*pmas(10+jl,1)/pmas(23,1))**2
925 eps=(2d0*pmas(24,1))**2/sh
926 epsp=(2d0*pmas(24,1)/pmas(23,1))**2
932 rln=
log((1d0+root)/(1d0-root))
936 phire=-0.25d0*(rln**2-paru(1)**2)
937 phiim=0.5d0*paru(1)*rln
939 psiim=-0.5d0*root*paru(1)
941 phire=(asin(1d0/sqrt(
eps)))**2
943 psire=sqrt(
eps-1d0)*asin(1d0/sqrt(
eps))
948 IF(epsp.GT.1d-4)
THEN
949 rln=
log((1d0+root)/(1d0-root))
951 rln=
log(4d0/epsp-2d0)
953 phirep=-0.25d0*(rln**2-paru(1)**2)
954 phiimp=0.5d0*paru(1)*rln
955 psirep=0.5d0*root*rln
956 psiimp=-0.5d0*root*paru(1)
958 phirep=(asin(1d0/sqrt(epsp)))**2
960 psirep=sqrt(epsp-1d0)*asin(1d0/sqrt(epsp))
963 fxyre=
eps*epsp/(8d0*(
eps-epsp))*(1d0+
eps*epsp/(
eps-epsp)*
964 & (phire-phirep)+2d0*
eps/(
eps-epsp)*(psire-psirep))
965 fxyim=
eps**2*epsp/(8d0*(
eps-epsp)**2)*
966 & (epsp*(phiim-phiimp)+2d0*(psiim-psiimp))
967 f1re=-
eps*epsp/(2d0*(
eps-epsp))*(phire-phirep)
968 f1im=-
eps*epsp/(2d0*(
eps-epsp))*(phiim-phiimp)
969 IF(
j.LE.3*
mstp(1))
THEN
971 IF(ihigg.EQ.3) fxyre=0d0
972 IF(ihigg.EQ.3) fxyim=0d0
973 IF(
j.LE.2*
mstp(1).AND.mod(
j,2).EQ.1)
THEN
975 ejh=paru(151+10*ihigg)
976 ELSEIF(
j.LE.2*
mstp(1))
THEN
978 ejh=paru(152+10*ihigg)
981 ejh=paru(153+10*ihigg)
983 IF(
mstp(4).EQ.0.AND.ihigg.EQ.1) ejh=1d0
984 etarej=ejc*ejh*(fxyre-0.25d0*f1re)
985 etaimj=ejc*ejh*(fxyim-0.25d0*f1im)
986 ELSEIF(
j.EQ.3*
mstp(1)+1)
THEN
988 heps=(1d0+2d0/
eps)*xw/xw1-(5d0+2d0/
eps)
989 etarej=-xw1*((3d0-xw/xw1)*f1re+heps*fxyre)
990 etaimj=-xw1*((3d0-xw/xw1)*f1im+heps*fxyim)
991 IF(
mstp(4).GE.1.OR.ihigg.GE.2)
THEN
992 etarej=etarej*paru(155+10*ihigg)
993 etaimj=etaimj*paru(155+10*ihigg)
997 fachhh=(pmas(24,1)/pmas(37,1))**2*(1d0-2d0*xw)*
998 & paru(158+10*ihigg+2*(ihigg/3))
1005 eta2=(etare**2+etaim**2)/(xw*xw1)
1006 wdtp(
i)=fac*(aem/paru(1))**2*(1d0-pmas(23,1)**2/sh)**3*eta2
1009 ELSEIF(
i.LE.17)
THEN
1011 pm1=pmas(iabs(kfdp(idc,1)),1)
1012 pg1=pmas(iabs(kfdp(idc,1)),2)
1013 IF(
mint(62).GE.1)
THEN
1014 IF(
mstp(42).EQ.0.OR.(4d0*(pm1+10d0*pg1)**2.LT.sh.AND.
1015 & ckin(46).LT.ckin(45).AND.ckin(48).LT.ckin(47).AND.
1016 &
max(ckin(45),ckin(47)).LT.pm1-10d0*pg1))
THEN
1018 widw=(1d0-4d0*rm1+12d0*rm1**2)*sqrt(
max(0d0,
1023 rmas=sqrt(
max(0d0,sh))
1024 CALL
pyofsh(1,kfla,kfdp(idc,1),kfdp(idc,2),rmas,widw,
1026 widwsv(ihigg,
i-15)=widw
1027 wid2sv(ihigg,
i-15)=wid2
1030 IF(mofsv(ihigg,
i-15).EQ.0)
THEN
1031 widw=(1d0-4d0*rm1+12d0*rm1**2)*sqrt(
max(0d0,
1035 widw=widwsv(ihigg,
i-15)
1036 wid2=wid2sv(ihigg,
i-15)
1039 wdtp(
i)=fac*widw/(2d0*(18-
i))
1040 IF(
mstp(49).NE.0) wdtp(
i)=wdtp(
i)*pmas(kfhigg,1)**2/shfs
1041 IF(
mstp(4).GE.1.OR.ihigg.GE.2) wdtp(
i)=wdtp(
i)*
1042 & paru(138+
i+10*ihigg)**2
1043 wid2=wid2*wids(7+
i,1)
1045 ELSEIF(
i.EQ.18.AND.ihigg.GE.2)
THEN
1047 wdtp(
i)=fac*0.5d0*sqrt(
max(0d0,
1048 & (1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1050 wdtp(
i)=wdtp(
i)*paru(179)**2
1051 ELSEIF(ihigg.EQ.3)
THEN
1052 wdtp(
i)=wdtp(
i)*paru(186)**2
1054 wid2=wids(23,2)*wids(25,2)
1056 ELSEIF(
i.EQ.19.AND.ihigg.GE.2)
THEN
1059 & pmas(23,1)**4/sh**2*sqrt(
max(0d0,1d0-4d0*rm1))
1061 wdtp(
i)=wdtp(
i)*paru(176)**2
1062 ELSEIF(ihigg.EQ.3)
THEN
1063 wdtp(
i)=wdtp(
i)*paru(169)**2
1066 ELSEIF((
i.EQ.20.OR.
i.EQ.21).AND.ihigg.GE.2)
THEN
1068 wdtp(
i)=fac*0.5d0*sqrt(
max(0d0,
1069 & (1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1070 & *paru(195+ihigg)**2
1072 wid2=wids(24,2)*wids(37,3)
1073 ELSEIF(
i.EQ.21)
THEN
1074 wid2=wids(24,3)*wids(37,2)
1077 ELSEIF(
i.EQ.22.AND.ihigg.EQ.2)
THEN
1079 wdtp(
i)=fac*0.5d0*paru(187)**2*sqrt(
max(0d0,
1080 & (1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1081 wid2=wids(36,2)*wids(23,2)
1083 ELSEIF(
i.EQ.23.AND.ihigg.EQ.2)
THEN
1085 wdtp(
i)=fac*0.5d0*paru(180)**2*
1086 & pmas(23,1)**4/sh**2*sqrt(
max(0d0,1d0-4d0*rm1))
1087 wid2=wids(25,2)*wids(36,2)
1089 ELSEIF(
i.EQ.24.AND.ihigg.EQ.2)
THEN
1091 wdtp(
i)=fac*0.25d0*paru(177)**2*
1092 & pmas(23,1)**4/sh**2*sqrt(
max(0d0,1d0-4d0*rm1))
1100 wfac0=1d0+rm10**2+rm20**2-2d0*(rm10+rm20+rm10*rm20)
1101 wfac=1d0+rm1**2+rm2**2-2d0*(rm1+rm2+rm1*rm2)
1102 IF(wfac.LE.0d0 .OR. wfac0.LE.0d0)
THEN
1107 wdtp(
i)=pmas(kfla,2)*brat(idc)*(shr/pmr)*sqrt(wfac)
1109 IF(kfc2.EQ.kfc1)
THEN
1113 IF(kfdp(idc,1).LT.0) ksgn1=3
1115 IF(kfdp(idc,2).LT.0) ksgn2=3
1116 wid2=wids(kfc1,ksgn1)*wids(kfc2,ksgn2)
1119 wdtp(
i)=fudge*wdtp(
i)
1120 wdtp(0)=wdtp(0)+wdtp(
i)
1121 IF(mdme(idc,1).GT.0)
THEN
1122 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1123 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
1124 wdte(
i,0)=wdte(
i,mdme(idc,1))
1125 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1129 ELSEIF(kfla.EQ.32)
THEN
1132 xwc=1d0/(16d0*xw*xw1)
1133 fac=(aem*xwc/3d0)*shr
1136 IF(
mint(61).GE.1.AND.icase.EQ.2)
THEN
1144 IF(
mint(61).EQ.1.AND.icase.EQ.2)
THEN
1150 IF(kfai.LE.10.AND.mod(kfai,2).EQ.0) kfaic=2
1151 IF(kfai.GT.10.AND.mod(kfai,2).NE.0) kfaic=3
1152 IF(kfai.GT.10.AND.mod(kfai,2).EQ.0) kfaic=4
1153 IF(kfai.LE.2.OR.kfai.EQ.11.OR.kfai.EQ.12)
THEN
1154 vpi=paru(119+2*kfaic)
1155 api=paru(120+2*kfaic)
1156 ELSEIF(kfai.LE.4.OR.kfai.EQ.13.OR.kfai.EQ.14)
THEN
1157 vpi=parj(178+2*kfaic)
1158 api=parj(179+2*kfaic)
1160 vpi=parj(186+2*kfaic)
1161 api=parj(187+2*kfaic)
1170 & 2d0*xwc*sh*(sh-sqmz)/((sh-sqmz)**2+hz**2)
1172 & 2d0*xwc*sh*(sh-sqmzp)/((sh-sqmzp)**2+hzp**2)
1174 &
mstp(44).EQ.7)
vint(114)=xwc**2*sh**2/((sh-sqmz)**2+hz**2)
1176 & 2d0*xwc**2*sh**2*((sh-sqmz)*(sh-sqmzp)+hz*hzp)/
1177 & (((sh-sqmz)**2+hz**2)*((sh-sqmzp)**2+hzp**2))
1179 &
mstp(44).EQ.7)
vint(116)=xwc**2*sh**2/((sh-sqmzp)**2+hzp**2)
1181 DO 290
i=1,mdcy(kc,3)
1183 IF(mdme(idc,1).LT.0) goto 290
1184 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
1185 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
1186 IF(sqrt(rm1)+sqrt(rm2).GT.1d0.OR.mdme(idc,1).LT.0) goto 290
1192 af=
sign(1d0,ef+0.1d0)
1195 vpf=paru(123-2*mod(
i,2))
1196 apf=paru(124-2*mod(
i,2))
1198 vpf=parj(182-2*mod(
i,2))
1199 apf=parj(183-2*mod(
i,2))
1201 vpf=parj(190-2*mod(
i,2))
1202 apf=parj(191-2*mod(
i,2))
1205 IF(
i.GE.6.AND.
mstp(35).GE.1) fcof=fcof*
1207 IF(
i.EQ.6) wid2=wids(6,1)
1208 IF((
i.EQ.7.OR.
i.EQ.8)) wid2=wids(
i,1)
1209 ELSEIF(
i.LE.16)
THEN
1212 af=
sign(1d0,ef+0.1d0)
1215 vpf=paru(127-2*mod(
i,2))
1216 apf=paru(128-2*mod(
i,2))
1217 ELSEIF(
i.LE.12)
THEN
1218 vpf=parj(186-2*mod(
i,2))
1219 apf=parj(187-2*mod(
i,2))
1221 vpf=parj(194-2*mod(
i,2))
1222 apf=parj(195-2*mod(
i,2))
1225 IF((
i.EQ.15.OR.
i.EQ.16)) wid2=wids(2+
i,1)
1227 be34=sqrt(
max(0d0,1d0-4d0*rm1))
1229 wdtpz=fcof*(vf**2*(1d0+2d0*rm1)+af**2*(1d0-4d0*rm1))*be34
1230 wdtp(
i)=fac*fcof*(vpf**2*(1d0+2d0*rm1)+
1231 & apf**2*(1d0-4d0*rm1))*be34
1232 ELSEIF(
mint(61).EQ.1.AND.icase.EQ.2)
THEN
1234 & ef*vf+
ei*vpi*
vint(113)*ef*vpf+(vi**2+ai**2)*
vint(114)*
1235 & vf**2+(vi*vpi+ai*api)*
vint(115)*vf*vpf+(vpi**2+api**2)*
1236 &
vint(116)*vpf**2)*(1d0+2d0*rm1)+((vi**2+ai**2)*
vint(114)*
1237 & af**2+(vi*vpi+ai*api)*
vint(115)*af*apf+(vpi**2+api**2)*
1238 &
vint(116)*apf**2)*(1d0-4d0*rm1))*be34
1239 ELSEIF(
mint(61).EQ.2)
THEN
1240 fggf=fcof*ef**2*(1d0+2d0*rm1)*be34
1241 fgzf=fcof*ef*vf*(1d0+2d0*rm1)*be34
1242 fgzpf=fcof*ef*vpf*(1d0+2d0*rm1)*be34
1243 fzzf=fcof*(vf**2*(1d0+2d0*rm1)+af**2*(1d0-4d0*rm1))*be34
1244 fzzpf=fcof*(vf*vpf*(1d0+2d0*rm1)+af*apf*(1d0-4d0*rm1))*
1246 fzpzpf=fcof*(vpf**2*(1d0+2d0*rm1)+apf**2*(1d0-4d0*rm1))*
1249 ELSEIF(
i.EQ.17)
THEN
1251 wdtpzp=paru(129)**2*xw1**2*
1252 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1253 & (1d0+10d0*rm1+10d0*rm2+rm1**2+rm2**2+10d0*rm1*rm2)
1257 ELSEIF(
mint(61).EQ.1.AND.icase.EQ.2)
THEN
1258 wdtp(
i)=fac*(vpi**2+api**2)*
vint(116)*wdtpzp
1259 ELSEIF(
mint(61).EQ.2)
THEN
1268 ELSEIF(
i.EQ.18)
THEN
1270 czc=2d0*(1d0-2d0*xw)
1271 be34c=(1d0-4d0*rm1)*sqrt(
max(0d0,1d0-4d0*rm1))
1273 wdtpz=0.25d0*paru(142)**2*czc**2*be34c
1274 wdtp(
i)=fac*0.25d0*paru(143)**2*czc**2*be34c
1275 ELSEIF(
mint(61).EQ.1.AND.icase.EQ.2)
THEN
1276 wdtp(
i)=fac*0.25d0*(
ei**2*
vint(111)+paru(142)*
ei*vi*
1277 &
vint(112)*czc+paru(143)*
ei*vpi*
vint(113)*czc+paru(142)**2*
1278 & (vi**2+ai**2)*
vint(114)*czc**2+paru(142)*paru(143)*
1279 & (vi*vpi+ai*api)*
vint(115)*czc**2+paru(143)**2*
1280 & (vpi**2+api**2)*
vint(116)*czc**2)*be34c
1281 ELSEIF(
mint(61).EQ.2)
THEN
1283 fgzf=0.25d0*paru(142)*czc*be34c
1284 fgzpf=0.25d0*paru(143)*czc*be34c
1285 fzzf=0.25d0*paru(142)**2*czc**2*be34c
1286 fzzpf=0.25d0*paru(142)*paru(143)*czc**2*be34c
1287 fzpzpf=0.25d0*paru(143)**2*czc**2*be34c
1290 ELSEIF(
i.EQ.19)
THEN
1292 ELSEIF(
i.EQ.20)
THEN
1294 flam=sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
1295 wdtpzp=paru(145)**2*4d0*abs(1d0-2d0*xw)*
1296 & (3d0*rm1+0.25d0*flam**2)*flam
1300 ELSEIF(
mint(61).EQ.1.AND.icase.EQ.2)
THEN
1301 wdtp(
i)=fac*(vpi**2+api**2)*
vint(116)*wdtpzp
1302 ELSEIF(
mint(61).EQ.2)
THEN
1310 wid2=wids(23,2)*wids(25,2)
1311 ELSEIF(
i.EQ.21.OR.
i.EQ.22)
THEN
1313 be34c=sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1323 wdtp(
i)=fac*czpah**2*be34c
1324 ELSEIF(
mint(61).EQ.1.AND.icase.EQ.2)
THEN
1325 wdtp(
i)=fac*(czah**2*(vi**2+ai**2)*
vint(114)+czah*czpah*
1326 & (vi*vpi+ai*api)*
vint(115)+czpah**2*(vpi**2+api**2)*
1328 ELSEIF(
mint(61).EQ.2)
THEN
1333 fzzpf=czah*czpah*be34c
1334 fzpzpf=czpah**2*be34c
1336 IF(
i.EQ.21) wid2=wids(25,2)*wids(36,2)
1337 IF(
i.EQ.22) wid2=wids(35,2)*wids(36,2)
1341 wdtp(
i)=fudge*wdtp(
i)
1342 wdtp(0)=wdtp(0)+wdtp(
i)
1344 IF(mdme(idc,1).GT.0)
THEN
1345 IF((icase.EQ.1.AND.
mint(61).NE.1).OR.
1346 & (icase.EQ.2.AND.
mint(61).EQ.1))
THEN
1347 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1348 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+
1349 & wdte(
i,mdme(idc,1))
1350 wdte(
i,0)=wdte(
i,mdme(idc,1))
1351 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1353 IF(
mint(61).EQ.2.AND.icase.EQ.2)
THEN
1369 IF(
mint(61).GE.1) icase=3-icase
1370 IF(icase.EQ.2) goto 280
1372 ELSEIF(kfla.EQ.34)
THEN
1374 fac=(aem/(24d0*xw))*shr
1375 DO 300
i=1,mdcy(kc,3)
1377 IF(mdme(idc,1).LT.0) goto 300
1378 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
1379 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
1380 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 300
1385 fcof=3d0*radc*(paru(131)**2+paru(132)**2)*
1386 & vckm((
i-1)/4+1,mod(
i-1,4)+1)
1388 IF(mod(
i,4).EQ.3) wid2=wids(6,2)
1389 IF(mod(
i,4).EQ.0) wid2=wids(8,2)
1390 IF(
i.GE.13) wid2=wid2*wids(7,3)
1392 IF(mod(
i,4).EQ.3) wid2=wids(6,3)
1393 IF(mod(
i,4).EQ.0) wid2=wids(8,3)
1394 IF(
i.GE.13) wid2=wid2*wids(7,2)
1396 ELSEIF(
i.LE.20)
THEN
1398 fcof=paru(133)**2+paru(134)**2
1400 IF(
i.EQ.20) wid2=wids(17,3)*wids(18,2)
1402 IF(
i.EQ.20) wid2=wids(17,2)*wids(18,3)
1405 wdtp(
i)=fac*fcof*0.5d0*(2d0-rm1-rm2-(rm1-rm2)**2)*
1406 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
1407 ELSEIF(
i.EQ.21)
THEN
1409 wdtp(
i)=fac*paru(135)**2*0.5d0*xw1*(rm1/rm2)*
1410 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1411 & (1d0+10d0*rm1+10d0*rm2+rm1**2+rm2**2+10d0*rm1*rm2)
1412 IF(kflr.GT.0) wid2=wids(24,2)*wids(23,2)
1413 IF(kflr.LT.0) wid2=wids(24,3)*wids(23,2)
1414 ELSEIF(
i.EQ.23)
THEN
1416 flam=sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
1417 wdtp(
i)=fac*paru(146)**2*2d0*(3d0*rm1+0.25d0*flam**2)*flam
1418 IF(kflr.GT.0) wid2=wids(24,2)*wids(25,2)
1419 IF(kflr.LT.0) wid2=wids(24,3)*wids(25,2)
1421 wdtp(
i)=fudge*wdtp(
i)
1422 wdtp(0)=wdtp(0)+wdtp(
i)
1423 IF(mdme(idc,1).GT.0)
THEN
1424 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1425 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
1426 wdte(
i,0)=wdte(
i,mdme(idc,1))
1427 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1431 ELSEIF(kfla.EQ.37)
THEN
1438 fac=(aem/(8d0*xw))*(shfs/pmas(24,1)**2)*shr
1439 DO 310
i=1,mdcy(kc,3)
1441 IF(mdme(idc,1).LT.0) goto 310
1444 rm1=pmas(kfc1,1)**2/sh
1445 rm2=pmas(kfc2,1)**2/sh
1446 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 310
1450 rm1r=
pymrun(kfdp(idc,1),sh)**2/sh
1451 rm2r=
pymrun(kfdp(idc,2),sh)**2/sh
1452 wdtp(
i)=fac*3d0*radc*
max(0d0,(rm1r*paru(141)**2+
1453 & rm2r/paru(141)**2)*(1d0-rm1r-rm2r)-4d0*rm1r*rm2r)*
1454 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*(sh/shfs)
1456 IF(
i.EQ.3) wid2=wids(6,2)
1457 IF(
i.EQ.4) wid2=wids(7,3)*wids(8,2)
1459 IF(
i.EQ.3) wid2=wids(6,3)
1460 IF(
i.EQ.4) wid2=wids(7,2)*wids(8,3)
1464 wdtp(
i)=fac*((rm1*paru(141)**2+rm2/paru(141)**2)*
1465 & (1d0-rm1-rm2)-4d0*rm1*rm2)*sqrt(
max(0d0,
1466 & (1d0-rm1-rm2)**2-4d0*rm1*rm2))*(sh/shfs)
1468 IF(
i.EQ.8) wid2=wids(17,3)*wids(18,2)
1470 IF(
i.EQ.8) wid2=wids(17,2)*wids(18,3)
1474 wdtp(
i)=fac*paru(195)**2*0.5d0*sqrt(
max(0d0,
1475 & (1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1476 IF(kflr.GT.0) wid2=wids(24,2)*wids(25,2)
1477 IF(kflr.LT.0) wid2=wids(24,3)*wids(25,2)
1484 wfac0=1d0+rm10**2+rm20**2-2d0*(rm10+rm20+rm10*rm20)
1485 wfac=1d0+rm1**2+rm2**2-2d0*(rm1+rm2+rm1*rm2)
1486 IF(wfac.LE.0d0 .OR. wfac0.LE.0d0)
THEN
1491 wdtp(
i)=pmas(kc,2)*brat(idc)*(shr/pmr)*sqrt(wfac)
1494 IF(kfls*kfdp(idc,1).LT.0.AND.kchg(kfc1,3).EQ.1) ksgn1=3
1496 IF(kfls*kfdp(idc,2).LT.0.AND.kchg(kfc2,3).EQ.1) ksgn2=3
1497 wid2=wids(kfc1,ksgn1)*wids(kfc2,ksgn2)
1499 wdtp(
i)=fudge*wdtp(
i)
1500 wdtp(0)=wdtp(0)+wdtp(
i)
1501 IF(mdme(idc,1).GT.0)
THEN
1502 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1503 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
1504 wdte(
i,0)=wdte(
i,mdme(idc,1))
1505 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1509 ELSEIF(kfla.EQ.41)
THEN
1511 fac=(aem/(12d0*xw))*shr
1512 DO 320
i=1,mdcy(kc,3)
1514 IF(mdme(idc,1).LT.0) goto 320
1515 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
1516 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
1517 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 320
1526 wdtp(
i)=fac*fcof*(2d0-rm1-rm2-(rm1-rm2)**2)*
1527 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
1529 IF(
i.EQ.4) wid2=wids(6,3)
1530 IF(
i.EQ.5) wid2=wids(7,3)
1531 IF(
i.EQ.6) wid2=wids(6,2)*wids(8,3)
1532 IF(
i.EQ.9) wid2=wids(17,3)
1534 IF(
i.EQ.4) wid2=wids(6,2)
1535 IF(
i.EQ.5) wid2=wids(7,2)
1536 IF(
i.EQ.6) wid2=wids(6,3)*wids(8,2)
1537 IF(
i.EQ.9) wid2=wids(17,2)
1539 wdtp(
i)=fudge*wdtp(
i)
1540 wdtp(0)=wdtp(0)+wdtp(
i)
1541 IF(mdme(idc,1).GT.0)
THEN
1542 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1543 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
1544 wdte(
i,0)=wdte(
i,mdme(idc,1))
1545 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1549 ELSEIF(kfla.EQ.42)
THEN
1551 fac=(aem/4d0)*paru(151)*shr
1552 DO 330
i=1,mdcy(kc,3)
1554 IF(mdme(idc,1).LT.0) goto 330
1555 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
1556 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
1557 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 330
1558 wdtp(
i)=fac*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1560 ilqq=kfdp(idc,1)*isign(1,kflr)
1561 IF(ilqq.GE.6) wid2=wids(ilqq,2)
1562 IF(ilqq.LE.-6) wid2=wids(-ilqq,3)
1563 ilql=kfdp(idc,2)*isign(1,kflr)
1564 IF(ilql.GE.17) wid2=wid2*wids(ilql,2)
1565 IF(ilql.LE.-17) wid2=wid2*wids(-ilql,3)
1566 wdtp(
i)=fudge*wdtp(
i)
1567 wdtp(0)=wdtp(0)+wdtp(
i)
1568 IF(mdme(idc,1).GT.0)
THEN
1569 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1570 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
1571 wdte(
i,0)=wdte(
i,mdme(idc,1))
1572 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1576 ELSEIF(kfla.EQ.ktechn+111.OR.kfla.EQ.ktechn+221)
THEN
1578 fac=(1d0/(32d0*paru(1)*rtcm(1)**2))*shr
1579 DO 340
i=1,mdcy(kc,3)
1581 IF(mdme(idc,1).LT.0) goto 340
1582 pm1=pmas(
pycomp(kfdp(idc,1)),1)
1583 pm2=pmas(
pycomp(kfdp(idc,2)),1)
1586 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 340
1590 facp=(as/(4d0*paru(1))*itcm(1)/rtcm(1))**2
1591 & /(8d0*paru(1))*sh*shr
1592 IF(kfla.EQ.ktechn+111)
THEN
1601 ika=iabs(kfdp(idc,1))
1602 IF(ika.LT.10) fcof=3d0*radc
1605 IF(ika.GE.4.AND.ika.LE.6)
THEN
1606 fcof=fcof*rtcm(1+ika)**2
1607 hm1=
pymrun(kfdp(idc,1),sh)
1608 hm2=
pymrun(kfdp(idc,2),sh)
1609 ELSEIF(ika.EQ.15)
THEN
1610 fcof=fcof*rtcm(8)**2
1612 wdtp(
i)=fac*fcof*(hm1+hm2)**2*
1613 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
1615 wdtp(
i)=fudge*wdtp(
i)
1616 wdtp(0)=wdtp(0)+wdtp(
i)
1617 IF(mdme(idc,1).GT.0)
THEN
1618 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1619 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
1620 wdte(
i,0)=wdte(
i,mdme(idc,1))
1621 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1625 ELSEIF(kfla.EQ.ktechn+211)
THEN
1627 fac=(1d0/(32d0*paru(1)*rtcm(1)**2))*shr
1628 DO 350
i=1,mdcy(kc,3)
1630 IF(mdme(idc,1).LT.0) goto 350
1631 pm1=pmas(
pycomp(kfdp(idc,1)),1)
1632 pm2=pmas(
pycomp(kfdp(idc,2)),1)
1634 IF(
i.EQ.5) pm3=pmas(
pycomp(kfdp(idc,3)),1)
1638 IF(sqrt(rm1)+sqrt(rm2)+sqrt(rm3).GT.1d0) goto 350
1642 IF(iabs(kfdp(idc,1)).LT.10) fcof=3d0*radc
1644 IF(
i.EQ.5.AND.shr.LT.pmas(6,1)+pmas(5,1))
THEN
1646 xmt2=pmas(6,1)**2/sh
1647 facp=fac/(4d0*paru(1))*fcof*xmt2*rtcm(7)**2
1649 check = sqrt(rm1)+sqrt(rm2)+sqrt(rm3)
1651 t0 = (1d0-check**2)*
1652 & (xmt2*(6d0*xmt2**2+3d0*xmt2*rm1-4d0*rm1**2)-
1653 & (5d0*xmt2**2+2d0*xmt2*rm1-8d0*rm1**2))/(4d0*xmt2**2)
1654 t1 = (1d0-xmt2)*(rm1-xmt2)*((xmt2**2+xmt2*rm1+4d0*rm1**2)
1655 & -3d0*xmt2**2*(xmt2+rm1))/(2d0*xmt2**3)
1656 t3 = rm1**2/xmt2**3*(3d0*xmt2-4d0*rm1+4d0*xmt2*rm1)
1657 wdtp(
i)=facp*(t0 +
t1*
log((xmt2-check**2)/(xmt2-1d0))
1666 ika=iabs(kfdp(idc,1))
1667 IF(ika.LT.10) fcof=3d0*radc
1670 IF(
i.GE.1.AND.
i.LE.5)
THEN
1672 fcof=fcof*rtcm(5)**2
1674 fcof=fcof*rtcm(6)**2
1676 fcof=fcof*rtcm(7)**2
1678 hm1=
pymrun(kfdp(idc,1),sh)
1679 hm2=
pymrun(kfdp(idc,2),sh)
1681 fcof=fcof*rtcm(8)**2
1683 wdtp(
i)=fac*fcof*(hm1+hm2)**2*
1684 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
1686 wdtp(
i)=fudge*wdtp(
i)
1687 wdtp(0)=wdtp(0)+wdtp(
i)
1688 IF(mdme(idc,1).GT.0)
THEN
1689 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1690 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
1691 wdte(
i,0)=wdte(
i,mdme(idc,1))
1692 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1696 ELSEIF(kfla.EQ.ktechn+331)
THEN
1698 fac=(sh/
parp(46)**2)*shr
1699 DO 360
i=1,mdcy(kc,3)
1701 IF(mdme(idc,1).LT.0) goto 360
1702 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
1703 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
1704 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 360
1707 wdtp(
i)=fac*rm1*sqrt(
max(0d0,1d0-4d0*rm1))/(4d0*paru(1))
1708 IF(
i.EQ.2) wid2=wids(6,1)
1710 wdtp(
i)=fac*5d0*as**2/(96d0*paru(1)**3)
1712 wdtp(
i)=fudge*wdtp(
i)
1713 wdtp(0)=wdtp(0)+wdtp(
i)
1714 IF(mdme(idc,1).GT.0)
THEN
1715 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1716 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
1717 wdte(
i,0)=wdte(
i,mdme(idc,1))
1718 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1722 ELSEIF(kfla.EQ.ktechn+113)
THEN
1724 alprht=2.91d0*(3d0/itcm(1))
1725 fac=(alprht/12d0)*shr
1726 facf=(1d0/6d0)*(aem**2/alprht)*shr
1730 CALL
pywidx(23,shp,wdtpp,wdtep)
1732 xwrht=(1d0-2d0*xw)/(4d0*xw*(1d0-xw))
1733 bwzr=xwrht*sh*(sh-sqmz)/((sh-sqmz)**2+gmmz**2)
1734 bwzi=xwrht*sh*gmmz/((sh-sqmz)**2+gmmz**2)
1735 DO 370
i=1,mdcy(kc,3)
1737 IF(mdme(idc,1).LT.0) goto 370
1738 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
1739 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
1740 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 370
1744 wdtp(
i)=fac*rtcm(3)**4*
1745 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1749 wdtp(
i)=fac*rtcm(3)**2*(1d0-rtcm(3)**2)*
1750 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3+
1751 & aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
1752 & ((1d0-rm1-rm2)**2-4d0*rm1*rm2 + 6d0*sqmw/sh)*
1753 & (1d0-rtcm(3)**2)/4d0/xw/24d0/rtcm(13)**2*shr**3
1754 wid2=wids(24,2)*wids(
pycomp(ktechn+211),3)
1757 wdtp(
i)=fac*rtcm(3)**2*(1d0-rtcm(3)**2)*
1758 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3+
1759 & aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
1760 & ((1d0-rm1-rm2)**2-4d0*rm1*rm2 + 6d0*sqmw/sh)*
1761 & (1d0-rtcm(3)**2)/4d0/xw/24d0/rtcm(13)**2*shr**3
1762 wid2=wids(
pycomp(ktechn+211),2)*wids(24,3)
1765 wdtp(
i)=fac*(1d0-rtcm(3)**2)**2*
1766 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1767 wid2=wids(
pycomp(ktechn+211),1)
1770 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1771 & (2d0*rtcm(2)-1d0)**2*(1d0-rtcm(3)**2)/24d0/rtcm(12)**2*
1773 wid2=wids(
pycomp(ktechn+111),2)
1776 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1777 & (1d0-rtcm(4)**2)/24d0/rtcm(12)**2*shr**3
1778 wid2=wids(
pycomp(ktechn+221),2)
1781 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1782 & (2d0*rtcm(2)-1d0)**2*(1d0-rtcm(3)**2)/24d0/rtcm(12)**2*
1784 wid2=wids(23,2)*wids(
pycomp(ktechn+111),2)
1787 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1788 & (1d0-rtcm(4)**2)/24d0/rtcm(12)**2*(1d0-2d0*xw)**2/4d0/
1790 wid2=wids(23,2)*wids(
pycomp(ktechn+221),2)
1797 IF(ia.GE.6.AND.ia.LE.8) wid2=wids(ia,1)
1801 IF(ia.GE.17) wid2=wids(ia,1)
1808 wdtp(
i)=facf*fcof*sqrt(
max(0d0,1d0-4d0*rm1))*((1d0-rm1)*
1809 & ((
ei+vali*bwzr)**2+(vali*bwzi)**2+
1810 & (
ei+vari*bwzr)**2+(vari*bwzi)**2)+6d0*rm1*(
1811 & (
ei+vali*bwzr)*(
ei+vari*bwzr)+vali*vari*bwzi**2))
1813 wdtp(
i)=fudge*wdtp(
i)
1814 wdtp(0)=wdtp(0)+wdtp(
i)
1815 IF(mdme(idc,1).GT.0)
THEN
1816 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1817 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
1818 wdte(
i,0)=wdte(
i,mdme(idc,1))
1819 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1823 ELSEIF(kfla.EQ.ktechn+213)
THEN
1825 alprht=2.91d0*(3d0/itcm(1))
1826 fac=(alprht/12d0)*shr
1830 CALL
pywidx(24,shp,wdtpp,wdtep)
1832 facf=(1d0/12d0)*(aem**2/alprht)*shr*
1833 & (0.125d0/xw**2)*sh**2/((sh-sqmw)**2+gmmw**2)
1834 DO 380
i=1,mdcy(kc,3)
1836 IF(mdme(idc,1).LT.0) goto 380
1837 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
1838 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
1839 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 380
1843 wdtp(
i)=fac*rtcm(3)**4*
1844 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1846 wid2=wids(24,2)*wids(23,2)
1848 wid2=wids(24,3)*wids(23,2)
1852 wdtp(
i)=fac*rtcm(3)**2*(1d0-rtcm(3)**2)*
1853 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3+
1854 & aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
1855 & ((1d0-rm1-rm2)**2-4d0*rm1*rm2 + 6d0*sqmw/sh)*
1856 & (1d0-rtcm(3)**2)/4d0/xw/24d0/rtcm(13)**2*shr**3
1858 wid2=wids(24,2)*wids(
pycomp(ktechn+111),2)
1860 wid2=wids(24,3)*wids(
pycomp(ktechn+111),2)
1864 wdtp(
i)=fac*rtcm(3)**2*(1d0-rtcm(3)**2)*
1865 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3+
1866 & aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))*
1867 & ((1d0-rm1-rm2)**2-4d0*rm1*rm2 + 6d0*sqmz/sh)*
1868 & (1d0-rtcm(3)**2)/4d0/xw/xw1/24d0/rtcm(13)**2*shr**3+
1869 & aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1870 & (2d0*rtcm(2)-1d0)**2*(1d0-rtcm(3)**2)/24d0/rtcm(12)**2*
1873 wid2=wids(
pycomp(ktechn+211),2)*wids(23,2)
1875 wid2=wids(
pycomp(ktechn+211),3)*wids(23,2)
1879 wdtp(
i)=fac*(1d0-rtcm(3)**2)**2*
1880 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1882 wid2=wids(
pycomp(ktechn+211),2)*wids(
pycomp(ktechn+111),2)
1884 wid2=wids(
pycomp(ktechn+211),3)*wids(
pycomp(ktechn+111),2)
1888 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1889 & (2d0*rtcm(2)-1d0)**2*(1d0-rtcm(3)**2)/24d0/rtcm(12)**2*
1892 wid2=wids(
pycomp(ktechn+211),2)
1894 wid2=wids(
pycomp(ktechn+211),3)
1898 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1899 & (1d0-rtcm(4)**2)/4d0/xw/24d0/rtcm(12)**2*shr**3
1901 wid2=wids(24,2)*wids(
pycomp(ktechn+221),2)
1903 wid2=wids(24,3)*wids(
pycomp(ktechn+221),2)
1910 fcof=3d0*radc*vckm((ia-1)/4+1,mod(ia-1,4)+1)
1912 IF(mod(ia,4).EQ.3) wid2=wids(6,2)
1913 IF(mod(ia,4).EQ.0) wid2=wids(8,2)
1914 IF(ia.GE.13) wid2=wid2*wids(7,3)
1916 IF(mod(ia,4).EQ.3) wid2=wids(6,3)
1917 IF(mod(ia,4).EQ.0) wid2=wids(8,3)
1918 IF(ia.GE.13) wid2=wid2*wids(7,2)
1923 IF(ia.EQ.20) wid2=wids(17,3)*wids(18,2)
1925 IF(ia.EQ.20) wid2=wids(17,2)*wids(18,3)
1928 wdtp(
i)=facf*fcof*(2d0-rm1-rm2-(rm1-rm2)**2)*
1929 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
1931 wdtp(
i)=fudge*wdtp(
i)
1932 wdtp(0)=wdtp(0)+wdtp(
i)
1933 IF(mdme(idc,1).GT.0)
THEN
1934 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
1935 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
1936 wdte(
i,0)=wdte(
i,mdme(idc,1))
1937 wdte(0,0)=wdte(0,0)+wdte(
i,0)
1941 ELSEIF(kfla.EQ.ktechn+223)
THEN
1943 alprht=2.91d0*(3d0/itcm(1))
1944 fac=(alprht/12d0)*shr
1945 facf=(1d0/6d0)*(aem**2/alprht)*shr*(2d0*rtcm(2)-1d0)**2
1948 CALL
pywidx(23,shp,wdtpp,wdtep)
1950 bwzr=(0.5d0/(1d0-xw))*sh*(sh-sqmz)/((sh-sqmz)**2+gmmz**2)
1951 bwzi=-(0.5d0/(1d0-xw))*sh*gmmz/((sh-sqmz)**2+gmmz**2)
1952 DO 390
i=1,mdcy(kc,3)
1954 IF(mdme(idc,1).LT.0) goto 390
1955 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
1956 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
1957 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 390
1961 wdtp(
i)=aem/24d0/rtcm(12)**2*(1d0-rtcm(3)**2)*
1962 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*shr**3
1963 wid2=wids(
pycomp(ktechn+111),2)
1966 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1967 & (1d0-rtcm(3)**2)/24d0/rtcm(12)**2*(1d0-2d0*xw)**2/4d0/
1969 wid2=wids(23,2)*wids(
pycomp(ktechn+111),2)
1972 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1973 & (2d0*rtcm(2)-1d0)**2*(1d0-rtcm(4)**2)/24d0/rtcm(12)**2*
1975 wid2=wids(
pycomp(ktechn+221),2)
1978 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1979 & (2d0*rtcm(2)-1d0)**2*(1d0-rtcm(4)**2)/24d0/rtcm(12)**2*
1981 wid2=wids(23,2)*wids(
pycomp(ktechn+221),2)
1984 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1985 & (1d0-rtcm(3)**2)/4d0/xw/24d0/rtcm(12)**2*shr**3+
1986 & fac*rtcm(3)**2*(1d0-rtcm(3)**2)*rtcm(11)**2*
1987 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1988 wid2=wids(24,2)*wids(
pycomp(ktechn+211),3)
1991 wdtp(
i)=aem*sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3*
1992 & (1d0-rtcm(3)**2)/4d0/xw/24d0/rtcm(12)**2*shr**3+
1993 & fac*rtcm(3)**2*(1d0-rtcm(3)**2)*rtcm(11)**2*
1994 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
1995 wid2=wids(24,3)*wids(
pycomp(ktechn+211),2)
1998 wdtp(
i)=fac*rtcm(3)**4*rtcm(11)**2*
1999 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
2003 wdtp(
i)=fac*(1d0-rtcm(3)**2)**2*rtcm(11)**2*
2004 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))**3
2005 wid2=wids(
pycomp(ktechn+211),1)
2012 IF(ia.GE.6.AND.ia.LE.8) wid2=wids(ia,1)
2016 IF(ia.GE.17) wid2=wids(ia,1)
2023 wdtp(
i)=facf*fcof*sqrt(
max(0d0,1d0-4d0*rm1))*((1d0-rm1)*
2024 & ((
ei+vali*bwzr)**2+(vali*bwzi)**2+
2025 & (
ei+vari*bwzr)**2+(vari*bwzi)**2)+6d0*rm1*(
2026 & (
ei+vali*bwzr)*(
ei+vari*bwzr)+vali*vari*bwzi**2))
2028 wdtp(
i)=fudge*wdtp(
i)
2029 wdtp(0)=wdtp(0)+wdtp(
i)
2030 IF(mdme(idc,1).GT.0)
THEN
2031 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2032 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2033 wdte(
i,0)=wdte(
i,mdme(idc,1))
2034 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2039 ELSEIF(kfla.EQ.ktechn+100021)
THEN
2042 IF(itcm(2).EQ.0)
THEN
2044 ELSEIF(itcm(2).EQ.1)
THEN
2047 DO 400
i=1,mdcy(kc,3)
2049 IF(mdme(idc,1).LT.0) goto 400
2050 pm1=pmas(
pycomp(kfdp(idc,1)),1)
2052 IF(rm1.GT.0.25d0) goto 400
2054 IF(
i.EQ.5.OR.
i.EQ.6.OR.imdl.EQ.2)
THEN
2059 wdtp(
i)=fac*(1d0+2d0*rm1)*sqrt(1d0-4d0*rm1)*fmix
2060 IF(
i.EQ.6) wid2=wids(6,1)
2061 wdtp(
i)=fudge*wdtp(
i)
2062 wdtp(0)=wdtp(0)+wdtp(
i)
2063 IF(mdme(idc,1).GT.0)
THEN
2064 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2065 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2066 wdte(
i,0)=wdte(
i,mdme(idc,1))
2067 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2071 ELSEIF(kfla.EQ.ktechn+100111.OR.kfla.EQ.ktechn+200111)
THEN
2072 fac=(1d0/(4d0*paru(1)*rtcm(1)**2))*shr
2074 DO 410
i=1,mdcy(kc,3)
2076 IF(mdme(idc,1).LT.0) goto 410
2077 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2078 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2079 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 410
2083 IF(kfla.EQ.ktechn+100111)
THEN
2088 facp=(as/(8d0*paru(1))*itcm(1)/rtcm(1))**2
2089 & /(2d0*paru(1))*sh*shr*clebg
2093 IF(
i.EQ.6) wid2=wids(6,1)
2095 ika=iabs(kfdp(idc,1))
2096 IF(ika.LT.10) fcof=3d0*radc
2097 hm1=
pymrun(kfdp(idc,1),sh)
2098 wdtp(
i)=fac*fcof*hm1**2*clebf*
2099 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
2101 wdtp(
i)=fudge*wdtp(
i)
2102 wdtp(0)=wdtp(0)+wdtp(
i)
2103 IF(mdme(idc,1).GT.0)
THEN
2104 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2105 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2106 wdte(
i,0)=wdte(
i,mdme(idc,1))
2107 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2111 ELSEIF(kfla.GE.ktechn+100113.AND.kfla.LE.ktechn+400113)
THEN
2113 alprht=2.91d0*(3d0/itcm(1))
2115 sin2t=2d0*tant3/(tant3**2+1d0)
2116 sint3=tant3/sqrt(tant3**2+1d0)
2119 x12=(rtcm(29)*sqrt(1d0-rtcm(29)**2)*cos(rtcm(30))+
2120 & rtcm(31)*sqrt(1d0-rtcm(31)**2)*cos(rtcm(32)))/sqrt(2d0)
2121 x21=(rtcm(29)*sqrt(1d0-rtcm(29)**2)*sin(rtcm(30))+
2122 & rtcm(31)*sqrt(1d0-rtcm(31)**2)*sin(rtcm(32)))/sqrt(2d0)
2123 x11=(.25d0*(rtcm(29)**2+rtcm(31)**2+2d0)-
2125 x22=(.25d0*(2d0-rtcm(29)**2-rtcm(31)**2)-
2127 CALL
pywidx(ktechn+100021,sh,wdtpp,wdtep)
2129 IF(wdtpp(0).GT.rtcm(33)*shr) wdtpp(0)=rtcm(33)*shr
2131 rmv8=pmas(
pycomp(ktechn+100021),1)
2132 fv8re=sh*(sh-rmv8**2)/((sh-rmv8**2)**2+gmv8**2)
2133 fv8im=sh*gmv8/((sh-rmv8**2)**2+gmv8**2)
2134 IF(itcm(2).EQ.0)
THEN
2139 DO 420
i=1,mdcy(kc,3)
2140 IF(
i.EQ.7.AND.(kfla.EQ.ktechn+200113.OR.
2141 & kfla.EQ.ktechn+300113)) goto 420
2143 IF(mdme(idc,1).LT.0) goto 420
2144 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2145 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2146 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 420
2149 IF(
i.EQ.6) wid2=wids(6,1)
2151 IF(kfla.EQ.ktechn+200113)
THEN
2154 ELSEIF(kfla.EQ.ktechn+300113)
THEN
2157 ELSEIF(kfla.EQ.ktechn+100113)
THEN
2162 IF(
i.EQ.5.OR.
i.EQ.6.OR.imdl.EQ.2)
THEN
2163 fmix=1d0/tant3/sin2t
2167 xfac=(xig+fmix*xij*fv8re)**2+(fmix*xij*fv8im)**2
2168 wdtp(
i)=fac*(1d0+2d0*rm1)*sqrt(1d0-4d0*rm1)*as/alprht*xfac
2170 wdtp(
i)=shr*as**2/(4d0*alprht)
2171 ELSEIF(kfla.EQ.ktechn+400113.AND.
i.LE.9)
THEN
2172 psh=shr*(1d0-rm1)/2d0
2173 wdtp(
i)=as/9d0*psh**3/rm82
2175 wdtp(
i)=2d0*wdtp(
i)*csxpp**2
2176 wid2=wids(
pycomp(kfdp(idc,1)),2)
2179 wid2=wids(
pycomp(kfdp(idc,1)),2)
2182 wdtp(
i)=fudge*wdtp(
i)
2183 wdtp(0)=wdtp(0)+wdtp(
i)
2184 IF(mdme(idc,1).GT.0)
THEN
2185 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2186 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2187 wdte(
i,0)=wdte(
i,mdme(idc,1))
2188 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2192 ELSEIF(kfla.EQ.kexcit+1)
THEN
2194 fac=(sh/rtcm(41)**2)*shr
2195 DO 430
i=1,mdcy(kc,3)
2197 IF(mdme(idc,1).LT.0) goto 430
2198 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2199 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2200 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 430
2204 wdtp(
i)=fac*as*rtcm(45)**2/3d0
2208 qf=-rtcm(43)/2d0+rtcm(44)/6d0
2209 wdtp(
i)=fac*aem*qf**2/4d0
2213 qf=-rtcm(43)*xw1/2d0-rtcm(44)*xw/6d0
2214 wdtp(
i)=fac*aem*qf**2/(8d0*xw*xw1)*
2215 & (1d0-rm1)**2*(2d0+rm1)
2219 wdtp(
i)=fac*aem*rtcm(43)**2/(16d0*xw)*
2220 & (1d0-rm1)**2*(2d0+rm1)
2221 IF(kflr.GT.0) wid2=wids(24,3)
2222 IF(kflr.LT.0) wid2=wids(24,2)
2224 wdtp(
i)=fudge*wdtp(
i)
2225 wdtp(0)=wdtp(0)+wdtp(
i)
2226 IF(mdme(idc,1).GT.0)
THEN
2227 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2228 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2229 wdte(
i,0)=wdte(
i,mdme(idc,1))
2230 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2234 ELSEIF(kfla.EQ.kexcit+2)
THEN
2236 fac=(sh/rtcm(41)**2)*shr
2237 DO 440
i=1,mdcy(kc,3)
2239 IF(mdme(idc,1).LT.0) goto 440
2240 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2241 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2242 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 440
2246 wdtp(
i)=fac*as*rtcm(45)**2/3d0
2250 qf=rtcm(43)/2d0+rtcm(44)/6d0
2251 wdtp(
i)=fac*aem*qf**2/4d0
2255 qf=rtcm(43)*xw1/2d0-rtcm(44)*xw/6d0
2256 wdtp(
i)=fac*aem*qf**2/(8d0*xw*xw1)*
2257 & (1d0-rm1)**2*(2d0+rm1)
2261 wdtp(
i)=fac*aem*rtcm(43)**2/(16d0*xw)*
2262 & (1d0-rm1)**2*(2d0+rm1)
2263 IF(kflr.GT.0) wid2=wids(24,2)
2264 IF(kflr.LT.0) wid2=wids(24,3)
2266 wdtp(
i)=fudge*wdtp(
i)
2267 wdtp(0)=wdtp(0)+wdtp(
i)
2268 IF(mdme(idc,1).GT.0)
THEN
2269 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2270 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2271 wdte(
i,0)=wdte(
i,mdme(idc,1))
2272 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2276 ELSEIF(kfla.EQ.kexcit+11)
THEN
2278 fac=(sh/rtcm(41)**2)*shr
2279 DO 450
i=1,mdcy(kc,3)
2281 IF(mdme(idc,1).LT.0) goto 450
2282 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2283 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2284 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 450
2288 qf=-rtcm(43)/2d0-rtcm(44)/2d0
2289 wdtp(
i)=fac*aem*qf**2/4d0
2293 qf=-rtcm(43)*xw1/2d0+rtcm(44)*xw/2d0
2294 wdtp(
i)=fac*aem*qf**2/(8d0*xw*xw1)*
2295 & (1d0-rm1)**2*(2d0+rm1)
2299 wdtp(
i)=fac*aem*rtcm(43)**2/(16d0*xw)*
2300 & (1d0-rm1)**2*(2d0+rm1)
2301 IF(kflr.GT.0) wid2=wids(24,3)
2302 IF(kflr.LT.0) wid2=wids(24,2)
2304 wdtp(
i)=fudge*wdtp(
i)
2305 wdtp(0)=wdtp(0)+wdtp(
i)
2306 IF(mdme(idc,1).GT.0)
THEN
2307 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2308 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2309 wdte(
i,0)=wdte(
i,mdme(idc,1))
2310 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2314 ELSEIF(kfla.EQ.kexcit+12)
THEN
2316 fac=(sh/rtcm(41)**2)*shr
2317 DO 460
i=1,mdcy(kc,3)
2319 IF(mdme(idc,1).LT.0) goto 460
2320 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2321 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2322 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 460
2326 qf=rtcm(43)*xw1/2d0+rtcm(44)*xw/2d0
2327 wdtp(
i)=fac*aem*qf**2/(8d0*xw*xw1)*
2328 & (1d0-rm1)**2*(2d0+rm1)
2332 wdtp(
i)=fac*aem*rtcm(43)**2/(16d0*xw)*
2333 & (1d0-rm1)**2*(2d0+rm1)
2334 IF(kflr.GT.0) wid2=wids(24,2)
2335 IF(kflr.LT.0) wid2=wids(24,3)
2337 wdtp(
i)=fudge*wdtp(
i)
2338 wdtp(0)=wdtp(0)+wdtp(
i)
2339 IF(mdme(idc,1).GT.0)
THEN
2340 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2341 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2342 wdte(
i,0)=wdte(
i,mdme(idc,1))
2343 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2347 ELSEIF(kfla.EQ.kdimen+39)
THEN
2349 fac=(
parp(50)**2/paru(1))*shr
2350 DO 470
i=1,mdcy(kc,3)
2352 IF(mdme(idc,1).LT.0) goto 470
2353 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2354 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2355 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 470
2360 IF(
i.GE.6.AND.
mstp(35).GE.1) fcof=fcof*
2362 wdtp(
i)=fac*fcof*sqrt(
max(0d0,1d0-4d0*rm1))**3*
2363 & (1d0+8d0*rm1/3d0)/320d0
2364 IF(
i.EQ.6) wid2=wids(6,1)
2365 IF(
i.EQ.7.OR.
i.EQ.8) wid2=wids(
i,1)
2366 ELSEIF(
i.LE.16)
THEN
2369 wdtp(
i)=fac*sqrt(
max(0d0,1d0-4d0*rm1))**3*
2370 & (1d0+8d0*rm1/3d0)/320d0
2371 IF(
i.EQ.15.OR.
i.EQ.16) wid2=wids(2+
i,1)
2372 ELSEIF(
i.EQ.17)
THEN
2375 ELSEIF(
i.EQ.18)
THEN
2378 ELSEIF(
i.EQ.19)
THEN
2380 wdtp(
i)=fac*sqrt(
max(0d0,1d0-4d0*rm1))*(13d0/12d0+
2381 & 14d0*rm1/3d0+4d0*rm1**2)/160d0
2383 ELSEIF(
i.EQ.20)
THEN
2385 wdtp(
i)=fac*sqrt(
max(0d0,1d0-4d0*rm1))*(13d0/12d0+
2386 & 14d0*rm1/3d0+4d0*rm1**2)/80d0
2389 wdtp(
i)=fudge*wdtp(
i)
2390 wdtp(0)=wdtp(0)+wdtp(
i)
2391 IF(mdme(idc,1).GT.0)
THEN
2392 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2393 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2394 wdte(
i,0)=wdte(
i,mdme(idc,1))
2395 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2399 ELSEIF(kfla.EQ.9900012.OR.kfla.EQ.9900014.OR.kfla.EQ.9900016)
THEN
2401 pmwr=
max(1.001d0*shr,pmas(
pycomp(9900024),1))
2402 fac=(aem**2/(768d0*paru(1)*xw**2))*shr**5/pmwr**4
2403 DO 480
i=1,mdcy(kc,3)
2405 IF(mdme(idc,1).LT.0) goto 480
2406 pm1=pmas(
pycomp(kfdp(idc,1)),1)
2407 pm2=pmas(
pycomp(kfdp(idc,2)),1)
2408 pm3=pmas(
pycomp(kfdp(idc,3)),1)
2409 IF(pm1+pm2+pm3.GE.shr) goto 480
2413 fcof=3d0*radc*vckm((
i-1)/3+1,mod(
i-1,3)+1)
2414 IF(mod(
i,3).EQ.0) wid2=wids(6,2)
2415 ELSEIF(
i.LE.18)
THEN
2417 fcof=3d0*radc*vckm((
i-10)/3+1,mod(
i-10,3)+1)
2418 IF(mod(
i-9,3).EQ.0) wid2=wids(6,3)
2422 wid2=wids(
pycomp(kfdp(idc,3)),2)
2425 fx=1d0-8d0*
x**2+8d0*
x**6-
x**8-24d0*
x**4*
log(
x)
2427 fy=(12d0*(1d0-
y)*
log(1d0-
y)+12d0*
y-6d0*
y**2-2d0*
y**3)/
y**4
2428 wdtp(
i)=fac*fcof*fx*fy
2429 wdtp(
i)=fudge*wdtp(
i)
2430 wdtp(0)=wdtp(0)+wdtp(
i)
2431 IF(mdme(idc,1).GT.0)
THEN
2432 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2433 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2434 wdte(
i,0)=wdte(
i,mdme(idc,1))
2435 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2439 ELSEIF(kfla.EQ.9900023)
THEN
2441 fac=(aem/(48d0*xw*xw1*(1d0-2d0*xw)))*shr
2442 DO 490
i=1,mdcy(kc,3)
2444 IF(mdme(idc,1).LT.0) goto 490
2445 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2446 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2447 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 490
2453 af=
sign(1d0,ef+0.1d0)*(1d0-2d0*xw)
2454 vf=
sign(1d0,ef+0.1d0)-4d0*ef*xw
2456 IF(
i.EQ.6) wid2=wids(6,1)
2457 ELSEIF(
i.EQ.7.OR.
i.EQ.10.OR.
i.EQ.13)
THEN
2462 ELSEIF(
i.EQ.8.OR.
i.EQ.11.OR.
i.EQ.14)
THEN
2468 ELSEIF(
i.LE.15)
THEN
2473 wid2=wids(
pycomp(kfdp(idc,1)),1)
2476 wdtp(
i)=fac*fcof*(vf**2*(1d0+2d0*rm1)+af**2*(1d0-4d0*rm1))*
2477 & sqrt(
max(0d0,1d0-4d0*rm1))*symmet
2478 wdtp(
i)=fudge*wdtp(
i)
2479 wdtp(0)=wdtp(0)+wdtp(
i)
2480 IF(mdme(idc,1).GT.0)
THEN
2481 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2482 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2483 wdte(
i,0)=wdte(
i,mdme(idc,1))
2484 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2488 ELSEIF(kfla.EQ.9900024)
THEN
2490 fac=(aem/(24d0*xw))*shr
2491 DO 500
i=1,mdcy(kc,3)
2493 IF(mdme(idc,1).LT.0) goto 500
2494 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2495 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2496 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 500
2500 fcof=3d0*radc*vckm((
i-1)/3+1,mod(
i-1,3)+1)
2502 IF(mod(
i,3).EQ.0) wid2=wids(6,2)
2504 IF(mod(
i,3).EQ.0) wid2=wids(6,3)
2506 ELSEIF(
i.LE.12)
THEN
2510 wdtp(
i)=fac*fcof*(2d0-rm1-rm2-(rm1-rm2)**2)*
2511 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
2512 wdtp(
i)=fudge*wdtp(
i)
2513 wdtp(0)=wdtp(0)+wdtp(
i)
2514 IF(mdme(idc,1).GT.0)
THEN
2515 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2516 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2517 wdte(
i,0)=wdte(
i,mdme(idc,1))
2518 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2522 ELSEIF(kfla.EQ.9900041)
THEN
2524 fac=(1d0/(8d0*paru(1)))*shr
2525 DO 510
i=1,mdcy(kc,3)
2527 IF(mdme(idc,1).LT.0) goto 510
2528 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2529 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2530 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 510
2534 fcof=
parp(180+3*((iabs(kfdp(idc,1))-11)/2)+
2535 & (iabs(kfdp(idc,2))-9)/2)**2
2536 IF(kfdp(idc,1).NE.kfdp(idc,2)) fcof=2d0*fcof
2539 fcof=0.5d0*
parp(190)**4*
parp(192)**2/pmas(24,1)**2*
2540 & (3d0*rm1+0.25d0/rm1-1d0)
2541 wid2=wids(24,4+(1-kfls)/2)
2544 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
2545 wdtp(
i)=fudge*wdtp(
i)
2546 wdtp(0)=wdtp(0)+wdtp(
i)
2547 IF(mdme(idc,1).GT.0)
THEN
2548 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2549 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2550 wdte(
i,0)=wdte(
i,mdme(idc,1))
2551 wdte(0,0)=wdte(0,0)+wdte(
i,0)
2555 ELSEIF(kfla.EQ.9900042)
THEN
2557 fac=(1d0/(8d0*paru(1)))*shr
2558 DO 520
i=1,mdcy(kc,3)
2560 IF(mdme(idc,1).LT.0) goto 520
2561 rm1=pmas(
pycomp(kfdp(idc,1)),1)**2/sh
2562 rm2=pmas(
pycomp(kfdp(idc,2)),1)**2/sh
2563 IF(sqrt(rm1)+sqrt(rm2).GT.1d0) goto 520
2567 fcof=
parp(180+3*((iabs(kfdp(idc,1))-11)/2)+
2568 & (iabs(kfdp(idc,2))-9)/2)**2
2569 IF(kfdp(idc,1).NE.kfdp(idc,2)) fcof=2d0*fcof
2572 fcof=
parp(191)**2*(3d0*rm1+0.25d0/rm1-1d0)
2573 wid2=wids(
pycomp(9900024),4+(1-kfls)/2)
2576 & sqrt(
max(0d0,(1d0-rm1-rm2)**2-4d0*rm1*rm2))
2577 wdtp(
i)=fudge*wdtp(
i)
2578 wdtp(0)=wdtp(0)+wdtp(
i)
2579 IF(mdme(idc,1).GT.0)
THEN
2580 wdte(
i,mdme(idc,1))=wdtp(
i)*wid2
2581 wdte(0,mdme(idc,1))=wdte(0,mdme(idc,1))+wdte(
i,mdme(idc,1))
2582 wdte(
i,0)=wdte(
i,mdme(idc,1))
2583 wdte(0,0)=wdte(0,0)+wdte(
i,0)