Changeset 36310 for branches/eam_branches/ipp-20130904/psphot/test/tap_psphot_galaxygrid.pro

Timestamp:

Nov 24, 2013, 3:28:13 PM (13 years ago)

Author:

eugene

Message:

add a radius check function to pcm fit source

File:

• branches/eam_branches/ipp-20130904/psphot/test/tap_psphot_galaxygrid.pro (modified) (4 diffs)

branches/eam_branches/ipp-20130904/psphot/test/tap_psphot_galaxygrid.pro

@@ -607 +607 @@
   $psphotConfig = $psphotConfig -Db PSPHOT:POISSON.ERRORS.PHOT.LMM F
   $psphotConfig = $psphotConfig -Db PSPHOT:EXTENDED_SOURCE_FITS_POISSON F
+   $psphotConfig = $psphotConfig -Di PSPHOT:EXT_FIT_ITER 40
 
   # ppImage / psphot on the output
@@ -614 +615 @@
   exec ppImage $ppImageConfig $psphotConfig -file $basename.fits $outname >>& $outname.log
   break -auto on
+end
+
+macro cmf.load.reset
+  $fields = X Y M T R r MT
+  foreach field $fields
+    foreach set in ot
+      delete -q $field\$set\_s
+    end
+  end
+
+  delete -q min_S Min_S
 end
 
@@ -1262 +1274 @@
   local radius Nrun 
 
+  if (1)
+    $Nrun = 0
+    foreach radius 1.0 1.5 2.0 2.5 3.0 4.0 6.0 8.0
+      # mkexp.devexp.single tests.20131120/testrad.ps1v1.exp.$Nrun EXP $radius 1.0
+      fitexp tests.20131120/testrad.ps1v1.exp.$Nrun tests.20131120/testrad.ps1v1.exp.$Nrun.t4s EXP_CONV
+      $Nrun ++
+    end
+  end
+
   $Nrun = 0
   foreach radius 1.0 1.5 2.0 2.5 3.0 4.0 6.0 8.0
-    mkexp.devexp.single tests.20131120/testrad.ps1v1.exp.$Nrun EXP $radius 1.0
-    fitexp tests.20131120/testrad.ps1v1.exp.$Nrun tests.20131120/testrad.ps1v1.exp.$Nrun.t4 EXP_CONV
-    $Nrun ++
-  end
-
-  $Nrun = 0
-  foreach radius 1.0 1.5 2.0 2.5 3.0 4.0 6.0 8.0
-    cmf.load.concat tests.20131120/testrad.ps1v1.exp.$Nrun.dat tests.20131120/testrad.ps1v1.exp.$Nrun.t4.cmf EXP
+    cmf.load.concat tests.20131120/testrad.ps1v1.exp.$Nrun.dat tests.20131120/testrad.ps1v1.exp.$Nrun.t4s.cmf EXP
     set dM_m = Mot_m - Min_m
     vstat -q dM_m
@@ -1284 +1298 @@
   grid.plots.devexp a
 end
+
+macro check.radius.loop
+  if ($0 != 2)
+    echo "USAGE: check.radius.loop (ext)"
+    break
+  end
+
+  # cmf.load.reset
+
+  local radius Nrun 
+
+  $Nrun = 0
+  foreach radius 1.0 1.5 2.0 2.5 3.0 4.0 6.0 8.0
+    cmf.load.concat tests.20131120/testrad.ps1v1.exp.$Nrun.dat tests.20131120/testrad.ps1v1.exp.$Nrun.$1.cmf EXP
+    set dM_m = Mot_m - Min_m
+    vstat -q dM_m
+    $mag_off = $MEDIAN
+    set dR_m = Rot_m - Rin_m
+    vstat -q dR_m
+    $rad_off = $MEDIAN
+    echo $radius $mag_off $rad_off
+    $Nrun ++
+  end
+
+  grid.plots.devexp a
+end
+
+macro check.fit
+  if ($0 != 3)
+    echo "USAGE: check.fit (dat) (cmf)"
+    break
+  end
+
+  cmf.load.concat $1 $2 EXP
+  set dM_m = Mot_m - Min_m
+  vstat -q dM_m
+  $mag_off = $MEDIAN
+  set dR_m = Rot_m - Rin_m
+  vstat -q dR_m
+  $rad_off = $MEDIAN
+  echo $mag_off $rad_off
+end