Changeset 42022 for branches/eam_branches/ipp-20211108/psphot/src/psphotFitSourcesLinear.c

Timestamp:

Feb 7, 2022, 11:36:08 AM (4 years ago)

Author:

eugene

Message:

add count of where linear-fit sources are skipped; in deteff code, model radius is too small: many faint sources are rejected because the model does not include enough of the profile: increase by a factor of 5

File:

• branches/eam_branches/ipp-20211108/psphot/src/psphotFitSourcesLinear.c (modified) (7 diffs)

branches/eam_branches/ipp-20211108/psphot/src/psphotFitSourcesLinear.c

@@ -217 +217 @@
     int Nsat = 0;
 
+    // track number kept at several stages of the analysis
+    int Nstep_0 = 0;
+    int Nstep_1 = 0;
+    int Nstep_2 = 0;
+    int Nstep_3 = 0;
+    int Nstep_4 = 0;
+    int Nstep_5 = 0;
+
     // select the sources which will be used for the fitting analysis
     for (int i = 0; i < sources->n; i++) {
@@ -246 +254 @@
         if (source->mode & PM_SOURCE_MODE_MOMENTS_FAILURE) continue;
 
+        Nstep_0++;
+
         // XXX count saturated stars
         if (source->mode & PM_SOURCE_MODE_SATSTAR) {
@@ -255 +265 @@
             if (!pmSourceCacheModel (source, maskVal)) continue;
         }
+        Nstep_1++;
 
         // save the original coords
@@ -265 +276 @@
         if (x > AnalysisRegion.x1) continue;
         if (y > AnalysisRegion.y1) continue;
+        Nstep_2++;
 
         // check the integral of the model : is it large enough?
@@ -283 +295 @@
             continue;
         }
+        Nstep_3++;
+
         bool isPSF = false;
         pmModel *model = pmSourceGetModel (&isPSF, source);
@@ -306 +320 @@
         } 
         if (modelSum  < cutModelSum * normFlux) continue;
+        Nstep_4++;
+
         if (maskedSum < cutMaskedSum * normFlux) continue;
+        Nstep_5++;
 
         // clear the 'mark' pixels and remask on the fit aperture 
@@ -322 +339 @@
     }
     psLogMsg ("psphot.ensemble", PS_LOG_MINUTIA, "built fitSources: %f sec (%ld objects)\n", psTimerMark ("psphot.linear"), sources->n);
+    psLogMsg ("psphot.ensemble", PS_LOG_INFO, "Number of sources kept at each stage: %d %d %d %d %d %d\n", Nstep_0, Nstep_1, Nstep_2, Nstep_3, Nstep_4, Nstep_5);
 
     if (fitSources->n == 0) {