Changeset 32616 for branches/eam_branches/ipp-20110906/psphot/src/psphotStackReadout.c

Timestamp:

Nov 5, 2011, 2:30:17 PM (15 years ago)

Author:

eugene

Message:

speed up psphotRadialApertures (mostly by avoiding the large mask reallocation); add profile wings to psphotStackReadout

File:

• branches/eam_branches/ipp-20110906/psphot/src/psphotStackReadout.c (modified) (4 diffs)

branches/eam_branches/ipp-20110906/psphot/src/psphotStackReadout.c

@@ -168 +168 @@
     psphotFitSourcesLinear (config, view, STACK_SRC, false);
     psphotStackVisualFilerule(config, view, STACK_SRC);
+
+    // measure the radial profiles to the sky
+    psphotRadialProfileWings (config, view, STACK_SRC);
 
     // re-measure the kron mags with models subtracted.  this pass starts with a circular
@@ -245 +248 @@
 pass1finish:
 
+    // measure the radial profiles to the sky (only measures new objects)
+    psphotRadialProfileWings (config, view, STACK_SRC);
+
     // re-measure the kron mags with models subtracted
     // psphotKronMasked(config, view, STACK_SRC);
@@ -254 +260 @@
 
     psMemDump("psfstats");
+
+    // XXX *** NOTE : this should happen above before the last psphotFitSourcesLinear in order
+    // to generate psf forced photometry for the undetected sources
 
     // generate the objects (objects unify the sources from the different images)
@@ -281 +290 @@
 
     // measure circular, radial apertures (objects sorted by S/N)
-    // XXX can we just use psphotRadialApertures
-    // XXX make sure the headers are consistent with this (which PSF convolutions, ie mark 'none')
-    // psphotRadialAperturesByObject (config, objectsRadial, view, STACK_SRC, nMatchedPSF); 
+    // this forces photometry on the undetected sources from other images
     psphotRadialApertures (config, view, STACK_SRC, 0); // XXX entry 0 == unmatched?
     psMemDump("extmeas");