2534{
2548
2549 double betalambda3 = 3.0*la1Q*la3Q + 3.0*la2Q*la3Q + 2.0*la3Q*la3Q + la1Q*la4Q + la2Q*la4Q + la4Q*la4Q + la5Q*la5Q
2550 + 3.0*la3Q*Yb1Q*Yb1Q + 3.0*la3Q*Yb2Q*Yb2Q + la3Q*Ytau1Q*Ytau1Q + la3Q*Ytau2Q*Ytau2Q + 3.0*la3Q*YtQ*YtQ
2551 - 6.0*Yb1Q*Yb1Q*(Yb2Q*Yb2Q + YtQ*YtQ) - 2.0*Ytau1Q*Ytau1Q*Ytau2Q*Ytau2Q;
2552 double betalambda5 = la5Q*(la1Q + la2Q + 4.0*la3Q + 6.0*la4Q)
2553 + 3.0*la5Q*Yb1Q*Yb1Q + la5Q*Ytau1Q*Ytau1Q + la5Q*(3.0*Yb2Q*Yb2Q + Ytau2Q*Ytau2Q + 3.0*YtQ*YtQ)
2554 - 6.0*Yb1Q*Yb1Q*Yb2Q*Yb2Q - 2.0*Ytau1Q*Ytau1Q*Ytau2Q*Ytau2Q;
2555
2556 double uniNLO12a = -la3Q/(16.0*M_PI);
2557 double uniNLO12b = 3.0*betalambda3/(256.0*M_PI*M_PI*M_PI);
2558 gslpp::complex uniNLO12c = (gslpp::complex::i()*M_PI-1.0)*(la3Q*la3Q+la5Q*la5Q)/(256.0*M_PI*M_PI*M_PI);
2559 gslpp::complex uniNLO12d = -la3Q/(512.0*M_PI*M_PI*M_PI) * WFRc3;
2560
2561 double uniNLO13a = -la5Q/(16.0*M_PI);
2562 double uniNLO13b = 3.0*betalambda5/(256.0*M_PI*M_PI*M_PI);
2563 gslpp::complex uniNLO13c = (gslpp::complex::i()*M_PI-1.0)*la3Q*la5Q/(128.0*M_PI*M_PI*M_PI);
2564 gslpp::complex uniNLO13d = -la4Q/(512.0*M_PI*M_PI*M_PI) * WFRc2;
2565
2566 gslpp::complex uniA=uniNLO12a+uniNLO12b+uniNLO12c+uniNLO12d;
2567 gslpp::complex uniC=uniNLO13a+uniNLO13b+uniNLO13c+uniNLO13d;
2568
2569 double Rp=0.01;
2570 double a0LO = uniNLO12a+uniNLO13a;
2571 if(fabs(a0LO)>Rpeps)
2572 {
2573 Rp = ((uniA+uniC)-a0LO).abs()/fabs(a0LO);
2574 }
2575
2576 return Rp;
2577}
THDMcache * getMyTHDMCache() const
A method get the THDMCache.