Index: trunk/doc/release.2015/ps1.analysis/analysis.tex
===================================================================
--- trunk/doc/release.2015/ps1.analysis/analysis.tex	(revision 39855)
+++ trunk/doc/release.2015/ps1.analysis/analysis.tex	(revision 39857)
@@ -1619,2 +1619,22 @@
 * table of mask image bit values
 * table of models?
+
+Work still needed:
+
+* Figures
+* Tables
+* refereces for other programs
+
+* moments & gauss sigma issue
+* words on the 2d maps
+* PS1 PSF profile discussion
+* PSF residual map
+* section 3.5.3 Model applied to detected objects needs to be reviewed
+
+* read for english & phrasing
+* convert XXX
+* background model description (see waters)
+* reduce / remove design goal wording?
+* too many bullets?
+* reduce coding description?
+* put engineering docs (psLib, psModules) on public website 
Index: trunk/doc/release.2015/ps1.analysis/plots.sh
===================================================================
--- trunk/doc/release.2015/ps1.analysis/plots.sh	(revision 39857)
+++ trunk/doc/release.2015/ps1.analysis/plots.sh	(revision 39857)
@@ -0,0 +1,49 @@
+
+macro peak.and.col
+
+  $scale = 50.0
+  $sigma = 3.0
+  $Ia = 500; $Xa = 100
+  $Ib = 100 ; $Xb = 50
+  $Ic = 30  ; $Xc = 180
+
+  create x 0 200
+  set ya = $Ia * exp (-1*(x - $Xa)^2/$scale^2)
+  set yb = $Ib * exp (-0.5*(x - $Xb)^2/$sigma^2)
+  set yc = $Ic * exp (-0.5*(x - $Xc)^2/$sigma^2)
+
+  set ys = ya + yb + yc
+  gaussdev ds x[] 0.0 1.0
+  set dy = ds*(sqrt(ys + 50))
+  set dY = (sqrt(ys + 20))
+  set yo = ys + dy
+  set ym = yo
+  vsmooth ym 2.0
+  
+  clear -s
+  section a 0.0 0.5 1.0 0.5
+  lim x yo; box; plot x yo -x hist
+
+  section b 0.0 0.0 1.0 0.5
+  lim x yo; box; plot x ym -x hist
+
+  set yd = 500 - ym
+
+  $dX = 5
+  peak x ym {$Xa - $dX} {$Xa + $dX}
+  line -c red $peakpos {$peakval + 10} to $peakpos $YMAX
+
+  peak x ym {$Xb - $dX} {$Xb + $dX}
+  line -c red $peakpos {$peakval + 10} to $peakpos $YMAX
+
+  peak x ym {$Xc - $dX} {$Xc + $dX}
+  line -c red $peakpos {$peakval + 10} to $peakpos $YMAX
+
+  $dX = 5
+  peak x yd $Xb {$Xb + 2*$dX}
+  line -c blue $peakpos {ym[$peakpos] - 10} to $peakpos $YMIN
+
+  peak x yd {$Xc - 2*$dX} $Xc 
+  line -c blue $peakpos {ym[$peakpos] - 10} to $peakpos $YMIN
+
+end