diff options
Diffstat (limited to 'src/find_peaks.c')
| -rw-r--r-- | src/find_peaks.c | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/src/find_peaks.c b/src/find_peaks.c index 7d350d5..6109896 100644 --- a/src/find_peaks.c +++ b/src/find_peaks.c @@ -197,8 +197,9 @@ void get_hough_transform(event *ev, double *pos, double *x, double *y, size_t n, /* Subtract off previous rings. */ for (i = 0; i < len; i++) { qhs_avg = get_qhs_avg(ev,pos,last+3*i); + for (j = 0; j < MAX_PMTS; j++) { - if (ev->pmt_hits[i].flags || pmts[i].pmt_type != PMT_NORMAL || !ev->pmt_hits[i].hit) continue; + if (ev->pmt_hits[j].flags || pmts[j].pmt_type != PMT_NORMAL || !ev->pmt_hits[j].hit) continue; SUB(pmt_dir,pmts[j].pos,pos); @@ -304,7 +305,7 @@ void find_peaks(event *ev, double *pos, size_t n, size_t m, double *peak_theta, /* Only add the latest peak to the results if it's more than `delta` * radians away from all other results so far. */ unique = 1; - for (j = 0; j < i; j++) { + for (j = 0; j < *npeaks; j++) { /* If the angle between the latest peak and a previous peak is * within `delta` radians, it's not unique. */ if (acos(DOT(dir+i*3,dir+j*3)) < delta) unique = 0; |