4939{
4940 bool flag_Quadratic = false;
4941
4942
4950
4951 gslpp::complex lambdaL = sqrt(2.) / ee * v * mt * (sbar * (
ewgc(
"CuWR",2,0) - gslpp::complex::i() *
ewgc(
"CuWI",2,0)) + cbar * (
ewgc(
"CuBR",2,0) - gslpp::complex::i() *
ewgc(
"CuBI",2,0)));
4952 gslpp::complex lambdaR = sqrt(2.) / ee * v * mt * (sbar * (
ewgc(
"CuWR",0,2) + gslpp::complex::i() *
ewgc(
"CuWI",0,2)) + cbar * (
ewgc(
"CuBR",0,2) + gslpp::complex::i() *
ewgc(
"CuBI",0,2)));
4953 double abslambdaL2 = (lambdaL * lambdaL.conjugate()).abs();
4954 double abslambdaR2 = (lambdaR * lambdaR.conjugate()).abs();
4955
4956
4957 if(flag_Quadratic){
4958
4959 return 0.;
4960
4961 }
4962 else{
4963
4964 double total = alpha / 2 * (abslambdaL2 + abslambdaR2) * mt / Gammatop;
4965
4966 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
4967
4968 return total;
4969
4970 }
4971
4972}
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 v() const
The Higgs vacuum expectation value.
const double getAle() const
A get method to retrieve the fine-structure constant .
const NPSMEFTd6General & GetmyNPSMEFTd6General() const