From 12b3ed330d74a0123150a2d5c670801baf782900 Mon Sep 17 00:00:00 2001
From: Radonirinaunimi <rrabeman@nikhef.nl>
Date: Thu, 23 Feb 2023 18:11:58 +0100
Subject: [PATCH] Fix threshold_ratio types in GLS

---
 src/nnusf/plot/fit.py        | 2 +-
 src/nnusf/sffit/check_gls.py | 3 ++-
 2 files changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/nnusf/plot/fit.py b/src/nnusf/plot/fit.py
index 0633712c..ee601d4f 100644
--- a/src/nnusf/plot/fit.py
+++ b/src/nnusf/plot/fit.py
@@ -235,7 +235,7 @@ def gls_sum_rules(**kwargs):
         marker="s",
         markersize=11,
         mfc="w",
-        label=r"$\rm{NN~Predictions}$",
+        label=r"$\rm{NNSF}\nu$",
         capsize=6,
         zorder=0,
     )
diff --git a/src/nnusf/sffit/check_gls.py b/src/nnusf/sffit/check_gls.py
index 62128d32..eee4c657 100644
--- a/src/nnusf/sffit/check_gls.py
+++ b/src/nnusf/sffit/check_gls.py
@@ -1,4 +1,5 @@
 # -*- coding: utf-8 -*-
+import math
 import numpy as np
 from eko.couplings import Couplings
 from eko.io import types as ekotypes
@@ -87,7 +88,7 @@ def alphas_eko(q2_value, order=3):
 
         # set heavy quark masses and their threshold ratios
         heavy_quark_masses = np.power([1.51, 4.92, 172.0], 2)
-        thresholds_ratios = ekotypes.MatchingScales(c=1.0, b=1.0, t=1.0)
+        thresholds_ratios = [1.0, 1.0, 1.0]
 
         # set (QCD,QED) perturbative order
         order = (order, 1)