Use exact expression for constant in pseudoscalar -> gg NNLO correction

based on that in NMSSMcalc
FlexibleSUSY · Feb 4, 2016 · f4335ab · f4335ab
1 parent d24935d
commit f4335ab
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/meta/EffectiveCouplings.m b/meta/EffectiveCouplings.m
@@ -95,7 +95,8 @@
            scalarFactor = Parameters`CreateLocalConstRefs[nloQCD + nnloQCD + nnnloQCD] <> "\n" <> scalarFactor;
            (* NLO, NNLO and NNNLO contributions to pseudoscalar coupling *)
            nloQCD = (97 / 4 - 7 / 6 Symbol["Nf"]) SARAH`strongCoupling^2 / (4 Pi^2);
-           nnloQCD = (171.544 + 5 Symbol["l"]) SARAH`strongCoupling^4 / (16 Pi^4);
+           nnloQCD = (237311 / 864 - 529 Zeta[2] / 24 - 445 Zeta[3] / 8 + 5 Symbol["l"]);
+           nnloQCD = nnloQCD SARAH`strongCoupling^4 / (16 Pi^4);
            nnnloQCD = 0;
            pseudoscalarFactor = pseudoscalarFactor <> "const double nlo_qcd = "
                                 <> CConversion`RValueToCFormString[nloQCD] <> ";\n";