Changeset 4115 for trunk/psphot/src/apply_psf_model.c

Timestamp:

Jun 4, 2005, 5:22:13 PM (21 years ago)

Author:

eugene

Message:

running sources in brightness order, using a mask

File:

• trunk/psphot/src/apply_psf_model.c (modified) (3 diffs)

trunk/psphot/src/apply_psf_model.c

@@ -22 +22 @@
     keep           = psRegionForImage (keep, image, keep);
 
-    // set the object surfact-brightness limit for fitted pixels
+    float shapeNsigma = psMetadataLookupF32 (&status, config, "PSF_SHAPE_NSIGMA");
+    float snFaint     = psMetadataLookupF32 (&status, config, "FAINT_SN_LIM");
+
+    // set the object surface-brightness limit for fitted pixels
     float FLUX_LIMIT  = FIT_NSIGMA * sky->sampleStdev;
     psLogMsg ("psphot.apply_psf_model", 3, "fitting pixels with at least %f object counts\n", FLUX_LIMIT);
@@ -31 +34 @@
     for (int i = 0; i < sources->n; i++) {
         psSource *source = sources->data[i];
+
+        // skip the existing PSF stars 
+        // XXX drop this -- refit?
         if (source->modelPSF != NULL) continue;
 
@@ -55 +61 @@
         pmSourceMaskSaturated (source, SATURATE);
 
-        // fit as PSF, not FLT (drop poor fits)
+        // fit as PSF, not FLT (skip poor fits)
         if (!pmSourceFitModel (source, model, true)) continue;
         source->modelPSF = model;
         Niter += model[0].nIter;
         Nfit ++;
+
+        mark_psf_source (source);
+        subtract_psf_source (source);
     }
     psLogMsg ("psphot", 3, "fit PSF models: %f sec for %d objects (%d total iterations)\n", psTimerMark ("psphot"), Nfit, Niter);