Not fitting at zphot if outside of confidence interval and BEST_AT_ZPHOT is not set

Author: cschreib

src/fast++-fitter.cpp

@@ -281,7 +281,7 @@ void fitter_t::fit_galaxies(const model_t& model, uint_t i0, uint_t i1) {
         uint_t is = i + i0;
 
         // Apply constraints on redshift
-        if ((idzp.safe[is] == npos || model.iz != idzp.safe[is]) &&
+        if ((!opts.best_at_zphot || idzp.safe[is] == npos || model.iz != idzp.safe[is]) &&
             (idz.safe[is] == npos || model.iz != idz.safe[is]) &&
             (model.iz < idzl.safe[is] || model.iz > idzu.safe[is])) {
             wsp.chi2.safe[i] = finf;