17 SUBROUTINE pycjdc(KFIN,XLAM,IDLAM,IKNT)
20 IMPLICIT DOUBLE PRECISION(a-
h, o-
z)
24 parameter(ksusy1=1000000,ksusy2=2000000,ktechn=3000000,
25 &kexcit=4000000,kdimen=5000000)
27 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
28 common/
pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
29 common/
pymssm/imss(0:99),rmss(0:99)
30 common/
pyssmt/zmix(4,4),umix(2,2),vmix(2,2),smz(4),smw(2),
31 &sfmix(16,4),zmixi(4,4),umixi(2,2),vmixi(2,2)
35 common/pyintc/xxc(10),cxc(8)
39 COMPLEX*16 zmixc(4,4),vmixc(2,2),umixc(2,2),olpp,orpp
40 COMPLEX*16 cal,cbl,car,cbr,ca,cb
42 DOUBLE PRECISION xmi,xmj,xmf,xmsf1,xmsf2,xmw,xmw2,
44 DOUBLE PRECISION s12min,s12max
45 DOUBLE PRECISION xmi2,xmi3,xmj2,xmh,xmh2,xmhp,xma2,xmb2,xmk
47 DOUBLE PRECISION tanw,xw,aem,c1,as,
ei,t3i,
beta,alfa
49 DOUBLE PRECISION xlam(0:400)
51 INTEGER lknt,ix,ih,
j,ij,
i,iknt
54 DOUBLE PRECISION etah(3),dh(3),eh(3)
56 DOUBLE PRECISION cbeta,sbeta,tanb
61 INTEGER kfnchi(4),kfcchi(2)
69 DATA etah/1d0,1d0,-1d0/
71 DATA pi/3.141592654d0/
73 DATA kfnchi/1000022,1000023,1000025,1000035/
74 DATA kfcchi/1000024,1000037/
88 IF(kfin.EQ.kfcchi(2)) ix=2
106 vmixc(
j,
i)=dcmplx(vmix(
j,
i),vmixi(
j,
i))
107 umixc(
j,
i)=dcmplx(umix(
j,
i),umixi(
j,
i))
113 IF(imss(11).EQ.1)
THEN
119 xfac=(xmi2/(xmp*xmgr))**2*axmi/48d0/
pi
120 IF(axmi.GT.xmgr+xmw)
THEN
126 & .5d0*(abs(vmixc(ix,1))**2+abs(umixc(ix,1))**2)+
127 & .5d0*((abs(vmixc(ix,2))*sbeta)**2+(abs(umixc(ix,2))*cbeta)**2))*
130 IF(axmi.GT.xmgr+pmas(37,1))
THEN
135 xlam(lknt)=xfac*(.5d0*((abs(vmixc(ix,2))*cbeta)**2+
136 & (abs(umixc(ix,2))*sbeta)**2))
137 & *(1d0-pmas(37,1)**2/xmi2)**4
148 IF(axmi.GE.axmj+xmz)
THEN
151 olpp=-vmixc(ij,1)*dconjg(vmixc(ix,1))-
152 & vmixc(ij,2)*dconjg(vmixc(ix,2))/2d0
153 orpp=-umixc(ix,1)*dconjg(umixc(ij,1))-
154 & umixc(ix,2)*dconjg(umixc(ij,2))/2d0
155 gx2=abs(olpp)**2+abs(orpp)**2
156 glr=dble(olpp*dconjg(orpp))
157 xlam(lknt)=
pyx2xg(c1/xmw2,xmi,xmj,xmz,gx2,glr)
158 idlam(lknt,1)=kfcchi(1)
163 ELSEIF(axmi.GE.axmj)
THEN
165 s12max=(axmi-axmj)**2
168 ei=kchg(iabs(ia),1)/3d0
174 xxc(5)=pmas(
pycomp(ksusy1+ja),1)
179 olpp=-vmixc(ij,1)*dconjg(vmixc(ix,1))-
180 & vmixc(ij,2)*dconjg(vmixc(ix,2))/2d0
181 orpp=-umixc(ix,1)*dconjg(umixc(ij,1))-
182 & umixc(ix,2)*dconjg(umixc(ij,2))/2d0
183 cxc(1)=dcmplx((t3i-xw*
ei)/xw/xw1)*orpp
184 cxc(2)=dcmplx(0d0,0d0)
185 cxc(3)=dcmplx((t3i-xw*
ei)/xw/xw1)*olpp
186 cxc(4)=-vmixc(ij,1)*dconjg(vmixc(ix,1))*dcmplx(t3i/xw)
187 cxc(5)=-dcmplx(
ei/xw1)*orpp
188 cxc(6)=dcmplx(0d0,0d0)
189 cxc(7)=-dcmplx(
ei/xw1)*olpp
190 cxc(8)=dcmplx(0d0,0d0)
191 IF( xxc(5).LT.axmi )
THEN
196 IF(axmi.GE.axmj+2d0*pmas(11,1))
THEN
198 xlam(lknt)=c1**2/xmi3/(16d0*
pi)*
200 idlam(lknt,1)=kfcchi(1)
203 IF(axmi.GE.axmj+2d0*pmas(13,1))
THEN
205 xlam(lknt)=xlam(lknt-1)
206 idlam(lknt,1)=kfcchi(1)
210 IF(axmi.GE.axmj+2d0*pmas(15,1))
THEN
212 xlam(lknt)=xlam(lknt-1)
213 idlam(lknt,1)=kfcchi(1)
223 ei=kchg(iabs(ia),1)/3d0
225 xxc(5)=pmas(
pycomp(ksusy1+ja),1)
227 cxc(1)=dcmplx((t3i-xw*
ei)/xw/xw1)*orpp
228 cxc(3)=dcmplx((t3i-xw*
ei)/xw/xw1)*olpp
229 cxc(4)=-umixc(ij,1)*dconjg(umixc(ix,1))*dcmplx(t3i/xw)
230 cxc(5)=-dcmplx(
ei/xw1)*orpp
231 cxc(7)=-dcmplx(
ei/xw1)*olpp
232 IF( xxc(5).LT.axmi )
THEN
237 IF(axmi.GE.axmj+2d0*pmas(12,1))
THEN
239 xlam(lknt)=c1**2/xmi3/(16d0*
pi)*
241 idlam(lknt,1)=kfcchi(1)
245 xlam(lknt)=xlam(lknt-1)
246 idlam(lknt,1)=kfcchi(1)
250 IF(axmi.GE.axmj+2d0*pmas(16,1))
THEN
251 IF(abs(sfmix(15,1)).GT.abs(sfmix(15,2)))
THEN
252 xxc(5)=pmas(
pycomp(ksusy1+15),1)
254 xxc(5)=pmas(
pycomp(ksusy2+15),1)
256 IF( xxc(5).LT.axmi )
THEN
261 xlam(lknt)=c1**2/xmi3/(16d0*
pi)*
263 idlam(lknt,1)=kfcchi(1)
272 ei=kchg(iabs(ia),1)/3d0
274 xxc(5)=pmas(
pycomp(ksusy1+ja),1)
276 cxc(1)=dcmplx((t3i-xw*
ei)/xw/xw1)*orpp
277 cxc(2)=dcmplx(0d0,0d0)
278 cxc(3)=dcmplx((t3i-xw*
ei)/xw/xw1)*olpp
279 cxc(4)=-vmixc(ij,1)*dconjg(vmixc(ix,1))*dcmplx(t3i/xw)
280 cxc(5)=-dcmplx(
ei/xw1)*orpp
281 cxc(6)=dcmplx(0d0,0d0)
282 cxc(7)=-dcmplx(
ei/xw1)*olpp
283 cxc(8)=dcmplx(0d0,0d0)
284 IF( xxc(5).LT.axmi )
THEN
289 IF(axmi.GE.axmj+2d0*pmas(1,1))
THEN
291 xlam(lknt)=3d0*c1**2/xmi3/(16d0*
pi)*
293 idlam(lknt,1)=kfcchi(1)
296 IF(axmi.GE.axmj+2d0*pmas(3,1))
THEN
298 xlam(lknt)=xlam(lknt-1)
299 idlam(lknt,1)=kfcchi(1)
304 IF(axmi.GE.axmj+2d0*pmas(5,1))
THEN
305 IF(abs(sfmix(5,1)).GT.abs(sfmix(5,2)))
THEN
306 xxc(5)=pmas(
pycomp(ksusy1+5),1)
308 xxc(5)=pmas(
pycomp(ksusy2+5),1)
310 IF( xxc(5).LT.axmi )
THEN
315 xlam(lknt)=3d0*c1**2/xmi3/(16d0*
pi)*
317 idlam(lknt,1)=kfcchi(1)
326 ei=kchg(iabs(ia),1)/3d0
328 xxc(5)=pmas(
pycomp(ksusy1+ja),1)
330 cxc(1)=dcmplx((t3i-xw*
ei)/xw/xw1)*orpp
331 cxc(2)=dcmplx(0d0,0d0)
332 cxc(3)=dcmplx((t3i-xw*
ei)/xw/xw1)*olpp
333 cxc(4)=-umixc(ij,1)*dconjg(umixc(ix,1))*dcmplx(t3i/xw)
334 cxc(5)=-dcmplx(
ei/xw1)*orpp
335 cxc(6)=dcmplx(0d0,0d0)
336 cxc(7)=-dcmplx(
ei/xw1)*olpp
337 cxc(8)=dcmplx(0d0,0d0)
338 IF( xxc(5).LT.axmi )
THEN
343 IF(axmi.GE.axmj+2d0*pmas(2,1))
THEN
345 xlam(lknt)=3d0*c1**2/xmi3/(16d0*
pi)*
347 idlam(lknt,1)=kfcchi(1)
350 IF(axmi.GE.axmj+2d0*pmas(4,1))
THEN
352 xlam(lknt)=xlam(lknt-1)
353 idlam(lknt,1)=kfcchi(1)
372 IF(axmi.GE.axmj+xmh)
THEN
375 olpp=(vmixc(2,1)*dconjg(umixc(1,2))*eh(ih) -
376 & vmixc(2,2)*dconjg(umixc(1,1))*dh(ih))/sr2
377 orpp=(dconjg(vmixc(1,1))*umixc(2,2)*eh(ih) -
378 & dconjg(vmixc(1,2))*umixc(2,1)*dh(ih))/sr2
380 gx2=abs(olpp)**2+abs(orpp)**2
381 glr=dble(olpp*dconjg(orpp))
382 xlam(lknt)=
pyx2xh(c1,xmi,xmk,xmh,gx2,glr)
383 idlam(lknt,1)=kfcchi(1)
384 idlam(lknt,2)=ith(ih)
397 IF(axmi.GE.axmj+xmw)
THEN
400 zmixc(ij,
i)=dcmplx(zmix(ij,
i),zmixi(ij,
i))
402 cxc(1)=(dconjg(zmixc(ij,2))*vmixc(ix,1)-
403 & dconjg(zmixc(ij,4))*vmixc(ix,2)/sr2)
404 cxc(3)=(zmixc(ij,2)*dconjg(umixc(ix,1))+
405 & zmixc(ij,3)*dconjg(umixc(ix,2))/sr2)
406 gx2=abs(cxc(1))**2+abs(cxc(3))**2
407 glr=dble(cxc(1)*dconjg(cxc(3)))
408 xlam(lknt)=
pyx2xg(c1/xmw2,xmi,xmj,xmw,gx2,glr)
409 idlam(lknt,1)=kfnchi(ij)
413 ELSEIF(axmi.GE.axmj)
THEN
415 s12max=(axmi-axmj)**2
417 zmixc(ij,
i)=dcmplx(zmix(ij,
i),zmixi(ij,
i))
419 cxc(1)=(dconjg(zmixc(ij,2))*vmixc(ix,1)-
420 & dconjg(zmixc(ij,4))*vmixc(ix,2)/sr2)/sr2
421 cxc(3)=(zmixc(ij,2)*dconjg(umixc(ix,1))+
422 & zmixc(ij,3)*dconjg(umixc(ix,2))/sr2)/sr2
423 cxc(5)=dcmplx(0d0,0d0)
424 cxc(7)=dcmplx(0d0,0d0)
430 t3j=
sign(1d0,ej+1d-6)/2d0
431 cxc(2)=vmixc(ix,1)*dconjg(zmixc(ij,1)*(ej-t3j)*
432 & tanw+zmixc(ij,2)*t3j)/sr2
433 cxc(4)=-dconjg(umixc(ix,1))*(
434 & zmixc(ij,1)*(
ei-t3i)*tanw+zmixc(ij,2)*t3i)/sr2
435 cxc(6)=dcmplx(0d0,0d0)
436 cxc(8)=dcmplx(0d0,0d0)
441 xxc(5)=pmas(
pycomp(ksusy1+ja),1)
442 xxc(6)=pmas(
pycomp(ksusy1+ia),1)
446 IF(xxc(5).LT.axmi)
THEN
448 ELSEIF(xxc(6).LT.axmi)
THEN
455 IF(axmi.GE.axmj+pmas(11,1)+pmas(12,1))
THEN
458 xlam(lknt)=c1**2/xmi3/(16d0*
pi)*temp
459 idlam(lknt,1)=kfnchi(ij)
463 IF( imss(12).NE. 0 ) goto 260
464 IF(axmi.GE.axmj+pmas(13,1)+pmas(14,1))
THEN
466 xlam(lknt)=xlam(lknt-1)
467 idlam(lknt,1)=kfnchi(ij)
472 IF(axmi.GE.axmj+pmas(15,1)+pmas(16,1))
THEN
474 IF(abs(sfmix(15,1)).GT.abs(sfmix(15,2)))
THEN
475 xxc(6)=pmas(
pycomp(ksusy1+15),1)
477 xxc(6)=pmas(
pycomp(ksusy2+15),1)
479 xxc(5)=pmas(
pycomp(ksusy1+16),1)
480 IF(xxc(5).LT.axmi)
THEN
482 ELSEIF(xxc(6).LT.axmi)
THEN
488 xlam(lknt)=c1**2/xmi3/(16d0*
pi)*temp
489 idlam(lknt,1)=kfnchi(ij)
501 t3j=
sign(1d0,ej+1d-6)/2d0
502 cxc(2)=vmixc(ix,1)*dconjg(zmixc(ij,1)*(ej-t3j)*
503 & tanw+zmixc(ij,2)*t3j)
504 cxc(4)=-dconjg(umixc(ix,1))*(
505 & zmixc(ij,1)*(
ei-t3i)*tanw+zmixc(ij,2)*t3i)
506 xxc(5)=pmas(
pycomp(ksusy1+ja),1)
507 xxc(6)=pmas(
pycomp(ksusy1+ia),1)
508 IF( xxc(5).LT.axmi .AND. xxc(6).LT.axmi ) goto 210
509 IF(xxc(5).LT.axmi)
THEN
512 IF(xxc(6).LT.axmi)
THEN
517 IF(axmi.GE.axmj+pmas(1,1)+pmas(2,1))
THEN
519 xlam(lknt)=3d0*c1**2/xmi3/(16d0*
pi)*
521 idlam(lknt,1)=kfnchi(ij)
524 IF(axmi.GE.axmj+pmas(3,1)+pmas(4,1))
THEN
526 xlam(lknt)=xlam(lknt-1)
527 idlam(lknt,1)=kfnchi(ij)
542 IF(axmi.GE.axmj+xmhp)
THEN
544 olpp=cbeta*(zmixc(ij,4)*dconjg(vmixc(ix,1))+(zmixc(ij,2)+
545 & zmixc(ij,1)*tanw)*dconjg(vmixc(ix,2))/sr2)
546 orpp=sbeta*(dconjg(zmixc(ij,3))*umixc(ix,1)-
547 & (dconjg(zmixc(ij,2))+dconjg(zmixc(ij,1))*tanw)*
549 gx2=abs(olpp)**2+abs(orpp)**2
550 glr=dble(olpp*dconjg(orpp))
551 xlam(lknt)=
pyx2xh(c1,xmi,xmj,xmhp,gx2,glr)
552 idlam(lknt,1)=kfnchi(ij)
562 IF(
j.GE.7.AND.
j.LE.10) goto 240
563 IF(mod(
j,2).EQ.0)
THEN
579 IF(mod(
j,2).EQ.0)
THEN
582 cbl=-xmf*vmixc(ix,2)/xmw/sbeta/sr2
583 car=-xmfp*umixc(ix,2)/xmw/cbeta/sr2
589 cbl=-xmf*umixc(ix,2)/xmw/cbeta/sr2
591 car=-xmfp*vmixc(ix,2)/xmw/sbeta/sr2
596 IF(axmi.GE.xmf+xmsf1)
THEN
601 ca=cal*sfmix(isf,1)+car*sfmix(isf,2)
602 cb=cbl*sfmix(isf,1)+cbr*sfmix(isf,2)
603 xlam(lknt)=fcol*c1/8d0/xmi3*sqrt(xl)*( (xmi2+xmb2-xma2)*
604 & (abs(ca)**2+abs(cb)**2)+4d0*dble(ca*dconjg(cb))*xmf*xmi)
606 IF(mod(
j,2).EQ.0)
THEN
616 IF(axmi.GE.xmf+xmsf2)
THEN
620 ca=cal*sfmix(isf,3)+car*sfmix(isf,4)
621 cb=cbl*sfmix(isf,3)+cbr*sfmix(isf,4)
623 xlam(lknt)=fcol*c1/8d0/xmi3*sqrt(xl)*( (xmi2+xmb2-xma2)*
624 & (abs(ca)**2+abs(cb)**2)+4d0*dble(ca*dconjg(cb))*xmf*xmi)
626 IF(mod(
j,2).EQ.0)
THEN
638 xmj=pmas(
pycomp(ksusy1+21),1)
642 s12max=(axmi-axmj)**2
647 xxc(5)=pmas(
pycomp(ksusy1+1),1)
648 xxc(6)=pmas(
pycomp(ksusy1+2),1)
651 olpp=dcmplx(cos(rmss(32)),sin(rmss(32)))
653 cxc(1)=dcmplx(0d0,0d0)
654 cxc(3)=dcmplx(0d0,0d0)
655 cxc(5)=dcmplx(0d0,0d0)
656 cxc(7)=dcmplx(0d0,0d0)
657 cxc(2)=umixc(ix,1)*olpp/sr2
658 cxc(4)=-dconjg(vmixc(ix,1))*orpp/sr2
659 cxc(6)=dcmplx(0d0,0d0)
660 cxc(8)=dcmplx(0d0,0d0)
661 IF(xxc(5).LT.axmi)
THEN
663 ELSEIF(xxc(6).LT.axmi)
THEN
668 IF( xxc(5).LT.axmi .OR. xxc(6).LT.axmi ) goto 250
669 IF(axmi.GE.axmj+pmas(1,1)+pmas(2,1))
THEN
671 xlam(lknt)=4d0*c1*as/xmi3/(16d0*
pi)*
673 idlam(lknt,1)=ksusy1+21
676 IF(axmi.GE.axmj+pmas(3,1)+pmas(4,1))
THEN
678 xlam(lknt)=xlam(lknt-1)
679 idlam(lknt,1)=ksusy1+21
688 CALL
pyrvch(kfin,xlam,idlam,lknt)
693 xlam(0)=xlam(0)+xlam(
i)
694 IF(xlam(
i).LT.0d0)
THEN
695 WRITE(mstu(11),*)
' XLAM(I) = ',xlam(
i),kcin,
700 IF(xlam(0).EQ.0d0)
THEN
702 WRITE(mstu(11),*)
' XLAM(0) = ',xlam(0)
703 WRITE(mstu(11),*) lknt
704 WRITE(mstu(11),*) (xlam(
j),
j=1,lknt)