Skip to content

Commit

Permalink
Init template for pDIS generators
Browse files Browse the repository at this point in the history
  • Loading branch information
Radonirinaunimi committed Oct 12, 2023
1 parent d6366b5 commit 621d32f
Show file tree
Hide file tree
Showing 5 changed files with 338 additions and 2 deletions.
289 changes: 289 additions & 0 deletions extras/data/E155/e155_g1f1.yaml
Original file line number Diff line number Diff line change
@@ -0,0 +1,289 @@
bins:
- x:
min: null
mid: 0.015
max: null
Q2:
min: null
mid: 1.22
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.5
max: null
Q2:
min: null
mid: 5.77
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.5
max: null
Q2:
min: null
mid: 14.02
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.5
max: null
Q2:
min: null
mid: 26.86
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.75
max: null
Q2:
min: null
mid: 15.7
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.75
max: null
Q2:
min: null
mid: 34.72
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.025
max: null
Q2:
min: null
mid: 1.59
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.035
max: null
Q2:
min: null
mid: 2.05
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.05
max: null
Q2:
min: null
mid: 2.58
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.05
max: null
Q2:
min: null
mid: 4.01
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.08
max: null
Q2:
min: null
mid: 3.24
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.08
max: null
Q2:
min: null
mid: 5.36
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.125
max: null
Q2:
min: null
mid: 4.03
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.125
max: null
Q2:
min: null
mid: 7.17
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.125
max: null
Q2:
min: null
mid: 10.99
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.175
max: null
Q2:
min: null
mid: 4.62
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.175
max: null
Q2:
min: null
mid: 8.9
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.175
max: null
Q2:
min: null
mid: 13.19
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.25
max: null
Q2:
min: null
mid: 5.06
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.25
max: null
Q2:
min: null
mid: 10.64
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.25
max: null
Q2:
min: null
mid: 17.21
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.35
max: null
Q2:
min: null
mid: 5.51
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.35
max: null
Q2:
min: null
mid: 12.6
max: null
y:
min: null
mid: 0.0
max: null
- x:
min: null
mid: 0.35
max: null
Q2:
min: null
mid: 22.73
max: null
y:
min: null
mid: 0.0
max: null
5 changes: 5 additions & 0 deletions extras/data/E155/metadata.yaml
Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
E155_NC_9GEV_EP_G1F1RATIO:
description: Measurements with Polarised electron-proton beams
nnpdf_id: E155
y_label: A1
y_label_tex: $g_1 / F_1$
38 changes: 38 additions & 0 deletions extras/data/machinery/generate/E155.py
Original file line number Diff line number Diff line change
@@ -0,0 +1,38 @@
from .utils import load, obs_template


def dump(src_path, _target):
"""Generate the input card for E155 measurements.
Parameters
----------
src_path : str
input path
Returns
-------
dict
observables dictionary, corresponding to the runcard
"""
obs = obs_template.copy()
data = load(src_path, 0, ["x", "Q2"])
dict_kins = [
dict(x=d["x"]["mid"], y=d["y"]["mid"], Q2=d["Q2"]["mid"])
for d in data
]

# Details regarding the observables
obs["prDIS"] = "NC"
obs["ProjectileDIS"] = "electron"
obs["PolarizationDIS"] = 1.0
obs["TargetDIS"] = "proton"
obs["observables"] = {"g1_total": dict_kins}

return obs


# renaming
new_names = {
"e155_g1f1": "E155_NC_9GEV_EP_G1F1RATIO",
}
4 changes: 2 additions & 2 deletions extras/data/machinery/generate/__init__.py
Original file line number Diff line number Diff line change
@@ -1,8 +1,8 @@
from . import BCDMS, CHORUS, HERA, NMC, NUTEV, POS, SLAC, utils
from . import BCDMS, CHORUS, E155, HERA, NMC, NUTEV, POS, SLAC, utils

exps = {
getattr(m, "__name__").rsplit(".", maxsplit=1)[-1]: m
for m in [BCDMS, CHORUS, HERA, NMC, NUTEV, POS, SLAC]
for m in [BCDMS, CHORUS, E155, HERA, NMC, NUTEV, POS, SLAC]
}


Expand Down
4 changes: 4 additions & 0 deletions extras/data/machinery/generate/utils.py
Original file line number Diff line number Diff line change
Expand Up @@ -28,6 +28,10 @@ def load(path, skiprows, fields):
list(dict)
list of datapoints
"""
if path.suffix == ".yaml":
infile = yaml.safe_load(path.read_text())
return infile["bins"]

data = np.loadtxt(path, skiprows=skiprows)
data = data[:, : len(fields)]
return [dict(zip(fields, d.tolist())) for d in data]
Expand Down

0 comments on commit 621d32f

Please sign in to comment.