Changeset 35661 for branches/eam_branches/ipp-20130509/psphot/test/tap_psphot_galaxies.pro

Timestamp:

Jun 14, 2013, 6:21:39 AM (13 years ago)

Author:

eugene

Message:

extend galaxy tests; add recipe selection option to refine the set of extended fits used; cleanup extended fits outputs and leaks

File:

• branches/eam_branches/ipp-20130509/psphot/test/tap_psphot_galaxies.pro (modified) (6 diffs)

branches/eam_branches/ipp-20130509/psphot/test/tap_psphot_galaxies.pro

@@ -75 +75 @@
 end
 
+macro go
+  mkexp test.exp 1.0 EXP
+  fitexp test.exp test.exp.fit EXP_CONV
+
+  mkexp test.ser 1.0 SERSIC
+  fitexp test.ser test.ser.fit SER_CONV
+
+  mkexp test.dev 1.0 DEV
+  fitexp test.dev test.dev.fit DEV_CONV
+
+  mkexp test.gau 1.0 GAUSS
+  fitexp test.gau test.gau.fit GAU_CONV
+
+  mkexp test.pg 1.0 PGAUSS
+  fitexp test.pg test.pg.fit PGA_CONV
+
+  mkexp test.qga 1.0 QGAUSS
+  fitexp test.qga test.qga.fit QGA_CONV
+
+  mkexp test.p1 1.0 PS1_V1
+  fitexp test.p1 test.p1.fit PS1_CONV
+end
+
+macro go.ckgalaxy
+  foreach type exp ser dev gau pg qga p1
+    ckgalaxy test.$type.dat test.$type.fit.cmf
+    wait $type
+  end
+end
+
 # create a reference database of fake stars to be used by ppSim below
 macro mkref
@@ -115 +145 @@
 end
 
-# compare two cmf files with extname Chip.psf 
-# things to compare:
-# * completeness (which sources in (1) are not detected in (2)
-# * positions (X_PSF, Y_PSF) 
-# * instrumental psf mags
-# * position errors (no input errors; use a model?)
-# * measured FWHM?
-# * kron mags (fluxes)
-# * etc, etc
+# create a realistic distribution of fake stars, GAUSS PSF
+macro fitexp
+  if ($0 != 4)
+    echo "USAGE: fitexp basename outname (fitModel)"
+    break
+  end
+
+  local basename fitModel outname
+  $basename = $1
+  $outname  = $2
+  $fitModel = $3
+
+  $psphotConfig = 
+  $psphotConfig = $psphotConfig -recipe PSPHOT GALAXY_TEST
+  $psphotConfig = $psphotConfig -threads 4
+  $psphotConfig = $psphotConfig -Db PSPHOT:LMM_FIT_CHISQ_CONVERGENCE F
+  $psphotConfig = $psphotConfig -Df PSPHOT:EXT_FIT_MIN_TOL 0.1
+  $psphotConfig = $psphotConfig -Di PSPHOT:LMM_FIT_GAIN_FACTOR_MODE 2
+  $psphotConfig = $psphotConfig -Db PSPHOT:SAVE.RESID T
+  $psphotConfig = $psphotConfig -D  PSPHOT:EXTENDED_SOURCE_MODELS_SELECTION $fitModel
+
+  # ppImage / psphot on the output
+  echo ppImage $ppImageConfig $psphotConfig -file $basename.fits $outname
+  exec ppImage $ppImageConfig $psphotConfig -file $basename.fits $outname
+end
+
 macro ckchip
   if ($0 != 5)
@@ -173 +220 @@
 
 macro ckgalaxy
-
-  data test.exp.dat
+  if ($0 != 3)
+    echo "USAGE: ckgalaxy (dat) (cmf)"
+    break
+  end
+
+  data $1
   read Xin_all 1 Yin_all 2 Type 4 Min_all 5 RmajIn_all 7 RminIn_all 8 ThetaIn_all 9
 
@@ -185 +236 @@
   subset RminIn = RminIn_all if (Type == 1)
 
-  data test.exp2.cmf 
+  data $2
   read -fits Chip.xfit X_EXT Y_EXT EXT_INST_MAG EXT_WIDTH_MAJ EXT_WIDTH_MIN EXT_THETA
   set EXT_THETA_ALT = EXT_THETA * (EXT_THETA >= 0.0) + (EXT_THETA + 3.14159265) * (EXT_THETA < 0.0)
@@ -200 +251 @@
   set dX = Xin_m - Xot_m  
   set dY = Yin_m - Yot_m  
-  lim -n 0 dX dY; clear; box; plot dX dY
+  # lim -n 0 dX dY; clear; box; plot dX dY
 
   reindex Mot_m = EXT_INST_MAG using index1
@@ -218 +269 @@
   set dM = Min_m - Mot_m  
   set dT = (Tin_m - Tot_m) * 180 / 3.14159265
+  set Tin_deg = Tin_m * 180 / 3.14159265
+
+  lim -n 0 Tin_deg -5 5; clear; box; plot Tin_deg dT; label -y "delta theta (deg)"
   lim -n 1 n dM; clear; box; plot n dM; label -y "delta mag"
-  lim -n 2 n dT; clear; box; plot n dT; label -y "delta theta (deg)"
-  lim -n 3 dT dM; clear; box; plot dT dM; label -x "delta theta (deg)" -y "delta mag"
+  lim -n 2 n -5 5; clear; box; plot n dT; label -y "delta theta (deg)"
+  lim -n 3 -5 5 dM; clear; box; plot dT dM; label -x "delta theta (deg)" -y "delta mag"
 
   set dR = Rin_m - Rot_m