nnpdf-data 0.0.1__py3-none-any.whl → 0.0.3__py3-none-any.whl

nnpdf_data/process_options.py

@@ -1,7 +1,7 @@
 """
-    Module to hold process dependent options
+Module to hold process dependent options
 
-    Only variables included in the `_Vars` enum and processes included in the ``Processes`` dictionary are allowed.
+Only variables included in the `_Vars` enum and processes included in the ``Processes`` dictionary are allowed.
 """
 
 import dataclasses
@@ -16,6 +16,7 @@ class _Vars:
     Q2 = "Q2"
     Q = "Q"
     y = "y"
+    abs_y = "abs_y"
     pT = "pT"
     ET = "ET"
     sqrts = "sqrts"
@@ -26,6 +27,7 @@ class _Vars:
     y_t = "y_t"
     y_ttBar = "y_ttBar"
     m_t2 = "m_t2"
+    m_t = "m_t"
     pT_t = "pT_t"
     m_ttBar = "m_ttBar"
     eta = "eta"
@@ -33,11 +35,16 @@ class _Vars:
     m_W2 = "m_W2"
     m_Z2 = "m_Z2"
     m_V2 = "m_V2"
+    m_W = "m_W"
+    m_Z = "m_Z"
+    M2 = "M2"
     abs_eta_1 = "abs_eta_1"
     abs_eta_2 = "abs_eta_2"
     eta_1 = "eta_1"
     eta_2 = "eta_2"
+    m_ll = "m_ll"
     m_ll2 = "m_ll2"
+    abs_y = "abs_y"
 
 
 class _KinematicsInformation:
@@ -92,7 +99,9 @@ class _Process:
             return True
         # We check if kin_cov is a subset of self.accepted_variables
         kin_cov = [v for v in kin_cov if not v.startswith("extra_")]
-        return set(self.accepted_variables).union(set(kin_cov)) == set(self.accepted_variables)
+        # sqrts is always accepted
+        accepted = set(self.accepted_variables).union({_Vars.sqrts})
+        return accepted.union(set(kin_cov)) == accepted
 
     def xq2map(self, kin_df, metadata):
         """Transform the kinematics dataframe into a x q dataframe
@@ -135,6 +144,9 @@ def _dis_xq2map(kin_info):
     x = kin_info.get_one_of("k1", _Vars.x)
     if "k2" in kin_info._kins:
         q2 = kin_info.get_one_of("k2") ** 2
+    elif _Vars.Q in kin_info._kins:
+        q = kin_info.get_one_of(_Vars.Q)
+        q2 = q**2
     else:
         q2 = kin_info.get_one_of(_Vars.Q2)
     return x, q2
@@ -164,6 +176,20 @@ def _shp_xq2map(kin_info):
     return np.clip(x, a_min=None, a_max=1, out=x), np.concatenate((q2, q2))
 
 
+def _pht_xq2map(kin_info):
+    # Then compute x, Q2
+    ET = kin_info.get_one_of(_Vars.ET)
+    ET2 = ET**2
+    eta = kin_info.get_one_of(_Vars.eta, _Vars.y)
+    sqrts = kin_info[_Vars.sqrts]
+
+    # eta = y for massless particles
+    x1 = ET / sqrts * np.exp(-eta)
+    x2 = ET / sqrts * np.exp(eta)
+    x = np.concatenate((x1, x2))
+    return np.clip(x, a_min=None, a_max=1, out=x), np.concatenate((ET2, ET2))
+
+
 def _dijets_xq2map(kin_info):
     # Here we can have either ystar or ydiff, but in either case we need to do the same
     ylab_1 = kin_info.get_one_of(_Vars.ystar, _Vars.ydiff, _Vars.eta_1, _Vars.abs_eta_1)
@@ -211,9 +237,11 @@ def _hqp_mqq_xq2map(kin_info):
 
 def _inc_xq2map(kin_info):
     # Compute x, Q2
-    # k2 necessary to take the mass for DY inclusive cross sections still not migrated
-    mass2 = kin_info.get_one_of(_Vars.m_W2, _Vars.m_Z2, _Vars.m_t2, "k2")
-
+    try:
+        mass = kin_info.get_one_of(_Vars.m_W, _Vars.m_Z, _Vars.m_t, _Vars.m_ll)
+        mass2 = mass**2
+    except KeyError:
+        mass2 = kin_info.get_one_of(_Vars.m_W2, _Vars.m_Z2, _Vars.m_t2, _Vars.m_ll)
     return np.sqrt(mass2) / kin_info[_Vars.sqrts], mass2
 
 
@@ -235,8 +263,8 @@ def _dyboson_xq2map(kin_info):
     Computes x and q2 mapping for pseudo rapidity observables
     originating from a W boson DY process.
     """
-    mass2 = kin_info.get_one_of(_Vars.m_W2, _Vars.m_Z2, _Vars.m_V2)
-    eta = kin_info.get_one_of(_Vars.eta, _Vars.y)
+    mass2 = kin_info.get_one_of(_Vars.m_W2, _Vars.m_Z2, _Vars.m_V2, _Vars.m_ll2, _Vars.M2)
+    eta = kin_info.get_one_of(_Vars.eta, _Vars.y, _Vars.abs_eta, _Vars.abs_y)
     sqrts = kin_info[_Vars.sqrts]
 
     # eta = y for massless particles
@@ -265,8 +293,8 @@ def _dybosonptrap_xq2map(kin_info):
     using the rapidity of the final lepton pair.
     """
     pT = kin_info[_Vars.pT]
-    eta = kin_info.get_one_of(_Vars.eta, _Vars.y)
-    m_ll2 = kin_info[_Vars.m_ll2]
+    eta = kin_info.get_one_of(_Vars.eta, _Vars.y, _Vars.abs_y)
+    m_ll2 = kin_info.get_one_of(_Vars.m_ll2, _Vars.m_Z2)
     sqrts = kin_info[_Vars.sqrts]
     ET2 = m_ll2 + pT * pT
     x1 = (np.sqrt(ET2) + pT) / sqrts * np.exp(-eta)
@@ -275,6 +303,34 @@ def _dybosonptrap_xq2map(kin_info):
     return x, np.concatenate((ET2, ET2))
 
 
+def _singletop_xq2map(kin_dict):
+
+    y_t = kin_dict[_Vars.y_t]
+    sqrts = kin_dict[_Vars.sqrts]
+    m_t2 = kin_dict[_Vars.m_t2]
+
+    q2 = m_t2
+    ratio = np.sqrt(q2) / sqrts
+    x1 = ratio * np.exp(y_t)
+    x2 = ratio * np.exp(-y_t)
+    x = np.concatenate((x1, x2))
+    return np.clip(x, a_min=None, a_max=1, out=x), np.concatenate((q2, q2))
+
+
+def _dymll_xq2map(kin_info):
+    """
+    Computes x and q2 mapping for DY Z -> 2 leptons mass.
+    x is approximated as x = sqrt(x1*x2) with m_ll^2 = x1*x2*s
+    """
+
+    m_ll = kin_info.get_one_of(_Vars.m_ll)
+    sqrts = kin_info.get_one_of(_Vars.sqrts)
+    m_ll2 = m_ll**2
+    x = m_ll / sqrts
+
+    return x, m_ll2
+
+
 DIS = _Process(
     "DIS",
     "Deep Inelastic Scattering",
@@ -341,7 +397,14 @@ HQP_YQ = _Process(
 HQP_PTQ = _Process(
     "HQP_PTQ",
     "Transverse momentum of top quark in top pair production",
-    accepted_variables=(_Vars.pT_t, _Vars.y_t, _Vars.y_ttBar, _Vars.sqrts, _Vars.m_t2),
+    accepted_variables=(
+        _Vars.pT_t,
+        _Vars.m_ttBar,
+        _Vars.y_t,
+        _Vars.y_ttBar,
+        _Vars.sqrts,
+        _Vars.m_t2,
+    ),
     xq2map_function=_hqp_ptq_xq2map,
 )
 
@@ -355,7 +418,7 @@ HQP_MQQ = _Process(
 INC = _Process(
     "INC",
     "Inclusive cross section",
-    accepted_variables=("zero", _Vars.sqrts, _Vars.m_W2, _Vars.m_Z2, _Vars.m_t2),
+    accepted_variables=("zero", _Vars.sqrts, _Vars.m_W2, _Vars.m_Z2, _Vars.m_t2, _Vars.m_ll),
     xq2map_function=_inc_xq2map,
 )
 
@@ -370,10 +433,29 @@ HERAJET = _Process(
 DY_2L = _Process(
     "DY_2L",
     "DY W or Z -> 2 leptons ",
-    accepted_variables=(_Vars.y, _Vars.eta, _Vars.m_W2, _Vars.m_Z2, _Vars.m_V2, _Vars.sqrts),
+    accepted_variables=(
+        _Vars.y,
+        _Vars.abs_y,
+        _Vars.eta,
+        _Vars.m_W2,
+        _Vars.m_Z2,
+        _Vars.m_V2,
+        _Vars.m_ll,
+        _Vars.m_ll2,
+        _Vars.sqrts,
+        _Vars.abs_eta,
+        _Vars.m_ll2,
+        _Vars.M2,
+    ),
     xq2map_function=_dyboson_xq2map,
 )
 
+DY_MLL = _Process(
+    "DY_MLL",
+    "DY Z -> ll mass of lepton pair",
+    accepted_variables=(_Vars.m_ll, _Vars.sqrts),
+    xq2map_function=_dymll_xq2map,
+)
 
 DY_PT = _Process(
     "DY_PT",
@@ -391,6 +473,7 @@ DY_PT_RAP = _Process(
         _Vars.m_Z2,
         _Vars.sqrts,
         _Vars.y,
+        _Vars.abs_y,
         _Vars.eta,
         _Vars.m_ll2,
     ),
@@ -404,12 +487,35 @@ POS_DIS = _Process(
     "POS_DIS", "Positivity of F2 structure functions", accepted_variables=(_Vars.x, _Vars.Q2)
 )
 
+POS_POL = _Process(
+    "POS_POL",
+    "Polarized DIS positivity boundary condition",
+    accepted_variables=(_Vars.x, _Vars.Q2, _Vars.y),
+)
+
+PHT = _Process(
+    "PHT",
+    "Photon production",
+    accepted_variables=(_Vars.eta, _Vars.ET, _Vars.sqrts),
+    xq2map_function=_pht_xq2map,
+)
+
+SINGLETOP = _Process(
+    "SINGLETOP",
+    "Single top production",
+    accepted_variables=(_Vars.m_t2, _Vars.sqrts, _Vars.y_t, _Vars.pT_t),
+    xq2map_function=_singletop_xq2map,
+)
+
+
 PROCESSES = {
     "DIS": DIS,
     "DIS_NC": dataclasses.replace(DIS, name="DIS_NC"),
     "DIS_CC": dataclasses.replace(DIS, name="DIS_CC"),
     "DIS_NCE": dataclasses.replace(DIS, name="DIS_NCE"),
     "DIS_POL": dataclasses.replace(DIS, name="DIS_POL"),
+    "DIS_NC_CHARM": dataclasses.replace(DIS, name="DIS_NC_CHARM"),
+    "DIS_NC_BOTTOM": dataclasses.replace(DIS, name="DIS_NC_BOTTOM"),
     "JET": JET,
     "DIJET": DIJET,
     "SHP_ASY": SHP,
@@ -422,7 +528,9 @@ PROCESSES = {
     "HERADIJET": dataclasses.replace(HERAJET, name="HERADIJET", description="DIS + jj production"),
     "JET_POL": JET_POL,
     "DIJET_POL": DIJET_POL,
-    "DY_Z_Y": dataclasses.replace(DY_2L, name="DY_Z_Y", description="DY Z -> ll rapidity"),
+    "DY_Z_Y": dataclasses.replace(DY_2L, name="DY_Z_Y", description="DY Z -> ll (pseudo)rapidity"),
+    "DYP_FT": dataclasses.replace(DY_2L, name="DYP_FT", description="DY Z -> ll (pseudo)rapidity"),
+    "DY_MLL": DY_MLL,
     "DY_W_ETA": dataclasses.replace(
         DY_2L, name="DY_W_ETA", description="DY W -> l nu pseudorapidity"
     ),
@@ -430,9 +538,15 @@ PROCESSES = {
         DY_2L, name="DY_VB_ETA", description="DY Z/W -> ll pseudorapidity"
     ),
     "DY_NC_PT": dataclasses.replace(DY_PT, name="DY_NC_PT", description="DY Z (ll) + j"),
+    "DY_CC_PT": dataclasses.replace(DY_PT, name="DY_CC_PT", description="DY W + j"),
     "DY_NC_PTRAP": dataclasses.replace(DY_PT_RAP, name="DY_NC_PTRAP", description="DY Z (ll) + j"),
     "POS_XPDF": POS_XPDF,
     "POS_DIS": POS_DIS,
+    "POS_POL": POS_POL,
+    "POS_DY": dataclasses.replace(POS_DIS, name="POS_DY"),
+    "INTEG": dataclasses.replace(POS_DIS, name="INTEG", description="Integrability dataset"),
+    "PHT": PHT,
+    "SINGLETOP": SINGLETOP,
 }
 
 

nnpdf_data/theory.py

@@ -4,12 +4,9 @@ that a theory card can contain.
 """
 
 import dataclasses
-from functools import lru_cache
 import logging
 from typing import Literal, Optional
 
-from .utils import parse_yaml_inp
-
 DEPRECATED_KEYS = ["MaxNfAs", "SxRes", "SxOrd" "EScaleVar", "Qedref", "global_nx"]
 
 log = logging.getLogger(__name__)
@@ -158,13 +155,3 @@ class TheoryCard:
 
     def asdict(self):
         return dataclasses.asdict(self)
-
-
-@lru_cache
-def parse_theory_card(theory_card):
-    """Read the theory card using validobj parsing
-    Returns the theory as a dictionary
-    """
-    if theory_card.exists():
-        return parse_yaml_inp(theory_card, TheoryCard)
-    raise TheoryNotFoundInDatabase(f"Theory card {theory_card} not found")

nnpdf_data/theory_cards/40005000.yaml

@@ -0,0 +1,48 @@
+ID: 40_005_000
+Comments: NNPDF4.0 N3LO alphas=0.118 and FHMRUVV splittings as Oct '24
+CKM:
+- 0.97428
+- 0.2253
+- 0.00347
+- 0.2252
+- 0.97345
+- 0.041
+- 0.00862
+- 0.0403
+- 0.999152
+DAMP: 0
+FNS: FONLL-E
+GF: 1.1663787e-05
+HQ: POLE
+IC: 1
+MP: 0.938
+MW: 80.398
+MZ: 91.1876
+MaxNfPdf: 5
+ModEv: TRN
+ModSV: null
+NfFF: 4
+PTO: 3
+Q0: 1.65
+QED: 0
+Qmb: 4.92
+Qmc: 1.51
+Qmt: 172.5
+Qref: 91.2
+SIN2TW: 0.23126
+TMC: 1
+XIF: 1.0
+XIR: 1.0
+alphaqed: 0.007496252
+alphas: 0.118
+kbThr: 1.0
+kcThr: 1.0
+ktThr: 1.0
+mb: 4.92
+mc: 1.51
+mt: 172.5
+nfref: 5
+nf0: 4
+n3lo_ad_variation: [0,0,0,0,0,0,0]
+n3lo_cf_variation: 0
+use_fhmruvv: true

nnpdf_data/theory_cards/40005001.yaml

@@ -0,0 +1,48 @@
+ID: 40_005_001
+Comments: NNPDF4.0 N3LO alphas=0.118 and FHMRUVV splittings as Oct '24, Pgg up variation
+CKM:
+- 0.97428
+- 0.2253
+- 0.00347
+- 0.2252
+- 0.97345
+- 0.041
+- 0.00862
+- 0.0403
+- 0.999152
+DAMP: 0
+FNS: FONLL-E
+GF: 1.1663787e-05
+HQ: POLE
+IC: 1
+MP: 0.938
+MW: 80.398
+MZ: 91.1876
+MaxNfPdf: 5
+ModEv: TRN
+ModSV: null
+NfFF: 4
+PTO: 3
+Q0: 1.65
+QED: 0
+Qmb: 4.92
+Qmc: 1.51
+Qmt: 172.5
+Qref: 91.2
+SIN2TW: 0.23126
+TMC: 1
+XIF: 1.0
+XIR: 1.0
+alphaqed: 0.007496252
+alphas: 0.118
+kbThr: 1.0
+kcThr: 1.0
+ktThr: 1.0
+mb: 4.92
+mc: 1.51
+mt: 172.5
+nfref: 5
+nf0: 4
+n3lo_ad_variation: [1,0,0,0,0,0,0]
+n3lo_cf_variation: 0
+use_fhmruvv: true

nnpdf_data/theory_cards/40005002.yaml

@@ -0,0 +1,48 @@
+ID: 40_005_002
+Comments: NNPDF4.0 N3LO alphas=0.118 and FHMRUVV splittings as Oct '24, Pgg down variation
+CKM:
+- 0.97428
+- 0.2253
+- 0.00347
+- 0.2252
+- 0.97345
+- 0.041
+- 0.00862
+- 0.0403
+- 0.999152
+DAMP: 0
+FNS: FONLL-E
+GF: 1.1663787e-05
+HQ: POLE
+IC: 1
+MP: 0.938
+MW: 80.398
+MZ: 91.1876
+MaxNfPdf: 5
+ModEv: TRN
+ModSV: null
+NfFF: 4
+PTO: 3
+Q0: 1.65
+QED: 0
+Qmb: 4.92
+Qmc: 1.51
+Qmt: 172.5
+Qref: 91.2
+SIN2TW: 0.23126
+TMC: 1
+XIF: 1.0
+XIR: 1.0
+alphaqed: 0.007496252
+alphas: 0.118
+kbThr: 1.0
+kcThr: 1.0
+ktThr: 1.0
+mb: 4.92
+mc: 1.51
+mt: 172.5
+nfref: 5
+nf0: 4
+n3lo_ad_variation: [2,0,0,0,0,0,0]
+n3lo_cf_variation: 0
+use_fhmruvv: true

nnpdf_data/theory_cards/40005003.yaml

@@ -0,0 +1,48 @@
+ID: 40_005_003
+Comments: NNPDF4.0 N3LO alphas=0.118 and FHMRUVV splittings as Oct '24, Pgq up variation
+CKM:
+- 0.97428
+- 0.2253
+- 0.00347
+- 0.2252
+- 0.97345
+- 0.041
+- 0.00862
+- 0.0403
+- 0.999152
+DAMP: 0
+FNS: FONLL-E
+GF: 1.1663787e-05
+HQ: POLE
+IC: 1
+MP: 0.938
+MW: 80.398
+MZ: 91.1876
+MaxNfPdf: 5
+ModEv: TRN
+ModSV: null
+NfFF: 4
+PTO: 3
+Q0: 1.65
+QED: 0
+Qmb: 4.92
+Qmc: 1.51
+Qmt: 172.5
+Qref: 91.2
+SIN2TW: 0.23126
+TMC: 1
+XIF: 1.0
+XIR: 1.0
+alphaqed: 0.007496252
+alphas: 0.118
+kbThr: 1.0
+kcThr: 1.0
+ktThr: 1.0
+mb: 4.92
+mc: 1.51
+mt: 172.5
+nfref: 5
+nf0: 4
+n3lo_ad_variation: [0,1,0,0,0,0,0]
+n3lo_cf_variation: 0
+use_fhmruvv: true

nnpdf_data/theory_cards/40005004.yaml

@@ -0,0 +1,48 @@
+ID: 40_005_004
+Comments: NNPDF4.0 N3LO alphas=0.118 and FHMRUVV splittings as Oct '24, Pgq down variation
+CKM:
+- 0.97428
+- 0.2253
+- 0.00347
+- 0.2252
+- 0.97345
+- 0.041
+- 0.00862
+- 0.0403
+- 0.999152
+DAMP: 0
+FNS: FONLL-E
+GF: 1.1663787e-05
+HQ: POLE
+IC: 1
+MP: 0.938
+MW: 80.398
+MZ: 91.1876
+MaxNfPdf: 5
+ModEv: TRN
+ModSV: null
+NfFF: 4
+PTO: 3
+Q0: 1.65
+QED: 0
+Qmb: 4.92
+Qmc: 1.51
+Qmt: 172.5
+Qref: 91.2
+SIN2TW: 0.23126
+TMC: 1
+XIF: 1.0
+XIR: 1.0
+alphaqed: 0.007496252
+alphas: 0.118
+kbThr: 1.0
+kcThr: 1.0
+ktThr: 1.0
+mb: 4.92
+mc: 1.51
+mt: 172.5
+nfref: 5
+nf0: 4
+n3lo_ad_variation: [0,2,0,0,0,0,0]
+n3lo_cf_variation: 0
+use_fhmruvv: true

nnpdf_data/theory_cards/40005005.yaml

@@ -0,0 +1,48 @@
+ID: 40_005_005
+Comments: NNPDF4.0 N3LO alphas=0.118 and FHMRUVV splittings as Oct '24, Pqg up variation
+CKM:
+- 0.97428
+- 0.2253
+- 0.00347
+- 0.2252
+- 0.97345
+- 0.041
+- 0.00862
+- 0.0403
+- 0.999152
+DAMP: 0
+FNS: FONLL-E
+GF: 1.1663787e-05
+HQ: POLE
+IC: 1
+MP: 0.938
+MW: 80.398
+MZ: 91.1876
+MaxNfPdf: 5
+ModEv: TRN
+ModSV: null
+NfFF: 4
+PTO: 3
+Q0: 1.65
+QED: 0
+Qmb: 4.92
+Qmc: 1.51
+Qmt: 172.5
+Qref: 91.2
+SIN2TW: 0.23126
+TMC: 1
+XIF: 1.0
+XIR: 1.0
+alphaqed: 0.007496252
+alphas: 0.118
+kbThr: 1.0
+kcThr: 1.0
+ktThr: 1.0
+mb: 4.92
+mc: 1.51
+mt: 172.5
+nfref: 5
+nf0: 4
+n3lo_ad_variation: [0,0,1,0,0,0,0]
+n3lo_cf_variation: 0
+use_fhmruvv: true

nnpdf_data/theory_cards/40005006.yaml

@@ -0,0 +1,48 @@
+ID: 40_005_006
+Comments: NNPDF4.0 N3LO alphas=0.118 and FHMRUVV splittings as Oct '24, Pqg down variation
+CKM:
+- 0.97428
+- 0.2253
+- 0.00347
+- 0.2252
+- 0.97345
+- 0.041
+- 0.00862
+- 0.0403
+- 0.999152
+DAMP: 0
+FNS: FONLL-E
+GF: 1.1663787e-05
+HQ: POLE
+IC: 1
+MP: 0.938
+MW: 80.398
+MZ: 91.1876
+MaxNfPdf: 5
+ModEv: TRN
+ModSV: null
+NfFF: 4
+PTO: 3
+Q0: 1.65
+QED: 0
+Qmb: 4.92
+Qmc: 1.51
+Qmt: 172.5
+Qref: 91.2
+SIN2TW: 0.23126
+TMC: 1
+XIF: 1.0
+XIR: 1.0
+alphaqed: 0.007496252
+alphas: 0.118
+kbThr: 1.0
+kcThr: 1.0
+ktThr: 1.0
+mb: 4.92
+mc: 1.51
+mt: 172.5
+nfref: 5
+nf0: 4
+n3lo_ad_variation: [0,0,2,0,0,0,0]
+n3lo_cf_variation: 0
+use_fhmruvv: true

nnpdf_data/theory_cards/40005007.yaml

@@ -0,0 +1,48 @@
+ID: 40_005_007
+Comments: NNPDF4.0 N3LO alphas=0.118 and FHMRUVV splittings as Oct '24, Pqq up variation
+CKM:
+- 0.97428
+- 0.2253
+- 0.00347
+- 0.2252
+- 0.97345
+- 0.041
+- 0.00862
+- 0.0403
+- 0.999152
+DAMP: 0
+FNS: FONLL-E
+GF: 1.1663787e-05
+HQ: POLE
+IC: 1
+MP: 0.938
+MW: 80.398
+MZ: 91.1876
+MaxNfPdf: 5
+ModEv: TRN
+ModSV: null
+NfFF: 4
+PTO: 3
+Q0: 1.65
+QED: 0
+Qmb: 4.92
+Qmc: 1.51
+Qmt: 172.5
+Qref: 91.2
+SIN2TW: 0.23126
+TMC: 1
+XIF: 1.0
+XIR: 1.0
+alphaqed: 0.007496252
+alphas: 0.118
+kbThr: 1.0
+kcThr: 1.0
+ktThr: 1.0
+mb: 4.92
+mc: 1.51
+mt: 172.5
+nfref: 5
+nf0: 4
+n3lo_ad_variation: [0,0,0,1,0,0,0]
+n3lo_cf_variation: 0
+use_fhmruvv: true