11 IMPLICIT DOUBLE PRECISION(a-
h, o-
z)
15 parameter(ksusy1=1000000,ksusy2=2000000,ktechn=3000000,
16 &kexcit=4000000,kdimen=5000000)
18 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
19 common/
pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
20 common/
pymssm/imss(0:99),rmss(0:99)
21 common/
pyssmt/zmix(4,4),umix(2,2),vmix(2,2),smz(4),smw(2),
22 &sfmix(16,4),zmixi(4,4),umixi(2,2),vmixi(2,2)
26 DOUBLE PRECISION xmw,xmz,xm(4)
27 DOUBLE PRECISION ar(4,4),wr(4),zr(4,4),zi(4,4),ai(4,4)
28 DOUBLE PRECISION wi(4),fv1(4),fv2(4),fv3(4)
29 DOUBLE PRECISION cosw,sinw
31 DOUBLE PRECISION tanb,cosb,sinb
32 DOUBLE PRECISION xm1,xm2,xm3,
beta
33 DOUBLE PRECISION q2,aem,a1,a2,aq,rm1,rm2
34 DOUBLE PRECISION arg,
x0,x1,ax0,ax1,at,bt
35 DOUBLE PRECISION y0,y1,amgx0,am1x0,amgx1,am1x1
36 DOUBLE PRECISION argx0,ar1x0,argx1,ar1x1
39 COMPLEX*16 car(4,4),cai(4,4),ca1,ca2
40 INTEGER ierr,
index(4),
i,
j,
k,iopt,ilr,kfnchi(4)
41 DATA kfnchi/1000022,1000023,1000025,1000035/
52 ELSEIF(iopt.GE.1)
THEN
56 a1=aem/(1d0-paru(102))
59 IF(imss(1).EQ.2) xm1=rmss(1)/rmss(20)*a1*5d0/3d0
63 ELSEIF(iopt.EQ.3)
THEN
70 WRITE(mstu(11),*)
' ERROR WITH M3 = ',xm3
77 pmas(
pycomp(ksusy1+21),1)=abs(xm3)
82 rm1=pmas(
pycomp(ilr*ksusy1+
i),1)**2/xm3**2
83 aq=aq+0.5d0*((2d0-rm1)*(rm1*
log(rm1)-1d0)
84 & +(1d0-rm1)**2*
log(abs(1d0-rm1)))
90 rm1=pmas(
pycomp(ilr*ksusy1+
i),1)**2/xm3**2
91 rm2=pmas(
i,1)**2/xm3**2
92 arg=(rm1-rm2-1d0)**2-4d0*rm2**2
94 x0=0.5d0*(1d0+rm2-rm1-sqrt(arg))
96 x1=0.5d0*(1d0+rm2-rm1+sqrt(arg))
101 ELSEIF(
x0.EQ.0d0)
THEN
105 at=0.5d0*
log(abs(1d0-
x0))*(1d0-
x0**2)+
106 & 0.5d0*
x0**2*
log(ax0)
112 ELSEIF(x1.EQ.0d0)
THEN
116 at=0.5d0*
log(abs(1d0-x1))*(1d0-x1**2)+0.5d0*
118 bt=(-1d0-2d0*x1)/4d0+bt
122 x0=0.5d0*(1d0+rm2-rm1)
124 amgx0=sqrt(
x0**2+y0**2)
125 am1x0=sqrt((1d0-
x0)**2+y0**2)
127 ar1x0=atan2(1d0-
x0,y0)
133 ar1x1=atan2(1d0-x1,y1)
134 at=0.5d0*
log(am1x0)*(1d0-
x0**2+3d0*y0**2)
135 & +0.5d0*(
x0**2-y0**2)*
log(amgx0)
136 bt=(-1d0-2d0*
x0)/4d0+
x0*y0*( ar1x0-argx0 )
137 at=at+0.5d0*
log(am1x1)*(1d0-x1**2+3d0*y1**2)
138 & +0.5d0*(x1**2-y1**2)*
log(amgx1)
139 bt=bt+(-1d0-2d0*x1)/4d0+x1*y1*( ar1x1-argx1 )
145 & /(2d0*paru(2))*(15d0+aq))
158 cosw=sqrt(1d0-paru(102))
167 ar(1,1) = xm1*cos(rmss(30))
168 ai(1,1) = xm1*sin(rmss(30))
169 ar(2,2) = xm2*cos(rmss(31))
170 ai(2,2) = xm2*sin(rmss(31))
175 ar(1,3) = -xmz*sinw*cosb
177 ar(1,4) = xmz*sinw*sinb
179 ar(2,3) = xmz*cosw*cosb
181 ar(2,4) = -xmz*cosw*sinb
183 ar(3,4) = -xmu*cos(rmss(33))
184 ai(3,4) = -xmu*sin(rmss(33))
185 ar(4,3) = -xmu*cos(rmss(33))
186 ai(4,3) = -xmu*sin(rmss(33))
188 CALL
pyeicg(4,4,ar,ai,wr,wi,1,zr,zi,fv1,fv2,fv3,ierr)
190 WRITE(mstu(11),*)
' PROBLEM WITH PYEICG IN PYINOM '
199 IF(xm(
k).LT.xm(
j))
THEN
217 pmas(
pycomp(kfnchi(
i)),1)=abs(smz(
i))
223 zmix(
i,
j)=zr(
j,
k)/sqrt(
s)
224 zmixi(
i,
j)=zi(
j,
k)/sqrt(
s)
225 IF(abs(zmix(
i,
j)).LT.1d-6) zmix(
i,
j)=0d0
226 IF(abs(zmixi(
i,
j)).LT.1d-6) zmixi(
i,
j)=0d0
234 ar(1,1) = xm2**2+2d0*xmw**2*sinb**2
235 ar(2,2) = xmu**2+2d0*xmw**2*cosb**2
236 ar(1,2) = sqrt(2d0)*xmw*(xm2*cos(rmss(31))*cosb+
237 &xmu*cos(rmss(33))*sinb)
238 ai(1,2) = sqrt(2d0)*xmw*(xm2*sin(rmss(31))*cosb-
239 &xmu*sin(rmss(33))*sinb)
240 ar(2,1) = sqrt(2d0)*xmw*(xm2*cos(rmss(31))*cosb+
241 &xmu*cos(rmss(33))*sinb)
242 ai(2,1) = sqrt(2d0)*xmw*(-xm2*sin(rmss(31))*cosb+
243 &xmu*sin(rmss(33))*sinb)
244 CALL
pyeicg(4,2,ar,ai,wr,wi,1,zr,zi,fv1,fv2,fv3,ierr)
246 WRITE(mstu(11),*)
' PROBLEM WITH PYEICG IN PYINOM '
250 IF(wr(2).LT.wr(1))
THEN
263 umix(
i,
j)=zr(
j,
k)/sqrt(
s)
264 umixi(
i,
j)=-zi(
j,
k)/sqrt(
s)
265 IF(abs(umix(
i,
j)).LT.1d-6) umix(
i,
j)=0d0
266 IF(abs(umixi(
i,
j)).LT.1d-6) umixi(
i,
j)=0d0
270 IF(abs(smw(1)).LT.abs(smz(1))+2d0*pmas(
pycomp(111),1))
THEN
271 CALL
pyerrm(18,
'(PYINOM:) '//
272 &
'forcing m(~chi+_1) > m(~chi0_1) + 2m(pi0)')
273 smw(1)=
sign(abs(smz(1))+2d0*pmas(
pycomp(111),1),smw(1))
275 pmas(
pycomp(ksusy1+24),1)=smw(1)
276 pmas(
pycomp(ksusy1+37),1)=smw(2)
281 ar(1,1) = xm2**2+2d0*xmw**2*cosb**2
282 ar(2,2) = xmu**2+2d0*xmw**2*sinb**2
283 ar(1,2) = sqrt(2d0)*xmw*(xm2*cos(rmss(31))*sinb+
284 &xmu*cos(rmss(33))*cosb)
285 ai(1,2) = sqrt(2d0)*xmw*(-xm2*sin(rmss(31))*sinb+
286 &xmu*sin(rmss(33))*cosb)
287 ar(2,1) = sqrt(2d0)*xmw*(xm2*cos(rmss(31))*sinb+
288 &xmu*cos(rmss(33))*cosb)
289 ai(2,1) = sqrt(2d0)*xmw*(xm2*sin(rmss(31))*sinb-
290 &xmu*sin(rmss(33))*cosb)
291 CALL
pyeicg(4,2,ar,ai,wr,wi,1,zr,zi,fv1,fv2,fv3,ierr)
293 WRITE(mstu(11),*)
' PROBLEM WITH PYEICG IN PYINOM '
297 IF(wr(2).LT.wr(1))
THEN
309 vmix(
i,
j)=zr(
j,
k)/sqrt(
s)
310 vmixi(
i,
j)=-zi(
j,
k)/sqrt(
s)
311 IF(abs(vmix(
i,
j)).LT.1d-6) vmix(
i,
j)=0d0
312 IF(abs(vmixi(
i,
j)).LT.1d-6) vmixi(
i,
j)=0d0