Changeset 27726

Timestamp:

Apr 22, 2010, 11:05:27 AM (16 years ago)

Author:

Paul Price

Message:

Based on simulation and real data, I now believe I have the correct expression for the limiting magnitude. Added a log message with the input parameters to the limiting magnitude.

File:

• branches/pap/psphot/src/psphotEfficiency.c (modified) (2 diffs)

branches/pap/psphot/src/psphotEfficiency.c

@@ -73 +73 @@
     psTrace("psphot.fake", 1, "Limiting flux: %f\n", fluxLim);
     psTrace("psphot.fake", 1, "Limiting mag: %f\n", *magLim);
+    psLogMsg("psphot", PS_LOG_INFO,
+             "Detection efficiency:\n"
+             "  Mean variance: %f * %f\n"
+             "  Threshold: %f\n"
+             "  Normalisation: %f\n"
+             "  Limiting magnitude: %f\n",
+             meanVar, *covarFactor, thresh, *norm, *magLim);
 
     *radius = smoothSigma * smoothNsigma;
@@ -166 +173 @@
     // loop over the available readouts
     for (int i = 0; i < num; i++) {
-        if (i == chisqNum) continue; // skip chisq image
-        if (!psphotEfficiencyReadout (config, view, "PSPHOT.INPUT", i, recipe)) {
+        if (i == chisqNum) continue; // skip chisq image
+        if (!psphotEfficiencyReadout (config, view, "PSPHOT.INPUT", i, recipe)) {
             psError (PSPHOT_ERR_CONFIG, false, "failed to measure detection efficiency for PSPHOT.INPUT entry %d", i);
-            return false;
-        }
+            return false;
+        }
     }
     return true;