Changeset 29311

Timestamp:

Oct 4, 2010, 2:47:45 PM (16 years ago)

Author:

eugene

Message:

start to add minor axis test

File:

• branches/eam_branches/ipp-20100823/psphot/src/psphotSourceSize.c (modified) (3 diffs)

branches/eam_branches/ipp-20100823/psphot/src/psphotSourceSize.c

@@ -333 +333 @@
         psF32 Mxx = source->moments->Mxx;
         psF32 Myy = source->moments->Myy;
+        psF32 Mxy = source->moments->Mxy;
 
         // replace object in image
@@ -365 +366 @@
         float nSigmaMXX = (Mxx - psfClump->X) / hypot(psfClump->dX, psfClump->X*psfClump->X*source->errMag);
         float nSigmaMYY = (Myy - psfClump->Y) / hypot(psfClump->dY, psfClump->Y*psfClump->Y*source->errMag);
+
+        // XXX should I change psfClump to use minor, major?
 
         // fprintf (stderr, "%f %f : Mxx: %f, Myy: %f, dx: %f, dy: %f, psfMag: %f, apMag: %f, dMag: %f, errMag: %f, nSigmaMag: %f, nSigmaMxx: %f, nSigmaMyy: %f\n", 
@@ -407 +410 @@
         // XXX only accept brightish detections as CRs
         // (nSigmaMAG < -options->nSigmaApResid) ||
+
+        // ** CRs are reliably flagged by a combination on Mminor < X && mag (or flux) > Y
+
+        Mminor = 0.5*(Mxx + Myy) - 0.5*sqrt(PS_SQR(Mxx - Myy) + 4.0*PS_SQR(Mxy));
+
+        // XXX do I need to find the Mminor, Mmajor distribution?
+
         bool isCR = isCR = (source->errMag < 1.0 / SIZE_SN_LIM) && ((Mxx < minMxx) || (Myy < minMyy));
         if (isCR) {