Changeset 4574 for trunk/psphot/src/mark_psf_source.c

Timestamp:

Jul 16, 2005, 6:02:29 AM (21 years ago)

Author:

eugene

Message:

improved galaxy model guesses, better handling of errors

File:

• trunk/psphot/src/mark_psf_source.c (modified) (3 diffs)

trunk/psphot/src/mark_psf_source.c

@@ -15 +15 @@
 // PS_SOURCE_BRIGHTSTAR 
 # define MIN_DS 0.01
-bool mark_psf_source (psSource *source, float shapeNsigma)
+bool mark_psf_source (psSource *source, float shapeNsigma, float SATURATE)
 { 
     float dSX, dSY, SX, SY, SN;
@@ -22 +22 @@
     if (source->modelPSF == NULL) return (false);
 
-    // SATSTAR objects are fitted with PSF regardless
-    if (source->type == PS_SOURCE_SATSTAR) return (true);
+    // if object has fitted peak above saturation, label as SATSTAR
+    // remember: fit does not use saturated pixels (masked)
+    if (source->modelPSF->params->data.F32[1] >= SATURATE) {
+        if (source->type == PS_SOURCE_PSFSTAR) {
+            psLogMsg ("psphot", 3, "PSFSTAR marked saturated\n");
+        }
+        source->type = PS_SOURCE_SATSTAR;
+        return (true);
+    } 
+    if (source->type == PS_SOURCE_SATSTAR) {
+        psLogMsg ("psphot", 4, "SATSTAR marked bright (fitted peak below saturation)\n");
+        source->type = PS_SOURCE_BRIGHTSTAR;
+    }
 
     SN  = source->modelPSF->params->data.F32[1]/source->modelPSF->dparams->data.F32[1];
@@ -35 +46 @@
     // sigma = 1 / SX
     // dsx = 1 / (SX * SN)
-    // dsx_o = hypot (1/SX*SN + MIN_DSX)
+    // dsx_o = hypot (1/(SX*SN), MIN_DSX)
 
     // assign PS_SOURCE_BRIGHTSTAR to bright objects within PSF region of dparams[]
     if ((fabs(nSx) < shapeNsigma) && (fabs(nSy) < shapeNsigma)) {
-      if (source->type == PS_SOURCE_PSFSTAR) return (true);
-      source->type = PS_SOURCE_BRIGHTSTAR;
-      return (true);
+        if (source->type == PS_SOURCE_PSFSTAR) return (true);
+        source->type = PS_SOURCE_BRIGHTSTAR;
+        return (true);
     }
     
     if (source->type == PS_SOURCE_PSFSTAR) {
-      psLogMsg ("psphot", 3, "PSFSTAR demoted based on dSx, dSy\n");
+        psLogMsg ("psphot", 3, "PSFSTAR demoted based on dSx, dSy\n");
     }
 
     if ((nSx >= shapeNsigma) || (nSy >= shapeNsigma)) {
-      source->type = PS_SOURCE_GALAXY;
-      return (false);
+        source->type = PS_SOURCE_GALAXY;
+        return (false);
     }
     // replace DEFECT with COSMIC?
     if ((nSx <= -shapeNsigma) || (nSy <= -shapeNsigma)) {
-      source->type = PS_SOURCE_DEFECT;
-      return (false);
+        source->type = PS_SOURCE_DEFECT;
+        return (false);
     }
     psTrace (".psphot.mark_psf_source", 2, "unexpected result: unmarked object\n");