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; }