Changeset 39865 for trunk/doc/release.2015/ps1.analysis/plots.sh

Timestamp:

Dec 14, 2016, 7:13:40 PM (10 years ago)

Author:

eugene

Message:

add figures, updates to analysis

File:

• trunk/doc/release.2015/ps1.analysis/plots.sh (modified) (4 diffs)

trunk/doc/release.2015/ps1.analysis/plots.sh

@@ +1 @@
+
+macro choose.seed
+
+  for i 0 100
+    rndseed $i
+    peak.and.col 
+    echo $i
+    cursor
+  end
+end
 
 macro peak.and.col
+
+  # using 3 gives a pretty look
+  rndseed 3
 
   $scale = 50.0
@@ -6 +19 @@
   $Ia = 500; $Xa = 100
   $Ib = 100 ; $Xb = 50
-  $Ic = 30  ; $Xc = 180
+  $Ic = 20  ; $Xc = 180
 
   create x 0 200
@@ -22 +35 @@
   
   clear -s
-  section a 0.0 0.5 1.0 0.5
-  lim x yo; box; plot x yo -x hist
+  resize 1200 600
+  label -fn helvetica 24
+  section a 0.0 0.5 1.0 0.45
+  lim x yo; box -lw 2 -xpad 0.5 -labels 0100 -ticks 1100; plot x yo -x hist -lw 2
+  label -y "Raw Counts"
 
-  section b 0.0 0.0 1.0 0.5
-  lim x yo; box; plot x ym -x hist
+  section b 0.0 0.0 1.0 0.55
+  lim x yo; box -lw 2 +xpad 0.5 -xpad 3.5 -labelpadx 3.0 -ticks 1100; plot x ym -x hist
+  label -y "Smoothed Counts" -x "Pixel Coordinate"
 
   set yd = 500 - ym
@@ -32 +49 @@
   $dX = 5
   peak x ym {$Xa - $dX} {$Xa + $dX}
-  line -c red $peakpos {$peakval + 10} to $peakpos $YMAX
+  line -c red $peakpos {$peakval - 10} to $peakpos 400; textline 90 350 "Primary Peak"
 
   peak x ym {$Xb - $dX} {$Xb + $dX}
-  line -c red $peakpos {$peakval + 10} to $peakpos $YMAX
+  line -c red $peakpos {$peakval + 10} to $peakpos 350; textline 30 400 "Significant Peak"
 
   peak x ym {$Xc - $dX} {$Xc + $dX}
-  line -c red $peakpos {$peakval + 10} to $peakpos $YMAX
+  line -c red $peakpos {$peakval + 10} to $peakpos 250; textline 160 300 "Insignificant Peak"
 
   $dX = 5
   peak x yd $Xb {$Xb + 2*$dX}
-  line -c blue $peakpos {ym[$peakpos] - 10} to $peakpos $YMIN
+  line -c blue $peakpos {ym[$peakpos] - 10} to $peakpos 150; textline 50 100 "Col"
 
   peak x yd {$Xc - 2*$dX} $Xc 
-  line -c blue $peakpos {ym[$peakpos] - 10} to $peakpos $YMIN
+  line -c blue $peakpos {ym[$peakpos] + 10} to $peakpos 150; textline 170 190 "Col"
+
+  png -name peaks.png
+  ps  -name peaks.ps
 
 end