Changeset 13225 for trunk/psphot/src/psphotChoosePSF.c

Timestamp:

May 3, 2007, 4:51:40 PM (19 years ago)

Author:

eugene

Message:

load PSF from external file, load sources file external file, move default values to master recipe

File:

• trunk/psphot/src/psphotChoosePSF.c (modified) (3 diffs)

trunk/psphot/src/psphotChoosePSF.c

@@ -12 +12 @@
     psTimerStart ("psphot");
 
-    // check if a PSF model is supplied by the user
-    psf = psMetadataLookupPtr (NULL, readout->analysis, "PSPHOT.PSF");
-    if (psf != NULL) {
-        if (!psphotPSFstats (readout, recipe, psf)) psAbort("cannot measure PSF shape terms");
-        return psf;
-    }
-
     // examine PSF sources in S/N order (brightest first)
     sources = psArraySort (sources, pmSourceSortBySN);
@@ -24 +17 @@
     // array to store candidate PSF stars
     int NSTARS = psMetadataLookupS32 (&status, recipe, "PSF_MAX_NSTARS");
-    if (!status) {
-        NSTARS = PS_MIN (sources->n, 200);
-        psWarning("PSF_MAX_NSTARS is not set in the recipe --- defaulting to %d\n", NSTARS);
-    }
+    PS_ASSERT (status, NULL);
 
     float PSF_MIN_DS = psMetadataLookupF32 (&status, recipe, "PSF_MIN_DS");
-    if (!status) {
-        PSF_MIN_DS = 0.01;
-        psWarning("PSF_MIN_DS is not set --- defaulting to %f\n", PSF_MIN_DS);
-    }
+    PS_ASSERT (status, NULL);
 
     // supply the measured sky variance for optional constant errors (non-poissonian)
     float SKY_SIG = psMetadataLookupF32 (&status, recipe, "SKY_SIG");
-    if (!status) {
-        SKY_SIG = 1.0;
-        psWarning("SKY_SIG is not set --- defaulting to %f\n", SKY_SIG);
-    }
+    PS_ASSERT (status, NULL);
+
     // use poissonian errors or local-sky errors
     bool POISSON_ERRORS = psMetadataLookupBool (&status, recipe, "POISSON_ERRORS");
-    if (!status) {
-        POISSON_ERRORS = true;
-        psWarning("POISSON_ERRORS is not set in the recipe --- defaulting to true.\n");
-    }
-    pmSourceFitModelInit (15, 0.01, PS_SQR(SKY_SIG), POISSON_ERRORS);
+    PS_ASSERT (status, NULL);
 
     // use poissonian errors or local-sky errors
     bool PSF_PARAM_WEIGHTS = psMetadataLookupBool (&status, recipe, "PSF_PARAM_WEIGHTS");
-    if (!status) {
-        PSF_PARAM_WEIGHTS = false;
-        psWarning("PSF_PARAM_WEIGHTS is not set in the recipe --- defaulting to false.\n");
-    }
+    PS_ASSERT (status, NULL);
 
     // how to model the PSF variations across the field
-    // XXX make a default value?  or not?
     psMetadata *md = psMetadataLookupMetadata (&status, recipe, "PSF.TREND.MASK");
-    psPolynomial2D *psfTrendMask;
-    if (!status || !md) {
-        psWarning("PSF.TREND.MASK is not set in the recipe --- defaulting to use zeroth order.\n");
-        psfTrendMask = psPolynomial2DAlloc(PS_POLYNOMIAL_ORD, 0, 0);
-    } else {
-        psfTrendMask = psPolynomial2DfromMetadata (md);
-        if (!psfTrendMask) {
-            psError(PSPHOT_ERR_PSF, true, "Unable to construct polynomial from PSF.TREND.MASK in the recipe");
-            return NULL;
-        }
+    PS_ASSERT (status, NULL);
+
+    // get the fixed PSF fit radius
+    // XXX check that PSF_FIT_RADIUS < SKY_OUTER_RADIUS
+    float RADIUS = psMetadataLookupF32 (&status, recipe, "PSF_FIT_RADIUS");
+    PS_ASSERT (status, NULL);
+
+    pmSourceFitModelInit (15, 0.01, PS_SQR(SKY_SIG), POISSON_ERRORS);
+
+    psPolynomial2D *psfTrendMask = psPolynomial2DfromMetadata (md);
+    if (!psfTrendMask) {
+        psError(PSPHOT_ERR_PSF, true, "Unable to construct polynomial from PSF.TREND.MASK in the recipe");
+        return NULL;
     }
 
@@ -90 +70 @@
         psError(PSPHOT_ERR_PSF, true, "Failed to find any PSF candidates");
         return NULL;
-    }
-
-    // get the fixed PSF fit radius
-    // XXX EAM : check that PSF_FIT_RADIUS < SKY_OUTER_RADIUS
-    float RADIUS = psMetadataLookupF32 (&status, recipe, "PSF_FIT_RADIUS");
-    if (!status) {
-        RADIUS = 20.0;
-        psWarning("PSF_FIT_RADIUS is not set in the recipe --- defaulting to %g\n", RADIUS);
     }