From c11b9949669846849c86bc412ed6d1fd1e91630f Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 27 Jul 2020 11:18:51 -0500 Subject: add 20 MeV cut to dc, dc-closure-test, and plot-dc --- utils/dc | 1 + utils/dc-closure-test | 1 + utils/plot-dc | 1 + 3 files changed, 3 insertions(+) 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) diff --git a/utils/dc-closure-test b/utils/dc-closure-test index e23f554..b688b0f 100755 --- a/utils/dc-closure-test +++ b/utils/dc-closure-test @@ -455,6 +455,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) diff --git a/utils/plot-dc b/utils/plot-dc index 1ecbbf9..d2e461d 100755 --- a/utils/plot-dc +++ b/utils/plot-dc @@ -123,6 +123,7 @@ if __name__ == '__main__': ev = pd.concat(evs) ev = ev[ev.prompt & ~np.isnan(ev.fmin)] + ev = ev[ev.ke > 20] with pd.option_context('display.max_rows', None, 'display.max_columns', None): print("Noise events") -- cgit