aboutsummaryrefslogtreecommitdiff
path: root/utils/plot-michels
diff options
context:
space:
mode:
authortlatorre <tlatorre@uchicago.edu>2021-01-04 13:03:03 -0600
committertlatorre <tlatorre@uchicago.edu>2021-01-04 13:03:03 -0600
commit3d13f8f409d412c6d2aa3447b81e29a1a3d7b2a7 (patch)
tree7c4a216353b7a149c25148a8b2435d20ca07b196 /utils/plot-michels
parent64e5d47398a428ae7eca9e3c1251b9d9a35a823d (diff)
downloadsddm-3d13f8f409d412c6d2aa3447b81e29a1a3d7b2a7.tar.gz
sddm-3d13f8f409d412c6d2aa3447b81e29a1a3d7b2a7.tar.bz2
sddm-3d13f8f409d412c6d2aa3447b81e29a1a3d7b2a7.zip
get rid of nhit_thresh
Diffstat (limited to 'utils/plot-michels')
-rwxr-xr-xutils/plot-michels5
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)