Changeset 29124 for branches/eam_branches/ipp-20100823/psphot/src/psphotExtendedSourceFits.c

Timestamp:

Sep 7, 2010, 9:16:10 AM (16 years ago)

Author:

eugene

Message:

updates from trunk

Location:

branches/eam_branches/ipp-20100823

Files:

• . (modified) (1 prop)
• psphot/src/psphotExtendedSourceFits.c (modified) (6 diffs)

branches/eam_branches/ipp-20100823/psphot/src/psphotExtendedSourceFits.c

@@ -200 +200 @@
             }
             psFree(job);
-            }
+        }
     }
     psFree (cellGroups);
@@ -270 +270 @@
         // set the radius based on the footprint (also sets the mask pixels)
         if (!psphotSetRadiusFootprint(&radius, readout, source, markVal, 1.0)) {
+            fprintf (stderr, "skipping (1) %f, %f\n", source->peak->xf, source->peak->yf);
             psFree (fitOptions)
             return false;
@@ -280 +281 @@
         // XXX save the psf-based moments for output
         if (!pmSourceMoments (source, radius, 0.0, 0.0, maskVal)) {
+            fprintf (stderr, "skipping (2) %f, %f\n", source->peak->xf, source->peak->yf);
             // subtract the best fit from the object, leave local sky
             pmSourceSub (source, PM_MODEL_OP_FULL, maskVal);
@@ -289 +291 @@
         psImage *modelFluxStart = psMemIncrRefCounter (source->modelFlux);
         if (!modelFluxStart) {
-          // XXX raise an error of some kind?
-          continue;
+            pmSourceCacheModel (source, maskVal);
+            modelFluxStart = psMemIncrRefCounter (source->modelFlux);
+            if (!modelFluxStart) {
+                fprintf (stderr, "skipping (3) %f, %f\n", source->peak->xf, source->peak->yf);
+                // XXX raise an error of some kind?
+                continue;
+            }
         }
-
+        
         if (savePics) {
           psphotSaveImage (NULL, readout->image, "image.xp.fits");
@@ -317 +324 @@
           float SNlim = psMetadataLookupF32 (&status, model, "SNLIM_VALUE");
           assert (status);
+
+          // fprintf (stderr, "xfit: %f, %f : %f\n", source->peak->xf, source->peak->yf, source->peak->SN); 
 
           // limit selection to some SN limit
@@ -449 +458 @@
     psFree (fitOptions);
 
+    // fprintf (stderr, "xfit : tried %ld objects\n", sources->n);
+
     // change the value of a scalar on the array (wrap this and put it in psArray.h)
     scalar = job->args->data[7];