From cf7f1c5fc243b9eaef562592ad62c3db8fc62367 Mon Sep 17 00:00:00 2001 From: Luca Mantani Date: Thu, 28 Nov 2024 11:58:09 +0100 Subject: [PATCH] moved cuts --- src/smefit/rge.py | 13 +++++++++---- 1 file changed, 9 insertions(+), 4 deletions(-) diff --git a/src/smefit/rge.py b/src/smefit/rge.py index c94b54ea..ee309c57 100644 --- a/src/smefit/rge.py +++ b/src/smefit/rge.py @@ -358,7 +358,7 @@ def RGEevolve(self, wcs, scale): return self.map_to_smefit(wc_final, scale) -def load_scales(datasets, theory_path, default_scale=1e3): +def load_scales(datasets, theory_path, default_scale=1e3, cutoff_scale=None): """ Load the energy scales for the datasets. @@ -370,6 +370,8 @@ def load_scales(datasets, theory_path, default_scale=1e3): path to the theory files default_scale: float default scale to use if the dataset does not have a scale + cutoff_scale: float + cutoff scale to use for the scales Returns ------- @@ -397,6 +399,9 @@ def load_scales(datasets, theory_path, default_scale=1e3): _logger.info(f"Loaded scales for dataset {dataset}: {dataset_scales}") + if cutoff_scale is not None: + scales = [scale for scale in scales if scale < cutoff_scale] + return scales @@ -436,9 +441,9 @@ def load_rge_matrix( return rgemat.values, operators_to_keep elif obs_scale == "dynamic": - scales = load_scales(datasets, theory_path, default_scale=init_scale) - if cutoff_scale is not None: - scales = [scale for scale in scales if scale < cutoff_scale] + scales = load_scales( + datasets, theory_path, default_scale=init_scale, cutoff_scale=cutoff_scale + ) operators_to_keep = {} rgemat = []