Index: branches/eam_branches/ipp-20211108/ippTests/ippCheck/camstats.pro
===================================================================
--- branches/eam_branches/ipp-20211108/ippTests/ippCheck/camstats.pro	(revision 41980)
+++ branches/eam_branches/ipp-20211108/ippTests/ippCheck/camstats.pro	(revision 41981)
@@ -1,4 +1,2 @@
-# 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,4 +22,9 @@
 
  list name -glob $1
+ if ($name:n == 0)
+   echo "ERROR: no files found matching $1"
+   break
+ end
+
  for i 0 $name:n
 
@@ -106,7 +109,7 @@
   $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