Index: branches/eam_branches/20090715/psphot/src/psphotReadout.c
===================================================================
--- branches/eam_branches/20090715/psphot/src/psphotReadout.c	(revision 25612)
+++ branches/eam_branches/20090715/psphot/src/psphotReadout.c	(revision 25613)
@@ -124,5 +124,4 @@
         return psphotReadoutCleanup (config, readout, recipe, detections, psf, sources);
     }
-
     psphotVisualShowPSFModel (readout, psf);
 
@@ -130,11 +129,8 @@
     psphotLoadExtSources (config, view, sources);
 
-    // construct an initial model for each object
+    // construct an initial model for each object, set the radius to fitRadius, set circular fit mask
     psphotGuessModels (config, readout, sources, psf);
 
-    // XXX test output of models
-    // psphotTestSourceOutput (readout, sources, recipe, psf);
-
-    // linear PSF fit to source peaks
+    // linear PSF fit to source peaks, subtract the models from the image (in PSF mask)
     psphotFitSourcesLinear (readout, sources, recipe, psf, FALSE);
 
@@ -151,4 +147,5 @@
 
     // non-linear PSF and EXT fit to brighter sources
+    // replace model flux, adjust mask as needed, fit, subtract the models (full stamp)
     psphotBlendFit (config, readout, sources, psf);
 
@@ -156,5 +153,5 @@
     psphotReplaceAllSources (sources, recipe);
 
-    // linear fit to include all sources
+    // linear fit to include all sources (subtract again)
     psphotFitSourcesLinear (readout, sources, recipe, psf, TRUE);
 
@@ -163,10 +160,5 @@
         goto pass1finish;
     }
-
-    // XXX for the moment, drop the re-calc of the background (prove this works)
-    // replace background in residual image
-    // psphotSkyReplace (config, view);
-    // re-measure background model (median, smoothed image)
-    // psphotImageMedian (config, view);
+    // NOTE: possibly re-measure background model here with objects subtracted
 
     // add noise for subtracted objects
@@ -188,5 +180,5 @@
     }
 
-    // create full input models
+    // create full input models, set the radius to fitRadius, set circular fit mask
     psphotGuessModels (config, readout, newSources, psf);
 
@@ -212,7 +204,4 @@
 finish:
 
-    // XXX drop this : test dev version
-    // psphotPetrosianAnalysis (readout, sources, recipe);
-
     // plot positive sources
     // psphotSourcePlots (readout, sources, recipe);