4896{
4897 bool flag_Quadratic = false;
4898
4899
4907
4908 gslpp::complex lambdaL = sqrt(2.) / ee * v * mt * (sbar * (
ewgc(
"CuWR",2,1) - gslpp::complex::i() *
ewgc(
"CuWI",2,1)) + cbar * (
ewgc(
"CuBR",2,1) - gslpp::complex::i() *
ewgc(
"CuBI",2,1)));
4909 gslpp::complex lambdaR = sqrt(2.) / ee * v * mt * (sbar * (
ewgc(
"CuWR",1,2) + gslpp::complex::i() *
ewgc(
"CuWI",1,2)) + cbar * (
ewgc(
"CuBR",1,2) + gslpp::complex::i() *
ewgc(
"CuBI",1,2)));
4910 double abslambdaL2 = (lambdaL * lambdaL.conjugate()).abs();
4911 double abslambdaR2 = (lambdaR * lambdaR.conjugate()).abs();
4912
4913
4914 if(flag_Quadratic){
4915
4916 return 0.;
4917
4918 }
4919 else{
4920
4921 double total = alpha / 2 * (abslambdaL2 + abslambdaR2) * mt / Gammatop;
4922
4923 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
4924
4925 return total;
4926
4927 }
4928
4929}
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