7 SUBROUTINE pyrghm(MCHI,MA,TANB,MQ,MUR,MD,MTOP,AU,AD,MU,
8 * mhp,hmp,mch,sa,ca,sab,cab,tanba,mglu,deltamt,deltamb)
9 IMPLICIT DOUBLE PRECISION(a-
h,l,m,o-
z)
10 dimension vh(2,2),m2(2,2),m2p(2,2)
13 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
26 alpha3 = alpha3z/(1d0 +(11d0 - 10d0/3d0)/4d0/
pi*alpha3z*
29 rmtop = mtop/(1d0+4d0*alpha3/3d0/
pi)
30 tq =
log((mq**2+mtop**2)/mtop**2)
31 tu =
log((mur**2 + mtop**2)/mtop**2)
32 td =
log((md**2 + mtop**2)/mtop**2)
38 tglu =
log(mglu**2/mtop**2)
39 sinb = tanb/dsqrt(1d0 + tanb**2)
42 *tanba = tanb*(1d0-3d0/32d0/
pi**2*
43 *(rmtop**2/
v**2/sinb**2-mb**2/
v**2/cosb**2)*
45 IF(ma.LT.mtop.OR.ma.EQ.mtop) tanbt = tanba
46 sinb = tanbt/sqrt(1d0 + tanbt**2)
47 cosb = 1d0/dsqrt(1d0 + tanbt**2)
48 g1 = sqrt(alpha1*4d0*
pi)
49 g2 = sqrt(alpha2*4d0*
pi)
50 g3 = sqrt(alpha3*4d0*
pi)
53 CALL
pygfxx(ma,tanba,mq,mur,md,mtop,au,ad,mu,mglu,vh,stop1,stop2,
54 *sbot1,sbot2,deltamt,deltamb)
55 IF(mq.GT.mur) tp = tq - tu
56 IF(mq.LT.mur.OR.mq.EQ.mur) tp = tu - tq
57 IF(mq.GT.mur) tdp = tu
58 IF(mq.LT.mur.OR.mq.EQ.mur) tdp = tq
59 IF(mq.GT.md) tpd = tq - td
60 IF(mq.LT.md.OR.mq.EQ.md) tpd = td - tq
61 IF(mq.GT.md) tdpd = td
62 IF(mq.LT.md.OR.mq.EQ.md) tdpd = tq
64 IF(mq.GT.md) dlambda1 = 6d0/96d0/
pi**2*g1**2*hd**2*tpd
65 IF(mq.LT.md.OR.mq.EQ.md) dlambda1 = 3d0/32d0/
pi**2*
66 * hd**2*(g1**2/3d0+g2**2)*tpd
68 IF(mq.GT.mur) dlambda2 =12d0/96d0/
pi**2*g1**2*hu**2*tp
69 IF(mq.LT.mur.OR.mq.EQ.mur) dlambda2 = 3d0/32d0/
pi**2*
70 * hu**2*(-g1**2/3d0+g2**2)*tp
83 IF(mglu.LT.mur.OR.mglu.LT.mq)
THEN
84 IF(mq.GT.mur.AND.mglu.GT.mur)
THEN
85 dlambdap2 = -4d0/(16d0*
pi**2)**2*hu**4*(tq**2-tglu**2)
88 IF(mq.GT.mur.AND.mglu.LT.mur)
THEN
89 dlambdap2 = -4d0/(16d0*
pi**2)**2*hu**4*(tq**2-tu**2)
92 IF(mq.GT.mur.AND.mglu.EQ.mur)
THEN
93 dlambdap2 = -4d0/(16d0*
pi**2)**2*hu**4*(tq**2-tu**2)
96 IF(mur.GT.mq.AND.mglu.GT.mq)
THEN
97 dlambdap2 = -4d0/(16d0*
pi**2)**2*hu**4*(tu**2-tglu**2)
100 IF(mur.GT.mq.AND.mglu.LT.mq)
THEN
101 dlambdap2 = -4d0/(16d0*
pi**2)**2*hu**4*(tu**2-tq**2)
104 IF(mur.GT.mq.AND.mglu.EQ.mq)
THEN
105 dlambdap2 = -4d0/(16d0*
pi**2)**2*hu**4*(tu**2-tq**2)
110 IF(mq.GT.md) dlambda3 = -1d0/32d0/
pi**2*g1**2*hd**2*tpd
111 IF(mq.LT.md.OR.mq.EQ.md) dlambda3 = 3d0/64d0/
pi**2*hd**2*
112 *(g2**2-g1**2/3d0)*tpd
113 IF(mq.GT.mur) dlambda3 = dlambda3 -
114 *1d0/16d0/
pi**2*g1**2*hu**2*tp
115 IF(mq.LT.mur.OR.mq.EQ.mur) dlambda3 = dlambda3 +
116 * 3d0/64d0/
pi**2*hu**2*(g2**2+g1**2/3d0)*tp
117 IF(mq.LT.mur) dlambda4 = -3d0/32d0/
pi**2*g2**2*hu**2*tp
118 IF(mq.LT.md) dlambda4 = dlambda4 - 3d0/32d0/
pi**2*g2**2*
120 lambda1 = ((g1**2 + g2**2)/4d0)*
121 * (1d0-3d0*hd**2*(tpd + tdpd)/8d0/
pi**2)
122 *+(3d0*hd**4d0/16d0/
pi**2) *tpd*(1d0
123 *+ (3d0*hd**2/2d0 + hu**2/2d0
124 *- 8d0*g3**2) * (tpd + 2d0*tdpd)/16d0/
pi**2)
125 *+(3d0*hd**4d0/8d0/
pi**2) *tdpd*(1d0 + (3d0*hd**2/2d0 + hu**2/2d0
126 *- 8d0*g3**2) * tdpd/16d0/
pi**2) + dlambda1
127 lambda2 = ((g1**2 + g2**2)/4d0)*(1d0-3d0*hu**2*
128 *(tp + tdp)/8d0/
pi**2)
129 *+(3d0*hu**4d0/16d0/
pi**2) *tp*(1d0
130 *+ (3d0*hu**2/2d0 + hd**2/2d0
131 *- 8d0*g3**2) * (tp + 2d0*tdp)/16d0/
pi**2)
132 *+(3d0*hu**4d0/8d0/
pi**2) *tdp*(1d0 + (3d0*hu**2/2d0 + hd**2/2d0
133 *- 8d0*g3**2) * tdp/16d0/
pi**2) + dlambda2 + dlambdap2
134 lambda3 = ((g2**2 - g1**2)/4d0)*(1d0-3d0*
135 *(hu**2)*(tp + tdp)/16d0/
pi**2 -3d0*
136 *(hd**2)*(tpd + tdpd)/16d0/
pi**2) +dlambda3
137 lambda4 = (- g2**2/2d0)*(1d0
138 *-3d0*(hu**2)*(tp + tdp)/16d0/
pi**2
139 *-3d0*(hd**2)*(tpd + tdpd)/16d0/
pi**2) +dlambda4
145 m2(1,1) = 2d0*
v**2*(lambda1*cosb**2+2d0*lambda6*
146 *cosb*sinb + lambda5*sinb**2) + ma**2*sinb**2
148 m2(2,2) = 2d0*
v**2*(lambda5*cosb**2+2d0*lambda7*
149 *cosb*sinb + lambda2*sinb**2) + ma**2*cosb**2
150 m2(1,2) = 2d0*
v**2*(lambda6*cosb**2+(lambda3+lambda4)*
151 *cosb*sinb + lambda7*sinb**2) - ma**2*sinb*cosb
158 mssusy=dsqrt(.5d0*(mq**2+mur**2)+mtop**2)
160 IF(mchi.GT.mssusy) goto 100
161 IF(mchi.LT.mtop) mchi=mtop
163 tchar=
log(mssusy**2/mchi**2)
165 deltal12=(9d0/64d0/
pi**2*g2**4+5d0/192d0/
pi**2*g1**4)*tchar
166 deltal3p4=(3d0/64d0/
pi**2*g2**4+7d0/192d0/
pi**2*g1**4
167 *+4d0/32d0/
pi**2*g1**2*g2**2)*tchar
169 deltam112=2d0*deltal12*
v**2*cosb**2
170 deltam222=2d0*deltal12*
v**2*sinb**2
171 deltam122=2d0*deltal3p4*
v**2*sinb*cosb
173 m2(1,1)=m2(1,1)+deltam112
174 m2(2,2)=m2(2,2)+deltam222
175 m2(1,2)=m2(1,2)+deltam122
176 m2(2,1)=m2(2,1)+deltam122
186 m2p(
i,
j) = m2(
i,
j) + vh(
i,
j)
189 trm2p = m2p(1,1) + m2p(2,2)
190 detm2p = m2p(1,1)*m2p(2,2) - m2p(1,2)*m2p(2,1)
191 mh2p = (trm2p - dsqrt(trm2p**2 - 4d0* detm2p))/2d0
192 hm2p = (trm2p + dsqrt(trm2p**2 - 4d0* detm2p))/2d0
194 mch2=ma**2+(lambda5-lambda4)*
v**2
196 IF(mh2p.LT.0.) goto 130
198 sin2alpha = 2d0*m2p(1,2)/sqrt(trm2p**2-4d0*detm2p)
199 cos2alpha = (m2p(1,1)-m2p(2,2))/sqrt(trm2p**2-4d0*detm2p)
200 IF(cos2alpha.GE.0.)
THEN
201 alpha = asin(sin2alpha)/2d0
203 alpha = -
pi/2d0-asin(sin2alpha)/2d0
215 sab = sa*(1d0-deltamb/(1d0+deltamb)*(1d0+ca/sa/tanb))
216 cab = ca*(1d0-deltamb/(1d0+deltamb)*(1d0-sa/ca/tanb))