From b11000942eb96c5b50a18030b0b8bb130a3701af Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 16 Nov 2020 08:00:44 -0600 Subject: loop over MC filenames --- utils/dc-closure-test | 10 +++++++++- 1 file changed, 9 insertions(+), 1 deletion(-) (limited to 'utils/dc-closure-test') 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)] -- cgit