rust: Activate NNLO OME

NNPDF · Sep 10, 2024 · 8f09ea0 · 8f09ea0
1 parent 9b498b7
commit 8f09ea0
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 2 deletions.
diff --git a/crates/ekore/src/operator_matrix_elements/unpolarized/spacelike.rs b/crates/ekore/src/operator_matrix_elements/unpolarized/spacelike.rs
@@ -35,6 +35,10 @@ pub fn A_singlet(
     if matching_order_qcd >= 1 {
         A_s[0] = as1::A_singlet(c, nf, L);
     }
+    if matching_order_qcd >= 2 {
+        // TODO recover MSbar mass
+        A_s[1] = as2::A_singlet(c, nf, L, false);
+    }
     A_s
 }
 
@@ -55,5 +59,8 @@ pub fn A_non_singlet(
     if matching_order_qcd >= 1 {
         A_ns[0] = as1::A_ns(c, nf, L);
     }
+    if matching_order_qcd >= 2 {
+        A_ns[1] = as2::A_ns(c, nf, L);
+    }
     A_ns
 }
diff --git a/crates/ekore/src/operator_matrix_elements/unpolarized/spacelike/as2.rs b/crates/ekore/src/operator_matrix_elements/unpolarized/spacelike/as2.rs
@@ -260,14 +260,14 @@ pub fn A_gg(c: &mut Cache, _nf: u8, L: f64) -> Complex<f64> {
 }
 
 /// |NNLO| singlet |OME|.
-pub fn A_singlet(c: &mut Cache, nf: u8, L: f64, is_msbar: bool) -> [[Complex<f64>; 3]; 3] {
+pub fn A_singlet(c: &mut Cache, nf: u8, L: f64, is_msbar_mass: bool) -> [[Complex<f64>; 3]; 3] {
     let A_hq_2 = A_hq_ps(c, nf, L);
     let A_qq_2 = A_qq_ns(c, nf, L);
     let mut A_hg_2 = A_hg(c, nf, L);
     let A_gq_2 = A_gq(c, nf, L);
     let mut A_gg_2 = A_gg(c, nf, L);
 
-    if is_msbar {
+    if is_msbar_mass {
         A_hg_2 -= 2.0 * 4.0 * CF * as1::A_hg(c, nf, 1.0);
         A_gg_2 -= 2.0 * 4.0 * CF * as1::A_gg(c, nf, 1.0);
     }