macro plot_wico
* 
* Plot parameterizations of Bcal resolution
* 10/31/06  ES
*
* get hitogram file created by geant3
*
hi/file 1 wico.hbook
*
* set options
*
* hi/create/title_global 'Winston Cone Transmission'
option ndate
option nbox
set stat 1111111
option stat
set fit 111
option fit
option grid
*
* plotting options
*
* set *
*set xmgl 4.
*set ymgl 4.
*set asiz 0.4
*set xlab 2.
*set ylab 1.
*set xsiz 20.
*set xmgl 3.
*set ymgl 3.
*set ysiz 20.
*set gsiz 0.4
*
* set font definitions to bold roman 
*
set CFON -21
set GFON -21
set LFON -21
set TFON -21
set VFON -21
set txfp -21
set SMGU 0.02
set SMGR 0.02
set CSIZ 0.33
set VSIZ 0.3
set TSIZ 0.35
set YHTI 0.9
set HWID 3.0
set BWID 3.0
set asiz 0.38
*
* get color map
*
exe color
zone 2 2 
*
* open metafile
*
for/file 66 plot_wico.ps
meta 66 -111
*
npts1 = 101
*
*
* plot 1
*
csize = 0.05
igset chhe [csize]
* hplot/null [emin] [emax] [ymin] [ymax]
zone 2 2
hi/plot 3 colz
hi/plot 2 colz
hi/operations/divide 2 3 101
option logz
* option linz
hi/plot 101 colz
igset chhe 0.3
exe window#push
itx 0.75 -0.18 'Ratio'
exe window#pop
hi/plot 10
wait

option linz
* plot 2
zone 2 2
hi/plot 4
hi/plot 5
hi/plot 6
hi/plot 9 colz
wait

* plot 3
* for/file 66 plot_wico_plot3.eps
* meta 66 -113
zone 2 2
hi/plot 7
hi/plot 17
hi/operations/divide 17 7 107
hi/plot 107
igset chhe 0.3
exe window#push
itx 1.25 0.8 'Distribution of [q]'
itx 1.25 0.7 'at Entrance and Exit'
itx 1.25 0.6  'of Winston Cone'
itx 0.75 -0.18 'Ratio'
exe window#pop
* close 66
wait

* plot 4
zone 2 2
hi/plot 8
hi/plot 18
hi/operations/divide 18 8 108
hi/plot 108
igset chhe 0.3
exe window#push
itx 0.75 -0.18 'Ratio'
exe window#pop
wait

* plot 5
zone 2 2
hi/plot 14
hi/plot 15
hi/plot 16
hi/plot 19 colz
wait
*
 close 66
*
exitm 
return