aboutsummaryrefslogtreecommitdiff
path: root/src/fit.c
diff options
context:
space:
mode:
authortlatorre <tlatorre@uchicago.edu>2020-04-13 12:00:30 -0500
committertlatorre <tlatorre@uchicago.edu>2020-04-13 12:00:30 -0500
commitaa2699809e9cd6cd48b57a011af278660e30d9ab (patch)
treeb8659df8188d828ab73d13ac4b216b20b1c2935c /src/fit.c
parent441871fe030ef8ba7ea9ad02fda4a1e9caf03f93 (diff)
downloadsddm-aa2699809e9cd6cd48b57a011af278660e30d9ab.tar.gz
sddm-aa2699809e9cd6cd48b57a011af278660e30d9ab.tar.bz2
sddm-aa2699809e9cd6cd48b57a011af278660e30d9ab.zip
fit for up to 5 peaks
This commit updates fit.c to start with 5 peaks for the direction seeds. I chose this number because I did some testing with the test-find-peaks program on the atmospheric MC and it looks like 5 peaks were necessary to capture the majority of the peaks.
Diffstat (limited to 'src/fit.c')
-rw-r--r--src/fit.c4
1 files changed, 2 insertions, 2 deletions
diff --git a/src/fit.c b/src/fit.c
index 2fb11b9..6e56f9b 100644
--- a/src/fit.c
+++ b/src/fit.c
@@ -5490,8 +5490,8 @@ int fit_event2(event *ev, double *xopt, double *fmin, int *id, size_t n, double
/* Find the peaks in the Hough transform of the event. */
find_peaks(ev, pos, t0, 100, 100, peak_theta, peak_phi, &npeaks, LEN(peak_theta), 0.1);
- /* Don't fit more than 3 peaks for now. */
- if (npeaks > 3) npeaks = 3;
+ /* Don't fit more than 5 peaks for now. */
+ if (npeaks > 5) npeaks = 5;
result = malloc(sizeof(size_t)*n*ipow(npeaks,n));