Changeset 26253 for branches/eam_branches/20091113/psastro/src/psastroZeroPoint.c

Timestamp:

Nov 22, 2009, 2:17:42 PM (17 years ago)

Author:

eugene

Message:

slightly more careful analysis of bootstrap errors

File:

• branches/eam_branches/20091113/psastro/src/psastroZeroPoint.c (modified) (3 diffs)

branches/eam_branches/20091113/psastro/src/psastroZeroPoint.c

@@ -94 +94 @@
     if (byExposure) {
         psMetadata *header = psMetadataLookupMetadata (&status, fpa->analysis, "PSASTRO.HEADER");
+        if (!header) {
+          header = psMetadataAlloc ();
+          psMetadataAddMetadata (fpa->analysis, PS_LIST_TAIL, "PSASTRO.HEADER",  PS_META_REPLACE, "psastro header stats", header);
+          psFree (header);
+        }
         psastroZeroPointAnalysis (header, dMag, zeropt, recipe);
         psFree (dMag);
@@ -184 +189 @@
 }
 
-#define MAG_RESOLUTION 0.001
+#define MAG_RESOLUTION 0.0002
 
 // set the bin closest to the corresponding value.  if USE_END is +/- 1,
@@ -304 +309 @@
     psTrace("psastro", 6, "subset stats: Sum: %f, S2: %f, Npts: %d\n", Sum, S2, edgeSample);
 
-    stats->clippedStdev = sqrt(S2 / edgeSample - PS_SQR(Sum/edgeSample));
+    stats->clippedStdev = PS_MAX (sqrt(S2 / edgeSample - PS_SQR(Sum/edgeSample)), MAG_RESOLUTION);
     psTrace("psastro", 5, "percentile stats %f +/- %f\n", stats->clippedMean, stats->clippedStdev);