diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2020-12-09 09:57:27 -0600 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2020-12-09 09:57:27 -0600 |
| commit | bf98750c36c0d5ebb0ed320ea20cad1dced8f02e (patch) | |
| tree | eb2bece62ae4768bbd99139c3138eed329adba80 | |
| parent | 8ca819c334cae0280118a2779041b0b7aaf4bbe2 (diff) | |
| download | sddm-bf98750c36c0d5ebb0ed320ea20cad1dced8f02e.tar.gz sddm-bf98750c36c0d5ebb0ed320ea20cad1dced8f02e.tar.bz2 sddm-bf98750c36c0d5ebb0ed320ea20cad1dced8f02e.zip | |
discard negative weights
| -rwxr-xr-x | utils/chi2 | 26 | ||||
| -rwxr-xr-x | utils/dm-search | 26 |
2 files changed, 52 insertions, 0 deletions
@@ -516,6 +516,32 @@ if __name__ == '__main__': muon_mc = get_events(args.muon_mc, merge_fits=True, nhit_thresh=args.nhit_thresh, mc=True) weights = pd.concat([read_hdf(filename, "weights") for filename in args.weights],ignore_index=True) + # There are a handful of weights which turn out to be slightly negative for + # some reason. For example: + # + # run evn universe weight + # 10970 25 597 -0.000055 + # 11389 87 729 -0.021397 + # 11701 204 2 -0.000268 + # 11919 120 82 -0.002245 + # 11976 163 48 -0.000306 + # 11976 163 710 -0.000022 + # 12131 76 175 -0.000513 + # 12207 70 255 -0.002925 + # 12207 70 282 -0.014856 + # 12207 70 368 -0.030593 + # 12207 70 453 -0.019011 + # 12207 70 520 -0.020748 + # 12207 70 834 -0.028754 + # 12207 70 942 -0.020309 + # 12233 230 567 -0.000143 + # 12618 168 235 -0.000020 + # 13428 128 42 -0.083639 + # 14264 23 995 -0.017637 + # 15034 69 624 -0.000143 + # 15752 154 957 -0.006827 + weights = weights[weights.weight > 0] + ev_mc = correct_energy_bias(ev_mc) muon_mc = correct_energy_bias(muon_mc) diff --git a/utils/dm-search b/utils/dm-search index 2679d53..34509ac 100755 --- a/utils/dm-search +++ b/utils/dm-search @@ -585,6 +585,32 @@ if __name__ == '__main__': muon_mc = get_events(args.muon_mc, merge_fits=True, nhit_thresh=args.nhit_thresh, mc=True) weights = pd.concat([read_hdf(filename, "weights") for filename in args.weights],ignore_index=True) + # There are a handful of weights which turn out to be slightly negative for + # some reason. For example: + # + # run evn universe weight + # 10970 25 597 -0.000055 + # 11389 87 729 -0.021397 + # 11701 204 2 -0.000268 + # 11919 120 82 -0.002245 + # 11976 163 48 -0.000306 + # 11976 163 710 -0.000022 + # 12131 76 175 -0.000513 + # 12207 70 255 -0.002925 + # 12207 70 282 -0.014856 + # 12207 70 368 -0.030593 + # 12207 70 453 -0.019011 + # 12207 70 520 -0.020748 + # 12207 70 834 -0.028754 + # 12207 70 942 -0.020309 + # 12233 230 567 -0.000143 + # 12618 168 235 -0.000020 + # 13428 128 42 -0.083639 + # 14264 23 995 -0.017637 + # 15034 69 624 -0.000143 + # 15752 154 957 -0.006827 + weights = weights[weights.weight > 0] + ev_mc = correct_energy_bias(ev_mc) muon_mc = correct_energy_bias(muon_mc) |