SIA datasets implementation for Pions, Kaons and Protons

nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-KA_91P2GEV_EE-KASUM/data.yaml

@@ -0,0 +1,30 @@
+data_central:
+- 12.4
+- 13.27
+- 15.33
+- 17.43
+- 18.33
+- 19.62
+- 20.02
+- 21.66
+- 25.84
+- 27.46
+- 27.63
+- 10.6
+- 9.53
+- 9.15
+- 8.41
+- 7.96
+- 7.26
+- 6.34
+- 5.63
+- 4.94
+- 4.39
+- 4.22
+- 3.63
+- 3.1
+- 2.245
+- 1.538
+- 0.841
+- 0.2936
+- 0.0596

nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-KA_91P2GEV_EE-KASUM/filter.py

@@ -0,0 +1,15 @@
+from filter_core import magic
+import yaml
+
+table = "rawdata/Table2.yaml"
+ndat = 29
+var_name = 'xp'
+
+data, kin, err = magic(table, ndat, var_name)
+
+with open('data.yaml', 'w') as f:
+    yaml.dump(data, f, sort_keys=False)
+with open('kinematics.yaml', 'w') as f:
+    yaml.dump(kin, f, sort_keys=False)
+with open('uncertainties.yaml', 'w') as f:
+    yaml.dump(err, f, sort_keys=False)

nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-KA_91P2GEV_EE-KASUM/filter_core.py

@@ -0,0 +1,63 @@
+import yaml
+
+from nnpdf_data.filter_utils.utils import percentage_to_absolute as pta
+from nnpdf_data.filter_utils.utils import prettify_float
+
+yaml.add_representer(float, prettify_float)
+
+
+def magic(table, ndat, var_name):
+    with open(table, 'r') as f:
+        input = yaml.safe_load(f)
+
+    data_central = []
+    kin = []
+    error = []
+
+    values = input['dependent_variables'][0]['values']
+
+    for i in range(ndat):
+        kin_min = input['independent_variables'][0]['values'][i]['low']
+        kin_max = input['independent_variables'][0]['values'][i]['high']
+        if 'value' in input['independent_variables'][0]['values'][i]:
+            kin_mid = input['independent_variables'][0]['values'][i]['value']
+        else:
+            kin_mid = (kin_min + kin_max) / 2
+
+        kin_value = {var_name: {'min': kin_min, 'mid': kin_mid, 'max': kin_max}}
+
+        data_central_value = values[i]['value']
+        error_value = {}
+        error_value['stat'] = values[i]['errors'][0]['symerror']
+        error_value['sys'] = values[i]['errors'][1]['symerror']
+        if i < 8:
+            error_value['sys_norm'] = pta(5, data_central_value)
+        else:
+            error_value['sys_norm'] = pta(3, data_central_value)
+
+        kin.append(kin_value)
+        data_central.append(data_central_value)
+        error.append(error_value)
+
+    error_definition = {}
+    error_definition['stat'] = {
+        'definition': 'statistical uncertainty',
+        'treatment': 'ADD',
+        'type': 'UNCORR',
+    }
+    error_definition['sys'] = {
+        'definition': 'systematic uncertainty',
+        'treatment': 'ADD',
+        'type': 'UNCORR',
+    }
+    error_definition['sys_norm'] = {
+        'definition': 'systematic uncertainty (normalization)',
+        'treatment': 'MULT',
+        'type': 'CORR',
+    }
+
+    data_central_yaml = {'data_central': data_central}
+    kin_yaml = {'bins': kin}
+    uncertainties_yaml = {'definitions': error_definition, 'bins': error}
+
+    return data_central_yaml, kin_yaml, uncertainties_yaml

nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-KA_91P2GEV_EE-KASUM/kinematics.yaml

@@ -0,0 +1,117 @@
+bins:
+- xp:
+    min: 0.0055
+    mid: 0.00574
+    max: 0.006
+- xp:
+    min: 0.006
+    mid: 0.00622
+    max: 0.0065
+- xp:
+    min: 0.0065
+    mid: 0.00673
+    max: 0.007
+- xp:
+    min: 0.007
+    mid: 0.00722
+    max: 0.0075
+- xp:
+    min: 0.0075
+    mid: 0.00773
+    max: 0.008
+- xp:
+    min: 0.008
+    mid: 0.00822
+    max: 0.0085
+- xp:
+    min: 0.0085
+    mid: 0.00871
+    max: 0.009
+- xp:
+    min: 0.009
+    mid: 0.00922
+    max: 0.0095
+- xp:
+    min: 0.013
+    mid: 0.0135
+    max: 0.014
+- xp:
+    min: 0.014
+    mid: 0.0147
+    max: 0.016
+- xp:
+    min: 0.016
+    mid: 0.0167
+    max: 0.018
+- xp:
+    min: 0.07
+    mid: 0.0719
+    max: 0.075
+- xp:
+    min: 0.075
+    mid: 0.0769
+    max: 0.08
+- xp:
+    min: 0.08
+    mid: 0.0819
+    max: 0.085
+- xp:
+    min: 0.085
+    mid: 0.087
+    max: 0.09
+- xp:
+    min: 0.09
+    mid: 0.0942
+    max: 0.1
+- xp:
+    min: 0.1
+    mid: 0.104
+    max: 0.11
+- xp:
+    min: 0.11
+    mid: 0.114
+    max: 0.12
+- xp:
+    min: 0.12
+    mid: 0.124
+    max: 0.13
+- xp:
+    min: 0.13
+    mid: 0.134
+    max: 0.14
+- xp:
+    min: 0.14
+    mid: 0.144
+    max: 0.15
+- xp:
+    min: 0.15
+    mid: 0.154
+    max: 0.16
+- xp:
+    min: 0.16
+    mid: 0.169
+    max: 0.18
+- xp:
+    min: 0.18
+    mid: 0.189
+    max: 0.2
+- xp:
+    min: 0.2
+    mid: 0.222
+    max: 0.25
+- xp:
+    min: 0.25
+    mid: 0.272
+    max: 0.3
+- xp:
+    min: 0.3
+    mid: 0.342
+    max: 0.4
+- xp:
+    min: 0.4
+    mid: 0.476
+    max: 0.6
+- xp:
+    min: 0.6
+    mid: 0.674
+    max: 0.8

nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-KA_91P2GEV_EE-KASUM/metadata.yaml

@@ -0,0 +1,41 @@
+setname: 'ALEPH_SIA-KA_91P2GEV_EE-KASUM'
+
+nnpdf_metadata:
+  nnpdf31_process: "SIA-KA"
+  experiment: "ALEPH"
+
+arXiv:
+  url: "N/A"
+iNSPIRE:
+  url: "https://inspirehep.net/literature/382179"
+hepdata:
+  url: "https://www.hepdata.net/record/ins382179"
+  version: 1
+
+version: 1
+version_comment: "Initial implementation"
+
+implemented_observables:
+
+- observable: {description: "normalized x-sec differential wrt $x_p$", label: '$\frac{1}{\sigma_{tot}}\frac{d\sigma^h}{dx_p}$', units: ''}
+  observable_name: XP
+  process_type: SIA_XP
+  ndata: 29
+  tables: [2]
+  kinematics:
+    variables:
+      xp: {description: "hadron 3-momentum scaled to beam 3-momentum", label: '$x_p$', units: ''}
+    file: kinematics.yaml
+  data_central: data.yaml
+  data_uncertainties:
+  - uncertainties.yaml
+  kinematic_coverage: [xp]
+  plotting:
+    dataset_label: 'ALEPH SIA KA- + KA+ : $\frac{1}{\sigma_{tot}}\frac{d\sigma^h}{dx_p}$'
+    x_scale: log
+    plot_x: xp
+    y_label: '$\frac{1}{\sigma_{tot}}\frac{d\sigma^h}{dx_p}$'
+  # theory:
+  #   FK_tables:
+  #   - - 
+  #   operation: 'null'

nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-KA_91P2GEV_EE-KASUM/rawdata

@@ -0,0 +1 @@
+../ALEPH_SIA-PI_91P2GEV_EE-PISUM/rawdata/

nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-KA_91P2GEV_EE-KASUM/uncertainties.yaml

@@ -0,0 +1,101 @@
+definitions:
+  stat:
+    definition: statistical uncertainty
+    treatment: ADD
+    type: UNCORR
+  sys:
+    definition: systematic uncertainty
+    treatment: ADD
+    type: UNCORR
+  sys_norm:
+    definition: systematic uncertainty (normalization)
+    treatment: MULT
+    type: CORR
+bins:
+- stat: 1.12
+  sys: 0.01
+  sys_norm: 0.62
+- stat: 0.91
+  sys: 0.01
+  sys_norm: 0.6635
+- stat: 0.9
+  sys: 0.01
+  sys_norm: 7.66500000e-01
+- stat: 0.92
+  sys: 0.02
+  sys_norm: 0.8715
+- stat: 0.88
+  sys: 0.02
+  sys_norm: 0.9165
+- stat: 0.9
+  sys: 0.02
+  sys_norm: 9.81000000e-01
+- stat: 0.86
+  sys: 0.05
+  sys_norm: 1.001
+- stat: 0.88
+  sys: 0.12
+  sys_norm: 1.083
+- stat: 0.66
+  sys: 0.5
+  sys_norm: 0.7752
+- stat: 0.47
+  sys: 0.68
+  sys_norm: 0.8238
+- stat: 0.53
+  sys: 2.2
+  sys_norm: 0.8289
+- stat: 0.3
+  sys: 1.28
+  sys_norm: 0.318
+- stat: 0.26
+  sys: 0.98
+  sys_norm: 0.2859
+- stat: 0.23
+  sys: 0.83
+  sys_norm: 0.2745
+- stat: 0.21
+  sys: 0.71
+  sys_norm: 0.2523
+- stat: 0.14
+  sys: 0.56
+  sys_norm: 2.38800000e-01
+- stat: 0.13
+  sys: 0.47
+  sys_norm: 2.17800000e-01
+- stat: 0.11
+  sys: 0.37
+  sys_norm: 0.1902
+- stat: 0.11
+  sys: 0.32
+  sys_norm: 1.68900000e-01
+- stat: 0.1
+  sys: 0.28
+  sys_norm: 0.1482
+- stat: 0.09
+  sys: 0.24
+  sys_norm: 1.31700000e-01
+- stat: 0.09
+  sys: 0.22
+  sys_norm: 0.1266
+- stat: 0.06
+  sys: 0.18
+  sys_norm: 1.08900000e-01
+- stat: 0.05
+  sys: 0.15
+  sys_norm: 9.30000000e-02
+- stat: 0.029
+  sys: 0.109
+  sys_norm: 6.73500000e-02
+- stat: 0.025
+  sys: 0.076
+  sys_norm: 0.04614
+- stat: 0.013
+  sys: 0.043
+  sys_norm: 0.02523
+- stat: 0.0053
+  sys: 0.0146
+  sys_norm: 0.008808
+- stat: 0.0022
+  sys: 0.0031
+  sys_norm: 0.001788