Skip to content

Commit

Permalink
init scratcing kin thr from asy
Browse files Browse the repository at this point in the history
  • Loading branch information
@giacomomagni
giacomomagni committed Sep 12, 2023
1 parent 5071e98 commit 65ab2ed
Show file tree
Hide file tree
Showing 2 changed files with 51 additions and 51 deletions.
68 changes: 34 additions & 34 deletions src/yadism/coefficient_functions/asy/f2_nc.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -7,17 +7,17 @@
class AsyLLGluon(pc.NeutralCurrentBaseAsy):
def NLO(self):
def cg_LL_NLO(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return raw_nc.c2g1am0_aq(z) * L

return RSL(cg_LL_NLO, args=[self.L])

def NNLO(self):
def cg_LL_NNLO(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return (raw_nc.c2g2am0_aq2(z) - raw_nc.c2g2am0_aqf(z)) * L**2

Expand All @@ -27,16 +27,16 @@ def cg_LL_NNLO(z, args):
class AsyNLLGluon(pc.NeutralCurrentBaseAsy):
def NLO(self):
def cg_NLL_NLO(z, _args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
return raw_nc.c2g1am0_a0(z)

return RSL(cg_NLL_NLO, args=[self.L])

def NNLO(self):
def cg_NLL_NNLO(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return (raw_nc.c2g2am0_aq(z) - raw_nc.c2g2am0_af(z)) * L

Expand All @@ -46,8 +46,8 @@ def cg_NLL_NNLO(z, args):
class AsyNNLLGluon(pc.NeutralCurrentBaseAsy):
def NNLO(self):
def cg_NNLL_NNLO(z, _args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
return raw_nc.c2g2am0_a0(z)

return RSL(cg_NNLL_NNLO)
Expand All @@ -60,8 +60,8 @@ class AsyNNNLLGluon(pc.NeutralCurrentBaseAsy):
class AsyLLSinglet(pc.NeutralCurrentBaseAsy):
def NNLO(self):
def cps_LL_NNLO(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return (raw_nc.c2ps2am0_aq2(z) - raw_nc.c2ps2am0_aqf(z)) * L**2

Expand All @@ -71,8 +71,8 @@ def cps_LL_NNLO(z, args):
class AsyNLLSinglet(pc.NeutralCurrentBaseAsy):
def NNLO(self):
def cps_NLL_NNLO(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return (raw_nc.c2ps2am0_aq(z) - raw_nc.c2ps2am0_af(z)) * L

Expand All @@ -82,8 +82,8 @@ def cps_NLL_NNLO(z, args):
class AsyNNLLSinglet(pc.NeutralCurrentBaseAsy):
def NNLO(self):
def cps_NNLL_NNLO(z, _args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
return raw_nc.c2ps2am0_a0(z)

return RSL(cps_NNLL_NNLO)
Expand All @@ -96,20 +96,20 @@ class AsyNNNLLSinglet(pc.NeutralCurrentBaseAsy):
class AsyLLNonSinglet(pc.NeutralCurrentBaseAsy):
def NNLO(self):
def cns_LL_NNLO_reg(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return raw_nc.c2ns2am0_aq2(z) * L**2

def cns_LL_NNLO_sing(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return raw_nc.c2ns2bm0_aq2(z) * L**2

def cns_LL_NNLO_loc(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return raw_nc.c2ns2cm0_aq2(z) * L**2

Expand All @@ -119,20 +119,20 @@ def cns_LL_NNLO_loc(z, args):
class AsyNLLNonSinglet(pc.NeutralCurrentBaseAsy):
def NNLO(self):
def cns_NLL_NNLO_reg(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return raw_nc.c2ns2am0_aq(z) * L

def cns_NLL_NNLO_sing(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return raw_nc.c2ns2bm0_aq(z) * L

def cns_NLL_NNLO_loc(z, args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
L = args[0]
return raw_nc.c2ns2cm0_aq(z) * L

Expand All @@ -142,18 +142,18 @@ def cns_NLL_NNLO_loc(z, args):
class AsyNNLLNonSinglet(pc.NeutralCurrentBaseAsy):
def NNLO(self):
def cns_NNLL_NNLO_reg(z, _args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
return raw_nc.c2ns2am0_a0(z)

def cns_NNLL_NNLO_sing(z, _args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
return raw_nc.c2ns2bm0_a0(z)

def cns_NNLL_NNLO_loc(z, _args):
if self.is_below_pair_threshold(z):
return 0.0
# if self.is_below_pair_threshold(z):
# return 0.0
return raw_nc.c2ns2cm0_a0(z)

return RSL(cns_NNLL_NNLO_reg, cns_NNLL_NNLO_sing, cns_NNLL_NNLO_loc)
Expand Down
34 changes: 17 additions & 17 deletions src/yadism/coefficient_functions/asy/partonic_channel.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -138,23 +138,23 @@ def NLO(self):


class NeutralCurrentBaseAsy(PartonicChannelAsy):
def decorator(self, f):
"""
Apply hadronic threshold
Parameters
----------
f : callable
input
Returns
-------
f : callable
output
"""
if self.is_below_pair_threshold(self.ESF.x):
return lambda: pc.RSL()
return f
# def decorator(self, f):
# """
# Apply hadronic threshold

# Parameters
# ----------
# f : callable
# input

# Returns
# -------
# f : callable
# output
# """
# if self.is_below_pair_threshold(self.ESF.x):
# return lambda: pc.RSL()
# return f

def is_below_pair_threshold(self, z):
"""
Expand Down

0 comments on commit 65ab2ed

Please sign in to comment.