From da14d3f2362af2d8f49dd805c0c762cb3276409f Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 9 Sep 2019 11:21:10 -0500 Subject: use EHS instead of QHS in is_muon() This commit updates is_muon() to use the EHS charge for the OWL tubes since based on looking at the charges for the OWL tubes it looks like the QHS values are completely random. --- src/dc.c | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) (limited to 'src') diff --git a/src/dc.c b/src/dc.c index 3701bf8..5f2d4c6 100644 --- a/src/dc.c +++ b/src/dc.c @@ -231,7 +231,7 @@ int is_muon(event *ev) size_t i; double ept[MAX_PMTS]; size_t nhit, nhit_owl, nhit_owl_early; - double early_time, max_owl_qhs; + double early_time, max_owl_ehs; nhit = 0; nhit_owl = 0; @@ -253,9 +253,9 @@ int is_muon(event *ev) switch (pmts[i].pmt_type) { case PMT_OWL: if (nhit_owl == 0) { - max_owl_qhs = ev->pmt_hits[i].qhs; - } else if (ev->pmt_hits[i].qhs > max_owl_qhs) { - max_owl_qhs = ev->pmt_hits[i].qhs; + max_owl_ehs = ev->pmt_hits[i].ehs; + } else if (ev->pmt_hits[i].ehs > max_owl_ehs) { + max_owl_ehs = ev->pmt_hits[i].ehs; } nhit_owl += 1; break; @@ -272,7 +272,7 @@ int is_muon(event *ev) if (nhit_owl < MUON_MIN_OWL_NHIT) return 0; /* Require at least 1 OWL tube to have a QHS >= 20.0. */ - if (max_owl_qhs < MUON_OWL_QHS) return 0; + if (max_owl_ehs < MUON_OWL_QHS) return 0; /* Sort the times. */ gsl_sort(ept,1,nhit); -- cgit