Index: /trunk/psphot/src/psphotApResid.c
===================================================================
--- /trunk/psphot/src/psphotApResid.c	(revision 6917)
+++ /trunk/psphot/src/psphotApResid.c	(revision 6918)
@@ -20,5 +20,5 @@
     // 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,5 +86,5 @@
     // 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,9 +215,11 @@
 
     // 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,5 +243,6 @@
 
     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],