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/plot-michels | 10 +++++++++- 1 file changed, 9 insertions(+), 1 deletion(-) (limited to 'utils/plot-michels') diff --git a/utils/plot-michels b/utils/plot-michels index ff88315..b2727ca 100755 --- a/utils/plot-michels +++ b/utils/plot-michels @@ -113,7 +113,15 @@ if __name__ == '__main__': evs.append(get_events(df.filename.values, merge_fits=True, nhit_thresh=args.nhit_thresh)) ev = pd.concat(evs) ev = correct_energy_bias(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 = correct_energy_bias(ev_mc) # Drop events without fits -- cgit