double NPbase::deltaGV_f(const Particle f) const
{
if (f.is("TOP")) return 0.;
/* SM values */
double alpha = trueSM.alphaMz();
double sW2SM = trueSM.sW2();
double cW2SM = trueSM.cW2();
double gVSM = trueSM.gV_f(f).real();
double gASM = trueSM.gA_f(f).real();
return ( gVSM * (alpha * obliqueT() - DeltaGF()) / 2.0
+ (gVSM - gASM) / 4.0 / sW2SM / (cW2SM - sW2SM)
*(alpha * (obliqueS() - 4.0 * cW2SM * sW2SM * obliqueT())
+ 4.0 * cW2SM * sW2SM * DeltaGF()));
}
double NPbase::deltaGA_f(const Particle f) const
{
if (f.is("TOP")) return 0.;
/* SM values */
double alpha = trueSM.alphaMz();
double gASM = trueSM.gA_f(f).real();
return ( gASM * (alpha * obliqueT() - DeltaGF()) / 2.0);
}
double NPSTUVWXY::GammaW() const
{
double Gamma_W = trueSM.GammaW();
double Wbar = (obliqueV() - obliqueW()) / alphaMz();
double alpha = StandardModel::alphaMz();
double c2 = trueSM.cW2();
double s2 = trueSM.sW2();
Gamma_W *= 1.0 - 3.0 * alpha / 4.0 / (c2 - s2)
*(obliqueS() - 2.0 * c2 * obliqueT()
- (c2 - s2) * obliqueU() / 2.0 / s2 - 2.0 * (c2 - s2) * Wbar)
- (1.0 + c2) / 2.0 / (c2 - s2) * DeltaGF();
return Gamma_W;
}
double NPEffectiveBS::GammaW() const
{
double Gamma_W = trueSM.GammaW();
double alpha = alphaMz();
double c2 = trueSM.cW2();
double s2 = trueSM.sW2();
double ratio = v() * v() / Lambda_NP / Lambda_NP;
Gamma_W *= 1.0 - 3.0 * alpha / 4.0 / (c2 - s2)
*(obliqueS() - 2.0 * c2 * obliqueT() - (c2 - s2) * obliqueU() / 2.0 / s2)
- (1.0 + c2) / 2.0 / (c2 - s2) * DeltaGF()
//+ (cHL1p + cHL2p + cHL3p + cHQ1p + cHQ2p) * ratio; // incorrect
+ 2.0 / 9.0 * (cHL1p + cHL2p + cHL3p + 3.0 * cHQ1p + 3.0 * cHQ2p) * ratio;
return Gamma_W;
}
double NPbase::GammaW() const
{
double Gamma_W = trueSM.GammaW();
double alpha = trueSM.alphaMz();
double c2 = trueSM.cW2();
double s2 = trueSM.sW2();
Gamma_W *= 1.0 - 3.0 * alpha / 4.0 / (c2 - s2)
*(obliqueS() - 2.0 * c2 * obliqueT() - (c2 - s2) * obliqueU() / 2.0 / s2)
- (1.0 + c2) / 2.0 / (c2 - s2) * DeltaGF();
//std::cout << "Gw: c_S=" << - 3.0*alpha/4.0/(c2-s2) << std::endl;
//std::cout << "Gw: c_T=" << - 3.0*alpha/4.0/(c2-s2)*(- 2.0*c2) << std::endl;
//std::cout << "Gw: c_U=" << - 3.0*alpha/4.0/(c2-s2)*(- (c2-s2)/2.0/s2) << std::endl;
return Gamma_W;
}
double NPbase::Mw() const
{
double myMw = trueSM.Mw();
double alpha = trueSM.alphaMz();
double c2 = trueSM.cW2();
double s2 = trueSM.sW2();
myMw *= 1.0 - alpha / 4.0 / (c2 - s2)
*(obliqueS() - 2.0 * c2 * obliqueT() - (c2 - s2) * obliqueU() / 2.0 / s2)
- s2 / 2.0 / (c2 - s2) * DeltaGF();
//std::cout << "Mw: c_S=" << - alpha/4.0/(c2-s2) << std::endl;
//std::cout << "Mw: c_T=" << - alpha/4.0/(c2-s2)*(- 2.0*c2) << std::endl;
//std::cout << "Mw: c_U=" << - alpha/4.0/(c2-s2)*(- (c2-s2)/2.0/s2) << std::endl;
return myMw;
}
double NPEffectiveBS::Mw_tree() const
{
double GF0 = GF * (1.0 - DeltaGF());
double tmp = 4.0 * M_PI * ale / sqrt(2.0) / GF0 / Mz / Mz;
return ( Mz / sqrt(2.0) * sqrt(1.0 + sqrt(1.0 - tmp)));
}
