4841 bool flag_Quadratic =
false;
4847 double sbar2 = sbar * sbar;
4848 double cbar2 = cbar * cbar;
4851 double mZ2 = mZ * mZ;
4852 double mt2 = mt * mt;
4853 double mt3 = mt2 * mt;
4858 gslpp::complex XL = 0.5 * v2 * (
ewgc(
"CHq3R",0,2) + gslpp::complex::i() *
ewgc(
"CHq3I",0,2) +
ewgc(
"CHq3R",2,0) - gslpp::complex::i() *
ewgc(
"CHq3I",2,0) -
ewgc(
"CHq1R",0,2) - gslpp::complex::i() *
ewgc(
"CHq1I",0,2) -
ewgc(
"CHq1R",2,0) + gslpp::complex::i() *
ewgc(
"CHq1I",2,0));
4859 double absXL2 = (XL * XL.conjugate()).abs();
4860 gslpp::complex XR = - 0.5 * v2 * (
ewgc(
"CHuR",0,2) + gslpp::complex::i() *
ewgc(
"CHuI",0,2) +
ewgc(
"CHuR",2,0) - gslpp::complex::i() *
ewgc(
"CHuI",2,0));
4861 double absXR2 = (XR * XR.conjugate()).abs();
4862 gslpp::complex kL = sqrt(2.) * v2 * (cbar * (
ewgc(
"CuWR",2,0) - gslpp::complex::i() *
ewgc(
"CuWI",2,0)) - sbar * (
ewgc(
"CuBR",2,0) - gslpp::complex::i() *
ewgc(
"CuBI",2,0)));
4863 double abskL2 = (kL * kL.conjugate()).abs();
4864 gslpp::complex kR = sqrt(2.) * v2 * (cbar * (
ewgc(
"CuWR",0,2) + gslpp::complex::i() *
ewgc(
"CuWI",0,2)) - sbar * (
ewgc(
"CuBR",0,2) + gslpp::complex::i() *
ewgc(
"CuBI",0,2)));
4865 double abskR2 = (kR * kR.conjugate()).abs();
4875 double Gammag = alpha / (32. * sbar2 * cbar2) * (absXL2 + absXR2) * mt3 / mZ2 * (1. - mZ2 / mt2) * (1. - mZ2 / mt2) * (1. + 2. * mZ2 / mt2);
4876 double Gammas = alpha / (16. * sbar2 * cbar2) * (abskL2 + abskR2) * mt * (1. - mZ2 / mt2) * (1. - mZ2 / mt2) * (2. + mZ2 / mt2);
4877 double Gammags = alpha / (32. * sbar2 * cbar2) * 12. * ((XL * kR.conjugate()).real() + (XR * kL.conjugate()).real()) * mt * (1. - mZ2 / mt2) * (1. - mZ2 / mt2);
4878 double total = (Gammag + Gammas + Gammags) / Gammatop;
4880 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
double ewgc(const std::string name) const
const double getOptionalParameter(std::string name) const
A method to get parameters that are specific to only one set of observables.
const double getMtpole() const
A get method to access the pole mass of the top quark.
const double getMz() const
A get method to access the mass of the boson .
const double v() const
The Higgs vacuum expectation value.
const double getAle() const
A get method to retrieve the fine-structure constant .
const NPSMEFTd6General & GetmyNPSMEFTd6General() const