Changeset 41981 for branches/eam_branches/ipp-20211108

Timestamp:

Jan 7, 2022, 4:38:07 PM (5 years ago)

Author:

eugene

Message:

gracefully handle certain errors; less verbosity on gauss fits

File:

• branches/eam_branches/ipp-20211108/ippTests/ippCheck/camstats.pro (modified) (3 diffs)

branches/eam_branches/ipp-20211108/ippTests/ippCheck/camstats.pro

@@ -1 @@
-# this 
-
 # load.camstats (glob) : load data from a set of stf files and calculate detection histograms
 # plot.camstats.resid  : plot photometric and astrometric residuals
@@ -24 +22 @@
 
  list name -glob $1
+ if ($name:n == 0)
+   echo "ERROR: no files found matching $1"
+   break
+ end
+
  for i 0 $name:n
 
@@ -106 +109 @@
   $C3 = 0.0
 
-  peak -q dx NdRs; $C0 = $peakpos; $C2 = $peaknum; $C1 = 2*$HISTBIN; vgauss dx NdRs con NdRf; $SigmaRA  = abs($C1)
-  peak -q dx NdDs; $C0 = $peakpos; $C2 = $peaknum; $C1 = 2*$HISTBIN; vgauss dx NdDs con NdDf; $SigmaDec = abs($C1)
-  peak -q dx NdMs; $C0 = $peakpos; $C2 = $peaknum; $C1 = 2*$HISTBIN; vgauss dx NdMs con NdMf; $SigmaMag = abs($C1)
+  peak -q dx NdRs; $C0 = $peakpos; $C2 = $peaknum; $C1 = 2*$HISTBIN; vgauss -q dx NdRs con NdRf; $SigmaRA  = abs($C1)
+  peak -q dx NdDs; $C0 = $peakpos; $C2 = $peaknum; $C1 = 2*$HISTBIN; vgauss -q dx NdDs con NdDf; $SigmaDec = abs($C1)
+  peak -q dx NdMs; $C0 = $peakpos; $C2 = $peaknum; $C1 = 2*$HISTBIN; vgauss -q dx NdMs con NdMf; $SigmaMag = abs($C1)
 
   section aM; lim 13 23 {-50*$HISTBIN} {50*$HISTBIN}; box $BST1