7 SUBROUTINE pyglui(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 common/pyintc/xxc(10),cxc(8)
29 COMPLEX*16 zmixc(4,4),vmixc(2,2),umixc(2,2),olpp,orpp,glij,grij
30 DOUBLE PRECISION xmi,xmj,xmf,axmj,axmi
31 DOUBLE PRECISION xmi2,xmi3,xma2,xmb2,xmfp
33 DOUBLE PRECISION tanw,xw,aem,c1,as,s12max,s12min
34 DOUBLE PRECISION ca,cb,al,ar,bl,br
35 DOUBLE PRECISION xlam(0:400)
37 INTEGER lknt,ix,ilr,
i,iknt,ifl
44 INTEGER kfnchi(4),kfcchi(2)
45 DATA pi/3.141592654d0/
48 DATA kfnchi/1000022,1000023,1000025,1000035/
49 DATA kfcchi/1000024,1000037/
53 IF(kfin.NE.ksusy1+21)
RETURN
57 tanw = sqrt(xw/(1d0-xw))
71 IF(imss(11).EQ.1)
THEN
75 xfac=(xmi2/(xmp*xmgr))**2*axmi/48d0/
pi
89 xmj=pmas(
pycomp(ilr*ksusy1+ifl),1)
92 IF(axmi.GE.axmj+xmf)
THEN
111 xlam(lknt)=4d0/8d0*as/4d0/xmi3*sqrt(xl)*((xmi2+xmb2-xma2)*
112 & (ca**2+cb**2)-4d0*ca*cb*xmi*xmf)
113 idlam(lknt,1)=ilr*ksusy1+ifl
117 xlam(lknt)=xlam(lknt-1)
118 idlam(lknt,1)=-idlam(lknt-1,1)
119 idlam(lknt,2)=-idlam(lknt-1,2)
130 IF(axmi.GE.axmj)
THEN
132 zmixc(ix,
i)=dcmplx(zmix(ix,
i),zmixi(ix,
i))
134 olpp=dcmplx(cos(rmss(32)),sin(rmss(32)))/sr2
141 xxc(5)=pmas(
pycomp(ksusy1+ia),1)
142 xxc(6)=pmas(
pycomp(ksusy2+ia),1)
149 glij=(t3i*zmixc(ix,2)-tanw*(t3i-
ei)*zmixc(ix,1))*olpp
150 grij=zmixc(ix,1)*(
ei*tanw)*orpp
160 s12max=(axmi-axmj)**2
161 IF( xxc(5).LT.axmi .OR. xxc(6).LT.axmi ) goto 130
162 IF(axmi.GE.axmj+2d0*pmas(1,1))
THEN
164 xlam(lknt)=c1*as/xmi3/(16d0*
pi)*
166 idlam(lknt,1)=kfnchi(ix)
170 IF(axmi.GE.axmj+2d0*pmas(3,1))
THEN
172 xlam(lknt)=xlam(lknt-1)
173 idlam(lknt,1)=kfnchi(ix)
178 IF(axmi.GE.axmj+2d0*pmas(5,1))
THEN
179 pmold=pmas(
pycomp(ksusy1+5),1)
180 IF(axmi.GT.pmas(
pycomp(ksusy2+5),1)+pmas(5,1))
THEN
182 ELSEIF(axmi.GT.pmas(
pycomp(ksusy1+5),1)+pmas(5,1))
THEN
183 pmas(
pycomp(ksusy1+5),1)=100d0*xmi
185 CALL
pytbbn(ix,100,-1d0/3d0,xmi,gam)
188 idlam(lknt,1)=kfnchi(ix)
191 pmas(
pycomp(ksusy1+5),1)=pmold
196 xxc(5)=pmas(
pycomp(ksusy1+ia),1)
197 xxc(6)=pmas(
pycomp(ksusy2+ia),1)
203 glij=(t3i*zmixc(ix,2)-tanw*(t3i-
ei)*zmixc(ix,1))*olpp
204 grij=zmixc(ix,1)*(
ei*tanw)*orpp
209 IF( xxc(5).LT.axmi .OR. xxc(6).LT.axmi ) goto 150
210 IF(axmi.GE.axmj+2d0*pmas(2,1))
THEN
212 xlam(lknt)=c1*as/xmi3/(16d0*
pi)*
214 idlam(lknt,1)=kfnchi(ix)
218 IF(axmi.GE.axmj+2d0*pmas(4,1))
THEN
220 xlam(lknt)=xlam(lknt-1)
221 idlam(lknt,1)=kfnchi(ix)
229 IF(axmi.GE.axmj+2d0*xmf)
THEN
230 pmold=pmas(
pycomp(ksusy1+6),1)
231 IF(axmi.GT.pmas(
pycomp(ksusy2+6),1)+xmf)
THEN
233 ELSEIF(axmi.GT.pmas(
pycomp(ksusy1+6),1)+xmf)
THEN
234 pmas(
pycomp(ksusy1+6),1)=100d0*xmi
236 CALL
pytbbn(ix,100,2d0/3d0,xmi,gam)
239 idlam(lknt,1)=kfnchi(ix)
242 pmas(
pycomp(ksusy1+6),1)=pmold
252 IF(axmi.GE.axmj)
THEN
254 vmixc(ix,
i)=dcmplx(vmix(ix,
i),vmixi(ix,
i))
255 umixc(ix,
i)=dcmplx(umix(ix,
i),umixi(ix,
i))
258 s12max=(axmi-axmj)**2
263 xxc(5)=pmas(
pycomp(ksusy1+1),1)
264 xxc(6)=pmas(
pycomp(ksusy1+2),1)
267 olpp=dcmplx(cos(rmss(32)),sin(rmss(32)))
269 cxc(1)=dcmplx(0d0,0d0)
270 cxc(3)=dcmplx(0d0,0d0)
271 cxc(5)=dcmplx(0d0,0d0)
272 cxc(7)=dcmplx(0d0,0d0)
273 cxc(2)=umixc(ix,1)*olpp/sr2
274 cxc(4)=-dconjg(vmixc(ix,1))*orpp/sr2
275 cxc(6)=dcmplx(0d0,0d0)
276 cxc(8)=dcmplx(0d0,0d0)
277 IF(xxc(5).LT.axmi)
THEN
279 ELSEIF(xxc(6).LT.axmi)
THEN
284 IF( xxc(5).LT.axmi .OR. xxc(6).LT.axmi ) goto 190
285 IF(axmi.GE.axmj+pmas(1,1)+pmas(2,1))
THEN
287 xlam(lknt)=0.5d0*c1*as/xmi3/(16d0*
pi)*
289 idlam(lknt,1)=kfcchi(ix)
293 xlam(lknt)=xlam(lknt-1)
294 idlam(lknt,1)=-idlam(lknt-1,1)
295 idlam(lknt,2)=-idlam(lknt-1,2)
296 idlam(lknt,3)=-idlam(lknt-1,3)
298 IF(axmi.GE.axmj+pmas(3,1)+pmas(4,1))
THEN
300 xlam(lknt)=xlam(lknt-1)
301 idlam(lknt,1)=kfcchi(ix)
305 xlam(lknt)=xlam(lknt-1)
306 idlam(lknt,1)=-idlam(lknt-1,1)
307 idlam(lknt,2)=-idlam(lknt-1,2)
308 idlam(lknt,3)=-idlam(lknt-1,3)
314 IF(axmi.GE.axmj+xmf+xmfp)
THEN
315 IF(xmi.GT.min(pmas(
pycomp(ksusy1+5),1)+xmfp,
316 $ pmas(
pycomp(ksusy2+6),1)+xmf)) goto 200
317 pmolt2=pmas(
pycomp(ksusy2+6),1)
318 pmolb2=pmas(
pycomp(ksusy2+5),1)
319 pmolt1=pmas(
pycomp(ksusy1+6),1)
320 pmolb1=pmas(
pycomp(ksusy1+5),1)
321 IF(xmi.GT.pmolt2+xmf) pmas(
pycomp(ksusy2+6),1)=100d0*axmi
322 IF(xmi.GT.pmolt1+xmf) pmas(
pycomp(ksusy1+6),1)=100d0*axmi
323 IF(xmi.GT.pmolb2+xmfp) pmas(
pycomp(ksusy2+5),1)=100d0*axmi
324 IF(xmi.GT.pmolb1+xmfp) pmas(
pycomp(ksusy1+5),1)=100d0*axmi
325 CALL
pytbbc(ix,100,xmi,gam)
328 idlam(lknt,1)=kfcchi(ix)
332 xlam(lknt)=xlam(lknt-1)
333 idlam(lknt,1)=-idlam(lknt-1,1)
334 idlam(lknt,2)=-idlam(lknt-1,2)
335 idlam(lknt,3)=-idlam(lknt-1,3)
336 pmas(
pycomp(ksusy2+6),1)=pmolt2
337 pmas(
pycomp(ksusy2+5),1)=pmolb2
338 pmas(
pycomp(ksusy1+6),1)=pmolt1
339 pmas(
pycomp(ksusy1+5),1)=pmolb1
346 CALL
pyrvgl(kfin,xlam,idlam,lknt)
351 IF(xlam(
i).LT.0d0) xlam(
i)=0d0
352 xlam(0)=xlam(0)+xlam(
i)
354 IF(xlam(0).EQ.0d0) xlam(0)=1d-6