Added warning when zphot falls out of confidence interval

Author: cschreib

src/fast++-fitter.cpp

@@ -69,8 +69,17 @@ fitter_t::fitter_t(const options_t& opt, const input_state_t& inp, const gridder
                 idzu.safe[is] = idz.safe[is];
             } else if (is_finite(input.zphot.safe(is,ilow)) && is_finite(input.zphot.safe(is,iup))) {
                 // Get range from zlow and zup
-                idzl.safe[is] = min_id(abs(output_z - input.zphot.safe(is,ilow)));
-                idzu.safe[is] = min_id(abs(output_z - input.zphot.safe(is,iup)));
+                idzl.safe[is] = min_id(abs(output.z - input.zphot.safe(is,ilow)));
+                idzu.safe[is] = min_id(abs(output.z - input.zphot.safe(is,iup)));
+
+                // Check that zphot falls inside confidence interval
+                if (opts.best_at_zphot && idzp.safe[is] != npos &&
+                    (idzp.safe[is] < idzl.safe[is] || idzp.safe[is] > idzu.safe[is])) {
+                    warning("for galaxy ", input.id.safe[is], " the photo-z (",
+                        input.zphot.safe(is,0), ") falls outside of the chosen confidence "
+                        "interval (", input.zphot.safe(is,ilow), " to ", input.zphot.safe(is,iup),
+                        ")");
+                }
             }
         }
     }