4982{
4983 bool flag_Quadratic = false;
4984
4985
4988 double sbar2 = sbar * sbar;
4990 double mt2 = mt * mt;
4992 double mH2 = mH * mH;
4994 double v2 = v * v;
4996
4997 gslpp::complex gL = 1.5 * v2 * (
ewgc(
"CuHR",2,1) - gslpp::complex::i() *
ewgc(
"CuHI",2,1));
4998 gslpp::complex gR = 1.5 * v2 * (
ewgc(
"CuHR",1,2) + gslpp::complex::i() *
ewgc(
"CuHI",1,2));
4999 double absgL2 = (gL * gL.conjugate()).abs();
5000 double absgR2 = (gR * gR.conjugate()).abs();
5001
5002
5003 if(flag_Quadratic){
5004
5005 return 0.;
5006
5007 }
5008 else{
5009
5010 double total = alpha / (32. * sbar2) * mt * (absgL2 + absgR2) * (1. - mH2 / mt2) * (1. - mH2 / mt2) / Gammatop;
5011
5012 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
5013
5014 return total;
5015
5016 }
5017
5018}
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