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/test-find-peaks.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'src/test-find-peaks.c') diff --git a/src/test-find-peaks.c b/src/test-find-peaks.c index a825f4b..dd3a9d8 100644 --- a/src/test-find-peaks.c +++ b/src/test-find-peaks.c @@ -363,7 +363,7 @@ int main(int argc, char **argv) get_hough_transform(&ev,pos,x,y,n,m,result,0,0); - find_peaks(&ev,pos,n,m,peak_theta,peak_phi,&npeaks,max_npeaks); + find_peaks(&ev,pos,n,m,peak_theta,peak_phi,&npeaks,max_npeaks,0.1); printf("gtid %i\n", ev.gtid); for (i = 0; i < npeaks; i++) { -- cgit