From 2f2cd1797190335468e956f4c286c6e4de27518e Mon Sep 17 00:00:00 2001 From: tlatorre Date: Sun, 2 Jun 2019 13:21:30 -0400 Subject: update find_peaks() to only return unique peaks This commit updates find_peaks() to only return peaks which are at least a certain number of degrees apart from each other. This is because I found that for many events the first few peaks would all be essentially the same direction and so the fit was taking a lot of time fitting essentially the same seed points. Since I now have to only try 3 peaks in order to get my grid jobs to run for less than a few hours it's necessary to make sure we aren't just fitting the same three directions for the "quick" minimization. I also updated the fit to only use a maximum of 3 seed directions. --- src/find_peaks.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'src/find_peaks.h') diff --git a/src/find_peaks.h b/src/find_peaks.h index 8511af6..0219561 100644 --- a/src/find_peaks.h +++ b/src/find_peaks.h @@ -22,6 +22,6 @@ void find_peaks_array(double *x, size_t n, size_t m, size_t *imax, size_t *jmax, size_t *npeaks, size_t max_peaks, double threshold); void get_hough_transform(event *ev, double *pos, double *x, double *y, size_t n, size_t m, double *result, double *last, size_t len); -void find_peaks(event *ev, double *pos, size_t n, size_t m, double *peak_theta, double *peak_phi, size_t *npeaks, size_t max_peaks); +void find_peaks(event *ev, double *pos, size_t n, size_t m, double *peak_theta, double *peak_phi, size_t *npeaks, size_t max_peaks, double delta); #endif -- cgit