27 double cos2b = cosb*cosb-sinb*sinb;
33 double lambdau11 = mUP / v2 *
sqrt(2.);
34 double mHusq = mA*mA * cosb*cosb + 0.5 * MZ*MZ * cos2b - MuH.
abs2();
36 return (lambdau11*lambdau11*(msQhat2(0,0).abs()+msUhat2(0,0).abs()+mHusq)-TUhat(0,0).abs2()) / fabs(msQhat2(0,0).abs()+msUhat2(0,0).abs()+mHusq);
51 double cos2b = cosb*cosb-sinb*sinb;
57 double lambdau22 = mCHARM / v2 *
sqrt(2.);
58 double mHusq = mA*mA * cosb*cosb + 0.5 * MZ*MZ * cos2b - MuH.
abs2();
60 return lambdau22*lambdau22-TUhat(1,1).abs2()/(msQhat2(1,1).abs()+msUhat2(1,1).abs()+mHusq);
75 double cos2b = cosb*cosb-sinb*sinb;
81 double lambdau33 = mTOP / v2 *
sqrt(2.);
82 double mHusq = mA*mA * cosb*cosb + 0.5 * MZ*MZ * cos2b - MuH.
abs2();
84 return lambdau33*lambdau33-TUhat(2,2).abs2()/(msQhat2(2,2).abs()+msUhat2(2,2).abs()+mHusq);
99 double cos2b = cosb*cosb-sinb*sinb;
105 double lambdau22 = mCHARM / v2 *
sqrt(2.);
106 double mHusq = mA*mA * cosb*cosb + 0.5 * MZ*MZ * cos2b - MuH.
abs2();
108 return lambdau22*lambdau22-TUhat(0,1).abs2()/(msQhat2(0,0).abs()+msUhat2(1,1).abs()+mHusq);
123 double cos2b = cosb*cosb-sinb*sinb;
129 double lambdau33 = mTOP / v2 *
sqrt(2.);
130 double mHusq = mA*mA * cosb*cosb + 0.5 * MZ*MZ * cos2b - MuH.
abs2();
132 return lambdau33*lambdau33-TUhat(0,2).abs2()/(msQhat2(0,0).abs()+msUhat2(2,2).abs()+mHusq);
147 double cos2b = cosb*cosb-sinb*sinb;
153 double lambdau33 = mTOP / v2 *
sqrt(2.);
154 double mHusq = mA*mA * cosb*cosb + 0.5 * MZ*MZ * cos2b - MuH.
abs2();
156 return lambdau33*lambdau33-TUhat(1,2).abs2()/(msQhat2(1,1).abs()+msUhat2(2,2).abs()+mHusq);
171 double cos2b = cosb*cosb-sinb*sinb;
177 double lambdad11 = mDOWN / v1 *
sqrt(2.);
178 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
180 return lambdad11*lambdad11-TDhat(0,0).abs2()/(msQhat2(0,0).abs()+msDhat2(0,0).abs()+mHdsq);
195 double cos2b = cosb*cosb-sinb*sinb;
201 double lambdad22 = mSTRANGE / v1 *
sqrt(2.);
202 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
204 return lambdad22*lambdad22-TDhat(1,1).abs2()/(msQhat2(1,1).abs()+msDhat2(1,1).abs()+mHdsq);
219 double cos2b = cosb*cosb-sinb*sinb;
225 double lambdad33 = mBOTTOM / v1 *
sqrt(2.);
226 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
228 return lambdad33*lambdad33-TDhat(2,2).abs2()/(msQhat2(2,2).abs()+msDhat2(2,2).abs()+mHdsq);
243 double cos2b = cosb*cosb-sinb*sinb;
249 double lambdad22 = mSTRANGE / v1 *
sqrt(2.);
250 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
252 return lambdad22*lambdad22-TDhat(0,1).abs2()/(msQhat2(0,0).abs()+msDhat2(1,1).abs()+mHdsq);
267 double cos2b = cosb*cosb-sinb*sinb;
273 double lambdad33 = mBOTTOM / v1 *
sqrt(2.);
274 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
276 return lambdad33*lambdad33-TDhat(0,2).abs2()/(msQhat2(0,0).abs()+msDhat2(2,2).abs()+mHdsq);
291 double cos2b = cosb*cosb-sinb*sinb;
297 double lambdad33 = mBOTTOM / v1 *
sqrt(2.);
298 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
300 return lambdad33*lambdad33-TDhat(1,2).abs2()/(msQhat2(1,1).abs()+msDhat2(2,2).abs()+mHdsq);
315 double cos2b = cosb*cosb-sinb*sinb;
321 double lambdal11 = mE / v1 *
sqrt(2.);
322 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
324 return lambdal11*lambdal11-TEhat(0,0).abs2()/(msLhat2(0,0).abs()+msEhat2(0,0).abs()+mHdsq);
339 double cos2b = cosb*cosb-sinb*sinb;
345 double lambdal22 = mMU / v1 *
sqrt(2.);
346 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
348 return lambdal22*lambdal22-TEhat(1,1).abs2()/(msLhat2(1,1).abs()+msEhat2(1,1).abs()+mHdsq);
363 double cos2b = cosb*cosb-sinb*sinb;
369 double lambdal33 = mTAU / v1 *
sqrt(2.);
370 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
372 return lambdal33*lambdal33-TEhat(2,2).abs2()/(msLhat2(2,2).abs()+msEhat2(2,2).abs()+mHdsq);
387 double cos2b = cosb*cosb-sinb*sinb;
393 double lambdal22 = mMU / v1 *
sqrt(2.);
394 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
396 return lambdal22*lambdal22-TEhat(0,1).abs2()/(msLhat2(0,0).abs()+msEhat2(1,1).abs()+mHdsq);
411 double cos2b = cosb*cosb-sinb*sinb;
417 double lambdal33 = mTAU / v1 *
sqrt(2.);
418 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
420 return lambdal33*lambdal33-TEhat(0,2).abs2()/(msLhat2(0,0).abs()+msEhat2(2,2).abs()+mHdsq);
435 double cos2b = cosb*cosb-sinb*sinb;
441 double lambdal33 = mTAU / v1 *
sqrt(2.);
442 double mHdsq = mA*mA * sinb*sinb - 0.5 * MZ*MZ * cos2b - MuH.
abs2();
444 return lambdal33*lambdal33-TEhat(1,2).abs2()/(msLhat2(1,1).abs()+msEhat2(2,2).abs()+mHdsq);