diff options
Diffstat (limited to 'utils/plot-michels')
| -rwxr-xr-x | utils/plot-michels | 5 |
1 files changed, 2 insertions, 3 deletions
diff --git a/utils/plot-michels b/utils/plot-michels index e0e4dd7..300562a 100755 --- a/utils/plot-michels +++ b/utils/plot-michels @@ -96,7 +96,6 @@ if __name__ == '__main__': parser.add_argument("filenames", nargs='+', help="input files") parser.add_argument("--save", action='store_true', default=False, help="save corner plots for backgrounds") parser.add_argument("--mc", nargs='+', required=True, help="atmospheric MC files") - parser.add_argument("--nhit-thresh", type=int, default=None, help="nhit threshold to apply to events before processing (should only be used for testing to speed things up)") parser.add_argument("--print-nll", action='store_true', default=False, help="print nll values") parser.add_argument("--steps", type=int, default=1000, help="number of steps in the MCMC chain") parser.add_argument("--walkers", type=int, default=100, help="number of walkers") @@ -110,7 +109,7 @@ if __name__ == '__main__': evs = [] rhdr = pd.concat([read_hdf(filename, "rhdr").assign(filename=filename) for filename in args.filenames],ignore_index=True) for run, df in rhdr.groupby('run'): - evs.append(get_events(df.filename.values, merge_fits=True, nhit_thresh=args.nhit_thresh)) + evs.append(get_events(df.filename.values, merge_fits=True)) ev = pd.concat(evs) ev = correct_energy_bias(ev) @@ -120,7 +119,7 @@ if __name__ == '__main__': # 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, mc=True)) + ev_mcs.append(get_events([filename], merge_fits=True, mc=True)) ev_mc = pd.concat(ev_mcs) ev_mc = correct_energy_bias(ev_mc) |