4785 bool flag_Quadratic =
false;
4791 double sbar2 = sbar * sbar;
4792 double cbar2 = cbar * cbar;
4795 double mZ2 = mZ * mZ;
4796 double mt2 = mt * mt;
4797 double mt3 = mt2 * mt;
4802 gslpp::complex XL = 0.5 * v2 * (
ewgc(
"CHq3R",1,2) + gslpp::complex::i() *
ewgc(
"CHq3I",1,2) +
ewgc(
"CHq3R",2,1) - gslpp::complex::i() *
ewgc(
"CHq3I",2,1) -
ewgc(
"CHq1R",1,2) - gslpp::complex::i() *
ewgc(
"CHq1I",1,2) -
ewgc(
"CHq1R",2,1) + gslpp::complex::i() *
ewgc(
"CHq1I",2,1));
4803 double absXL2 = (XL * XL.conjugate()).abs();
4804 gslpp::complex XR = - 0.5 * v2 * (
ewgc(
"CHuR",1,2) + gslpp::complex::i() *
ewgc(
"CHuI",1,2) +
ewgc(
"CHuR",2,1) - gslpp::complex::i() *
ewgc(
"CHuI",2,1));
4805 double absXR2 = (XR * XR.conjugate()).abs();
4806 gslpp::complex kL = sqrt(2.) * v2 * (cbar * (
ewgc(
"CuWR",2,1) - gslpp::complex::i() *
ewgc(
"CuWI",2,1)) - sbar * (
ewgc(
"CuBR",2,1) - gslpp::complex::i() *
ewgc(
"CuBI",2,1)));
4807 double abskL2 = (kL * kL.conjugate()).abs();
4808 gslpp::complex kR = sqrt(2.) * v2 * (cbar * (
ewgc(
"CuWR",1,2) + gslpp::complex::i() *
ewgc(
"CuWI",1,2)) - sbar * (
ewgc(
"CuBR",1,2) + gslpp::complex::i() *
ewgc(
"CuBI",1,2)));
4809 double abskR2 = (kR * kR.conjugate()).abs();
4819 double Gammag = alpha / (32. * sbar2 * cbar2) * (absXL2 + absXR2) * mt3 / mZ2 * (1. - mZ2 / mt2) * (1. - mZ2 / mt2) * (1. + 2. * mZ2 / mt2);
4820 double Gammas = alpha / (16. * sbar2 * cbar2) * (abskL2 + abskR2) * mt * (1. - mZ2 / mt2) * (1. - mZ2 / mt2) * (2. + mZ2 / mt2);
4821 double Gammags = alpha / (32. * sbar2 * cbar2) * 12. * ((XL * kR.conjugate()).real() + (XR * kL.conjugate()).real()) * mt * (1. - mZ2 / mt2) * (1. - mZ2 / mt2);
4822 double total = (Gammag + Gammas + Gammags) / Gammatop;
4824 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