1388{
1401
1402 double betalambda1 = 6.0*la1Q*la1Q + 2.0*la3Q*la3Q + 2.0*la3Q*la4Q + la4Q*la4Q + la5Q*la5Q
1403 + 6.0*la1Q*Yb1Q*Yb1Q + 2.0*la1Q*Ytau1Q*Ytau1Q
1404 - 6.0*Yb1Q*Yb1Q*Yb1Q*Yb1Q - 2.0*Ytau1Q*Ytau1Q*Ytau1Q*Ytau1Q;
1405 double betalambda2 = 6.0*la2Q*la2Q + 2.0*la3Q*la3Q + 2.0*la3Q*la4Q + la4Q*la4Q + la5Q*la5Q
1406 + 6.0*la2Q*Yb2Q*Yb2Q + 2.0*la2Q*Ytau2Q*Ytau2Q + 6.0*la2Q*YtQ*YtQ
1407 - 6.0*Yb2Q*Yb2Q*Yb2Q*Yb2Q - 2.0*Ytau2Q*Ytau2Q*Ytau2Q*Ytau2Q - 6.0*YtQ*YtQ*YtQ*YtQ;
1408 double betalambda4 = la1Q*la4Q + la2Q*la4Q + 4.0*la3Q*la4Q + 2.0*la4Q*la4Q + 4.0*la5Q*la5Q
1409 + 3.0*la4Q*Yb1Q*Yb1Q + 3.0*la4Q*Yb2Q*Yb2Q + la4Q*Ytau1Q*Ytau1Q + la4Q*Ytau2Q*Ytau2Q + 3.0*la4Q*YtQ*YtQ
1410 - 2.0*Ytau1Q*Ytau1Q*Ytau2Q*Ytau2Q - 6.0*Yb1Q*Yb1Q*(Yb2Q*Yb2Q - YtQ*YtQ);
1411
1412 double uniNLO6a = -la1Q/(16.0*M_PI);
1413 double uniNLO6b = 3.0*betalambda1/(256.0*M_PI*M_PI*M_PI);
1414 gslpp::complex uniNLO6c = (gslpp::complex::i()*M_PI-1.0)*(la1Q*la1Q+la4Q*la4Q)/(256.0*M_PI*M_PI*M_PI);
1415 gslpp::complex uniNLO6d = -la1Q/(512.0*M_PI*M_PI*M_PI) * WFRc1;
1416
1417 double uniNLO7a = -la2Q/(16.0*M_PI);
1418 double uniNLO7b = 3.0*betalambda2/(256.0*M_PI*M_PI*M_PI);
1419 gslpp::complex uniNLO7c = (gslpp::complex::i()*M_PI-1.0)*(la2Q*la2Q+la4Q*la4Q)/(256.0*M_PI*M_PI*M_PI);
1420 gslpp::complex uniNLO7d = -la2Q/(512.0*M_PI*M_PI*M_PI) * (-WFRc1+2.0*WFRc2);
1421
1422 double uniNLO9a = -la4Q/(16.0*M_PI);
1423 double uniNLO9b = 3.0*betalambda4/(256.0*M_PI*M_PI*M_PI);
1424 gslpp::complex uniNLO9c = (gslpp::complex::i()*M_PI-1.0)*(la1Q+la2Q)*la4Q/(256.0*M_PI*M_PI*M_PI);
1425 gslpp::complex uniNLO9d = -la5Q/(512.0*M_PI*M_PI*M_PI) * WFRc2;
1426
1427 gslpp::complex uniA=uniNLO6a+uniNLO6b+uniNLO6c+uniNLO6d;
1428 gslpp::complex uniB=uniNLO7a+uniNLO7b+uniNLO7c+uniNLO7d;
1429 gslpp::complex uniC=uniNLO9a+uniNLO9b+uniNLO9c+uniNLO9d;
1430
1431 return (0.5*(uniA+uniB+sqrt((uniA-uniB)*(uniA-uniB)+4.0*uniC*uniC)) - 0.5*gslpp::complex::i()).abs();
1432}
THDMcache * getMyTHDMCache() const
A method get the THDMCache.