Changeset 31313 for branches/eam_branches/ipp-20110404/psphot/src/psphotFitSourcesLinearStack.c

Timestamp:

Apr 18, 2011, 8:55:26 AM (15 years ago)

Author:

eugene

Message:

raise bit for PASS1 detections; set MODE_PSFMODEL in linear fit as well as nonlinear fit; move flags for fitting results into pmSourceFits not the calling function (more explicitly sets the bits in the right cases); adding test MaskedKron function

File:

• branches/eam_branches/ipp-20110404/psphot/src/psphotFitSourcesLinearStack.c (modified) (1 diff)

branches/eam_branches/ipp-20110404/psphot/src/psphotFitSourcesLinearStack.c

@@ -69 +69 @@
             if (!pmSourceCacheModel (source, maskVal)) continue;
           }
+
+          // check the integral of the model : is it large enough?
+          float modelSum = 0.0;
+          for (int iy = 0; iy < source->modelFlux->numRows; iy++) {
+              for (int ix = 0; ix < source->modelFlux->numCols; ix++) {
+                  modelSum += source->modelFlux->data.F32[iy][ix];
+              }
+          }
+          if (modelSum < 0.5) continue; // skip sources with no model constraint (somewhat arbitrary limit)
+          if (modelSum < 0.8) {
+              fprintf (stderr, "low-sig model @ %f, %f (%f sum, %f peak)\n",
+                       source->peak->xf, source->peak->yf, modelSum, source->peak->rawFlux);
+          }
 
           source->mode |= PM_SOURCE_MODE_LINEAR_FIT;