5028{
5029 bool flag_Quadratic = false;
5030
5031
5034 double sbar2 = sbar * sbar;
5036 double mt2 = mt * mt;
5038 double mH2 = mH * mH;
5040 double v2 = v * v;
5042
5043 gslpp::complex gL = 1.5 * v2 * (
ewgc(
"CuHR",2,0) - gslpp::complex::i() *
ewgc(
"CuHI",2,0));
5044 gslpp::complex gR = 1.5 * v2 * (
ewgc(
"CuHR",0,2) + gslpp::complex::i() *
ewgc(
"CuHI",0,2));
5045 double absgL2 = (gL * gL.conjugate()).abs();
5046 double absgR2 = (gR * gR.conjugate()).abs();
5047
5048
5049 if(flag_Quadratic){
5050
5051 return 0.;
5052
5053 }
5054 else{
5055
5056 double total = alpha / (32. * sbar2) * mt * (absgL2 + absgR2) * (1. - mH2 / mt2) * (1. - mH2 / mt2) / Gammatop;
5057
5058 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
5059
5060 return total;
5061
5062 }
5063
5064}
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.
virtual const double getMHl() const
A get method to retrieve the Higgs mass .
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