Index: /branches/eam_branches/ipp-20110505/psphot/src/psphotMagnitudes.c
===================================================================
--- /branches/eam_branches/ipp-20110505/psphot/src/psphotMagnitudes.c	(revision 31521)
+++ /branches/eam_branches/ipp-20110505/psphot/src/psphotMagnitudes.c	(revision 31522)
@@ -212,7 +212,5 @@
         if (status && isfinite(source->apFlux)) {
 	    Nap ++;
-	} else {
-	    fprintf (stderr, "failed to measure mag for source @ %f,%f\n", source->peak->xf, source->peak->yf);
-	}
+	} 
 
         // clear the mask bit
Index: /branches/eam_branches/ipp-20110505/psphot/src/psphotRoughClass.c
===================================================================
--- /branches/eam_branches/ipp-20110505/psphot/src/psphotRoughClass.c	(revision 31521)
+++ /branches/eam_branches/ipp-20110505/psphot/src/psphotRoughClass.c	(revision 31522)
@@ -180,6 +180,20 @@
             psAssert (status, "missing MOMENTS_SY_MAX");
         }
-
-        psfClump = pmSourcePSFClump (NULL, region, sources, PSF_SN_LIM, PSF_CLUMP_GRID_SCALE, MOMENTS_SX_MAX, MOMENTS_SY_MAX, MOMENTS_AR_MAX);
+        float MOMENTS_SX_MIN = psMetadataLookupF32(&status, analysis, "MOMENTS_SX_MIN");
+        if (!status) {
+            MOMENTS_SX_MIN = psMetadataLookupF32(&status, recipe, "MOMENTS_SX_MIN");
+            if (!status) {
+		MOMENTS_SX_MIN = 0.5;
+	    }
+        }
+        float MOMENTS_SY_MIN = psMetadataLookupF32(&status, analysis, "MOMENTS_SY_MIN");
+        if (!status) {
+            MOMENTS_SY_MIN = psMetadataLookupF32(&status, recipe, "MOMENTS_SY_MIN");
+            if (!status) {
+		MOMENTS_SY_MIN = 0.5;
+	    }
+        }
+
+        psfClump = pmSourcePSFClump (NULL, region, sources, PSF_SN_LIM, PSF_CLUMP_GRID_SCALE, MOMENTS_SX_MAX, MOMENTS_SY_MAX, MOMENTS_SX_MIN, MOMENTS_SY_MIN, MOMENTS_AR_MAX);
 
         psMetadataAddF32 (regionMD, PS_LIST_TAIL, "PSF.CLUMP.X",  PS_META_REPLACE, "psf clump center", psfClump.X);
Index: /branches/eam_branches/ipp-20110505/psphot/src/psphotSourceStats.c
===================================================================
--- /branches/eam_branches/ipp-20110505/psphot/src/psphotSourceStats.c	(revision 31521)
+++ /branches/eam_branches/ipp-20110505/psphot/src/psphotSourceStats.c	(revision 31522)
@@ -505,4 +505,13 @@
     psF32 MOMENTS_AR_MAX = psMetadataLookupF32(&status, recipe, "MOMENTS_AR_MAX"); psAssert (status, "missing MOMENTS_AR_MAX");
 
+    float MOMENTS_SX_MIN = psMetadataLookupF32(&status, recipe, "MOMENTS_SX_MIN");
+    if (!status) {
+	MOMENTS_SX_MIN = 0.5;
+    }
+    float MOMENTS_SY_MIN = psMetadataLookupF32(&status, recipe, "MOMENTS_SY_MIN");
+    if (!status) {
+	MOMENTS_SY_MIN = 0.5;
+    }
+
     // when we set the window, we are not attempting to measure spatial variations; we can use a somewhat higher S/N limit
     // since we are using all sources (true?)
@@ -544,5 +553,5 @@
 
         // determine the PSF parameters from the source moment values
-        pmPSFClump psfClump = pmSourcePSFClump (NULL, NULL, sources, PSF_SN_LIM, PSF_CLUMP_GRID_SCALE, MOMENTS_SX_MAX, MOMENTS_SY_MAX, MOMENTS_AR_MAX);
+        pmPSFClump psfClump = pmSourcePSFClump (NULL, NULL, sources, PSF_SN_LIM, PSF_CLUMP_GRID_SCALE, MOMENTS_SX_MAX, MOMENTS_SY_MAX, MOMENTS_SX_MIN, MOMENTS_SY_MIN, MOMENTS_AR_MAX);
         psLogMsg ("psphot", 3, "radius %.1f, nStars: %d of %d in clump, nSigma: %5.2f, X,  Y: %f, %f (%f, %f)\n", sigma[i], psfClump.nStars, psfClump.nTotal, psfClump.nSigma, psfClump.X, psfClump.Y, sqrt(psfClump.X) / sigma[i], sqrt(psfClump.Y) / sigma[i]);