From 0f46f5851f99b703a24f8e5f696c1f88f9827dc7 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Niccol=C3=B2=20Laurenti?= Date: Tue, 9 Jan 2024 21:10:19 +0100 Subject: [PATCH] Fix test sv --- src/eko/evolution_operator/__init__.py | 8 ++++---- src/eko/scale_variations/exponentiated.py | 4 ++-- tests/eko/scale_variations/test_exponentiated.py | 8 ++++---- 3 files changed, 10 insertions(+), 10 deletions(-) diff --git a/src/eko/evolution_operator/__init__.py b/src/eko/evolution_operator/__init__.py index 8092d31a0..92508f164 100644 --- a/src/eko/evolution_operator/__init__.py +++ b/src/eko/evolution_operator/__init__.py @@ -18,7 +18,7 @@ import ekore.anomalous_dimensions.unpolarized.time_like as ad_ut from .. import basis_rotation as br -from .. import interpolation, mellin +from .. import constants, interpolation, mellin from .. import scale_variations as sv from ..kernels import non_singlet as ns from ..kernels import non_singlet_qed as qed_ns @@ -505,7 +505,7 @@ def quad_ker_qed( # scale var exponentiated is directly applied on gamma if sv_mode == sv.Modes.exponentiated: gamma_s = sv.exponentiated.gamma_variation_qed( - gamma_s, order, nf, L, alphaem_running + gamma_s, order, nf, constants.lepton_number(mu2_to), L, alphaem_running ) ker = qed_s.dispatcher( order, @@ -533,7 +533,7 @@ def quad_ker_qed( # scale var exponentiated is directly applied on gamma if sv_mode == sv.Modes.exponentiated: gamma_v = sv.exponentiated.gamma_variation_qed( - gamma_v, order, nf, L, alphaem_running + gamma_v, order, nf, constants.lepton_number(mu2_to), L, alphaem_running ) ker = qed_v.dispatcher( order, @@ -558,7 +558,7 @@ def quad_ker_qed( # scale var exponentiated is directly applied on gamma if sv_mode == sv.Modes.exponentiated: gamma_ns = sv.exponentiated.gamma_variation_qed( - gamma_ns, order, nf, L, alphaem_running + gamma_ns, order, nf, constants.lepton_number(mu2_to), L, alphaem_running ) ker = qed_ns.dispatcher( order, diff --git a/src/eko/scale_variations/exponentiated.py b/src/eko/scale_variations/exponentiated.py index 13d21c34e..491c71395 100644 --- a/src/eko/scale_variations/exponentiated.py +++ b/src/eko/scale_variations/exponentiated.py @@ -46,7 +46,7 @@ def gamma_variation(gamma, order, nf, L): @nb.njit(cache=True) -def gamma_variation_qed(gamma, order, nf, L, alphaem_running): +def gamma_variation_qed(gamma, order, nf, nl, L, alphaem_running): """Adjust the anomalous dimensions with the scale variations. Parameters @@ -71,6 +71,6 @@ def gamma_variation_qed(gamma, order, nf, L, alphaem_running): gamma[1:, 0] = gamma_variation(gamma[1:, 0], order, nf, L) if alphaem_running: if order[1] >= 2: - beta0qed = beta.beta_qed((0, 2), nf) + beta0qed = beta.beta_qed((0, 2), nf, nl) gamma[0, 2] += beta0qed * gamma[0, 1] * L return gamma diff --git a/tests/eko/scale_variations/test_exponentiated.py b/tests/eko/scale_variations/test_exponentiated.py index 79cc6b386..ee76f9ed9 100644 --- a/tests/eko/scale_variations/test_exponentiated.py +++ b/tests/eko/scale_variations/test_exponentiated.py @@ -39,13 +39,13 @@ def test_gamma_ns_qed_fact(): [0.15, 0.0, 0.0, 0.0, 0.0], ] ).transpose() - gamma_ns_as1aem1_0 = gamma_variation_qed(gamma_ns.copy(), (1, 1), 3, 0, True) + gamma_ns_as1aem1_0 = gamma_variation_qed(gamma_ns.copy(), (1, 1), 3, 3, 0, True) np.testing.assert_allclose(gamma_ns_as1aem1_0, gamma_ns) - gamma_ns_as1aem1_1 = gamma_variation_qed(gamma_ns.copy(), (1, 1), 3, -1, True) + gamma_ns_as1aem1_1 = gamma_variation_qed(gamma_ns.copy(), (1, 1), 3, 3, -1, True) np.testing.assert_allclose(gamma_ns_as1aem1_1, gamma_ns) - gamma_ns_as2aem2_1 = gamma_variation_qed(gamma_ns.copy(), (2, 2), 3, -1, True) + gamma_ns_as2aem2_1 = gamma_variation_qed(gamma_ns.copy(), (2, 2), 3, 3, -1, True) assert gamma_ns_as2aem2_1[2, 0] < gamma_ns[2, 0] assert gamma_ns_as2aem2_1[0, 2] > gamma_ns[0, 2] # beta0qed < 0 - gamma_ns_as4aem2_1 = gamma_variation_qed(gamma_ns.copy(), (4, 2), 3, -1, True) + gamma_ns_as4aem2_1 = gamma_variation_qed(gamma_ns.copy(), (4, 2), 3, 3, -1, True) gamma_ns_N3LO_1 = gamma_variation(gamma_ns.copy()[1:, 0], (4, 0), 3, -1) np.testing.assert_allclose(gamma_ns_as4aem2_1[1:, 0], gamma_ns_N3LO_1)