intorduce PTO_match

extras/lh_bench_23/cfg.py

@@ -31,6 +31,7 @@
     ),
     xif=1.0,
     n3lo_ad_variation=(0, 0, 0, 0),
+    matching_order=[2, 0],
 )
 
 

src/eko/evolution_operator/grid.py

@@ -54,6 +54,7 @@ def __init__(
         intrinsic_flavors: list,
         xif: float,
         n3lo_ad_variation: tuple,
+        matching_order: Order,
         configs: Configs,
         debug: Debug,
         atlas: Atlas,
@@ -81,6 +82,7 @@ def __init__(
         config["debug_skip_non_singlet"] = debug.skip_non_singlet
         config["polarized"] = configs.polarized
         config["time_like"] = configs.time_like
+        config["matching_order"] = matching_order
 
         if order == (1, 0) and method != "iterate-exact":
             logger.warning("Evolution: In LO we use the exact solution always!")

src/eko/evolution_operator/operator_matrix_element.py

@@ -220,7 +220,7 @@ def __init__(self, config, managers, nf, q2, is_backward, L, is_msbar):
         self.L = L
         self.is_msbar = is_msbar
         # Note for the moment only QCD matching is implemented
-        self.order = (self.order[0] - 1, self.order[1])
+        self.order = tuple(config["matching_order"])
 
     @property
     def labels(self):

src/eko/io/runcards.py

@@ -49,6 +49,8 @@ class TheoryCard(DictLike):
     """Ratio between factorization scale and process scale."""
     n3lo_ad_variation: N3LOAdVariation
     """|N3LO| anomalous dimension variation: ``(gg_var, gq_var, qg_var, qq_var)``."""
+    matching_order: Order
+    """Matching conditions perturbative order tuple, ``(QCD, QED)``."""
 
 
 @dataclass
@@ -210,6 +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"]])
 
         return TheoryCard.from_dict(new)
 

src/eko/runner/legacy.py

@@ -77,6 +77,7 @@ def __init__(
             couplings=cs,
             interpol_dispatcher=bfd,
             n3lo_ad_variation=new_theory.n3lo_ad_variation,
+            matching_order=new_theory.matching_order,
         )
 
         with EKO.create(path) as builder:

src/eko/runner/parts.py

@@ -88,6 +88,7 @@ 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,
     )
 
 

src/ekobox/cards.py

@@ -27,6 +27,7 @@
     ),
     xif=1.0,
     n3lo_ad_variation=(0, 0, 0, 0),
+    matching_order=[0, 0],
 )
 
 _operator = dict(