Changeset 41722

Timestamp:

Jul 15, 2021, 3:54:04 PM (5 years ago)

Author:

eugene

Message:

clean up / improve plots

File:

• trunk/tools/eam/teleff/camstats.sh (modified) (7 diffs)

trunk/tools/eam/teleff/camstats.sh

@@ -83 +83 @@
 
   # allow mag and pos histograms to have different bin sizes?
-  histogram dRs NdRs -0.5 0.5 $HISTBIN -range dx
-  histogram dDs NdDs -0.5 0.5 $HISTBIN -range dx
-  histogram dMs NdMs -0.5 0.5 $HISTBIN -range dx
+  histogram dRs NdRs {-50*$HISTBIN} {50*$HISTBIN} $HISTBIN -range dx
+  histogram dDs NdDs {-50*$HISTBIN} {50*$HISTBIN} $HISTBIN -range dx
+  histogram dMs NdMs {-50*$HISTBIN} {50*$HISTBIN} $HISTBIN -range dx
 
   NdRs[0] = 0; NdRs[-1] = 0
@@ -93 +93 @@
   $C3 = 0.0
 
-  peak -q dx NdRs; $C0 = $peakpos; $C2 = $peaknum; $C1 = 3*$HISTBIN; vgauss dx NdRs con NdRf
-  peak -q dx NdDs; $C0 = $peakpos; $C2 = $peaknum; $C1 = 3*$HISTBIN; vgauss dx NdDs con NdDf
-  peak -q dx NdMs; $C0 = $peakpos; $C2 = $peaknum; $C1 = 3*$HISTBIN; vgauss dx NdMs con NdMf
-
-  section aM; lim 13 23 -0.5 0.5; box $BST1
+  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)
+
+  section aM; lim 13 23 {-50*$HISTBIN} {50*$HISTBIN}; box $BST1
   plot CAL_PSF_MAG dM -pt box -sz 0.5 -op 0.2 where keepGood
   label -x mag -y "&sd&h mag"
   
-  section aR; lim 13 23 -0.5 0.5; box $BST1
+  section aR; lim 13 23 {-50*$HISTBIN} {50*$HISTBIN}; box $BST1
   plot CAL_PSF_MAG dR -pt box -sz 0.5 -op 0.2 where keepGood
   label -x mag -y "&sd&h R.A."
   
-  section aD; lim 13 23 -0.5 0.5; box $BST1
+  section aD; lim 13 23 {-50*$HISTBIN} {50*$HISTBIN}; box $BST1
   plot CAL_PSF_MAG dD -pt box -sz 0.5 -op 0.2 where keepGood
   label -x mag -y "&sd&h Dec"
@@ -113 +113 @@
   plot dx NdMf -x line -c red -lw 2
   label -x "&sd&h mag" +y "N (15 < mag < 17)"
+  sprintf line "&ss&h (mag) = %.0f mmag" {1000*$SigmaMag}; textline -frac 0.02 0.95 -justify 5 "$line"
   
   section bR; lim dx NdRs; box $BST2
@@ -118 +119 @@
   plot dx NdRf -x line -c red -lw 2
   label -x "&sd&h R.A" +y "N (15 < mag < 17)"
+  sprintf line "&ss&h (mag) = %.0f mas" {1000*$SigmaRA}; textline -frac 0.02 0.95 -justify 5 "$line"
   
   section bD; lim dx NdDs; box $BST2
@@ -123 +125 @@
   plot dx NdDf -x line -c red -lw 2
   label -x "&sd&h Dec" +y "N (15 < mag < 17)"
+  sprintf line "&ss&h (Dec) = %.0f mas" {1000*$SigmaDec}; textline -frac 0.02 0.95 -justify 5 "$line"
 end
 
@@ -153 +156 @@
  lim magBin -0.02 1.02; box $BSTY -ticks 1110;
 
+ plot -x hist magBin fracPoss  -lw 3 -c blue70
  plot -x hist magBin fracFound -lw 3 -c red70
- plot -x hist magBin fracPoss  -lw 3 -c blue70
  plot -x hist magBin fracGood  -lw 3 -c red
 
- subset tmp = fracGood  where (magBin > 15) && (magBin < 20); vstat -q tmp; $FracBrightGood  = $MEDIAN
  subset tmp = fracPoss  where (magBin > 15) && (magBin < 20); vstat -q tmp; $FracBrightPoss  = $MEDIAN
  subset tmp = fracFound where (magBin > 15) && (magBin < 20); vstat -q tmp; $FracBrightFound = $MEDIAN
-
- line -c red  -lt dot -lw 2 100 $FracBrightGood to 0 $FracBrightGood
- line -c blue -lt dot -lw 2 100 $FracBrightPoss to 0 $FracBrightPoss
+ subset tmp = fracGood  where (magBin > 15) && (magBin < 20); vstat -q tmp; $FracBrightGood  = $MEDIAN
+
+ line -c blue70 -lt dot -lw 2 100 $FracBrightPoss  to 0 $FracBrightPoss
+ line -c red70  -lt dot -lw 2 100 $FracBrightFound to 0 $FracBrightFound
+ line -c red    -lt dot -lw 2 100 $FracBrightGood  to 0 $FracBrightGood
 
  fprintf "%.1f-pct of all bright sources were detected unmasked" {100*$FracBrightGood} 
  fprintf "%.1f-pct of unmasked bright sources were found" {100*$FracBrightFound / $FracBrightPoss}
  fprintf "%.1f-pct of unmasked bright sources were detected unmasked" {100*$FracBrightGood / $FracBrightPoss}
+
  fprintf "mag limit %.1f" $MagLim
+ sprintf line "mag limit: %.1f" $MagLim
+ textline -frac 0.95 0.95 -justify 3 "$line"
 
  label -x mag -y "detected fraction"
@@ -181 +188 @@
  label +y "detected number"
 
+ vstat -q numTotal; $Ntotal = $TOTAL
  vstat -q numFound; $Nfound = $TOTAL
  vstat -q numGood;  $Ngood  = $TOTAL
  
- fprintf "%d total sources detected, %d unmasked" $Nfound $Ngood
-
- $Xlegend = 0.02; $dXlegend = 0.03
- $Ylegend = 0.95; $dYlegend = 0.02
- line -frac $Xlegend $Ylegend to {$Xlegend + $dXlegend} $Ylegend -c blue  -lw 3; textline -frac -justify 5 {$Xlegend + $dXlegend + 0.01} $Ylegend "unmasked source positions (psfqf > 0.85)"; $Ylegend -= $dYlegend
- line -frac $Xlegend $Ylegend to {$Xlegend + $dXlegend} $Ylegend -c red70 -lw 3; textline -frac -justify 5 {$Xlegend + $dXlegend + 0.01} $Ylegend "all sources detected"                    ; $Ylegend -= $dYlegend
- line -frac $Xlegend $Ylegend to {$Xlegend + $dXlegend} $Ylegend -c red   -lw 3; textline -frac -justify 5 {$Xlegend + $dXlegend + 0.01} $Ylegend "sources detected with psfqf > 0.85"
+ fprintf "%d of %d total sources detected, %d unmasked" $Nfound $Ntotal $Ngood
+
+ legend.init 0.02 0.95 0.03 0.02
+ sprintf line "unmasked source positions (psfqf > 0.85) [%.1f pct]" {100*$FracBrightPoss};  legend.line blue  3 solid "$line"
+ sprintf line "all sources detected [%.1f pct]"                     {100*$FracBrightFound}; legend.line red70 3 solid "$line"
+ sprintf line "sources detected with psfqf > 0.85 [%.1f pct]"       {100*$FracBrightGood};  legend.line red   3 solid "$line"
+end
+
+macro legend.init
+ if ($0 != 5)
+  echo "USAGE: legend.init (X) (Y) (dX) (dY)"
+  break
+ end
+ 
+ $legend:X = $1; $legend:dX = $3
+ $legend:Y = $2; $legend:dY = $4
+end
+
+macro legend.line
+ if ($0 != 5)
+  echo "USAGE: legend.init (color) (lweight) (ltype) (label)"
+  break
+ end
+
+ local color myline mylt mylw
+ 
+ $color  = $1
+ $mylt   = $2
+ $mylw   = $3
+ $myline = $4
+
+ line -frac $legend:X $legend:Y to {$legend:X + $legend:dX} $legend:Y -c $color -lw $mylw -lt $mylt
+ textline -frac -justify 5 {$legend:X + $legend:dX + 0.01} $legend:Y "$myline"
+ $legend:Y -= $legend:dY
 end