hack to get rid of flasher and muon events in breakdown sample

3 files changed, 9 insertions, 0 deletions

diff --git a/utils/dc b/utils/dc
index 39ca008..ee9a1c5 100755
--- a/utils/dc
+++ b/utils/dc
@@ -282,6 +282,9 @@ if __name__ == '__main__':
         evs.append(get_events(df.filename.values, merge_fits=True))
     ev = pd.concat(evs)
 
+    # Hack to get rid of flasher and muon events in the breakdown sample.
+    ev.breakdown &= ev.r_psup < 0.9
+
     ev = ev[ev.prompt]
     ev = ev[ev.nhit_cal > 100]
 
diff --git a/utils/dc-closure-test b/utils/dc-closure-test
index 82ca005..2628edb 100755
--- a/utils/dc-closure-test
+++ b/utils/dc-closure-test
@@ -439,6 +439,9 @@ if __name__ == '__main__':
         evs.append(get_events(df.filename.values, merge_fits=True, nhit_thresh=args.nhit_thresh))
     ev = pd.concat(evs)
 
+    # Hack to get rid of flasher and muon events in the breakdown sample.
+    ev.breakdown &= ev.r_psup < 0.9
+
     ev = ev[ev.prompt]
     ev = ev[ev.nhit_cal > 100]
 
diff --git a/utils/plot-dc b/utils/plot-dc
index acd7e7b..bd520d6 100755
--- a/utils/plot-dc
+++ b/utils/plot-dc
@@ -127,6 +127,9 @@ if __name__ == '__main__':
             evs.append(get_events([filename], merge_fits=True, mc=True))
     ev = pd.concat([ev for ev in evs if len(ev) > 0])
 
+    # Hack to get rid of flasher and muon events in the breakdown sample.
+    ev.breakdown &= ev.r_psup < 0.9
+
     ev = ev[ev.prompt & ~np.isnan(ev.fmin)]
     ev = ev[ev.ke > 20]