Changeset 6918 for trunk/psphot/src/psphotApResid.c

Timestamp:

Apr 19, 2006, 5:03:11 PM (20 years ago)

Author:

eugene

Message:

added count of excluded sources

File:

• trunk/psphot/src/psphotApResid.c (modified) (4 diffs)

trunk/psphot/src/psphotApResid.c

@@ -20 +20 @@
     // measure the aperture loss as a function of radius for PSF
     float REF_RADIUS = psMetadataLookupF32 (&status, recipe, "PSF_REF_RADIUS");
-    psf->growth = pmGrowthCurveAlloc (3.0, 50.0, REF_RADIUS);
+    psf->growth = pmGrowthCurveAlloc (3.0, 100.0, REF_RADIUS);
     psphotGrowthCurve (readout, psf);
     
@@ -86 +86 @@
     // 3hi/1lo sigma clipping on the rflux vs metric fit
     psStats *stats = psStatsAlloc (PS_STAT_SAMPLE_MEDIAN | PS_STAT_SAMPLE_STDEV);
-    stats->min = 3.0;
+    stats->min = 2.0;
     stats->max = 3.0;
 
@@ -215 +215 @@
 
     // measure scatter for sources with dMag < 0.01 (S/N = 100)
+    int Nkeep = 0;
     psStats *residStats = psStatsAlloc (PS_STAT_SAMPLE_STDEV);
     for (int i = 0; i < dMag->n; i++) {
         if (dMag->data.F64[i] > 0.01) {
             mask->data.U8[i] |= 0x02;
-        }
+        } 
+        if (! mask->data.U8[i]) Nkeep ++;
     }
     residStats  = psVectorStats (residStats, resid, NULL, mask, 0x03);
@@ -241 +243 @@
 
     psLogMsg ("psphot.apresid", 3, "measure full-frame aperture residual: %f sec\n", psTimerMark ("psphot"));
-    psLogMsg ("psphot.apresid", 4, "aperture residual: %f +/- %f : %f bias, %f skysat\n", psf->ApResid, psf->dApResid, psf->skyBias, psf->skySat);
+    psLogMsg ("psphot.apresid", 4, "aperture residual: %f +/- %f : %f bias, %f skysat (%d of %d used)\n", 
+              psf->ApResid, psf->dApResid, psf->skyBias, psf->skySat, Nkeep, Npsf);
     psLogMsg ("psphot.apresid", 4, "apresid trends: %f %f %f %f %f\n", 
               1e3*psf->ApTrend->coeff[1][0][0][0],