Skip to content

Commit d886cf4

Browse files
johannbrehmerjohannbrehmer
authored
Merge pull request #324 from johannbrehmer/develop
Better test of the ratio estimation
2 parents 08620cc + 2fbdbfb commit d886cf4

File tree

3 files changed

+338
-294
lines changed

3 files changed

+338
-294
lines changed

‎examples/tutorial_particle_physics/4b_fisher_information.ipynb‎

Lines changed: 34 additions & 56 deletions
Original file line numberDiff line numberDiff line change
@@ -79,14 +79,14 @@
7979
"name": "stderr",
8080
"output_type": "stream",
8181
"text": [
82-
"11:05 madminer.analysis INFO Loading data from data/lhe_data_shuffled.h5\n",
83-
"11:05 madminer.analysis INFO Found 2 parameters\n",
84-
"11:05 madminer.analysis INFO Did not find nuisance parameters\n",
85-
"11:05 madminer.analysis INFO Found 6 benchmarks, of which 6 physical\n",
86-
"11:05 madminer.analysis INFO Found 3 observables\n",
87-
"11:05 madminer.analysis INFO Found 15030 events\n",
88-
"11:05 madminer.analysis INFO Found morphing setup with 6 components\n",
89-
"11:05 madminer.analysis INFO Did not find nuisance morphing setup\n"
82+
"11:10 madminer.analysis INFO Loading data from data/lhe_data_shuffled.h5\n",
83+
"11:10 madminer.analysis INFO Found 2 parameters\n",
84+
"11:10 madminer.analysis INFO Did not find nuisance parameters\n",
85+
"11:10 madminer.analysis INFO Found 6 benchmarks, of which 6 physical\n",
86+
"11:10 madminer.analysis INFO Found 3 observables\n",
87+
"11:10 madminer.analysis INFO Found 15030 events\n",
88+
"11:10 madminer.analysis INFO Found morphing setup with 6 components\n",
89+
"11:10 madminer.analysis INFO Did not find nuisance morphing setup\n"
9090
]
9191
}
9292
],
@@ -118,16 +118,16 @@
118118
"name": "stderr",
119119
"output_type": "stream",
120120
"text": [
121-
"11:05 madminer.fisherinfor INFO Found 2 parameters in SALLY model, matching 2 physical parameters in MadMiner file\n",
122-
"11:05 madminer.fisherinfor INFO Evaluating rate Fisher information\n",
123-
"11:05 madminer.fisherinfor INFO Evaluating kinematic Fisher information on batch 1 / 1\n"
121+
"11:10 madminer.fisherinfor INFO Found 2 parameters in SALLY model, matching 2 physical parameters in MadMiner file\n",
122+
"11:10 madminer.fisherinfor INFO Evaluating rate Fisher information\n",
123+
"11:10 madminer.fisherinfor INFO Evaluating kinematic Fisher information on batch 1 / 1\n"
124124
]
125125
},
126126
{
127127
"name": "stdout",
128128
"output_type": "stream",
129129
"text": [
130-
"Kinematic Fisher information after 30 ifb:\n",
130+
"Fisher information after 30 ifb:\n",
131131
"[[9.64206944 0.68312743]\n",
132132
" [0.68312743 8.25388798]]\n"
133133
]
@@ -152,19 +152,19 @@
152152
},
153153
{
154154
"cell_type": "code",
155-
"execution_count": 13,
155+
"execution_count": 5,
156156
"metadata": {},
157157
"outputs": [
158158
{
159159
"name": "stderr",
160160
"output_type": "stream",
161161
"text": [
162-
"11:09 madminer.fisherinfor INFO Bins with largest statistical uncertainties on rates:\n",
163-
"11:09 madminer.fisherinfor INFO Bin 1: (0.0085 +/- 0.0031) fb (36.4223 %)\n",
164-
"11:09 madminer.fisherinfor INFO Bin 3: (0.3612 +/- 0.0471) fb (13.0365 %)\n",
165-
"11:09 madminer.fisherinfor INFO Bin 5: (0.0030 +/- 0.0003) fb (10.7680 %)\n",
166-
"11:09 madminer.fisherinfor INFO Bin 2: (0.6017 +/- 0.0508) fb (8.4371 %)\n",
167-
"11:09 madminer.fisherinfor INFO Bin 4: (0.0520 +/- 0.0028) fb (5.3161 %)\n"
162+
"11:10 madminer.fisherinfor INFO Bins with largest statistical uncertainties on rates:\n",
163+
"11:10 madminer.fisherinfor INFO Bin 1: (0.0085 +/- 0.0031) fb (36.4223 %)\n",
164+
"11:10 madminer.fisherinfor INFO Bin 3: (0.3612 +/- 0.0471) fb (13.0365 %)\n",
165+
"11:10 madminer.fisherinfor INFO Bin 5: (0.0030 +/- 0.0003) fb (10.7680 %)\n",
166+
"11:10 madminer.fisherinfor INFO Bin 2: (0.6017 +/- 0.0508) fb (8.4371 %)\n",
167+
"11:10 madminer.fisherinfor INFO Bin 4: (0.0520 +/- 0.0028) fb (5.3161 %)\n"
168168
]
169169
},
170170
{
@@ -198,49 +198,19 @@
198198
},
199199
{
200200
"cell_type": "code",
201-
"execution_count": 14,
201+
"execution_count": 6,
202202
"metadata": {},
203203
"outputs": [
204204
{
205205
"name": "stderr",
206206
"output_type": "stream",
207207
"text": [
208-
"11:09 madminer.fisherinfor INFO Bins with largest statistical uncertainties on rates:\n",
209-
"11:09 madminer.fisherinfor INFO Bin (1, 2): (0.0033 +/- 0.0026) fb (78.0906 %)\n",
210-
"11:09 madminer.fisherinfor INFO Bin (6, 5): (0.0005 +/- 0.0002) fb (50.9767 %)\n",
211-
"11:09 madminer.fisherinfor INFO Bin (4, 3): (0.0831 +/- 0.0380) fb (45.7958 %)\n",
212-
"11:09 madminer.fisherinfor INFO Bin (1, 1): (0.0041 +/- 0.0017) fb (40.8387 %)\n",
213-
"11:09 madminer.fisherinfor INFO Bin (1, 4): (0.0003 +/- 0.0001) fb (26.7261 %)\n",
214-
"11:09 madminer.fisherinfor INFO Bin (7, 2): (0.0004 +/- 0.0001) fb (23.9490 %)\n",
215-
"11:09 madminer.fisherinfor INFO Bin (4, 1): (0.0501 +/- 0.0105) fb (20.8863 %)\n",
216-
"11:09 madminer.fisherinfor INFO Bin (5, 5): (0.0089 +/- 0.0017) fb (19.6110 %)\n",
217-
"11:09 madminer.fisherinfor INFO Bin (3, 1): (0.0592 +/- 0.0115) fb (19.4720 %)\n",
218-
"11:09 madminer.fisherinfor INFO Bin (1, 3): (0.0007 +/- 0.0001) fb (16.6667 %)\n",
219-
"11:09 madminer.fisherinfor INFO Bin (3, 5): (0.1642 +/- 0.0243) fb (14.7882 %)\n",
220-
"11:09 madminer.fisherinfor INFO Bin (2, 3): (0.3278 +/- 0.0451) fb (13.7567 %)\n",
221-
"11:09 madminer.fisherinfor INFO Bin (7, 1): (0.0005 +/- 0.0001) fb (13.7479 %)\n",
222-
"11:09 madminer.fisherinfor INFO Bin (2, 5): (0.1260 +/- 0.0171) fb (13.6011 %)\n",
223-
"11:09 madminer.fisherinfor INFO Bin (4, 2): (0.0638 +/- 0.0084) fb (13.1747 %)\n",
224-
"11:09 madminer.fisherinfor INFO Bin (2, 4): (0.0887 +/- 0.0108) fb (12.2208 %)\n",
225-
"11:09 madminer.fisherinfor INFO Bin (5, 3): (0.0064 +/- 0.0007) fb (11.1804 %)\n",
226-
"11:09 madminer.fisherinfor INFO Bin (6, 4): (0.0016 +/- 0.0002) fb (10.5590 %)\n",
227-
"11:09 madminer.fisherinfor INFO Bin (5, 4): (0.0106 +/- 0.0008) fb (7.6726 %)\n",
228-
"11:09 madminer.fisherinfor INFO Bin (5, 2): (0.0261 +/- 0.0018) fb (7.0872 %)\n",
229-
"11:09 madminer.fisherinfor INFO Bin (7, 3): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
230-
"11:09 madminer.fisherinfor INFO Bin (6, 3): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
231-
"11:09 madminer.fisherinfor INFO Bin (1, 5): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
232-
"11:09 madminer.fisherinfor INFO Bin (2, 1): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
233-
"11:09 madminer.fisherinfor INFO Bin (2, 2): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
234-
"11:09 madminer.fisherinfor INFO Bin (6, 2): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
235-
"11:09 madminer.fisherinfor INFO Bin (6, 1): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
236-
"11:09 madminer.fisherinfor INFO Bin (5, 1): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
237-
"11:09 madminer.fisherinfor INFO Bin (3, 2): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
238-
"11:09 madminer.fisherinfor INFO Bin (3, 3): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
239-
"11:09 madminer.fisherinfor INFO Bin (3, 4): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
240-
"11:09 madminer.fisherinfor INFO Bin (7, 4): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
241-
"11:09 madminer.fisherinfor INFO Bin (4, 4): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
242-
"11:09 madminer.fisherinfor INFO Bin (4, 5): (0.0000 +/- 0.0000) fb (0.0000 %)\n",
243-
"11:09 madminer.fisherinfor INFO Bin (7, 5): (0.0000 +/- 0.0000) fb (0.0000 %)\n"
208+
"11:10 madminer.fisherinfor INFO Bins with largest statistical uncertainties on rates:\n",
209+
"11:10 madminer.fisherinfor INFO Bin (1, 2): (0.0033 +/- 0.0026) fb (78.0906 %)\n",
210+
"11:10 madminer.fisherinfor INFO Bin (6, 5): (0.0005 +/- 0.0002) fb (50.9767 %)\n",
211+
"11:10 madminer.fisherinfor INFO Bin (4, 3): (0.0831 +/- 0.0380) fb (45.7958 %)\n",
212+
"11:10 madminer.fisherinfor INFO Bin (1, 1): (0.0041 +/- 0.0017) fb (40.8387 %)\n",
213+
"11:10 madminer.fisherinfor INFO Bin (1, 4): (0.0003 +/- 0.0001) fb (26.7261 %)\n"
244214
]
245215
},
246216
{
@@ -251,6 +221,14 @@
251221
"[[7.20570573 0.33434335]\n",
252222
" [0.33434335 6.21207077]]\n"
253223
]
224+
},
225+
{
226+
"name": "stderr",
227+
"output_type": "stream",
228+
"text": [
229+
"/Users/johannbrehmer/work/projects/madminer/madminer/madminer/fisherinformation.py:1063: RuntimeWarning: divide by zero encountered in true_divide\n",
230+
" inv_sigma = sanitize_array(1.0 / sigma) # Shape (n_events,)\n"
231+
]
254232
}
255233
],
256234
"source": [

0 commit comments

Comments
 (0)