Skip to content

Commit

Permalink
Merge branch 'numerical_fonll' into massive_n3lo
Browse files Browse the repository at this point in the history
  • Loading branch information
giacomomagni committed Sep 26, 2023
2 parents 5790ca1 + eef2355 commit fe4c7be
Show file tree
Hide file tree
Showing 2 changed files with 0 additions and 42 deletions.
1 change: 0 additions & 1 deletion src/yadism/coefficient_functions/asy/partonic_channel.py
Original file line number Diff line number Diff line change
Expand Up @@ -12,7 +12,6 @@ class PartonicChannelAsy(pc.PartonicChannel):
"""Massless limit of a coeficient function."""

def __init__(self, *args, m2hq, n3lo_cf_variation=0):
self.m2hq = m2hq
super().__init__(*args)
self.L = np.log(self.ESF.Q2 / m2hq)
self.n3lo_cf_variation = n3lo_cf_variation
Expand Down
41 changes: 0 additions & 41 deletions tests/yadism/cf/test_asy.py
Original file line number Diff line number Diff line change
Expand Up @@ -49,47 +49,6 @@ def test_h_g():
)


MISSING_NC_MAP = {
h_f2_nc.NonSinglet: [
f_f2_nc.AsyLLNonSinglet,
f_f2_nc.AsyNLLNonSinglet,
f_f2_nc.AsyNNLLNonSinglet,
],
h_fl_nc.NonSinglet: [
f_fl_nc.AsyLLNonSinglet,
f_fl_nc.AsyNLLNonSinglet,
f_fl_nc.AsyNNLLNonSinglet,
],
}


def test_missing_nnlo():
Q2 = 2e4
esf = MockESF("F2_charm", 0.1, Q2)
interp = interpolation.InterpolatorDispatcher(
interpolation.lambertgrid(50), 4, False
)
m2hq = 2
for nf in [3, 4]:
for z in [1e-4, 1e-3]:
for icls, iasyclss in MISSING_NC_MAP.items():
iobj = icls(esf, nf, m2hq=m2hq)
a = (
conv.convolute_vector(iobj.NNLO(), interp, z)[0]
if iobj.NNLO()
else 0.0
)
b = 0.0
if iasyclss:
for iasycls in iasyclss:
iasyobj = iasycls(esf, nf, m2hq=m2hq)
if iasyobj.NNLO():
b += conv.convolute_vector(iasyobj.NNLO(), interp, z)[0]
np.testing.assert_allclose(
a, b, rtol=7e-3, atol=3e-2, err_msg=f"{nf=},{z=}"
)


INTRINSIC_NC_MAP = {
i_f2_nc.Splus: [
f_f2_nc.AsyLLIntrinsic,
Expand Down

0 comments on commit fe4c7be

Please sign in to comment.