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);
CCBd11(const StandardModel &SM_i)
CCBd11 constructor.
CCBd12(const StandardModel &SM_i)
CCBd12 constructor.
CCBd13(const StandardModel &SM_i)
CCBd13 constructor.
CCBd22(const StandardModel &SM_i)
CCBd22 constructor.
CCBd23(const StandardModel &SM_i)
CCBd23 constructor.
CCBd33(const StandardModel &SM_i)
CCBd33 constructor.
CCBe11(const StandardModel &SM_i)
CCBe11 constructor.
CCBe12(const StandardModel &SM_i)
CCBe12 constructor.
CCBe13(const StandardModel &SM_i)
CCBe13 constructor.
CCBe22(const StandardModel &SM_i)
CCBe22 constructor.
CCBe23(const StandardModel &SM_i)
CCBe23 constructor.
CCBe33(const StandardModel &SM_i)
CCBe33 constructor.
CCBu11(const StandardModel &SM_i)
CCBu11 constructor.
CCBu12(const StandardModel &SM_i)
CCBu12 constructor.
CCBu13(const StandardModel &SM_i)
CCBu13 constructor.
CCBu22(const StandardModel &SM_i)
CCBu22 constructor.
CCBu23(const StandardModel &SM_i)
CCBu23 constructor.
CCBu33(const StandardModel &SM_i)
CCBu33 constructor.
const double & getMass() const
A get method to access the particle mass.
const Particle & getQuarks(const QCD::quark q) const
A get method to access a quark as an object of the type Particle.
A base class for SUSY models.
gslpp::matrix< gslpp::complex > getMsLhat2() const
gslpp::complex getMuH() const
Gets the parameter in the superpotential.
gslpp::matrix< gslpp::complex > getTDhat() const
Gets the trilinear-coupling matrix for down-type squarks.
const double getMHptree() const
Gets the charged-Higgs mass at tree-level.
gslpp::matrix< gslpp::complex > getMsDhat2() const
gslpp::matrix< gslpp::complex > getTEhat() const
Gets the trilinear-coupling matrix for charged sleptons.
gslpp::matrix< gslpp::complex > getMsUhat2() const
const double getCosb() const
Gets .
gslpp::matrix< gslpp::complex > getMsQhat2() const
const double getSinb() const
Gets .
gslpp::matrix< gslpp::complex > getMsEhat2() const
gslpp::matrix< gslpp::complex > getTUhat() const
Gets the trilinear-coupling matrix for up-type squarks.
A model class for the Standard Model.
const Particle & getLeptons(const QCD::lepton p) const
A get method to retrieve the member object of a lepton.
const double getMz() const
A get method to access the mass of the boson .
A class for a model prediction of an observable.