Error bands broken (Bethe Block and Mass_Model_Proton)
daughter package used to notice the bug: CCQENu1DME
location of ntuples used: /pnfs/minerva/persistent/users/drut1186/HopefullyFinalMateusTuples/
location of histograms that contain the bug: /pnfs/minerva/persistent/users/mateusc/CCQENu_v21r1p1_Pub_May2019_NEWMuonEnergySys_NOSysFluxUniverses_CV/
version of CCQENuUtils.cxx: 1.60
Bethe Block and Mass_Model_Proton error bands have points with infinite values. This was noticed after the implementation of the new number of universes for the Muon variables (2 vs the usual 100) and subsequent fixes in the lateral errors filling. The bug is present in the histograms with full statistics but not in every playlist, me1L does not present any problem.
#1 Updated by Daniel Ruterbories over 1 year ago
- Assignee set to Daniel Ruterbories
Okay looking this.
I think the issue is once an array element is set not passing the rest of the events are always 0. Because of how the code ran before the shifts could be reset by the other functions. This is wrong or incredibly clever and not transparent. The recent changes make the behavior explicit, but missed this last step. The solution is to change the if(!passRecoil) call.
if(passRecoi) set shifts to 0
else set shifts to non-physical
So no matter what we get the values we want.