diff options
Diffstat (limited to 'src/dc.c')
| -rw-r--r-- | src/dc.c | 6 |
1 files changed, 3 insertions, 3 deletions
@@ -254,7 +254,7 @@ int is_muon(event *ev) * in this case for some reason the ECA time still looks reasonable. * Therefore, since we only use the ECA time here we instead just make * sure that the ECA time is something reasonable. */ - if (!ev->pmt_hits[i].hit || pmts[i].pmt_type != PMT_OWL || ev->pmt_hits[i].ept <= -100) continue; + if (!ev->pmt_hits[i].hit || pmts[i].pmt_type != PMT_OWL || ev->pmt_hits[i].ept <= -100 || ev->pmt_hits[i].best_uncal_q < -100) continue; n_nearby = 0; for (j = 0; j < MAX_PMTS; j++) { @@ -265,7 +265,7 @@ int is_muon(event *ev) if (distance < MUON_OWL_NEARBY_DISTANCE) { t_nearby[n_nearby] = ev->pmt_hits[j].ept; - q_nearby[n_nearby] = ev->pmt_hits[j].ehs; + q_nearby[n_nearby] = ev->pmt_hits[j].best_uncal_q; n_nearby += 1; } } @@ -278,7 +278,7 @@ int is_muon(event *ev) t_median = gsl_stats_median_from_sorted_data(t_nearby,1,n_nearby); q_median = gsl_stats_median_from_sorted_data(q_nearby,1,n_nearby); - if (ev->pmt_hits[i].ehs/q_median > 1.0 && ev->pmt_hits[i].ept < t_median) + if (ev->pmt_hits[i].best_uncal_q/q_median > 1.0 && ev->pmt_hits[i].ept < t_median) nhit_owl_early += 1; nhit_owl += 1; |