diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2020-11-16 08:00:44 -0600 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2020-11-16 08:00:44 -0600 |
| commit | b11000942eb96c5b50a18030b0b8bb130a3701af (patch) | |
| tree | 8f3ea0f21d717db8827bea3643473ef7889bbab9 /utils/dc-closure-test | |
| parent | fbcdf03f205840c5cd550689e3f60f5168f8ecca (diff) | |
| download | sddm-b11000942eb96c5b50a18030b0b8bb130a3701af.tar.gz sddm-b11000942eb96c5b50a18030b0b8bb130a3701af.tar.bz2 sddm-b11000942eb96c5b50a18030b0b8bb130a3701af.zip | |
loop over MC filenames
Diffstat (limited to 'utils/dc-closure-test')
| -rwxr-xr-x | utils/dc-closure-test | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/utils/dc-closure-test b/utils/dc-closure-test index 3c53a1b..738e4f3 100755 --- a/utils/dc-closure-test +++ b/utils/dc-closure-test @@ -466,7 +466,15 @@ if __name__ == '__main__': ev = z_cut(ev) ev = udotr_cut(ev) - ev_mc = get_events(args.mc, merge_fits=True, apply_nhit_trigger=False) + # Note: We loop over the MC filenames here instead of just passing the + # whole list to get_events() because I had to rerun some of the MC events + # using SNOMAN and so most of the runs actually have two different files + # and otherwise the GTIDs will clash + ev_mcs = [] + for filename in args.mc: + ev_mcs.append(get_events([filename], merge_fits=True, nhit_thresh=args.nhit_thresh, apply_nhit_trigger=False)) + ev_mc = pd.concat(ev_mcs) + ev_mc = ev_mc[ev_mc.prompt] ev_mc = ev_mc[ev_mc.nhit_cal > 100] ev_mc = ev_mc[~np.isnan(ev_mc.fmin)] |