Test only alpha_em

NNPDF · Sep 21, 2023 · d2ecc69 · d2ecc69
1 parent 0e62fa3
commit d2ecc69
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 8 deletions.
diff --git a/validphys2/src/validphys/photon/compute.py b/validphys2/src/validphys/photon/compute.py
@@ -502,6 +502,8 @@ def rge(_t, alpha, beta_qcd_vec, beta_qcd_mix, beta_qed_vec, beta_qed_mix):
     we need alpha at very low scale, below the Landau pole of alpha_s.
     This makes the alpha evolution crash. For this reason we evolve alpha
     without the mixed terms, i.e. decoupling it from alpha_s.
+    rge_qcd is set to zero since we don't want an alpha_s value that is crap
+    even if it is not used.
     We left the betaQCD and beta_mix part implemented in the case we find a
     solution. Anyway, it is not slowing down the code.
 
@@ -514,7 +516,7 @@ def rge(_t, alpha, beta_qcd_vec, beta_qcd_mix, beta_qed_vec, beta_qed_mix):
             + np.sum([alpha[0] ** (k + 1) * b for k, b in enumerate(beta_qcd_vec[1:])])
             # + alpha[1] * beta_qcd_mix
         )
-    )
+    ) * 0.0
     rge_qed = (
         -(alpha[1] ** 2)
         * beta_qed_vec[0]

diff --git a/validphys2/src/validphys/tests/photon/test_compute.py b/validphys2/src/validphys/tests/photon/test_compute.py
@@ -123,27 +123,27 @@ def test_couplings_exa():
         coupl_ref = [theory["alphas"], theory["alphaqed"]]
         for q in [5, 10, 20, 50, 80, 100, 200]:
             np.testing.assert_allclose(
-                alpha.couplings_fixed_flavor(q, coupl_ref, theory["Qref"], 5),
-                alpha.couplings_variable_flavor(q),
+                alpha.couplings_fixed_flavor(q, coupl_ref, theory["Qref"], 5)[1],
+                alpha.couplings_variable_flavor(q)[1],
                 rtol=1e-10,
             )
             np.testing.assert_allclose(
-                alpha.couplings_fixed_flavor(q, coupl_ref, theory["Qref"], 5),
+                alpha.couplings_fixed_flavor(q, coupl_ref, theory["Qref"], 5)[1],
                 eko_alpha.compute_exact_alphaem_running(
                     np.array(coupl_ref) / (4 * np.pi), 5, theory["Qref"] ** 2, q**2
-                )
+                )[1]
                 * 4
                 * np.pi,
                 rtol=1e-7,
             )
         for q in [1, 2, 3, 4]:
             np.testing.assert_allclose(
-                alpha.couplings_variable_flavor(q), eko_alpha.a_em(q**2) * 4 * np.pi, rtol=1e-7
+                alpha.couplings_variable_flavor(q)[1], eko_alpha.a_em(q**2) * 4 * np.pi, rtol=1e-7
             )
         for nf in range(3, theory["MaxNfAs"]):
             np.testing.assert_allclose(
-                alpha.couplings_thresh[nf],
-                eko_alpha.a(mass_list[nf - 3] ** 2, nf) * 4 * np.pi,
+                alpha.couplings_thresh[nf][1],
+                eko_alpha.a_em(mass_list[nf - 3] ** 2, nf) * 4 * np.pi,
                 rtol=3e-7,
             )