diff --git a/src/ekomark/benchmark/external/pegasus_utils.py b/src/ekomark/benchmark/external/pegasus_utils.py
index e227222c7..9ac6c6273 100644
--- a/src/ekomark/benchmark/external/pegasus_utils.py
+++ b/src/ekomark/benchmark/external/pegasus_utils.py
@@ -82,7 +82,7 @@ def compute_pegasus_data(theory, operators, skip_pdfs, rotate_to_evolution_basis
 
     # run pegaus
     out_tabs = {}
-    for mu2 in np.array(operators["mugrid"]) ** 2:
+    for mu2 in operators["Q2grid"]:
         tab = {}
         for x in target_xgrid:
             # last two numbers are the min and max pid to calculate,