From 5c2653e210f5685350c20be0a551749c39994b84 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 6 Jul 2020 09:48:59 -0500 Subject: only look at Michel electrons where the muon had < 2500 nhit --- utils/plot-muons | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) (limited to 'utils/plot-muons') 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) -- cgit