diff options
Diffstat (limited to 'src')
| -rw-r--r-- | src/dc.c | 18 |
1 files changed, 9 insertions, 9 deletions
@@ -451,7 +451,12 @@ int crate_isotropy(event *ev) * throw away the 10% of the tubes that have the largest charge and then * divide this by 0.9* NHIT. Tag the event if this ratio is below 0.25. * - * I copied the logic from the SNOMAN file flt_q_nhit_cut.for. */ + * I copied the logic from the SNOMAN file flt_q_nhit_cut.for. + * + * Update: I changed the logic from the SNOMAN code to not require good + * calibrations. The reason is that for electrical pickup events, the charges + * on most of the hits can be very small or negative and so we wouldn't tag + * them with this cut. */ int qvnhit(event *ev) { size_t i; @@ -461,13 +466,10 @@ int qvnhit(event *ev) nhit = 0; for (i = 0; i < MAX_PMTS; i++) { - if (ev->pmt_hits[i].flags || pmts[i].pmt_type != PMT_NORMAL) continue; - - /* Require good calibrations. */ - if (ev->pmt_hits[i].pf & (KPF_NO_CAL | KPF_BAD_CAL)) continue; + if (!ev->pmt_hits[i].hit || pmts[i].pmt_type != PMT_NORMAL) continue; /* Remove unphysical charge hits. */ - if (ev->pmt_hits[i].qhl <= -100) continue; + if (ev->pmt_hits[i].ehl <= -100) continue; /* FIXME: SNOMAN code checks for pmt terminator. I don't know where * that info is stored. */ @@ -478,9 +480,7 @@ int qvnhit(event *ev) * we use QHL since it is a better measure of pickup (long integrate) * */ - if (ev->pmt_hits[i].hit) { - ehl[nhit++] = ev->pmt_hits[i].ehl/32.3; - } + ehl[nhit++] = ev->pmt_hits[i].ehl/32.3; } gsl_sort(ehl,1,nhit); |