From f92ebc3df77a24f67eeb7b78c45d03b7806bdfa3 Mon Sep 17 00:00:00 2001
From: tlatorre <tlatorre@uchicago.edu>
Date: Mon, 16 Nov 2020 08:31:20 -0600
Subject: fix integrate_mc_hist() by dividing by bin width

---
 utils/dm-search | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/utils/dm-search b/utils/dm-search
index 4e8772e..a486593 100755
--- a/utils/dm-search
+++ b/utils/dm-search
@@ -399,7 +399,7 @@ def integrate_mc_hist(mc, muons, id, x, a, b, atmo_scale_factor, muon_scale_fact
         i = np.digitize(x,bins)
         if i == len(bins):
             i = len(bins) - 1
-        return mc_hists[id][i-1] + muon_hists[id][i-1]
+        return (mc_hists[id][i-1] + muon_hists[id][i-1])/(bins[i]-bins[i-1])
 
     # Yes, this is a stupid simple integration that I don't need quad for, but
     # I don't want to have to code all the corner cases including when a and b
-- 
cgit