Changeset 42187

Timestamp:

Apr 30, 2022, 4:31:42 PM (4 years ago)

Author:

eugene

Message:

more log info; add check for stars selected for PSF model, use S/N 0.1 to measure the psf stats if none exist (e.g., in Forced situation)

Location:

branches/eam_branches/ipp-20220316/psphot/src

Files:

• psphotChoosePSF.c (modified) (1 diff)
• psphotSkyReplace.c (modified) (1 diff)

branches/eam_branches/ipp-20220316/psphot/src/psphotChoosePSF.c

@@ -594 +594 @@
     psVector *psfExtra2 = psVectorAllocEmpty (100, PS_DATA_F32);
 
+    // count how many PSF stars are defined
+    int nPSFstars = 0; // count the number of PSF stars
     for (int i = 0; i < sources->n; i++) {
-
         pmSource *source = sources->data[i];
-
         if (!(source->mode & PM_SOURCE_MODE_PSFSTAR)) continue;
+    }
+
+    for (int i = 0; i < sources->n; i++) {
+
+        pmSource *source = sources->data[i];
+
+        // If no sources are PSFSTARs, then let's just use all with S/N > 0.1.
+        // (This can happen if all are supplied)
+        if (nPSFstars) {
+            if (!(source->mode & PM_SOURCE_MODE_PSFSTAR)) continue;
+        } else {
+            if (source->psfMagErr > 0.1) continue;
+        }
 
         float xc = source->peak->xf;

branches/eam_branches/ipp-20220316/psphot/src/psphotSkyReplace.c

@@ -4 +4 @@
 {
     bool status = true;
+
+    fprintf (stdout, "\n");
+    psLogMsg ("psphot", PS_LOG_INFO, "--- psphot Replace Sky ---");
 
     int num = psphotFileruleCount(config, filerule);