analysis/test_sysvar.py
=== stdout of analysis/tests/test_sysvar.py ===
[INFO] Steering file: /data/b2soft/buildbot/light/build/analysis/tests/test_sysvar.py
{(11, 11): p_min p_max ... variable threshold
0 0.5 1.0 ... electronID 0.8
1 0.5 1.0 ... electronID 0.8
2 0.5 1.0 ... electronID 0.8
3 0.5 1.0 ... electronID 0.8
4 0.5 1.0 ... electronID 0.8
.. ... ... ... ... ...
79 5.0 6.0 ... electronID 0.8
80 5.0 6.0 ... electronID 0.8
81 5.0 6.0 ... electronID 0.8
82 5.0 6.0 ... electronID 0.8
83 5.0 6.0 ... electronID 0.8
[84 rows x 13 columns], (11, 211): p_min p_max ... variable threshold
0 0.5 1.0 ... electronID 0.8
1 0.5 1.0 ... electronID 0.8
2 0.5 1.0 ... electronID 0.8
3 0.5 1.0 ... electronID 0.8
4 0.5 1.0 ... electronID 0.8
.. ... ... ... ... ...
65 4.0 5.0 ... electronID 0.8
66 4.0 5.0 ... electronID 0.8
67 4.0 5.0 ... electronID 0.8
68 4.0 5.0 ... electronID 0.8
69 4.0 5.0 ... electronID 0.8
[70 rows x 13 columns]}
cosTheta p charge electronID ... mcPDG B0_dec_mode B0_PDG B0_sigProb
0 0.250 1.5 -1 0.850 ... 11 0 511 0.02
1 0.550 2.5 1 0.900 ... 211 1 -511 0.03
2 0.850 3.5 1 0.950 ... -11 2 511 0.02
3 0.925 4.5 1 0.975 ... 211 3 -511 0.90
[4 rows x 9 columns]
{(11, 11): p_min p_max ... variable threshold
0 0.5 1.0 ... electronID 0.8
1 0.5 1.0 ... electronID 0.8
2 0.5 1.0 ... electronID 0.8
3 0.5 1.0 ... electronID 0.8
4 0.5 1.0 ... electronID 0.8
.. ... ... ... ... ...
79 5.0 6.0 ... electronID 0.8
80 5.0 6.0 ... electronID 0.8
81 5.0 6.0 ... electronID 0.8
82 5.0 6.0 ... electronID 0.8
83 5.0 6.0 ... electronID 0.8
[84 rows x 13 columns], (11, 211): p_min p_max ... variable threshold
0 0.5 1.0 ... electronID 0.8
1 0.5 1.0 ... electronID 0.8
2 0.5 1.0 ... electronID 0.8
3 0.5 1.0 ... electronID 0.8
4 0.5 1.0 ... electronID 0.8
.. ... ... ... ... ...
65 4.0 5.0 ... electronID 0.8
66 4.0 5.0 ... electronID 0.8
67 4.0 5.0 ... electronID 0.8
68 4.0 5.0 ... electronID 0.8
69 4.0 5.0 ... electronID 0.8
[70 rows x 13 columns]}
cosTheta p charge electronID ... mcPDG B0_dec_mode B0_PDG B0_sigProb
0 0.250 1.5 -1 0.850 ... 11 0 511 0.02
1 0.550 2.5 1 0.900 ... 211 1 -511 0.03
2 0.850 3.5 1 0.950 ... -11 2 511 0.02
3 0.925 4.5 1 0.975 ... 211 3 -511 0.90
[4 rows x 9 columns]
{(11, 11): p_min p_max ... variable threshold
0 0.5 1.0 ... electronID 0.8
1 0.5 1.0 ... electronID 0.8
2 0.5 1.0 ... electronID 0.8
3 0.5 1.0 ... electronID 0.8
4 0.5 1.0 ... electronID 0.8
.. ... ... ... ... ...
79 5.0 6.0 ... electronID 0.8
80 5.0 6.0 ... electronID 0.8
81 5.0 6.0 ... electronID 0.8
82 5.0 6.0 ... electronID 0.8
83 5.0 6.0 ... electronID 0.8
[84 rows x 13 columns], (11, 211): p_min p_max ... variable threshold
0 0.5 1.0 ... electronID 0.8
1 0.5 1.0 ... electronID 0.8
2 0.5 1.0 ... electronID 0.8
3 0.5 1.0 ... electronID 0.8
4 0.5 1.0 ... electronID 0.8
.. ... ... ... ... ...
65 4.0 5.0 ... electronID 0.8
66 4.0 5.0 ... electronID 0.8
67 4.0 5.0 ... electronID 0.8
68 4.0 5.0 ... electronID 0.8
69 4.0 5.0 ... electronID 0.8
[70 rows x 13 columns]}
cosTheta p charge electronID ... mcPDG B0_dec_mode B0_PDG B0_sigProb
0 0.250 1.5 -1 0.850 ... 11 0 511 0.02
1 0.550 2.5 1 0.900 ... 211 1 -511 0.03
2 0.850 3.5 1 0.950 ... -11 2 511 0.02
3 0.925 4.5 1 0.975 ... 211 3 -511 0.90
[4 rows x 9 columns]
{(11, 11): p_min p_max ... variable threshold
0 0.5 1.0 ... electronID 0.8
1 0.5 1.0 ... electronID 0.8
2 0.5 1.0 ... electronID 0.8
3 0.5 1.0 ... electronID 0.8
4 0.5 1.0 ... electronID 0.8
.. ... ... ... ... ...
79 5.0 6.0 ... electronID 0.8
80 5.0 6.0 ... electronID 0.8
81 5.0 6.0 ... electronID 0.8
82 5.0 6.0 ... electronID 0.8
83 5.0 6.0 ... electronID 0.8
[84 rows x 13 columns], (11, 211): p_min p_max ... variable threshold
0 0.5 1.0 ... electronID 0.8
1 0.5 1.0 ... electronID 0.8
2 0.5 1.0 ... electronID 0.8
3 0.5 1.0 ... electronID 0.8
4 0.5 1.0 ... electronID 0.8
.. ... ... ... ... ...
65 4.0 5.0 ... electronID 0.8
66 4.0 5.0 ... electronID 0.8
67 4.0 5.0 ... electronID 0.8
68 4.0 5.0 ... electronID 0.8
69 4.0 5.0 ... electronID 0.8
[70 rows x 13 columns]}
cosTheta p charge electronID ... mcPDG B0_dec_mode B0_PDG B0_sigProb
0 0.250 1.5 -1 0.850 ... 11 0 511 0.02
1 0.550 2.5 1 0.900 ... 211 1 -511 0.03
2 0.850 3.5 1 0.950 ... -11 2 511 0.02
3 0.925 4.5 1 0.975 ... 211 3 -511 0.90
[4 rows x 9 columns]
{(11, 11): p_min p_max ... variable threshold
0 0.5 1.0 ... electronID 0.8
1 0.5 1.0 ... electronID 0.8
2 0.5 1.0 ... electronID 0.8
3 0.5 1.0 ... electronID 0.8
4 0.5 1.0 ... electronID 0.8
.. ... ... ... ... ...
79 5.0 6.0 ... electronID 0.8
80 5.0 6.0 ... electronID 0.8
81 5.0 6.0 ... electronID 0.8
82 5.0 6.0 ... electronID 0.8
83 5.0 6.0 ... electronID 0.8
[84 rows x 13 columns], (11, 211): p_min p_max ... variable threshold
0 0.5 1.0 ... electronID 0.8
1 0.5 1.0 ... electronID 0.8
2 0.5 1.0 ... electronID 0.8
3 0.5 1.0 ... electronID 0.8
4 0.5 1.0 ... electronID 0.8
.. ... ... ... ... ...
65 4.0 5.0 ... electronID 0.8
66 4.0 5.0 ... electronID 0.8
67 4.0 5.0 ... electronID 0.8
68 4.0 5.0 ... electronID 0.8
69 4.0 5.0 ... electronID 0.8
[70 rows x 13 columns]}
cosTheta p charge electronID ... mcPDG B0_dec_mode B0_PDG B0_sigProb
0 0.250 1.5 -1 0.850 ... 11 0 511 0.02
1 0.550 2.5 1 0.900 ... 211 1 -511 0.03
2 0.850 3.5 1 0.950 ... -11 2 511 0.02
3 0.925 4.5 1 0.975 ... 211 3 -511 0.90
[4 rows x 9 columns]
==================================================
------------------
ReweighterParticle
Type: FEI Prefix:
Merged table:
data_MC_ratio ... data_MC_uncertainty_sys_up
count 4.000000 ... 4.00000
mean 0.805000 ... 0.05000
std 0.042032 ... 0.03559
min 0.750000 ... 0.02000
25% 0.787500 ... 0.02750
50% 0.810000 ... 0.04000
75% 0.827500 ... 0.06250
max 0.850000 ... 0.10000
[8 rows x 5 columns]
PDG binning:
(511, 511): {'dec_mode': ['mode0', 'mode2']}
------------------
------------------
ReweighterParticle
Type: PID Prefix:
Merged table:
data_MC_ratio ... data_MC_uncertainty_sys_up
count 154.000000 ... 154.000000
mean 1.433314 ... 0.009324
std 0.501422 ... 0.005378
min 0.744701 ... 0.001103
25% 0.982177 ... 0.004378
50% 1.148172 ... 0.009014
75% 1.943355 ... 0.013712
max 2.329627 ... 0.019873
[8 rows x 5 columns]
PDG binning:
(11, 11): {'charge': array([-2., 0., 2.]), 'p': array([0.5, 1. , 2. , 3. , 4. , 5. , 6. ]), 'cosTheta': array([-0.7 , -0.5 , 0. , 0.2 , 0.5 , 0.8 , 0.9 , 0.95])}
(11, 211): {'charge': array([-2., 0., 2.]), 'p': array([0.5, 1. , 2. , 3. , 4. , 5. ]), 'cosTheta': array([-0.7 , -0.5 , 0. , 0.2 , 0.5 , 0.8 , 0.9 , 0.95])}
------------------
{'charge': array([-2., 0., 2.]), 'p': array([0.5, 1. , 2. , 3. , 4. , 5. , 6. ]), 'cosTheta': array([-0.7 , -0.5 , 0. , 0.2 , 0.5 , 0.8 , 0.9 , 0.95])}
Required variables: ['p', 'cosTheta', 'charge', 'PDG', 'mcPDG']
Skip charge for plotting!
Required variables: ['B0_dec_mode', 'B0_PDG']
=== stderr of analysis/tests/test_sysvar.py ===
/data/b2soft/buildbot/light/build/analysis/tests/test_sysvar.py:15: FutureWarning: This module will soon be deprecated and eventually removed in a future release. Its functionality is being taken over by the standalone SysVar package: https://gitlab.desy.de/belle2/software/sysvar
from sysvar import Reweighter
..../data/b2soft/buildbot/light/build/lib/Linux_x86_64/opt/sysvar.py:567: FutureWarning: Setting an item of incompatible dtype is deprecated and will raise an error in a future version of pandas. Value 'Rest' has dtype incompatible with float64, please explicitly cast to a compatible dtype first.
binning_df.loc[binning_df['PDG'] == reco_pdg, _fei_mode_col] = particle.merged_table.query(
/data/b2soft/buildbot/light/build/lib/Linux_x86_64/opt/sysvar.py:586: PerformanceWarning: DataFrame is highly fragmented. This is usually the result of calling `frame.insert` many times, which has poor performance. Consider joining all columns at once using pd.concat(axis=1) instead. To get a de-fragmented frame, use `newframe = frame.copy()`
ntuple_df[f'{particle.get_varname(col)}'] = binning_df[col]
/data/b2soft/buildbot/light/build/lib/Linux_x86_64/opt/sysvar.py:586: PerformanceWarning: DataFrame is highly fragmented. This is usually the result of calling `frame.insert` many times, which has poor performance. Consider joining all columns at once using pd.concat(axis=1) instead. To get a de-fragmented frame, use `newframe = frame.copy()`
ntuple_df[f'{particle.get_varname(col)}'] = binning_df[col]
/data/b2soft/buildbot/light/build/lib/Linux_x86_64/opt/sysvar.py:586: PerformanceWarning: DataFrame is highly fragmented. This is usually the result of calling `frame.insert` many times, which has poor performance. Consider joining all columns at once using pd.concat(axis=1) instead. To get a de-fragmented frame, use `newframe = frame.copy()`
ntuple_df[f'{particle.get_varname(col)}'] = binning_df[col]
/data/b2soft/buildbot/light/build/lib/Linux_x86_64/opt/sysvar.py:586: PerformanceWarning: DataFrame is highly fragmented. This is usually the result of calling `frame.insert` many times, which has poor performance. Consider joining all columns at once using pd.concat(axis=1) instead. To get a de-fragmented frame, use `newframe = frame.copy()`
ntuple_df[f'{particle.get_varname(col)}'] = binning_df[col]
F
======================================================================
FAIL: test_reweight (__main__.TestSysVar.test_reweight)
Tests reweighting of PID and FEI particles
----------------------------------------------------------------------
Traceback (most recent call last):
File "/data/b2soft/buildbot/light/build/analysis/tests/test_sysvar.py", line 153, in test_reweight
self.assertTrue((local_data.query('abs(mcPDG) != 11')[['Weight']+cols].std(axis=1) > 0.015).all())
AssertionError: False is not true
----------------------------------------------------------------------
Ran 5 tests in 0.270s
FAILED (failures=1)
=== end analysis/tests/test_sysvar.py ===