Changeset 17118

Timestamp:

Mar 24, 2008, 6:17:54 PM (18 years ago)

Author:

eugene

Message:

working on better ext detections

File:

• branches/eam_branch_20080324/psphot/src/psphotSourceSize.c (modified) (1 diff)

branches/eam_branch_20080324/psphot/src/psphotSourceSize.c

@@ -47 +47 @@
         // XXX for now, just skip any sources with masked pixels
         source->crNsigma  = -5.0;
+        bool keep = true;
+        for (int iy = -1; (iy <= +1) && keep; iy++) {
+            for (int ix = -1; (ix <= +1) && keep; ix++) {
+                if (mask[yPeak+iy][xPeak+ix]) { keep &= false; }
+            }
+        }
+        if (!keep) continue;
+
+        // need a reasonably fast way to follow the psf R = 1 annulus
+
+        // save the PSF model from the Ensemble fit
+        pmModel *PSF = source->modelPSF;
+
+        // convert model to polar terms
+        psEllipseShape shape;
+        psF32 *PAR = PSF->params->data.F32;
+        shape.sx   = PAR[PM_PAR_SXX] / M_SQRT2;
+        shape.sy   = PAR[PM_PAR_SYY] / M_SQRT2;
+        shape.sxy  = PAR[PM_PAR_SXY];
+
+        psEllipseShapeToAxes (shape, 20.0);
+        
+
+        for (theta = 0.0; theta < 360.0; theta += 10.0) {
+          
+
+
+        // XXX for now, just skip any masked pixels
+        extSum = 0.0;
         bool keep = true;
         for (int iy = -1; (iy <= +1) && keep; iy++) {