diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2020-12-21 02:22:55 -0600 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2020-12-21 02:22:55 -0600 |
| commit | 8353ec89afffc568f7160cd0c0e97cf458f89368 (patch) | |
| tree | f78850f5c532a5a5fdd9d5a59e6a7b573f0acf6f | |
| parent | 51e7b24254ae8261c455c5d8fe47544a29576c2c (diff) | |
| download | sddm-8353ec89afffc568f7160cd0c0e97cf458f89368.tar.gz sddm-8353ec89afffc568f7160cd0c0e97cf458f89368.tar.bz2 sddm-8353ec89afffc568f7160cd0c0e97cf458f89368.zip | |
use correct bins when calculation threshold
| -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 63f6438..823874b 100755 --- a/utils/dm-search +++ b/utils/dm-search @@ -514,7 +514,7 @@ def get_limits(dm_masses,data,muon,data_mc,atmo_scale_factor,muon_scale_factor,b # approximately equal to the number of trials so we need to # increase our discovery threshold to account for the look # elsewhere effect. - threshold = DISCOVERY_P_VALUE/(len(bins)-1) + threshold = DISCOVERY_P_VALUE/(len(bins[dm_particle_id])-1) discovery = poisson.ppf(1-threshold,n) + 1 - n # Here, we scale the discovery threshold by the fraction of the # dark matter hist in the histogram range. The idea is that if only |