From f8c99fdd302b6d2190ae619681fb1cea5f74f391 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Thu, 13 Jun 2019 11:08:04 -0500 Subject: add a data cleaning cut to tag incoming muons This commit adds a data cleaning cut to tag incoming muons by looking for early OWL hits. It also significantly updates the flasher cut to catch more flashers. In particular, the flasher cut now does the following: - loops over *all* paddle cards with at least 4 hits instead of just the paddle cards with the most hits - uses QLX to look for charge outliers in the paddle card - fixes a few bugs (for example, uninitialized values in the charge array) - adds a check to to see if the given slot is early with respect to all PMTs within 4 meters to catch the case where the flashing channel is missing from the event --- src/fit.c | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) (limited to 'src/fit.c') diff --git a/src/fit.c b/src/fit.c index 40f67f5..e4e054f 100644 --- a/src/fit.c +++ b/src/fit.c @@ -6172,7 +6172,7 @@ skip_mc: last_run = ev.run; } - rv = get_event(f,&ev,&b); + if (get_event(f,&ev,&b)) goto skip_event; nhit = get_nhit(&ev); @@ -6181,6 +6181,11 @@ skip_mc: fprintf(fout, " gtid: %i\n", ev.gtid); fprintf(fout, " nhit: %zu\n", nhit); fprintf(fout, " is_flasher: %i\n", is_flasher(&ev)); + fprintf(fout, " is_muon: %i\n", is_muon(&ev)); + fprintf(fout, " is_junk: %i\n", junk_cut(f, &bmast, &b)); + fprintf(fout, " is_crate_isotropy: %i\n", crate_isotropy(&ev)); + fprintf(fout, " is_qvnhit: %i\n", qvnhit(&ev)); + fprintf(fout, " is_neck_event: %i\n", is_neck_event(&ev)); } if (nhit < min_nhit) goto skip_event; -- cgit