135 LOGL_DEBUG(
"TRACE2.ip.diffim.BuildSingleKernelVisitor.processCandidate",
136 "Failed to cast SpatialCellCandidate to KernelCandidate %d",
139 "Failed to cast SpatialCellCandidate to KernelCandidate");
146 LOGL_DEBUG(
"TRACE1.ip.diffim.BuildSingleKernelVisitor.processCandidate",
147 "Processing candidate %d",
kCandidate->getId());
148 LOGL_DEBUG(
"TRACE4.ip.diffim.BuildSingleKernelVisitor.processCandidate",
155 if (_useRegularization)
161 kCandidate->setStatus(afwMath::SpatialCellCandidate::BAD);
162 LOGL_DEBUG(
"TRACE3.ip.diffim.BuildSingleKernelVisitor.processCandidate",
163 "Unable to process candidate %d; exception caught (%s)",
170 if (
kCandidate->getStatus() == afwMath::SpatialCellCandidate::BAD) {
171 LOGL_DEBUG(
"TRACE3.ip.diffim.BuildSingleKernelVisitor.processCandidate",
172 "Candidate %d Returned BAD upon build, exiting",
186 _imstats.apply(diffim, _coreRadius);
188 _imstats.apply(diffim);
190 LOGL_DEBUG(
"TRACE2.ip.diffim.BuildSingleKernelVisitor.processCandidate",
191 "Unable to calculate imstats for Candidate %d",
kCandidate->getId());
192 kCandidate->setStatus(afwMath::SpatialCellCandidate::BAD);
204 LOGL_DEBUG(
"TRACE4.ip.diffim.BuildSingleKernelVisitor.processCandidate",
206 LOGL_DEBUG(
"TRACE4.ip.diffim.BuildSingleKernelVisitor.processCandidate",
207 "Kernel Sum = %.3f",
kSum);
208 LOGL_DEBUG(
"TRACE4.ip.diffim.BuildSingleKernelVisitor.processCandidate",
209 "Background = %.3f", background);
210 LOGL_DEBUG(
"TRACE2.ip.diffim.BuildSingleKernelVisitor.processCandidate",
211 "Candidate %d resids = %.3f +/- %.3f sigma (%d pix)",
220 kCandidate->setStatus(afwMath::SpatialCellCandidate::BAD);
221 LOGL_DEBUG(
"TRACE3.ip.diffim.BuildSingleKernelVisitor.processCandidate",
222 "Rejecting candidate %d, encountered NaN",
228 if (_ps->getAsBool(
"singleKernelClipping")) {
229 if (fabs(_imstats.getMean()) > _ps->getAsDouble(
"candidateResidualMeanMax")) {
230 kCandidate->setStatus(afwMath::SpatialCellCandidate::BAD);
231 LOGL_DEBUG(
"TRACE3.ip.diffim.BuildSingleKernelVisitor.processCandidate",
232 "Rejecting candidate %d; bad mean residual : |%.3f| > %.3f",
235 _ps->getAsDouble(
"candidateResidualMeanMax"));
238 else if (_imstats.getRms() > _ps->getAsDouble(
"candidateResidualStdMax")) {
239 kCandidate->setStatus(afwMath::SpatialCellCandidate::BAD);
240 LOGL_DEBUG(
"TRACE3.ip.diffim.BuildSingleKernelVisitor.processCandidate",
241 "Rejecting candidate %d; bad residual rms : %.3f > %.3f",
244 _ps->getAsDouble(
"candidateResidualStdMax"));
248 kCandidate->setStatus(afwMath::SpatialCellCandidate::GOOD);
249 LOGL_DEBUG(
"TRACE3.ip.diffim.BuildSingleKernelVisitor.processCandidate",
254 kCandidate->setStatus(afwMath::SpatialCellCandidate::GOOD);
255 LOGL_DEBUG(
"TRACE5.ip.diffim.BuildSingleKernelVisitor.processCandidate",
256 "Sigma clipping not enabled");
260 if (!(_useCoreStats)) {
262 _imstats.apply(diffim, _coreRadius);
264 LOGL_DEBUG(
"TRACE2.ip.diffim.BuildSingleKernelVisitor.processCandidate",
265 "Unable to calculate core imstats for Candidate %d",
267 kCandidate->setStatus(afwMath::SpatialCellCandidate::BAD);
270 LOGL_DEBUG(
"TRACE3.ip.diffim.BuildSingleKernelVisitor.processCandidate",
271 "Candidate %d core resids = %.3f +/- %.3f sigma (%d pix)",