add 20 MeV cut to dc, dc-closure-test, and plot-dc

1 files changed, 1 insertions, 0 deletions

diff --git a/utils/dc b/utils/dc
index 40ae13e..0cc4ba3 100755
--- a/utils/dc
+++ b/utils/dc
@@ -303,6 +303,7 @@ if __name__ == '__main__':
             fitted_fraction[bg] = 1.0
 
     ev = ev[~np.isnan(ev.fmin)]
+    ev = ev[ev.ke > 20]
 
     # figure out bins for high level variables
     ev = radius_cut(ev)