Changeset 15400 for branches/eam_branch_20071023/psphot/src/psphotMakeResiduals.c

Timestamp:

Oct 28, 2007, 3:41:13 PM (19 years ago)

Author:

eugene

Message:

adding extended source measurements

File:

• branches/eam_branch_20071023/psphot/src/psphotMakeResiduals.c (modified) (3 diffs)

branches/eam_branch_20071023/psphot/src/psphotMakeResiduals.c

@@ -28 +28 @@
 
     float nSigma = psMetadataLookupF32(&status, recipe, "PSF.RESIDUALS.NSIGMA");
+    PS_ASSERT (status, false);
+
+    float pixelSN = psMetadataLookupF32(&status, recipe, "PSF.RESIDUALS.PIX.SN");
     PS_ASSERT (status, false);
 
@@ -195 +198 @@
                 //resid->weight->data.F32[oy][ox] = fluxStats->sampleStdev;
 
-                // if value < NRESID_SIGMA*sigma, mask pixel in resid map
-                if (resid->Ro->data.F32[oy][ox] < NRESID_SIGMA*fluxStats->sampleStdev) {
-                  resid->mask->data.U8[oy][ox] = RESID_SIG_MASK;
+                if (resid->Ro->data.F32[oy][ox] < pixelSN*fluxStats->sampleStdev) {
+                  resid->mask->data.U8[oy][ox] = 1;
                 }
 
@@ -234 +236 @@
                 resid->Ry->data.F32[oy][ox] = B->data.F64[2];
 
-                dRo = sqrt(A->data.F32[0][0]);
-                if (resid->Ro->data.F32[oy][ox] < NRESID_SIGMA*dRo) {
-                  resid->mask->data.U8[oy][ox] = RESID_SIG_MASK;
+                float dRo = sqrt(A->data.F32[0][0]);
+                if (resid->Ro->data.F32[oy][ox] < pixelSN*dRo) {
+                  resid->mask->data.U8[oy][ox] = 1;
                 }
                 //resid->weight->data.F32[oy][ox] = XXX;