diff options
author | tlatorre <tlatorre@uchicago.edu> | 2018-09-25 09:40:02 -0500 |
---|---|---|
committer | tlatorre <tlatorre@uchicago.edu> | 2018-09-25 09:40:02 -0500 |
commit | d50b813868ed5e1cb57cb7adae0a00f2e8f42e52 (patch) | |
tree | 1c9ece7c34a92f6b2ba5a301117299a8469779a6 | |
parent | 58c286a6f2d9d593d46888681bdf33066a21c70a (diff) | |
download | sddm-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.c | 2 |
1 files changed, 1 insertions, 1 deletions
@@ -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)); |