diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2020-11-16 08:31:20 -0600 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2020-11-16 08:31:20 -0600 |
| commit | f92ebc3df77a24f67eeb7b78c45d03b7806bdfa3 (patch) | |
| tree | 310508e443ce950bae711ba259b9c308d60a081e /utils/dm-search | |
| parent | 67dd8aa0737bd362dab382500149ee9c3e8637ab (diff) | |
| download | sddm-f92ebc3df77a24f67eeb7b78c45d03b7806bdfa3.tar.gz sddm-f92ebc3df77a24f67eeb7b78c45d03b7806bdfa3.tar.bz2 sddm-f92ebc3df77a24f67eeb7b78c45d03b7806bdfa3.zip | |
fix integrate_mc_hist() by dividing by bin width
Diffstat (limited to 'utils/dm-search')
| -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 |