1813{
1826
1827 double betalambda1 = 6.0*la1Q*la1Q + 2.0*la3Q*la3Q + 2.0*la3Q*la4Q + la4Q*la4Q + la5Q*la5Q
1828 + 6.0*la1Q*Yb1Q*Yb1Q + 2.0*la1Q*Ytau1Q*Ytau1Q
1829 - 6.0*Yb1Q*Yb1Q*Yb1Q*Yb1Q - 2.0*Ytau1Q*Ytau1Q*Ytau1Q*Ytau1Q;
1830 double betalambda2 = 6.0*la2Q*la2Q + 2.0*la3Q*la3Q + 2.0*la3Q*la4Q + la4Q*la4Q + la5Q*la5Q
1831 + 6.0*la2Q*Yb2Q*Yb2Q + 2.0*la2Q*Ytau2Q*Ytau2Q + 6.0*la2Q*YtQ*YtQ
1832 - 6.0*Yb2Q*Yb2Q*Yb2Q*Yb2Q - 2.0*Ytau2Q*Ytau2Q*Ytau2Q*Ytau2Q - 6.0*YtQ*YtQ*YtQ*YtQ;
1833 double betalambda5 = la5Q*(la1Q + la2Q + 4.0*la3Q + 6.0*la4Q)
1834 + 3.0*la5Q*Yb1Q*Yb1Q + la5Q*Ytau1Q*Ytau1Q + la5Q*(3.0*Yb2Q*Yb2Q + Ytau2Q*Ytau2Q + 3.0*YtQ*YtQ)
1835 - 6.0*Yb1Q*Yb1Q*Yb2Q*Yb2Q - 2.0*Ytau1Q*Ytau1Q*Ytau2Q*Ytau2Q;
1836
1837 double uniNLO18a = -la1Q/(16.0*M_PI);
1838 double uniNLO18b = 3.0*betalambda1/(256.0*M_PI*M_PI*M_PI);
1839 gslpp::complex uniNLO18c = (gslpp::complex::i()*M_PI-1.0)*(la1Q*la1Q+la5Q*la5Q)/(256.0*M_PI*M_PI*M_PI);
1840 gslpp::complex uniNLO18d = -la1Q/(512.0*M_PI*M_PI*M_PI) * WFRc1;
1841
1842 double uniNLO19a = -la2Q/(16.0*M_PI);
1843 double uniNLO19b = 3.0*betalambda2/(256.0*M_PI*M_PI*M_PI);
1844 gslpp::complex uniNLO19c = (gslpp::complex::i()*M_PI-1.0)*(la2Q*la2Q+la5Q*la5Q)/(256.0*M_PI*M_PI*M_PI);
1845 gslpp::complex uniNLO19d = -la2Q/(512.0*M_PI*M_PI*M_PI) * (-WFRc1+2.0*WFRc2);
1846
1847 double uniNLO20a = -la5Q/(16.0*M_PI);
1848 double uniNLO20b = 3.0*betalambda5/(256.0*M_PI*M_PI*M_PI);
1849 gslpp::complex uniNLO20c = (gslpp::complex::i()*M_PI-1.0)*(la1Q+la2Q)*la5Q/(256.0*M_PI*M_PI*M_PI);
1850 gslpp::complex uniNLO20d = -la4Q/(512.0*M_PI*M_PI*M_PI) * WFRc2;
1851
1852 gslpp::complex uniA=uniNLO18a+uniNLO18b+uniNLO18c+uniNLO18d;
1853 gslpp::complex uniB=uniNLO19a+uniNLO19b+uniNLO19c+uniNLO19d;
1854 gslpp::complex uniC=uniNLO20a+uniNLO20b+uniNLO20c+uniNLO20d;
1855
1856 return (0.5*(uniA+uniB-sqrt((uniA-uniB)*(uniA-uniB)+4.0*uniC*uniC)) - 0.5*gslpp::complex::i()).abs();
1857}
THDMcache * getMyTHDMCache() const
A method get the THDMCache.