Changeset 25261 for branches/pap/psphot/src/psphotFake.c

Timestamp:

Sep 2, 2009, 5:00:06 PM (17 years ago)

Author:

Paul Price

Message:

FWHM_X and FWHM_Y don't appear to get set any more. Instead, FWHM_MAJ and FWHM_MIN do. Does this explain the much-reduced limiting magnitudes from the difference images??? Either way, it's required to get the detection efficiency measurements working.

File:

• branches/pap/psphot/src/psphotFake.c (modified) (4 diffs)

branches/pap/psphot/src/psphotFake.c

@@ -17 +17 @@
                       const pmReadout *ro,     // Readout of interest
                       float thresh,            // Threshold for source identification
-                      float xFWHM, float yFWHM, // Size of PSF
+                      float fwhmMajor, float fwhmMinor, // Size of PSF
                       float smoothNsigma,       // Smoothing limit
                       psImageMaskType maskVal   // Value to mask
@@ -25 +25 @@
     PM_ASSERT_READOUT_VARIANCE(ro, false);
 
-    float smoothSigma  = 0.5*(xFWHM + yFWHM) / (2.0*sqrtf(2.0*log(2.0)));
+    float smoothSigma  = 0.5*(fwhmMajor + fwhmMinor) / (2.0*sqrtf(2.0*log(2.0)));
     psKernel *kernel = psImageSmoothKernel(smoothSigma, smoothNsigma); // Kernel used for smoothing
     psKernel *newCovar = psImageCovarianceCalculate(kernel, ro->covariance); // Covariance matrix
@@ -131 +131 @@
     // Collect recipe information
     bool mdok;                                                 // Status of MD lookup
-    float xFWHM = psMetadataLookupF32(NULL, recipe, "FWHM_X"); // PSF size in x
-    float yFWHM = psMetadataLookupF32(NULL, recipe, "FWHM_Y"); // PSF size in y
-    if (!isfinite(xFWHM) || !isfinite(yFWHM)) {
-        psError(PSPHOT_ERR_CONFIG, false, "Unable to find FWHM_X and FWHM_Y in recipe");
+    float fwhmMajor = psMetadataLookupF32(NULL, recipe, "FWHM_MAJ"); // PSF size in x
+    float fwhmMinor = psMetadataLookupF32(NULL, recipe, "FWHM_MIN"); // PSF size in y
+    if (!isfinite(fwhmMajor) || !isfinite(fwhmMinor) || fwhmMajor == 0.0 || fwhmMinor == 0.0) {
+        psError(PSPHOT_ERR_CONFIG, false, "Unable to find FWHM_MAJ and FWHM_MIN in recipe");
         return false;
     }
@@ -170 +170 @@
     float magLim;                       // Guess at limiting magnitude
     float minFlux;                      // Minimum flux for fake image
-    if (!fakeLimit(&magLim, &minFlux, readout, thresh, xFWHM, yFWHM, smoothNsigma, maskVal)) {
+    if (!fakeLimit(&magLim, &minFlux, readout, thresh, fwhmMajor, fwhmMinor, smoothNsigma, maskVal)) {
         psError(PS_ERR_UNKNOWN, false, "Unable to determine limits for image");
         return false;