Skip to content

Commit 61488fa

Browse files
t7phyt7phy
committed
aleph pro impl
1 parent 39c92f3 commit 61488fa

6 files changed

Lines changed: 343 additions & 0 deletions

File tree

‎nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-PRO_91P2GEV_EE-PROSUM/data.yaml‎

Lines changed: 27 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,27 @@
1+
data_central:
2+
- 8.32
3+
- 8.95
4+
- 9.8
5+
- 10.3
6+
- 10.7
7+
- 11.58
8+
- 12.37
9+
- 12.46
10+
- 5.315
11+
- 5.008
12+
- 4.445
13+
- 4.555
14+
- 3.742
15+
- 3.355
16+
- 2.905
17+
- 2.653
18+
- 2.371
19+
- 2.137
20+
- 1.878
21+
- 1.696
22+
- 1.299
23+
- 0.966
24+
- 0.614
25+
- 0.305
26+
- 0.0784
27+
- 0.0054

‎nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-PRO_91P2GEV_EE-PROSUM/filter.py‎

Lines changed: 15 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,15 @@
1+
from filter_core import magic
2+
import yaml
3+
4+
table = "rawdata/Table3.yaml"
5+
ndat = 26
6+
var_name = 'xp'
7+
8+
data, kin, err = magic(table, ndat, var_name)
9+
10+
with open('data.yaml', 'w') as f:
11+
yaml.dump(data, f, sort_keys=False)
12+
with open('kinematics.yaml', 'w') as f:
13+
yaml.dump(kin, f, sort_keys=False)
14+
with open('uncertainties.yaml', 'w') as f:
15+
yaml.dump(err, f, sort_keys=False)

‎nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-PRO_91P2GEV_EE-PROSUM/filter_core.py‎

Lines changed: 63 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,63 @@
1+
import yaml
2+
3+
from nnpdf_data.filter_utils.utils import percentage_to_absolute as pta
4+
from nnpdf_data.filter_utils.utils import prettify_float
5+
6+
yaml.add_representer(float, prettify_float)
7+
8+
9+
def magic(table, ndat, var_name):
10+
with open(table, 'r') as f:
11+
input = yaml.safe_load(f)
12+
13+
data_central = []
14+
kin = []
15+
error = []
16+
17+
values = input['dependent_variables'][0]['values']
18+
19+
for i in range(ndat):
20+
kin_min = input['independent_variables'][0]['values'][i]['low']
21+
kin_max = input['independent_variables'][0]['values'][i]['high']
22+
if 'value' in input['independent_variables'][0]['values'][i]:
23+
kin_mid = input['independent_variables'][0]['values'][i]['value']
24+
else:
25+
kin_mid = (kin_min + kin_max) / 2
26+
27+
kin_value = {var_name: {'min': kin_min, 'mid': kin_mid, 'max': kin_max}}
28+
29+
data_central_value = values[i]['value']
30+
error_value = {}
31+
error_value['stat'] = values[i]['errors'][0]['symerror']
32+
error_value['sys_1'] = values[i]['errors'][1]['symerror']
33+
if i < 6:
34+
error_value['sys_2'] = pta(5, data_central_value)
35+
else:
36+
error_value['sys_2'] = pta(3, data_central_value)
37+
38+
kin.append(kin_value)
39+
data_central.append(data_central_value)
40+
error.append(error_value)
41+
42+
error_definition = {}
43+
error_definition['stat'] = {
44+
'definition': 'statistical uncertainty',
45+
'treatment': 'ADD',
46+
'type': 'UNCORR',
47+
}
48+
error_definition['sys_1'] = {
49+
'definition': 'systematic uncertainty 1',
50+
'treatment': 'MULT',
51+
'type': 'CORR',
52+
}
53+
error_definition['sys_2'] = {
54+
'definition': 'systematic uncertainty 2',
55+
'treatment': 'MULT',
56+
'type': 'CORR',
57+
}
58+
59+
data_central_yaml = {'data_central': data_central}
60+
kin_yaml = {'bins': kin}
61+
uncertainties_yaml = {'definitions': error_definition, 'bins': error}
62+
63+
return data_central_yaml, kin_yaml, uncertainties_yaml

‎nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-PRO_91P2GEV_EE-PROSUM/kinematics.yaml‎

Lines changed: 105 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,105 @@
1+
bins:
2+
- xp:
3+
min: 0.01
4+
mid: 1.05000000e-02
5+
max: 0.011
6+
- xp:
7+
min: 0.011
8+
mid: 0.0115
9+
max: 0.012
10+
- xp:
11+
min: 0.012
12+
mid: 0.0125
13+
max: 0.013
14+
- xp:
15+
min: 0.013
16+
mid: 0.0135
17+
max: 0.014
18+
- xp:
19+
min: 0.014
20+
mid: 0.0147
21+
max: 0.016
22+
- xp:
23+
min: 0.016
24+
mid: 0.0167
25+
max: 0.018
26+
- xp:
27+
min: 0.024
28+
mid: 0.0247
29+
max: 0.026
30+
- xp:
31+
min: 0.026
32+
mid: 0.0268
33+
max: 0.028
34+
- xp:
35+
min: 0.07
36+
mid: 0.0719
37+
max: 0.075
38+
- xp:
39+
min: 0.075
40+
mid: 0.0769
41+
max: 0.08
42+
- xp:
43+
min: 0.08
44+
mid: 0.0819
45+
max: 0.085
46+
- xp:
47+
min: 0.085
48+
mid: 0.087
49+
max: 0.09
50+
- xp:
51+
min: 0.09
52+
mid: 0.0942
53+
max: 0.1
54+
- xp:
55+
min: 0.1
56+
mid: 0.104
57+
max: 0.11
58+
- xp:
59+
min: 0.11
60+
mid: 0.114
61+
max: 0.12
62+
- xp:
63+
min: 0.12
64+
mid: 0.124
65+
max: 0.13
66+
- xp:
67+
min: 0.13
68+
mid: 0.134
69+
max: 0.14
70+
- xp:
71+
min: 0.14
72+
mid: 0.144
73+
max: 0.15
74+
- xp:
75+
min: 0.15
76+
mid: 0.154
77+
max: 0.16
78+
- xp:
79+
min: 0.16
80+
mid: 0.169
81+
max: 0.18
82+
- xp:
83+
min: 0.18
84+
mid: 0.189
85+
max: 0.2
86+
- xp:
87+
min: 0.2
88+
mid: 0.222
89+
max: 0.25
90+
- xp:
91+
min: 0.25
92+
mid: 0.272
93+
max: 0.3
94+
- xp:
95+
min: 0.3
96+
mid: 0.342
97+
max: 0.4
98+
- xp:
99+
min: 0.4
100+
mid: 0.476
101+
max: 0.6
102+
- xp:
103+
min: 0.6
104+
mid: 0.674
105+
max: 0.8

‎nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-PRO_91P2GEV_EE-PROSUM/metadata.yaml‎

Lines changed: 41 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,41 @@
1+
setname: 'ALEPH SIA-PRO_91P2GEV_EE-PROSUM'
2+
3+
nnpdf_metadata:
4+
nnpdf31_process: "SIA-PRO"
5+
experiment: "ALEPH"
6+
7+
arXiv:
8+
url: "N/A"
9+
iNSPIRE:
10+
url: "https://inspirehep.net/literature/382179"
11+
hepdata:
12+
url: "https://www.hepdata.net/record/ins382179"
13+
version: 1
14+
15+
version: 1
16+
version_comment: "Initial implementation"
17+
18+
implemented_observables:
19+
20+
- observable: {description: "normalized x-sec differential wrt $x_p$", label: '$\frac{1}{\sigma_{tot}}\frac{d\sigma^h}{dx_p}$', units: ''}
21+
observable_name: XP
22+
process_type: SIA_XP
23+
ndata: 39
24+
tables: [1]
25+
kinematics:
26+
variables:
27+
xp: {description: "hadron 3-momentum scaled to beam 3-momentum", label: '$x_p$', units: ''}
28+
file: kinematics.yaml
29+
data_central: data.yaml
30+
data_uncertainties:
31+
- uncertainties.yaml
32+
kinematic_coverage: [xp]
33+
plotting:
34+
dataset_label: 'ALEPH SIA PI- + PI+ : $\frac{1}{\sigma_{tot}}\frac{d\sigma^h}{dx_p}$'
35+
x_scale: log
36+
plot_x: xp
37+
y_label: '$\frac{1}{\sigma_{tot}}\frac{d\sigma^h}{dx_p}$'
38+
# theory:
39+
# FK_tables:
40+
# - -
41+
# operation: 'null'

‎nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-PRO_91P2GEV_EE-PROSUM/uncertainties.yaml‎

Lines changed: 92 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,92 @@
1+
definitions:
2+
stat:
3+
definition: statistical uncertainty
4+
treatment: ADD
5+
type: UNCORR
6+
sys_1:
7+
definition: systematic uncertainty 1
8+
treatment: MULT
9+
type: CORR
10+
sys_2:
11+
definition: systematic uncertainty 2
12+
treatment: MULT
13+
type: CORR
14+
bins:
15+
- stat: 0.35
16+
sys_1: 0.0
17+
sys_2: 4.16000000e-01
18+
- stat: 0.36
19+
sys_1: 0.0
20+
sys_2: 0.4475
21+
- stat: 0.36
22+
sys_1: 0.01
23+
sys_2: 0.49
24+
- stat: 0.38
25+
sys_1: 0.01
26+
sys_2: 0.515
27+
- stat: 0.26
28+
sys_1: 0.01
29+
sys_2: 0.535
30+
- stat: 0.27
31+
sys_1: 0.04
32+
sys_2: 0.579
33+
- stat: 0.18
34+
sys_1: 0.23
35+
sys_2: 0.3711
36+
- stat: 0.18
37+
sys_1: 0.44
38+
sys_2: 0.3738
39+
- stat: 0.216
40+
sys_1: 0.876
41+
sys_2: 0.15945
42+
- stat: 0.183
43+
sys_1: 0.639
44+
sys_2: 0.15024
45+
- stat: 0.162
46+
sys_1: 0.549
47+
sys_2: 1.33350000e-01
48+
- stat: 0.154
49+
sys_1: 0.474
50+
sys_2: 0.13665
51+
- stat: 0.092
52+
sys_1: 0.355
53+
sys_2: 0.11226
54+
- stat: 0.084
55+
sys_1: 0.292
56+
sys_2: 0.10065
57+
- stat: 0.077
58+
sys_1: 0.232
59+
sys_2: 0.08715
60+
- stat: 0.072
61+
sys_1: 0.205
62+
sys_2: 0.07959
63+
- stat: 0.068
64+
sys_1: 0.178
65+
sys_2: 0.07113
66+
- stat: 0.064
67+
sys_1: 0.162
68+
sys_2: 0.06411
69+
- stat: 0.061
70+
sys_1: 0.146
71+
sys_2: 5.63400000e-02
72+
- stat: 0.041
73+
sys_1: 0.118
74+
sys_2: 0.05088
75+
- stat: 0.036
76+
sys_1: 0.099
77+
sys_2: 0.03897
78+
- stat: 0.02
79+
sys_1: 0.073
80+
sys_2: 0.02898
81+
- stat: 0.017
82+
sys_1: 0.054
83+
sys_2: 1.84200000e-02
84+
- stat: 0.009
85+
sys_1: 0.031
86+
sys_2: 0.00915
87+
- stat: 0.0034
88+
sys_1: 0.011
89+
sys_2: 0.002352
90+
- stat: 0.0011
91+
sys_1: 0.0022
92+
sys_2: 0.000162

0 commit comments

Comments
 (0)