5081{
5096
5097 double betalambda1 = 6.0*la1Q*la1Q + 2.0*la3Q*la3Q + 2.0*la3Q*la4Q + la4Q*la4Q + la5Q*la5Q
5098 + 6.0*la1Q*Yb1Q*Yb1Q + 2.0*la1Q*Ytau1Q*Ytau1Q
5099 - 6.0*Yb1Q*Yb1Q*Yb1Q*Yb1Q - 2.0*Ytau1Q*Ytau1Q*Ytau1Q*Ytau1Q;
5100 double betalambda2 = 6.0*la2Q*la2Q + 2.0*la3Q*la3Q + 2.0*la3Q*la4Q + la4Q*la4Q + la5Q*la5Q
5101 + 6.0*la2Q*Yb2Q*Yb2Q + 2.0*la2Q*Ytau2Q*Ytau2Q + 6.0*la2Q*YtQ*YtQ
5102 - 6.0*Yb2Q*Yb2Q*Yb2Q*Yb2Q - 2.0*Ytau2Q*Ytau2Q*Ytau2Q*Ytau2Q - 6.0*YtQ*YtQ*YtQ*YtQ;
5103 double betalambda5 = la5Q*(la1Q + la2Q + 4.0*la3Q + 6.0*la4Q)
5104 + 3.0*la5Q*Yb1Q*Yb1Q + la5Q*Ytau1Q*Ytau1Q + la5Q*(3.0*Yb2Q*Yb2Q + Ytau2Q*Ytau2Q + 3.0*YtQ*YtQ)
5105 - 6.0*Yb1Q*Yb1Q*Yb2Q*Yb2Q - 2.0*Ytau1Q*Ytau1Q*Ytau2Q*Ytau2Q;
5106
5107 double uniNLO15a = -la1Q/(16.0*M_PI);
5108 double uniNLO15b = 3.0*betalambda1/(256.0*M_PI*M_PI*M_PI);
5109 gslpp::complex uniNLO15c = (gslpp::complex::i()*M_PI-1.0)*(la1Q*la1Q+la5Q*la5Q)/(256.0*M_PI*M_PI*M_PI);
5110 gslpp::complex uniNLO15d = -la1Q/(512.0*M_PI*M_PI*M_PI) * (WFRc1-2.0*WFRc2+WFRc3+2.0*WFRc4);
5111
5112 double uniNLO16a = -la2Q/(16.0*M_PI);
5113 double uniNLO16b = 3.0*betalambda2/(256.0*M_PI*M_PI*M_PI);
5114 gslpp::complex uniNLO16c = (gslpp::complex::i()*M_PI-1.0)*(la2Q*la2Q+la5Q*la5Q)/(256.0*M_PI*M_PI*M_PI);
5115 gslpp::complex uniNLO16d = -la2Q/(512.0*M_PI*M_PI*M_PI) * (-WFRc1+2.0*WFRc2-WFRc3+2.0*WFRc4);
5116
5117 double uniNLO17a = -la5Q/(16.0*M_PI);
5118 double uniNLO17b = 3.0*betalambda5/(256.0*M_PI*M_PI*M_PI);
5119 gslpp::complex uniNLO17c = (gslpp::complex::i()*M_PI-1.0)*(la1Q+la2Q)*la5Q/(256.0*M_PI*M_PI*M_PI);
5120 gslpp::complex uniNLO17d = -2.0*la5Q/(512.0*M_PI*M_PI*M_PI) * WFRc4;
5121
5122 gslpp::complex uniA=uniNLO15a+uniNLO15b+uniNLO15c+uniNLO15d;
5123 gslpp::complex uniB=uniNLO16a+uniNLO16b+uniNLO16c+uniNLO16d;
5124 gslpp::complex uniC=uniNLO17a+uniNLO17b+uniNLO17c+uniNLO17d;
5125
5126 return (0.5*(uniA+uniB+sqrt((uniA-uniB)*(uniA-uniB)+4.0*uniC*uniC))).imag();
5127}
THDMcache * getMyTHDMCache() const
A method get the THDMCache.