3550{
3563
3564 double betalambda3 = 3.0*la1Q*la3Q + 3.0*la2Q*la3Q + 2.0*la3Q*la3Q + la1Q*la4Q + la2Q*la4Q + la4Q*la4Q + la5Q*la5Q
3565 + 3.0*la3Q*Yb1Q*Yb1Q + 3.0*la3Q*Yb2Q*Yb2Q + la3Q*Ytau1Q*Ytau1Q + la3Q*Ytau2Q*Ytau2Q + 3.0*la3Q*YtQ*YtQ
3566 - 6.0*Yb1Q*Yb1Q*(Yb2Q*Yb2Q + YtQ*YtQ) - 2.0*Ytau1Q*Ytau1Q*Ytau2Q*Ytau2Q;
3567 double betalambda5 = la5Q*(la1Q + la2Q + 4.0*la3Q + 6.0*la4Q)
3568 + 3.0*la5Q*Yb1Q*Yb1Q + la5Q*Ytau1Q*Ytau1Q + la5Q*(3.0*Yb2Q*Yb2Q + Ytau2Q*Ytau2Q + 3.0*YtQ*YtQ)
3569 - 6.0*Yb1Q*Yb1Q*Yb2Q*Yb2Q - 2.0*Ytau1Q*Ytau1Q*Ytau2Q*Ytau2Q;
3570
3571 double uniNLO10a = -la3Q/(16.0*M_PI);
3572 double uniNLO10b = 3.0*betalambda3/(256.0*M_PI*M_PI*M_PI);
3573 gslpp::complex uniNLO10c = (gslpp::complex::i()*M_PI-1.0)*(la3Q*la3Q+la5Q*la5Q)/(256.0*M_PI*M_PI*M_PI);
3574 gslpp::complex uniNLO10d = -(la3Q+la4Q-la5Q)/(512.0*M_PI*M_PI*M_PI) * WFRc2;
3575
3576 double uniNLO11a = -la5Q/(16.0*M_PI);
3577 double uniNLO11b = 3.0*betalambda5/(256.0*M_PI*M_PI*M_PI);
3578 gslpp::complex uniNLO11c = (gslpp::complex::i()*M_PI-1.0)*la3Q*la5Q/(128.0*M_PI*M_PI*M_PI);
3579 gslpp::complex uniNLO11d = -(la4Q-2.0*la5Q)/(512.0*M_PI*M_PI*M_PI) * WFRc2;
3580
3581 gslpp::complex uniA=uniNLO10a+uniNLO10b+uniNLO10c+uniNLO10d;
3582 gslpp::complex uniC=uniNLO11a+uniNLO11b+uniNLO11c+uniNLO11d;
3583
3584 double R=0.01;
3585 double a0LO = uniNLO10a+uniNLO11a;
3586 gslpp::complex a0NLO = (uniA+uniC);
3587 if(a0NLO.abs()>Rpeps)
3588 {
3589 R = (a0NLO-a0LO).abs()/(a0NLO.abs());
3590 }
3591
3592 return R;
3593}
THDMcache * getMyTHDMCache() const
A method get the THDMCache.