Changeset 13165 for trunk/psastro/src/psastroOneChip.c

Timestamp:

May 2, 2007, 3:00:32 PM (19 years ago)

Author:

eugene

Message:

testing alternate stats measurements

File:

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

trunk/psastro/src/psastroOneChip.c

@@ -88 +88 @@
     // XXX allow statitic to be set by the user
     psStats *fitStats = psStatsAlloc (PS_STAT_CLIPPED_MEAN | PS_STAT_CLIPPED_STDEV);
+    // psStats *fitStats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
     fitStats->clipSigma = psMetadataLookupF32 (&status, recipe, "PSASTRO.CHIP.NSIGMA");
     fitStats->clipIter = psMetadataLookupS32 (&status, recipe, "PSASTRO.CHIP.NITER");
@@ -107 +108 @@
     // pixError is the average 1D scatter in pixels ('results' are in FPA units = microns)
     float pixError = 0.5*(results->xStats->clippedStdev + results->yStats->clippedStdev) / pixelScale;
+    // float pixError = 0.5*(results->xStats->robustStdev + results->yStats->robustStdev) / pixelScale;
 
     // astError is the average 1D scatter in arcsec ('results' are in FPA units = microns)
     float astError = 0.5*(results->xStats->clippedStdev + results->yStats->clippedStdev) * plateScale;
+    // float astError = 0.5*(results->xStats->robustStdev + results->yStats->robustStdev) * plateScale;
     int astNstar = results->yStats->clippedNvalues;
 
@@ -144 +147 @@
         psastroDumpRawstars (rawstars, fpa, chip);
     }
+    if (psTraceGetLevel("psastro.dump") > 0) {
+        psastroDumpStars (refstars, "refstars.cal.dat");
+    }
 
     if (psTraceGetLevel("psastro.plot") > 0) {