From 32935df894537497cc53ea6a2aa1b09a37621e23 Mon Sep 17 00:00:00 2001
From: Felix Hekhorn <felix.hekhorn@uni-tuebingen.de>
Date: Thu, 11 Jan 2024 15:55:28 +0200
Subject: [PATCH 1/2] Update couplings docs
---
src/eko/couplings.py | 19 +++++++++++--------
1 file changed, 11 insertions(+), 8 deletions(-)
diff --git a/src/eko/couplings.py b/src/eko/couplings.py
index b0d39a070..06b4a5b33 100644
--- a/src/eko/couplings.py
+++ b/src/eko/couplings.py
@@ -388,10 +388,9 @@ class Couplings:
- although, we only provide methods for
:math:`a_i = \frac{\alpha_i(\mu^2)}{4\pi}` the reference value has to be
given in terms of :math:`\alpha_i(\mu_0^2)` due to legacy reasons
- - the ``order`` refers to the perturbative order of the beta function, thus
- ``order=(0,0)`` means leading order beta function, means evolution with :math:`\beta_qcd_as2`,
- means running at 1-loop - so there is a natural mismatch between ``order`` and the
- number of loops by one unit
+ - the ``order`` refers to the perturbative order of the beta functions, i.e.
+ the number of loops for QCD and QED respectively. QCD is always running with
+ at least 1 loop, while QED might not run at all (so 0 loop).
Normalization is given by :cite:`Herzog:2017ohr`:
@@ -407,9 +406,9 @@ class Couplings:
couplings :
reference configuration
order :
- Evaluated order of the beta function: ``0`` = LO, ...
+ Number of loops in beta functions (QCD, QED)
method :
- Applied method to solve the beta function
+ Applied method to solve the beta functions
masses :
list with quark masses squared
hqm_scheme :
@@ -436,9 +435,13 @@ def assert_positive(name, var):
assert_positive("alpha_em_ref", couplings.alphaem)
assert_positive("scale_ref", couplings.scale)
if order[0] not in [1, 2, 3, 4]:
- raise NotImplementedError("a_s beyond N3LO is not implemented")
+ raise NotImplementedError(
+ "QCD order has to be an integer between 1 (LO) and 4 (N3LO)"
+ )
if order[1] not in [0, 1, 2]:
- raise NotImplementedError("a_em beyond NLO is not implemented")
+ raise NotImplementedError(
+ "QED order has to be an integer between 0 (no running) and 2 (NLO)"
+ )
self.order = tuple(order)
if method.value not in ["expanded", "exact"]:
raise ValueError(f"Unknown method {method.value}")
From f980362c9801c7281dfd3c62e84b7631b070688e Mon Sep 17 00:00:00 2001
From: niccolo <niclaurenti@gmail.com>
Date: Thu, 11 Jan 2024 15:39:17 +0100
Subject: [PATCH 2/2] Fix leftover
---
src/eko/couplings.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/eko/couplings.py b/src/eko/couplings.py
index 06b4a5b33..f0849aaf8 100644
--- a/src/eko/couplings.py
+++ b/src/eko/couplings.py
@@ -529,7 +529,7 @@ def compute_exact_alphaem_running(self, a_ref, nf, scale_from, scale_to):
Parameters
----------
- as_ref : numpy.ndarray
+ a_ref : numpy.ndarray
reference alpha_s and alpha
nf : int
value of nf for computing alpha_i
@@ -615,7 +615,7 @@ def compute_exact_fixed_alphaem(self, a_ref, nf, scale_from, scale_to):
Parameters
----------
- as_ref : numpy.ndarray
+ a_ref : numpy.ndarray
reference alpha_s and alpha
nf : int
value of nf for computing alpha_i