aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authortlatorre <tlatorre@uchicago.edu>2018-09-25 09:40:02 -0500
committertlatorre <tlatorre@uchicago.edu>2018-09-25 09:40:02 -0500
commitd50b813868ed5e1cb57cb7adae0a00f2e8f42e52 (patch)
tree1c9ece7c34a92f6b2ba5a301117299a8469779a6
parent58c286a6f2d9d593d46888681bdf33066a21c70a (diff)
downloadsddm-d50b813868ed5e1cb57cb7adae0a00f2e8f42e52.tar.gz
sddm-d50b813868ed5e1cb57cb7adae0a00f2e8f42e52.tar.bz2
sddm-d50b813868ed5e1cb57cb7adae0a00f2e8f42e52.zip
increase maxeval to 20 for the "quick" minimization phase
When testing out the fitter on 500 MeV muons, there was at least one event which started out at a position very far from the true position. This event had a secondary electron like ring which is what I suspect caused the fit to start out in a position far from the true position. This fix correctly starts the minimization close to the true position. In the future I should look at updating get_direction() so that it finds the largest ring direction instead of just doing a weighted average of all the vectors from the position to the PMTs.
-rw-r--r--src/fit.c2
1 files changed, 1 insertions, 1 deletions
diff --git a/src/fit.c b/src/fit.c
index 8ebaeec..9bcc580 100644
--- a/src/fit.c
+++ b/src/fit.c
@@ -5189,7 +5189,7 @@ int fit_event(event *ev, double *xopt, double *fmin)
fpars.fast = 1;
fpars.print = 0;
nlopt_set_ftol_abs(opt, 1.0);
- nlopt_set_maxeval(opt, 10);
+ nlopt_set_maxeval(opt, 20);
for (i = 0; i < sizeof(startingParameters)/sizeof(startingParameters[0]); i++) {
memcpy(x,x0,sizeof(x));