Changeset 33643 for trunk/psastro/src/psastroZeroPoint.c

Timestamp:

Mar 30, 2012, 3:01:51 PM (14 years ago)

Author:

eugene

Message:

updates related to zero points from new reference catalogs & w-band mags

File:

• trunk/psastro/src/psastroZeroPoint.c (modified) (3 diffs)

trunk/psastro/src/psastroZeroPoint.c

@@ -189 +189 @@
       pmAstromObj *ref = refstars->data[match->ref];
 
+      // XXX I should be applying the airmass term here or in psastroLoadRefstars
+      // (ie, ref->Mag would be a prediction for a given star.
       float value = ref->Mag - raw->Mag - 2.5*log10(exptime);
       psVectorAppend (dMag, value);
@@ -214 +216 @@
 
     psStats *stats = NULL;
+    float zeroPtObs, zeroPtErr;
     if (useMean) {
         // stats structure and mask for use in measuring the clipping statistic
         // this analysis has too few data points to use the robust median method
-        stats = psStatsAlloc (PS_STAT_CLIPPED_MEAN | PS_STAT_CLIPPED_STDEV);
+        stats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
         if (!psVectorStats (stats, dMag, NULL, NULL, 0)) {
             psError(PS_ERR_UNKNOWN, false, "failure to measure zero point by mean");
             return false;
         }
+        zeroPtObs = stats->robustMedian;
+        zeroPtErr = stats->robustStdev;
     } else {
         stats = psastroStatsPercentile (dMag, recipe);
@@ -228 +233 @@
             return false;
         }
-    }
-    fprintf (stderr, "zero point %f +/- %f using %ld stars; transparency %f\n", stats->clippedMean, stats->clippedStdev, dMag->n, zeropt - stats->clippedMean);
-
-    psMetadataAddF32 (header, PS_LIST_TAIL, "ZPT_OBS", PS_META_REPLACE, "measured zero point", stats->clippedMean);
+        zeroPtObs = stats->clippedMean;
+        zeroPtErr = stats->clippedStdev;
+    }
+    fprintf (stderr, "zero point %f +/- %f using %ld stars; transparency %f\n", zeroPtObs, zeroPtErr, dMag->n, zeropt - zeroPtObs);
+
+    psMetadataAddF32 (header, PS_LIST_TAIL, "ZPT_OBS", PS_META_REPLACE, "measured zero point", zeroPtObs);
     psMetadataAddF32 (header, PS_LIST_TAIL, "ZPT_REF", PS_META_REPLACE, "reference zero point", zeropt);
-    psMetadataAddF32 (header, PS_LIST_TAIL, "ZPT_ERR", PS_META_REPLACE, "error on zero point",  stats->clippedStdev);
-    psMetadataAddF32 (header, PS_LIST_TAIL, "ZPT_OFF", PS_META_REPLACE, "zero point offset",    zeropt - stats->clippedMean);
+    psMetadataAddF32 (header, PS_LIST_TAIL, "ZPT_ERR", PS_META_REPLACE, "error on zero point",  zeroPtErr);
+    psMetadataAddF32 (header, PS_LIST_TAIL, "ZPT_OFF", PS_META_REPLACE, "zero point offset",    zeropt - zeroPtObs);
 
     psFree (stats);