From 782977045c55beb80268dc7211fcc657eb1ad5e0 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 17 Aug 2020 11:53:40 -0500 Subject: update stoppig muon criteria This commit updates the criteria for selecting stopping muons from: - calibrated nhit < 4000 - udotr < -0.5 to - reconstructed kinetic energy < 10 GeV The previous criteria were intended to remove through going atmospheric events but produced a strong bias in the comparison due to the nhit cut and an energy bias in the data relative to the Monte Carlo. The new cut does a good job of cutting through going muons but doesn't produce the same bias. --- utils/plot-muons | 7 ++----- 1 file changed, 2 insertions(+), 5 deletions(-) diff --git a/utils/plot-muons b/utils/plot-muons index fc8427c..d14f8b6 100755 --- a/utils/plot-muons +++ b/utils/plot-muons @@ -153,11 +153,8 @@ if __name__ == '__main__': stopping_muons = pd.merge(ev[ev.muon & ev.stopping_muon],michel,left_on=['run','gtid'],right_on=['run','muon_gtid'],suffixes=('','_michel')) stopping_muons_mc = pd.merge(ev_mc[ev_mc.muon & ev_mc.stopping_muon],michel_mc,left_on=['run','gtid'],right_on=['run','muon_gtid'],suffixes=('','_michel')) - stopping_muons = stopping_muons[stopping_muons.nhit_cal < 4000] - stopping_muons_mc = stopping_muons_mc[stopping_muons_mc.nhit_cal < 4000] - - stopping_muons = stopping_muons[stopping_muons.udotr < -0.5] - stopping_muons_mc = stopping_muons_mc[stopping_muons_mc.udotr < -0.5] + stopping_muons = stopping_muons[stopping_muons.ke < 10e3] + stopping_muons_mc = stopping_muons_mc[stopping_muons_mc.ke < 10e3] stopping_muons = stopping_muons[stopping_muons.cos_theta < -0.5] stopping_muons_mc = stopping_muons_mc[stopping_muons_mc.cos_theta < -0.5] -- cgit