Changeset 16800 for trunk/psastro/src/psastroConvert.c

Timestamp:

Mar 4, 2008, 12:24:29 PM (18 years ago)

Author:

eugene

Message:

report number and reason for skipped detections

File:

• trunk/psastro/src/psastroConvert.c (modified) (3 diffs)

trunk/psastro/src/psastroConvert.c

@@ -79 +79 @@
     float mMin = +100.0;
     float mMax = -100.0;
+    int nModeSkip = 0;
+    int nFaintSkip = 0;
+    int nBrightSkip = 0;
+    int nInfSkip = 0;
 
     for (int i = 0; (i < inStars->n) && (j < rawStars->n); i++) {
@@ -84 +88 @@
         pmSource *source = sources->data[n];
         // XXX this needs to be flexible
-        if (source->mode & skip) continue;
-        if ((iMagMax != 0.0) && (source->psfMag > iMagMax)) continue;
-        if ((iMagMin != 0.0) && (source->psfMag < iMagMin)) continue;
-        if (!isfinite(source->psfMag)) continue;
+        if (source->mode & skip) { 
+          nModeSkip ++;
+          continue;
+        }
+        if ((iMagMax != 0.0) && (source->psfMag > iMagMax)) {
+          nFaintSkip ++;
+          continue;
+        }
+        if ((iMagMin != 0.0) && (source->psfMag < iMagMin)) {
+          nBrightSkip ++;
+          continue;
+        }
+        if (!isfinite(source->psfMag)) {
+          nInfSkip ++;
+          continue;
+        }
         mMin = PS_MIN (mMin, source->psfMag);
         mMax = PS_MAX (mMax, source->psfMag);
@@ -97 +113 @@
     psMetadataAdd (readout->analysis, PS_LIST_TAIL, "PSASTRO.RAWSTARS", PS_DATA_ARRAY, "astrometry objects", rawStars);
     psLogMsg ("psastro", 4, "loaded %ld sources, using %ld of %ld good sources (inst mag: %f to %f)\n", sources->n, rawStars->n, inStars->n, mMin, mMax);
+    psLogMsg ("psastro", 4, "skip reasons: mode: %d, faint: %d, bright: %d, inf: %d\n", nModeSkip, nFaintSkip, nBrightSkip, nInfSkip);
 
     psFree (index);