7 SUBROUTINE pysfdc(KFIN,XLAM,IDLAM,IKNT)
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/
pymssm/imss(0:99),rmss(0:99)
20 common/
pyssmt/zmix(4,4),umix(2,2),vmix(2,2),smz(4),smw(2),
21 &sfmix(16,4),zmixi(4,4),umixi(2,2),vmixi(2,2)
25 COMPLEX*16 zmixc(4,4),vmixc(2,2),umixc(2,2)
26 COMPLEX*16 cal,car,cbl,cbr,calp,carp,cblp,cbrp,ca,cb
28 DOUBLE PRECISION xmi,xmj,xmf,xmsf1,xmsf2,xmw,xmw2,xmz,axmj
29 DOUBLE PRECISION xmi2,xmi3,xma2,xmb2,xmfp
31 DOUBLE PRECISION tanw,xw,aem,c1,as
32 DOUBLE PRECISION al,ar,bl,br
33 DOUBLE PRECISION ch1,ch2,ch3,ch4
34 DOUBLE PRECISION xmbot,xmtop
35 DOUBLE PRECISION xlam(0:400)
37 INTEGER lknt,ix,ilr,idu,
j,
i,iknt,ifl,ii
39 DOUBLE PRECISION cbeta,sbeta
41 DOUBLE PRECISION beta,alfa,xmu,at,ab,atrit,atrib,atril
42 DOUBLE PRECISION cosa,sina,tanb
44 DOUBLE PRECISION ghrr,ghll,ghlr,xmb,blr
46 INTEGER igg(4),kfnchi(4),kfcchi(2)
48 DATA pi/3.141592654d0/
50 DATA kfnchi/1000022,1000023,1000025,1000035/
51 DATA kfcchi/1000024,1000037/
57 IF(kfin.EQ.ksusy2+12.OR.kfin.EQ.ksusy2+14.OR.
58 &kfin.EQ.ksusy2+16)
RETURN
64 tanw = sqrt(xw/(1d0-xw))
69 zmixc(
j,
i)=dcmplx(zmix(
j,
i),zmixi(
j,
i))
74 vmixc(
j,
i)=dcmplx(vmix(
j,
i),vmixi(
j,
i))
75 umixc(
j,
i)=dcmplx(umix(
j,
i),umixi(
j,
i))
113 IF(imss(11).EQ.1)
THEN
117 xfac=(xmi2/(xmp*xmgr))**2*xmi/48d0/
pi
120 ELSEIF(ifl.EQ.6)
THEN
125 IF(xmi.GT.xmgr+xmf)
THEN
130 xlam(lknt)=xfac*(1d0-xmf**2/xmi2)**4
149 IF(xmi.GE.axmj+xmfp)
THEN
156 ELSEIF(ifl.LT.6)
THEN
161 cal=-xmfp*umixc(ix,2)/sr2/xmw/cbeta
162 cbr=-xmf*vmixc(ix,2)/sr2/xmw/sbeta
168 ELSEIF(ifl.LT.5)
THEN
173 cal=-xmfp*vmixc(ix,2)/sr2/xmw/sbeta
174 cbr=-xmf*umixc(ix,2)/sr2/xmw/cbeta
178 calp=sfmix(ifl,1)*cal + sfmix(ifl,2)*car
179 cblp=sfmix(ifl,1)*cbl + sfmix(ifl,2)*cbr
180 carp=sfmix(ifl,4)*car + sfmix(ifl,3)*cal
181 cbrp=sfmix(ifl,4)*cbr + sfmix(ifl,3)*cbl
199 xlam(lknt)=2d0*c1/8d0/xmi3*sqrt(xl)*((xmi2-xmb2-xma2)*
200 & (abs(ca)**2+abs(cb)**2)-4d0*dble(ca*dconjg(cb))*xmj*xmfp)
203 idlam(lknt,1)=-kfcchi(ix)
206 idlam(lknt,1)=kfcchi(ix)
218 IF(xmi.GE.axmj+xmf)
THEN
224 ELSEIF(ifl.LT.5)
THEN
227 cbl=-zmixc(ix,2)+tanw*zmixc(ix,1)*(2d0*
ei+1)
228 cal=xmf*zmixc(ix,3)/xmw/cbeta
229 car=-2d0*
ei*tanw*zmixc(ix,1)
234 ELSEIF(ifl.LT.5)
THEN
237 cbl=zmixc(ix,2)+tanw*zmixc(ix,1)*(2d0*
ei-1)
238 cal=xmf*zmixc(ix,4)/xmw/sbeta
239 car=-2d0*
ei*tanw*zmixc(ix,1)
243 calp=sfmix(ifl,1)*cal + sfmix(ifl,2)*car
244 cblp=sfmix(ifl,1)*cbl + sfmix(ifl,2)*cbr
245 carp=sfmix(ifl,4)*car + sfmix(ifl,3)*cal
246 cbrp=sfmix(ifl,4)*cbr + sfmix(ifl,3)*cbl
264 xlam(lknt)=c1/8d0/xmi3*sqrt(xl)*((xmi2-xmb2-xma2)*
265 & (abs(ca)**2+abs(cb)**2)-4d0*dble(ca*dconjg(cb))*xmj*xmf)
266 idlam(lknt,1)=kfnchi(ix)
276 IF(ilr.EQ.1) goto 160
278 xmsf1=pmas(
pycomp(kfin-ksusy1),1)
279 IF(xmi.LT.xmsf1+xmb) goto 160
284 ELSEIF(ig.EQ.25)
THEN
287 ELSEIF(ifl.EQ.6)
THEN
289 ELSEIF(ifl.LT.5)
THEN
295 ghll=xmz/cw*(0.5d0-
ei*xw)*(-sin(alfa+
beta))+
296 & xmf**2/xmw*cosa/sbeta
297 ghrr=xmz/cw*(
ei*xw)*(-sin(alfa+
beta))+
298 & xmf**2/xmw*cosa/sbeta
300 ghll=xmz/cw*(0.5d0-
ei*xw)*(-sin(alfa+
beta))+
301 & xmf**2/xmw*(-sina)/cbeta
302 ghrr=xmz/cw*(
ei*xw)*(-sin(alfa+
beta))+
303 & xmf**2/xmw*(-sina)/cbeta
307 ELSEIF(ifl.EQ.6)
THEN
309 ELSEIF(ifl.EQ.15)
THEN
316 ghlr=xmf/2d0/xmw/sbeta*(-xmu*sina+
319 ghlr=xmf/2d0/xmw/cbeta*(xmu*cosa-
325 ELSEIF(ig.EQ.35)
THEN
328 ELSEIF(ifl.EQ.6)
THEN
330 ELSEIF(ifl.LT.5)
THEN
336 ghll=xmz/cw*(0.5d0-
ei*xw)*cos(alfa+
beta)+
337 & xmf**2/xmw*sina/sbeta
338 ghrr=xmz/cw*(
ei*xw)*cos(alfa+
beta)+
339 & xmf**2/xmw*sina/sbeta
341 ghll=xmz/cw*(0.5d0-
ei*xw)*cos(alfa+
beta)+
342 & xmf**2/xmw*cosa/cbeta
343 ghrr=xmz/cw*(
ei*xw)*cos(alfa+
beta)+
344 & xmf**2/xmw*cosa/cbeta
348 ELSEIF(ifl.EQ.6)
THEN
350 ELSEIF(ifl.EQ.15)
THEN
357 ghlr=xmf/2d0/xmw/sbeta*(xmu*cosa+
360 ghlr=xmf/2d0/xmw/cbeta*(xmu*sina+
366 ELSEIF(ig.EQ.36)
THEN
371 ELSEIF(ifl.EQ.6)
THEN
373 ELSEIF(ifl.LT.5)
THEN
380 ELSEIF(ifl.EQ.6)
THEN
382 ELSEIF(ifl.EQ.15)
THEN
389 ghlr=xmf/2d0/xmw*(-xmu+at/tanb)
391 ghlr=xmf/2d0/xmw/(-xmu+at*tanb)
397 al=sfmix(ifl,1)*sfmix(ifl,3)*bl+
398 & sfmix(ifl,2)*sfmix(ifl,4)*br+
399 & (sfmix(ifl,1)*sfmix(ifl,4)+sfmix(ifl,3)*sfmix(ifl,2))*blr
400 xl=
pylamf(xmi2,xmsf1**2,xmb**2)
403 xlam(lknt)=c1/4d0/xmi3*xl**1.5d0/xmb**2*al**2
405 xlam(lknt)=c1/4d0/xmi3*sqrt(xl)*al**2
408 idlam(lknt,1)=kfin-ksusy1
414 IF(mod(ifl,2).EQ.0)
THEN
422 IF(xmi.GT.xmb+xmsf1)
THEN
423 IF(mod(ifl,2).EQ.0)
THEN
425 al=1d0/sr2*sfmix(ifl,1)*sfmix(ifl-1,1)
427 al=1d0/sr2*sfmix(ifl,3)*sfmix(ifl-1,1)
431 al=1d0/sr2*sfmix(ifl,1)*sfmix(ifl+1,1)
433 al=1d0/sr2*sfmix(ifl,3)*sfmix(ifl+1,1)
436 xl=
pylamf(xmi2,xmsf1**2,xmb**2)
438 xlam(lknt)=c1/4d0/xmi3*xl**1.5d0/xmb**2*al**2
441 idlam(lknt,2)=
sign(24,kchg(ifl,1))
443 IF(xmi.GT.xmb+xmsf2)
THEN
444 IF(mod(ifl,2).EQ.0)
THEN
446 al=1d0/sr2*sfmix(ifl,1)*sfmix(ifl-1,3)
448 al=1d0/sr2*sfmix(ifl,3)*sfmix(ifl-1,3)
452 al=1d0/sr2*sfmix(ifl,1)*sfmix(ifl+1,3)
454 al=1d0/sr2*sfmix(ifl,3)*sfmix(ifl+1,3)
457 xl=
pylamf(xmi2,xmsf2**2,xmb**2)
459 xlam(lknt)=c1/4d0/xmi3*xl**1.5d0/xmb**2*al**2
462 idlam(lknt,2)=
sign(24,kchg(ifl,1))
467 IF(mod(ifl,2).EQ.0)
THEN
475 IF(xmi.GT.xmb+xmsf1)
THEN
480 IF(mod(ifl,2).EQ.0)
THEN
483 ch1=-sfmix(ifl,1)*sfmix(ifl-1,1)
484 ch2= sfmix(ifl,2)*sfmix(ifl-1,2)
485 ch3=-sfmix(ifl,1)*sfmix(ifl-1,2)
486 ch4=-sfmix(ifl,2)*sfmix(ifl-1,1)
489 ch1= sfmix(ifl,3)*sfmix(ifl-1,1)
490 ch2=-sfmix(ifl,4)*sfmix(ifl-1,2)
491 ch3= sfmix(ifl,3)*sfmix(ifl-1,2)
492 ch4= sfmix(ifl,4)*sfmix(ifl-1,1)
503 ch1=-sfmix(ifl+1,1)*sfmix(ifl,1)
504 ch2= sfmix(ifl+1,2)*sfmix(ifl,2)
505 ch3=-sfmix(ifl+1,1)*sfmix(ifl,2)
506 ch4=-sfmix(ifl+1,2)*sfmix(ifl,1)
509 ch1= sfmix(ifl,3)*sfmix(ifl+1,1)
510 ch2=-sfmix(ifl,4)*sfmix(ifl+1,2)
511 ch3= sfmix(ifl,4)*sfmix(ifl+1,1)
512 ch4= sfmix(ifl,3)*sfmix(ifl+1,2)
521 xl=
pylamf(xmi2,xmsf1**2,xmb**2)
524 al=ch1*(xmw2*2d0*cbeta*sbeta-xmfp**2*tanb-xmf**2/tanb)+
525 & ch2*2d0*xmf*xmfp/(2d0*cbeta*sbeta)+
526 & ch3*xmfp*(-xmu+ab*tanb)+ch4*xmf*(-xmu+at/tanb)
527 xlam(lknt)=c1/8d0/xmi3*sqrt(xl)/xmw2*al**2
530 idlam(lknt,2)=
sign(37,kchg(ifl,1))
532 IF(xmi.GT.xmb+xmsf2)
THEN
537 IF(mod(ifl,2).EQ.0)
THEN
540 ch1= sfmix(ifl-1,3)*sfmix(ifl,1)
541 ch2=-sfmix(ifl-1,4)*sfmix(ifl,2)
542 ch3= sfmix(ifl-1,4)*sfmix(ifl,1)
543 ch4= sfmix(ifl-1,3)*sfmix(ifl,2)
546 ch1= -sfmix(ifl,3)*sfmix(ifl-1,3)
547 ch2= sfmix(ifl,4)*sfmix(ifl-1,4)
548 ch3= -sfmix(ifl,3)*sfmix(ifl-1,4)
549 ch4= -sfmix(ifl,4)*sfmix(ifl-1,3)
560 ch1= sfmix(ifl+1,3)*sfmix(ifl,1)
561 ch2=-sfmix(ifl+1,4)*sfmix(ifl,2)
562 ch3= sfmix(ifl+1,3)*sfmix(ifl,2)
563 ch4= sfmix(ifl+1,4)*sfmix(ifl,1)
566 ch1= -sfmix(ifl+1,3)*sfmix(ifl,3)
567 ch2= sfmix(ifl+1,4)*sfmix(ifl,4)
568 ch3= -sfmix(ifl+1,3)*sfmix(ifl,4)
569 ch4= -sfmix(ifl+1,4)*sfmix(ifl,3)
578 xl=
pylamf(xmi2,xmsf1**2,xmb**2)
581 al=ch1*(xmw2*2d0*cbeta*sbeta-xmfp**2*tanb-xmf**2/tanb)+
582 & ch2*2d0*xmf*xmfp/(2d0*cbeta*sbeta)+
583 & ch3*xmfp*(-xmu+ab*tanb)+ch4*xmf*(-xmu+at/tanb)
584 xlam(lknt)=c1/8d0/xmi3*sqrt(xl)/xmw2*al**2
587 idlam(lknt,2)=
sign(37,kchg(ifl,1))
595 IF(ifl.EQ.6) xmf=pmas(6,1)
596 IF(ifl.EQ.5) xmf=pmas(5,1)
597 xmj=pmas(
pycomp(ksusy1+21),1)
599 IF(xmi.GE.axmj+xmf)
THEN
617 xlam(lknt)=4d0/3d0*as/2d0/xmi3*sqrt(xl)*((xmi2-xmb2-xma2)*
618 & (xca**2+xcb**2)+4d0*xca*xcb*xmj*xmf)
619 idlam(lknt,1)=ksusy1+21
626 IF(kfin.EQ.ksusy1+6.AND.pmas(kcin,1).GT.
627 &pmas(
pycomp(ksusy1+22),1)+pmas(4,1))
THEN
634 xlam(lknt)=c1**3/64d0/
pi**2/xmi3*sqrt(xl)
635 IF(xlam(lknt).EQ.0) xlam(lknt)=1d-3
636 idlam(lknt,1)=ksusy1+22
642 CALL
pyrvsf(kfin,xlam,idlam,lknt)
647 IF(xlam(
i).LT.0d0) xlam(
i)=0d0
648 xlam(0)=xlam(0)+xlam(
i)
650 IF(xlam(0).EQ.0d0) xlam(0)=1d-3