Changeset 31452 for trunk/psphot/src/psphotFitSourcesLinearStack.c

Timestamp:

May 5, 2011, 11:09:38 AM (15 years ago)

Author:

eugene

Message:

merge updates from eam branch 20110404

Location:

trunk/psphot

Files:

• . (modified) (1 prop)
• src/psphotFitSourcesLinearStack.c (modified) (1 diff)

trunk/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;