Index: trunk/psphot/src/psphotRadialAperturesByObject.c
===================================================================
--- trunk/psphot/src/psphotRadialAperturesByObject.c	(revision 30624)
+++ trunk/psphot/src/psphotRadialAperturesByObject.c	(revision 31154)
@@ -43,7 +43,4 @@
     assert (maskVal);
 
-    // XXX temporary user-supplied systematic sky noise measurement (derive from background model)
-    float skynoise = psMetadataLookupF32 (&status, recipe, "SKY.NOISE");
-
     // S/N limit to perform full non-linear fits
     float SN_LIM = psMetadataLookupF32 (&status, recipe, "RADIAL_APERTURES_SN_LIM");
@@ -63,5 +60,5 @@
     
     // source analysis is done in S/N order (brightest first)
-    objects = psArraySort (objects, pmPhotObjSortBySN);
+    objects = psArraySort (objects, pmPhotObjSortByFlux);
 
     // generate look-up arrays for readouts
@@ -109,5 +106,5 @@
 	    // limit selection to some SN limit
 	    assert (source->peak); // how can a source not have a peak?
-	    if (source->peak->SN < SN_LIM) continue;
+	    if (sqrt(source->peak->detValue) < SN_LIM) continue;
 
 	    int index = source->imageID;
@@ -147,14 +144,13 @@
 
 	    // force source image to be a bit larger...
-	    // float radius = source->peak->xf - source->pixels->col0;
-	    // radius = PS_MAX (radius, source->peak->yf - source->pixels->row0);
-	    // radius = PS_MAX (radius, source->pixels->numRows - source->peak->yf + source->pixels->row0);
-	    // radius = PS_MAX (radius, source->pixels->numCols - source->peak->xf + source->pixels->col0);
 	    pmSourceRedefinePixels (source, readout, source->peak->xf, source->peak->yf, outerRadius + 2);
 
-	    if (!psphotRadialApertureSource (source, recipe, skynoise, maskVal, radMax, nMatchedPSF)) {
+	    if (!psphotRadialApertureSource (source, recipe, maskVal, radMax, nMatchedPSF)) {
 		psTrace ("psphot", 5, "failed to extract radial profile for source at %7.1f, %7.1f", source->moments->Mx, source->moments->My);
 	    } else {
 		source->mode |= PM_SOURCE_MODE_RADIAL_FLUX;
+		if (source->parent) {
+		    source->parent->mode |= PM_SOURCE_MODE_RADIAL_FLUX;
+		}
 	    }