Allow pineko to check if the grid contains SV

src/pineko/check.py

@@ -41,3 +41,29 @@ def check_grid_and_eko_compatible(pineappl_grid, operators):
     # x-grid
     if not np.all(in1d(np.unique(operators["targetgrid"]), np.array(x_grid))):
         raise ValueError("x grid in pineappl grid and eko operator are NOT compatible!")
+
+
+def check_grid_contains_sv(pineappl_grid, theory_card):
+    """
+    Raises a `ValueError if the theory_card asks for scale-variations but they are not
+    available in the pineappl grid.
+
+    Parameters
+    ----------
+        pineappl_grid : pineappl.grid.Grid
+            grid
+        theory_card : dict
+            theory card
+    """
+    xir = theory_card["XIR"]
+    xif = theory_card["XIF"]
+    ftr = theory_card["fact_to_ren_scale_ratio"]
+    if xir == 1 and xif == 1 and ftr == 1:
+        return 0
+    order_list = [order.as_tuple() for order in pineappl_grid.orders()]
+    for order in order_list:
+        if order[-1] != 0 or order[-2] != 0:
+            return 0
+    raise ValueError(
+        "Theory card is requesting scale variations but they are not available for this grid!"
+    )

src/pineko/cli/check_scalevar.py

@@ -0,0 +1,17 @@
+# -*- coding: utf-8 -*-
+import click
+import pineappl
+import rich
+
+from .. import check, theory_card
+from ._base import command
+
+
+@command.command("check_scalevar")
+@click.argument("pineappl_path", metavar="PINEAPPL", type=click.Path(exists=True))
+@click.argument("theory_ID", metavar="ID", type=click.Path(exists=True))
+def subcommand(pineappl_path, theory_ID):
+    """Check if PineAPPL grid contains scale variations if theory card needs them"""
+    t_card = theory_card.load(theory_ID)
+    pineappl_grid = pineappl.grid.Grid.read(pineappl_path)
+    check.check_grid_contains_sv(pineappl_grid, t_card)