Changeset 42088 for trunk/psphot/src/psphotSourceFits.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/psphotSourceFits.c (modified) (2 diffs)

trunk/psphot/src/psphotSourceFits.c

@@ -499 +499 @@
 
 # if (PS_TRACE_ON) 
-    if (psTraceGetLevel ("psphot") >= 6) {
+    if ( (source->peak->xf > 5500) &&
+         (source->peak->yf < 500) &&
+         psTraceGetLevel ("psphot") >= 6) {
 
       // Moments-based shapes parameters
@@ -546 +548 @@
 
 # if (PS_TRACE_ON) 
-    if (psTraceGetLevel ("psphot") >= 5) {
-      if (psTraceGetLevel ("psphot") >= 6) {
-        fprintf (stderr, "chisq: %f, nIter: %d, radius: %f, npix: %d\n", model->chisqNorm, model->nIter, model->fitRadius, model->nPix);
-      }
+    if ((source->peak->xf > 5500) &&
+        (source->peak->yf < 500) &&
+        psTraceGetLevel ("psphot") >= 6) {
+      fprintf (stderr, "chisq: %f, nIter: %d, radius: %f, npix: %d\n", model->chisqNorm, model->nIter, model->fitRadius, model->nPix);
       fprintf (stderr, "--- fitted values ---\n");
       for (int i = 0; i < model->params->n; i++) {
         fprintf (stderr, "par %d: %f\n", i, model->params->data.F32[i]);
       }
-      psTraceSetLevel("psLib.math.psMinimizeLMChi2", 0);
-    }
+    }
+    psTraceSetLevel("psLib.math.psMinimizeLMChi2", 0);
 # endif