Changeset 30707 for branches/eam_branches/ipp-20110213/psphot/src/psphotRadialAperturesByObject.c

Timestamp:

Feb 19, 2011, 10:40:50 AM (15 years ago)

Author:

eugene

Message:

remove invalid references to isophot and kron from extended source analysis; use the measured sky stdev instead of the user guess; optionally allow stars as well as galaxies to get petrosians; fix target vs input seeing in psphotStack; allow radial apertures for single-image analysis as well as stack; calculate the radial aperture flux errors; calculate the petrosian parameter errors; turn on threads for ext source fits; save covar matrix for ext source fits; calculate petrosian fill factor

Location:

branches/eam_branches/ipp-20110213/psphot

Files:

• . (modified) (1 prop)
• src/psphotRadialAperturesByObject.c (modified) (2 diffs)

branches/eam_branches/ipp-20110213/psphot/src/psphotRadialAperturesByObject.c

@@ -42 +42 @@
     psImageMaskType maskVal = psMetadataLookupImageMask(&status, recipe, "MASK.PSPHOT"); // Mask value for bad pixels
     assert (maskVal);
-
-    // XXX temporary user-supplied systematic sky noise measurement (derive from background model)
-    float skynoise = psMetadataLookupF32 (&status, recipe, "SKY.NOISE");
 
     // S/N limit to perform full non-linear fits
@@ -147 +144 @@
 
             // force source image to be a bit larger...
-            // float radius = source->peak->xf - source->pixels->col0;
-            // radius = PS_MAX (radius, source->peak->yf - source->pixels->row0);
-            // radius = PS_MAX (radius, source->pixels->numRows - source->peak->yf + source->pixels->row0);
-            // radius = PS_MAX (radius, source->pixels->numCols - source->peak->xf + source->pixels->col0);
             pmSourceRedefinePixels (source, readout, source->peak->xf, source->peak->yf, outerRadius + 2);
 
-            if (!psphotRadialApertureSource (source, recipe, skynoise, maskVal, radMax, nMatchedPSF)) {
+            if (!psphotRadialApertureSource (source, recipe, maskVal, radMax, nMatchedPSF)) {
                 psTrace ("psphot", 5, "failed to extract radial profile for source at %7.1f, %7.1f", source->moments->Mx, source->moments->My);
             } else {