diff options
| -rwxr-xr-x | utils/dm-search | 2 |
1 files changed, 1 insertions, 1 deletions
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 |