Skip to content

Commit 39c92f3

Browse files
t7phyt7phy
committed
aleph ka impl
1 parent 96677f4 commit 39c92f3

6 files changed

Lines changed: 367 additions & 0 deletions

File tree

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

Lines changed: 30 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,30 @@
1+
data_central:
2+
- 12.4
3+
- 13.27
4+
- 15.33
5+
- 17.43
6+
- 18.33
7+
- 19.62
8+
- 20.02
9+
- 21.66
10+
- 25.84
11+
- 27.46
12+
- 27.63
13+
- 10.6
14+
- 9.53
15+
- 9.15
16+
- 8.41
17+
- 7.96
18+
- 7.26
19+
- 6.34
20+
- 5.63
21+
- 4.94
22+
- 4.39
23+
- 4.22
24+
- 3.63
25+
- 3.1
26+
- 2.245
27+
- 1.538
28+
- 0.841
29+
- 0.2936
30+
- 0.0596

‎nnpdf_data/nnpdf_data/commondata/ALEPH_SIA-KA_91P2GEV_EE-KASUM/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/Table2.yaml"
5+
ndat = 29
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-KA_91P2GEV_EE-KASUM/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 < 8:
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-KA_91P2GEV_EE-KASUM/kinematics.yaml‎

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

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

Lines changed: 41 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,41 @@
1+
setname: 'ALEPH SIA-KA_91P2GEV_EE-KASUM'
2+
3+
nnpdf_metadata:
4+
nnpdf31_process: "SIA-KA"
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-KA_91P2GEV_EE-KASUM/uncertainties.yaml‎

Lines changed: 101 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,101 @@
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: 1.12
16+
sys_1: 0.01
17+
sys_2: 0.62
18+
- stat: 0.91
19+
sys_1: 0.01
20+
sys_2: 0.6635
21+
- stat: 0.9
22+
sys_1: 0.01
23+
sys_2: 7.66500000e-01
24+
- stat: 0.92
25+
sys_1: 0.02
26+
sys_2: 0.8715
27+
- stat: 0.88
28+
sys_1: 0.02
29+
sys_2: 0.9165
30+
- stat: 0.9
31+
sys_1: 0.02
32+
sys_2: 9.81000000e-01
33+
- stat: 0.86
34+
sys_1: 0.05
35+
sys_2: 1.001
36+
- stat: 0.88
37+
sys_1: 0.12
38+
sys_2: 1.083
39+
- stat: 0.66
40+
sys_1: 0.5
41+
sys_2: 0.7752
42+
- stat: 0.47
43+
sys_1: 0.68
44+
sys_2: 0.8238
45+
- stat: 0.53
46+
sys_1: 2.2
47+
sys_2: 0.8289
48+
- stat: 0.3
49+
sys_1: 1.28
50+
sys_2: 0.318
51+
- stat: 0.26
52+
sys_1: 0.98
53+
sys_2: 0.2859
54+
- stat: 0.23
55+
sys_1: 0.83
56+
sys_2: 0.2745
57+
- stat: 0.21
58+
sys_1: 0.71
59+
sys_2: 0.2523
60+
- stat: 0.14
61+
sys_1: 0.56
62+
sys_2: 2.38800000e-01
63+
- stat: 0.13
64+
sys_1: 0.47
65+
sys_2: 2.17800000e-01
66+
- stat: 0.11
67+
sys_1: 0.37
68+
sys_2: 0.1902
69+
- stat: 0.11
70+
sys_1: 0.32
71+
sys_2: 1.68900000e-01
72+
- stat: 0.1
73+
sys_1: 0.28
74+
sys_2: 0.1482
75+
- stat: 0.09
76+
sys_1: 0.24
77+
sys_2: 1.31700000e-01
78+
- stat: 0.09
79+
sys_1: 0.22
80+
sys_2: 0.1266
81+
- stat: 0.06
82+
sys_1: 0.18
83+
sys_2: 1.08900000e-01
84+
- stat: 0.05
85+
sys_1: 0.15
86+
sys_2: 9.30000000e-02
87+
- stat: 0.029
88+
sys_1: 0.109
89+
sys_2: 6.73500000e-02
90+
- stat: 0.025
91+
sys_1: 0.076
92+
sys_2: 0.04614
93+
- stat: 0.013
94+
sys_1: 0.043
95+
sys_2: 0.02523
96+
- stat: 0.0053
97+
sys_1: 0.0146
98+
sys_2: 0.008808
99+
- stat: 0.0022
100+
sys_1: 0.0031
101+
sys_2: 0.001788

0 commit comments

Comments
 (0)