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.

1 files 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]