Changeset 19908

Timestamp:

Oct 6, 2008, 3:08:55 AM (18 years ago)

Author:

eugene

Message:

added test function to count psf stars

File:

• trunk/psphot/src/psphotChoosePSF.c (modified) (6 diffs)

trunk/psphot/src/psphotChoosePSF.c

@@ -1 +1 @@
 # include "psphotInternal.h"
+
+void psphotCountPSFStars (psArray *sources) {
+
+    int nPSF = 0;
+
+    // select the candidate PSF stars (pointers to original sources)
+    for (int i = 0; i < sources->n; i++) {
+        pmSource *source = sources->data[i];
+        if (source->mode & PM_SOURCE_MODE_PSFSTAR) {
+            nPSF ++;
+        }
+    }
+
+    fprintf (stderr, "N PSF: %d\n", nPSF);
+}
 
 // try PSF models and select best option
@@ -87 +102 @@
 
     psArray *stars = psArrayAllocEmpty (sources->n);
+
+    // psphotCountPSFStars (sources);
 
     // select the candidate PSF stars (pointers to original sources)
@@ -100 +117 @@
         }
     }
+
+    // psphotCountPSFStars (sources);
+
     // check that the identified psf stars sufficiently cover the region; if not, extend the
     // limits somewhat
     psphotCheckStarDistribution (stars, sources, options);
+
+    // psphotCountPSFStars (sources);
 
     psLogMsg ("psphot.pspsf", PS_LOG_DETAIL, "selected candidate %ld PSF objects\n", stars->n);
@@ -295 +317 @@
     psMetadataAddF32 (recipe, PS_LIST_TAIL, "DSY_STDV", PS_META_REPLACE, "stdev of mean-corrected dSYY", nSy);
 
+    // psphotCountPSFStars (sources);
+
     // unset the PSFSTAR flag for stars not used for PSF model
     for (int i = 0; i < try->sources->n; i++) {
@@ -302 +326 @@
         }
     }
+
+    // psphotCountPSFStars (sources);
 
     // build a PSF residual image
@@ -402 +428 @@
     }
 
+    // psphotCountPSFStars (sources);
+
     char *modelName = pmModelClassGetName (psf->type);
     psLogMsg ("psphot.pspsf", PS_LOG_INFO, "select psf model: %f sec\n", psTimerMark ("psphot"));