Changeset 13426 for trunk/psModules/src/objects/pmPSF.c

Timestamp:

May 18, 2007, 4:08:34 AM (19 years ago)

Author:

rhl

Message:

Check {modelFromPSFFunc,pmSourcePhotometryAper,psEllipsePolToAxes}'s return code;

File:

• trunk/psModules/src/objects/pmPSF.c (modified) (8 diffs)

trunk/psModules/src/objects/pmPSF.c

@@ -6 +6 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.22 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2007-04-27 22:14:08 $
+ *  @version $Revision: 1.23 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2007-05-18 14:08:34 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -36 +36 @@
 #include "pmFPAMaskWeight.h"
 #include "psVectorBracket.h"
+#include "pmErrorCodes.h"
 
 /*****************************************************************************/
@@ -178 +179 @@
 
     // set model parameters for this source based on PSF information
-    modelFromPSFFunc (modelPSF, modelEXT, psf);
+    if (!modelFromPSFFunc (modelPSF, modelEXT, psf)) {
+        psError(PM_ERR_PSF, true, "Failed to set model params from PSF");
+        psFree(modelPSF);
+        return NULL;
+    }
+    // XXX note that model->residuals is just a reference
+    modelPSF->residuals = psf->residuals;
 
     return (modelPSF);
@@ -248 +255 @@
     pol.e2 = fittedPar[PM_PAR_E2];
 
-    psEllipseAxes axes = psEllipsePolToAxes (pol, minMinorAxis);
+    psEllipseAxes axes;
+    if (psEllipsePolToAxes (pol, minMinorAxis, &axes) != PS_ERR_NONE) {
+        psError(PM_ERR_PSF, false, "Failed to convert e[012] (%g,%g,%g) to axes",
+                pol.e0, pol.e1, pol.e2);
+        return false;
+    }
     psEllipseShape shape = psEllipseAxesToShape (axes);
 
@@ -519 +531 @@
     // no need to interpolate since we have forced the object center
     // to 0.5, 0.5 above
+    bool completeGrowthCurve = true;            // do we have a complete curve of growth?
     for (int i = 0; i < psf->growth->radius->n; i++) {
 
@@ -525 +538 @@
         // mask the given aperture and measure the apMag
         psImageKeepCircle (mask, xc, yc, radius, "OR", PM_MASK_MARK);
-        pmSourcePhotometryAper (&apMag, model, image, mask);
+        if (!pmSourcePhotometryAper (&apMag, model, image, mask)) {
+            psError(PM_ERR_PHOTOM, false, "Measuring apMag for radius == %g", radius);
+            completeGrowthCurve = false;
+            break;
+        }
 
         // XXX since we re-mask on each pass, this could be dropped.
@@ -538 +555 @@
     }
 
-    psf->growth->apRef = psVectorInterpolate (psf->growth->radius, psf->growth->apMag, psf->growth->refRadius);
-    psf->growth->apLoss = psf->growth->fitMag - psf->growth->apRef;
+    if (completeGrowthCurve) {
+        psf->growth->apRef = psVectorInterpolate (psf->growth->radius, psf->growth->apMag, psf->growth->refRadius);
+        psf->growth->apLoss = psf->growth->fitMag - psf->growth->apRef;
+    } else {
+        psf->growth->apRef = NAN;
+        psf->growth->apLoss = 0;
+    }
 
     psLogMsg ("psphot.growth", 4, "GrowthCurve : apLoss : %f\n", psf->growth->apLoss);
@@ -549 +571 @@
     psFree (modelRef);
 
-    return true;
-}
+    return completeGrowthCurve ? true : false;
+}