Changeset 30721

Timestamp:

Feb 21, 2011, 9:02:51 AM (15 years ago)

Author:

eugene

Message:

add tests for radial apertures

File:

• branches/eam_branches/ipp-20110213/psphot/test/tap_psphot_stackphot.pro (modified) (1 diff)

branches/eam_branches/ipp-20110213/psphot/test/tap_psphot_stackphot.pro

@@ -763 +763 @@
 end
 
+macro ckradialflux
+ data $1
+ read -fits Chip.xrad X_APER Y_APER PSF_FWHM APER_FLUX APER_FLUX_ERR APER_FLUX_STDEV APER_FILL
+ read -fits Chip.psf X_PSF Y_PSF PSF_INST_MAG
+ set mag = PSF_INST_MAG
+ # XXX include nradii in header
+ for i 0 11
+   set mA$i = -2.5*log(APER_FLUX:$i)
+   set dM$i = mA$i - PSF_INST_MAG
+   vstat dM$i
+   $DM$i = $MEDIAN
+ end
+ for i 1 11
+   $j = $i - 1
+   echo {$DM$i - $DM$j}
+ end
+
+ lim mag -2 2; clear; box
+ for i 0 11
+   plot mag dM$i -c black
+ end
+end
+
 if ($SCRIPT)
   fulltest 4