diff --git a/crates/ekore/src/operator_matrix_elements/unpolarized/spacelike.rs b/crates/ekore/src/operator_matrix_elements/unpolarized/spacelike.rs index fc17da638..3c57d33ef 100644 --- 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 index c01f54932..40cded5cb 100644 --- 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 { } /// |NNLO| singlet |OME|. -pub fn A_singlet(c: &mut Cache, nf: u8, L: f64, is_msbar: bool) -> [[Complex; 3]; 3] { +pub fn A_singlet(c: &mut Cache, nf: u8, L: f64, is_msbar_mass: bool) -> [[Complex; 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); }