Changeset 42125 for trunk/ippTests

Timestamp:

Mar 17, 2022, 1:54:21 PM (4 years ago)

Author:

eugene

Message:

update ippCheck to handle empty output files

Location:

trunk/ippTests

Files:

• . (modified) (1 prop)
• ippCheck (modified) (1 prop)
• ippCheck/cammask.pro (modified) (2 diffs)
• ippCheck/camstats.pro (modified) (4 diffs)
• ippCheck/diffmask.pro (modified) (1 diff)
• ippCheck/ippCheck (modified) (2 diffs)
• ippCheck/warpmask.pro (modified) (1 diff)
• ippCheck/warpstats.pro (modified) (3 diffs)

trunk/ippTests/ippCheck/cammask.pro

@@ -78 +78 @@
  if ($table:Ny == 0)
   echo "ERROR: no reference sources found in $refsrcs"
+
+  # generate an empty output file
+  $OUT_FIELDS = 
+  foreach field $SMF_FIELDS_CAM $REF_FIELDS_CAM $XTR_FIELDS_CAM
+    $OUT_FIELDS = $OUT_FIELDS $field\_out
+    create -empty $field\_out
+  end
+  
+  sprintf output "%s/%s.%s.cam.stf" $OUTDIR $expname $OUTVERSION
+  write -fits DATA $output $OUT_FIELDS
+  
+  # save this name for warp and diff analysis, if desired
+  $CAMSOURCE = $output
+
   break
  end
@@ -271 +285 @@
    if ($COPYMASKS) exec neb-copy $camMask:$i $OUTDIR
    $j ++
- end
- $camMask:n = $chipName:n
+   if ($j >= $tmpName:n) last
+ end
+ $camMask:n = $tmpName:n
 end
 

trunk/ippTests/ippCheck/camstats.pro

@@ -38 +38 @@
    # results in magFound, fracFound, numFound
    calc.camstats
-   if ($i == 0)
+   if (numTotal_sum?[] == 0)
      set numTotal_sum = numTotal
      set numFound_sum = numFound
@@ -57 +57 @@
  end
    
- set numTotal = numTotal_sum
- set numFound = numFound_sum
- set numPoss  = numPoss_sum
- set numGood  = numGood_sum
- delete numTotal_sum numFound_sum numPoss_sum numGood_sum
+ if (numTotal_sum?[] == 0)
+   calc.camstats 
+ else 
+   set numTotal = numTotal_sum
+   set numFound = numFound_sum
+   set numPoss  = numPoss_sum
+   set numGood  = numGood_sum
+   delete numTotal_sum numFound_sum numPoss_sum numGood_sum
+ end
 
  set fracFound = (numTotal > 0) ? numFound / numTotal : zero(numTotal)
@@ -79 +83 @@
   $BST2 = -tickpad 0.2 -xpad 3.5 -labelpadx 3.0 -ypad 0.5 -labelpady 0.0 +ypad 4.5 +labelpady 4.0 -labels 1001
   
-  section aM 0.00 0.66 0.65 0.33
-  section aR 0.00 0.33 0.65 0.33
-  section aD 0.00 0.00 0.65 0.33
-
-  section bM 0.65 0.66 0.35 0.33
-  section bR 0.65 0.33 0.35 0.33
-  section bD 0.65 0.00 0.35 0.33
+  section aM 0.00 0.66 0.65 0.33;  label -x mag -y "&sd&h mag"
+  section aR 0.00 0.33 0.65 0.33;  label -x mag -y "&sd&h R.A."
+  section aD 0.00 0.00 0.65 0.33;  label -x mag -y "&sd&h Dec" 
+
+  section bM 0.65 0.66 0.35 0.33;  label -x mag -y "&sd&h mag" 
+  section bR 0.65 0.33 0.35 0.33;  label -x mag -y "&sd&h R.A."
+  section bD 0.65 0.00 0.35 0.33;  label -x mag -y "&sd&h Dec" 
+
+  if (RA_PSF?[] == 0)
+    echo "no sources in datafile"
+    return
+  end
 
   set dR = 3600*(RA_PSF - Rref)*dcos(Dref)
@@ -256 +265 @@
  $Mdel = 0.25
  
+ if (Mref_out?[] == 0)
+   create magBin $Mmin $Mmax $Mdel
+   set numTotal = zero(magBin)
+   set numFound = zero(magBin)
+   set numPoss  = zero(magBin)
+   set numGood  = zero(magBin)
+   return
+ end
+
  delete -q  magBin numTotal numPoss numFound numGood
  for mag $Mmin $Mmax $Mdel

trunk/ippTests/ippCheck/diffmask.pro

@@ -51 +51 @@
   # load the reference source table
   # how to choose the name of the filter?
-  data $refsrcs; 
+  data $refsrcs
+
+  # first check if the FITS table has any entries
+  read -fits DATA -sizes
+  if ($table:Ny == 0) 
+    echo "WARNING: empty reference table $refsrcs"
+
+    $OUT_FIELDS = 
+    foreach field $SMF_FIELDS_DIF $REF_FIELDS_DIF $XTR_FIELDS_DIF
+      $OUT_FIELDS = $OUT_FIELDS $field\_out
+      create -empty $field\_out
+    end
+    
+    sprintf output "%s/%s.%s.dif.stf" $OUTDIR $expname $OUTVERSION
+    write -fits DATA $output $OUT_FIELDS
+    return
+  end
+
   if ($USE_CAMSOURCE) 
    read -fits DATA Rref_out Dref_out Mref_out

trunk/ippTests/ippCheck/ippCheck

@@ -59 +59 @@
   getargs -no-warp  BOOL  -var SKIP_WARP
   getargs -no-diff  BOOL  -var SKIP_DIFF
-  getargs -do-plots BOOL  -var DO_PLOTS;   
+  getargs -no-plots BOOL  -var SKIP_PLOTS;   
   getargs -use-refs BOOL  -var USE_REFSOURCE
   getargs -camera   VALUE -var CAMERA;
@@ -97 +97 @@
   end
 
-  if ($DO_PLOTS) plot.exposure.real $expname
+  if (not($SKIP_PLOTS)) plot.exposure.real $expname
 end
 

trunk/ippTests/ippCheck/warpmask.pro

@@ -44 +44 @@
   # NOTE: Mref is chosen to match filter for this exposure
   data $refsrcs                 ;
+
+  # first check if the FITS table has any entries
+  read -fits DATA -sizes
+  if ($table:Ny == 0) 
+    echo "WARNING: empty reference table $refsrcs"
+
+    # generate empty output vectors for output file
+    $OUT_FIELDS = 
+    foreach field $REF_FIELDS_WRP $XTR_FIELDS_WRP
+      $OUT_FIELDS = $OUT_FIELDS $field\_out
+      create -empty $field\_out
+    end
+    
+    sprintf output "%s/%s.%s.wrp.stf" $OUTDIR $expname $OUTVERSION
+    write -fits DATA $output $OUT_FIELDS
+
+    return
+  end
+
   if ($USE_CAMSOURCE) 
    read -fits DATA Rref_out Dref_out Mref_out

trunk/ippTests/ippCheck/warpstats.pro

@@ -24 +24 @@
    # results in magFound, fracFound, numFound
    calc.warpstats
-   if ($i == 0)
+   if (numTotal_sum?[] == 0)
      set numTotal_sum = numTotal
      set numPoss_sum  = numPoss
@@ -39 +39 @@
  end
    
- set numTotal = numTotal_sum
- set numPoss  = numPoss_sum
- delete numTotal_sum numPoss_sum
+ if (numTotal_sum?[] == 0)
+    # generate vectors for empty inputs
+   calc.warpstats
+ else   
+   set numTotal = numTotal_sum
+   set numPoss  = numPoss_sum
+   delete numTotal_sum numPoss_sum
+ end
 
  set fracPoss  = (numTotal > 0) ? numPoss  / numTotal : zero(numTotal)
@@ -89 +94 @@
  $Mdel = 0.25
  
+ if (Mref_out?[] == 0)
+   create magBin $Mmin $Mmax $Mdel
+   set numTotal = zero(magBin)
+   set numPoss  = zero(magBin)
+   return
+ end
+
  delete -q  magBin numTotal numPoss
  for mag $Mmin $Mmax $Mdel