Changeset 42293 for trunk/psastro/src/psastroConvert.c

Timestamp:

Sep 27, 2022, 5:08:10 PM (4 years ago)

Author:

tdeboer

Message:

adding CTE detection functionality and new CTE source flag

File:

• trunk/psastro/src/psastroConvert.c (modified) (5 diffs)

trunk/psastro/src/psastroConvert.c

@@ -18 +18 @@
 # define PHOT_SOURCE_MASK (PM_SOURCE_MODE_FAIL | PM_SOURCE_MODE_BLEND | PM_SOURCE_MODE_BADPSF | \
                            PM_SOURCE_MODE_DEFECT | PM_SOURCE_MODE_CR_LIMIT | PM_SOURCE_MODE_EXT_LIMIT | \
-                           PM_SOURCE_MODE_POOR) // Mask to apply to sources for rejection
+                           PM_SOURCE_MODE_POOR ) // Mask to apply to sources for rejection
 
 static psArray *chooseStars(psArray *inStars, char *listName, psArray *sources, psVector *index, int nMax, float iMagMin, float iMagMax, pmSourceMode skip);
@@ -193 +193 @@
     int nBrightSkip = 0;
     int nInfSkip = 0;
+    int nCTESkip = 0;
 
     for (int i = 0; (i < inStars->n) && (j < rawStars->n); i++) {
@@ -217 +218 @@
           continue;
         }
+        //Also kick out stars that touch the CTE region
+        if (source->mode2 & PM_SOURCE_MODE2_ON_CTE) {
+            nCTESkip ++;
+            continue;
+        }
+        
         mMin = PS_MIN (mMin, source->psfMag);
         mMax = PS_MAX (mMax, source->psfMag);
@@ -225 +232 @@
 
     psLogMsg ("psastro", 4, "loaded %ld %ssources, using %ld of %ld good sources (inst mag: %f to %f)\n", sources->n, listName, rawStars->n, inStars->n, mMin, mMax);
-    psLogMsg ("psastro", 4, "skip reasons: mode: %d, faint: %d, bright: %d, inf: %d\n", nModeSkip, nFaintSkip, nBrightSkip, nInfSkip);
+    psLogMsg ("psastro", 4, "skip reasons: mode: %d, faint: %d, bright: %d, inf: %d, CTE: %d\n", nModeSkip, nFaintSkip, nBrightSkip, nInfSkip,nCTESkip);
 
     return rawStars;
@@ -509 +516 @@
   float myPltScale = fabs(pltScale[chipID]);
 
+  psLogMsg ("psastro.correctKH", PS_LOG_INFO, "applying KH correction to %s (%d)\n", chipName, chipID);
+
   // apply the correction to the detections
   for (int i = 0; i < inStars->n; i++) {