diff options
Diffstat (limited to 'utils')
| -rwxr-xr-x | utils/chi2 | 8 | ||||
| -rwxr-xr-x | utils/dm-search | 8 |
2 files changed, 8 insertions, 8 deletions
@@ -660,8 +660,8 @@ if __name__ == '__main__': for i in range(args.coverage): # Calculate expected number of events - N = len(data_mc)*xtrue[0]/atmo_scale_factor - N_atm = len(data_atm_mc)*xtrue[0]/atmo_scale_factor + N = data_mc.flux_weight.sum()*xtrue[0]/atmo_scale_factor + N_atm = data_atm_mc.flux_weight.sum()*xtrue[0]/atmo_scale_factor N_muon = len(muon)*xtrue[5]/muon_scale_factor N_muon_atm = len(muon_atm)*xtrue[5]/muon_scale_factor @@ -761,8 +761,8 @@ if __name__ == '__main__': xtrue = truncnorm_scaled(PRIORS_LOW,PRIORS_HIGH,PRIORS,PRIOR_UNCERTAINTIES) # Calculate expected number of events - N = len(data_mc)*xtrue[0]/atmo_scale_factor - N_atm = len(data_atm_mc)*xtrue[0]/atmo_scale_factor + N = data_mc.flux_weight.sum()*xtrue[0]/atmo_scale_factor + N_atm = data_atm_mc.flux_weight.sum()*xtrue[0]/atmo_scale_factor N_muon = len(muon)*xtrue[5]/muon_scale_factor N_muon_atm = len(muon_atm)*xtrue[5]/muon_scale_factor diff --git a/utils/dm-search b/utils/dm-search index f8a8ff1..a99f210 100755 --- a/utils/dm-search +++ b/utils/dm-search @@ -730,8 +730,8 @@ if __name__ == '__main__': xtrue = sample_priors() # Calculate expected number of events - N = len(data_mc)*xtrue[0]/atmo_scale_factor - N_atm = len(data_atm_mc)*xtrue[0]/atmo_scale_factor + N = data_mc.flux_weight.sum()*xtrue[0]/atmo_scale_factor + N_atm = data_atm_mc.flux_weight.sum()*xtrue[0]/atmo_scale_factor N_muon = len(muon)*xtrue[5]/muon_scale_factor N_muon_atm = len(muon_atm)*xtrue[5]/muon_scale_factor N_dm = xtrue[6] @@ -811,8 +811,8 @@ if __name__ == '__main__': xtrue = sample_priors() # Calculate expected number of events - N = len(data_mc)*xtrue[0]/atmo_scale_factor - N_atm = len(data_atm_mc)*xtrue[0]/atmo_scale_factor + N = data_mc.flux_weight.sum()*xtrue[0]/atmo_scale_factor + N_atm = data_atm_mc.flux_weight.sum()*xtrue[0]/atmo_scale_factor N_muon = len(muon)*xtrue[5]/muon_scale_factor N_muon_atm = len(muon_atm)*xtrue[5]/muon_scale_factor |