Changeset 42088 for trunk/psphot/src/psphotFitSourcesLinear.c

Timestamp:

Feb 28, 2022, 2:40:32 PM (4 years ago)

Author:

eugene

Message:

for bright stars in stacks, the peak in the S/N image can have a central divot if the variance of the central pixels is relatively large compared to the ring of pixels further from the center. the detected peaks are then not well centered. this can be avoided if we detect on the signal image itself (after first determining an appropriate flux value for the desired S/N limit); added an option to use the signal image for the first pass peak detetion step

Location:

trunk/psphot

Files:

• . (modified) (1 prop)
• src/psphotFitSourcesLinear.c (modified) (7 diffs)

trunk/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) {