From 547911a3182af1c998055ecfde36744bbcee7396 Mon Sep 17 00:00:00 2001 From: giacomomagni Date: Wed, 8 Nov 2023 09:05:37 +0100 Subject: [PATCH] hard code QED matching to 0 --- src/eko/io/runcards.py | 4 ++-- src/eko/runner/parts.py | 6 +++--- tests/eko/runner/test_parts.py | 2 +- 3 files changed, 6 insertions(+), 6 deletions(-) diff --git a/src/eko/io/runcards.py b/src/eko/io/runcards.py index e9260a6eb..7dbbe2b34 100644 --- a/src/eko/io/runcards.py +++ b/src/eko/io/runcards.py @@ -212,8 +212,8 @@ def new_theory(self): new["xif"] = old["XIF"] new["n3lo_ad_variation"] = old.get("n3lo_ad_variation", (0, 0, 0, 0)) - # here PTO: 0 means truly LO - new["matching_order"] = old.get("PTO_matching", [old["PTO"], old["QED"]]) + # here PTO: 0 means truly LO, no QED matching is available so far. + new["matching_order"] = old.get("PTO_matching", [old["PTO"], 0]) return TheoryCard.from_dict(new) diff --git a/src/eko/runner/parts.py b/src/eko/runner/parts.py index 128a2553b..bbbf36408 100644 --- a/src/eko/runner/parts.py +++ b/src/eko/runner/parts.py @@ -88,9 +88,9 @@ def evolve_configs(eko: EKO) -> dict: n_integration_cores=ocard.configs.n_integration_cores, ModSV=ocard.configs.scvar_method, n3lo_ad_variation=tcard.n3lo_ad_variation, - matching_order=tcard.matching_order - if tcard.matching_order is not None - else (tcard.order[0] - 1, tcard.order[1] - 1), + matching_order=tcard.matching_order if tcard.matching_order is not None + # Here order is shifted by one, no QED matching is available so far. + else (tcard.order[0] - 1, 0), ) diff --git a/tests/eko/runner/test_parts.py b/tests/eko/runner/test_parts.py index 39c3bc843..2ec217dd7 100644 --- a/tests/eko/runner/test_parts.py +++ b/tests/eko/runner/test_parts.py @@ -13,7 +13,7 @@ def test_evolve_configs(eko_factory): e43 = eko_factory.get({}) assert e43.theory_card.order == (4, 3) p43 = parts.evolve_configs(e43) - assert p43["matching_order"] == (3, 2) + assert p43["matching_order"] == (3, 0) # QCD@N3LO + QED@N2LO w/ matching_order eko_factory.theory.matching_order = (3, 0) e43b = eko_factory.get({})