Merge pull request #368 from NNPDF/fix-pol-as3

Fix NNLO polarized
NNPDF · May 15, 2024 · c0ce0f0 · c0ce0f0
2 parents 8292f30 + 6d989f2
commit c0ce0f0
Show file tree

Hide file tree

Showing 3 changed files with 17 additions and 18 deletions.
diff --git a/src/ekore/anomalous_dimensions/polarized/space_like/as3.py b/src/ekore/anomalous_dimensions/polarized/space_like/as3.py
@@ -177,9 +177,9 @@ def gamma_gq(N, nf, cache):
 
     """
     S1 = c.get(c.S1, cache, N)
-    S2 = c.get(c.S1, cache, N)
-    S3 = c.get(c.S1, cache, N)
-    S4 = c.get(c.S1, cache, N)
+    S2 = c.get(c.S2, cache, N)
+    S3 = c.get(c.S3, cache, N)
+    S4 = c.get(c.S4, cache, N)
     B3 = (-(S1**3) - 3 * S1 * S2 - 2 * S3) / N
     B4 = (S1**4 + 6 * S1**2 * S2 + 3 * S2**2 + 8 * S1 * S3 + 6 * S4) / N
     E1 = S1 / N**2 + (-zeta2 + S2) / N
@@ -256,8 +256,8 @@ def gamma_ps(N, nf, cache):
 
     """
     S1 = c.get(c.S1, cache, N)
-    S2 = c.get(c.S1, cache, N)
-    S3 = c.get(c.S1, cache, N)
+    S2 = c.get(c.S2, cache, N)
+    S3 = c.get(c.S3, cache, N)
     B3 = (-(S1**3) - 3 * S1 * S2 - 2 * S3) / N
     B31 = (
         -((1 / (1 + N) + S1) ** 3)

diff --git a/src/ekore/operator_matrix_elements/polarized/space_like/as2.py b/src/ekore/operator_matrix_elements/polarized/space_like/as2.py
@@ -172,7 +172,7 @@ def A_hg(n, cache, L):
             -2
             * CF
             * (
-                3
+                6
                 * (
                     4
                     + 8 * n
@@ -187,7 +187,7 @@ def A_hg(n, cache, L):
                 * n**2
                 * (1 + n) ** 2
                 * (2 + n)
-                * (2 + 3 * n + 3 * n**3)
+                * (2 + 3 * n + 3 * n**2)
                 * np.pi**2
                 + 12 * n**2 * (1 + n) ** 3 * (-36 - 22 * n - 2 * n**2 + n**3) * S1
                 + 12 * n**2 * (1 + n) ** 3 * (-2 + 3 * n + 3 * n**2) * S1**2
@@ -206,7 +206,6 @@ def A_hg(n, cache, L):
             + CA
             * (
                 -24
-                * (1 + n) ** 3
                 * (
                     4
                     + 12 * n

diff --git a/tests/ekore/operator_matrix_elements/polarized/space_like/test_nnlo.py b/tests/ekore/operator_matrix_elements/polarized/space_like/test_nnlo.py
@@ -27,18 +27,18 @@ def test_quark_number_conservation():
 def test_hg():
     refs = {
         0: [
-            -242.3869306886845,
-            -480.6122566782359,
-            -814.6150038529145,
-            -1244.9445223204148,
-            -1771.5995558739958,
+            -14.120648885210752,
+            -15.093091588201201,
+            -14.014668171450229,
+            -12.815826533022285,
+            -11.761139164575077,
         ],
         10: [
-            -21.21409118251165,
-            -295.28287396218656,
-            -664.653836052461,
-            -1121.4595082763883,
-            -1667.8336593381266,
+            207.0521906209621,
+            170.2362911278482,
+            135.9464996290033,
+            110.6691875110041,
+            92.0047573712939,
         ],
     }
     for L, vals in refs.items():