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.

1 files changed, 5 insertions, 5 deletions

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);