Changeset 11177

Timestamp:

Jan 18, 2007, 6:56:45 PM (19 years ago)

Author:

eugene

Message:

allow faint sources to skip the moments

File:

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

trunk/psphot/src/psphotSourceStats.c

@@ -13 +13 @@
     float OUTER    = psMetadataLookupF32 (&status, recipe, "SKY_OUTER_RADIUS");
     float RADIUS   = psMetadataLookupF32 (&status, recipe, "PSF_MOMENTS_RADIUS");
+    float MIN_SN   = psMetadataLookupF32 (&status, recipe, "MOMENTS_SN_MIN");
     char *breakPt  = psMetadataLookupStr (&status, recipe, "BREAK_POINT");
 
     sources = psArrayAllocEmpty (peaks->n);
 
+    int Nfail = 0;
+    int Nmoments = 0;
     for (int i = 0; i < peaks->n; i++) {
 
@@ -31 +34 @@
         }
 
-        // XXX skip faint sources?
+        // skip faint sources
+        if (source->peak->SN < MIN_SN) {
+            psArrayAdd (sources, 100, source);
+            psFree (source);
+            continue;
+        }
 
         // measure a local sky value
@@ -38 +46 @@
         if (!status) {
           psFree (source);
+          Nfail ++;
           continue;
         }
@@ -46 +55 @@
         if (!status) {
           psFree (source);
+          Nfail ++;
           continue;
         }
@@ -55 +65 @@
             psArrayAdd (sources, 100, source);
             psFree (source);
+            Nmoments ++;
             continue;
         }
@@ -67 +78 @@
             psArrayAdd (sources, 100, source);
             psFree (source);
+            Nmoments ++;
             continue;
         }
 
         psFree (source);
+        Nfail ++;
         continue;
     }
 
-    psLogMsg ("psphot", PS_LOG_INFO, "%ld moments: %f sec\n", sources->n, psTimerMark ("psphot"));
+    psLogMsg ("psphot", PS_LOG_INFO, "%ld sources, %d moments, %d failed: %f sec\n", sources->n, Nmoments, Nfail, psTimerMark ("psphot"));
 
     return (sources);