only look at Michel electrons where the muon had < 2500 nhit

1 files changed, 6 insertions, 2 deletions

diff --git a/utils/plot-muons b/utils/plot-muons
index a55c1ba..11d2d2e 100755
--- a/utils/plot-muons
+++ b/utils/plot-muons
@@ -255,8 +255,12 @@ if __name__ == '__main__':
             plt.title("Stopping Muon Energy Resolution")
 
     # For the Michel energy plot, we only look at the single particle electron
-    # fit
-    michel = michel[michel.id == 20]
+    # fit and events where the corresponding muon had less than 2500 nhit.  The
+    # reason for only looking at Michel electrons from muons with less than
+    # 2500 nhit is because there is significant ringing and afterpulsing after
+    # a large muon which can cause the reconstruction to overestimate the
+    # energy.
+    michel = michel[(michel.muon_nhit < 2500) & (michel.id == 20)]
 
     fig = plt.figure()
     bins = np.linspace(0,100,41)