comment 51743fc for IC asy

NNPDF · Sep 13, 2023 · a2213e4 · a2213e4
1 parent 65ab2ed
commit a2213e4
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 15 deletions.
diff --git a/src/yadism/coefficient_functions/asy/kernels.py b/src/yadism/coefficient_functions/asy/kernels.py
@@ -147,7 +147,6 @@ def generate_intrinsic_asy(esf, nf, pto_evol, ihq):
     cfs = import_pc_module(kind, esf.process)
     m2hq = esf.info.m2hq[ihq - 4]
     if esf.process == "CC":
-        m1sq, m2sq = 0.0, m2hq
         w = kernels.cc_weights(
             esf.info.coupling_constants,
             esf.Q2,
@@ -157,7 +156,6 @@ def generate_intrinsic_asy(esf, nf, pto_evol, ihq):
         )
         weights = {k: v for k, v in w["ns"].items() if abs(k) == ihq}
     else:  # NC
-        m1sq, m2sq = m2hq, m2hq
         if is_pv:
             wVA = esf.info.coupling_constants.get_weight(ihq, esf.Q2, "VA")
             wAV = esf.info.coupling_constants.get_weight(ihq, esf.Q2, "AV")
@@ -174,19 +172,19 @@ def generate_intrinsic_asy(esf, nf, pto_evol, ihq):
     asys = [
         kernels.Kernel(
             weights,
-            cfs.AsyLLIntrinsic(esf, nf, m1sq=m1sq, m2sq=m2sq),
+            cfs.AsyLLIntrinsic(esf, nf, m2hq=m2hq),
         )
     ]
     if pto_evol > 0:
         asys.extend(
             [
                 kernels.Kernel(
                     weights,
-                    cfs.AsyNLLIntrinsicMatching(esf, nf, m1sq=m1sq, m2sq=m2sq),
+                    cfs.AsyNLLIntrinsicMatching(esf, nf, m2hq=m2hq),
                 ),
                 kernels.Kernel(
                     weights,
-                    cfs.AsyNLLIntrinsicLight(esf, nf, m1sq=m1sq, m2sq=m2sq),
+                    cfs.AsyNLLIntrinsicLight(esf, nf, m2hq=m2hq),
                 ),
             ]
         )

diff --git a/src/yadism/coefficient_functions/asy/partonic_channel.py b/src/yadism/coefficient_functions/asy/partonic_channel.py
@@ -5,7 +5,7 @@
 
 from .. import partonic_channel as pc
 from .. import splitting_functions as split
-from ..intrinsic.partonic_channel import NeutralCurrentBase as incb
+# from ..intrinsic.partonic_channel import NeutralCurrentBase as incb
 from ..partonic_channel import RSL
 
 
@@ -23,16 +23,16 @@ class PartonicChannelAsyLLIntrinsic(PartonicChannelAsy):
 
     light_cls = lambda _esf, _nf: None
 
-    def __init__(self, *args, m1sq, m2sq):
-        super().__init__(*args, m2hq=m2sq)
-        self.m1sq = m1sq
-        self.m2sq = m2sq
+    # def __init__(self, *args, m1sq, m2sq):
+    #     super().__init__(*args, m2hq=m2sq)
+    #     self.m1sq = m1sq
+    #     self.m2sq = m2sq
 
-    def convolution_point(self):
-        sigma_pm = self.ESF.Q2 + self.m2sq - self.m1sq
-        delta = incb.kinematic_delta(self.m1sq, self.m2sq, -self.ESF.Q2)
-        eta = 2.0 * self.ESF.Q2 / (sigma_pm + delta)
-        return self.ESF.x / eta
+    # def convolution_point(self):
+    #     sigma_pm = self.ESF.Q2 + self.m2sq - self.m1sq
+    #     delta = incb.kinematic_delta(self.m1sq, self.m2sq, -self.ESF.Q2)
+    #     eta = 2.0 * self.ESF.Q2 / (sigma_pm + delta)
+    #     return self.ESF.x / eta
 
     def LO(self):
         """Return |LO| from light."""