macro plot_resolutions
* 
* Read resolution parameters obtained from momres and plot
* 1/27/07
*
* set options
*
option ndate
option nbox
*set stat 1111111
set stat 111
option nstat
set fit 111
next
zone
*
* plotting options
*
set xmgl 4.
*set xmgl 2.
set ymgl 4.
*set ymgl 3.
set asiz 0.4
set xlab 2.
set ylab 1.
set xsiz 20.
set ysiz 20.
*
* 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 gsiz 0.35
set ygti 1.0
title 'Hall D CDC at 90 degrees. [s]=150[m], scatt chamber, st'

*
*
npts = 5
*
* resolution parameters for deltap/p, phi (bending angle), and vertex
*
vector/create layers([npts]) R [npts]*0.
vector/create dpsig1([npts]) R [npts]*0.
vector/create dpsig2([npts]) R [npts]*0.
vector/create dphisig1([npts]) R [npts]*0.
vector/create dphisig2([npts]) R [npts]*0.
vector/create dxsig1([npts]) R [npts]*0.
vector/create dxsig2([npts]) R [npts]*0.
*
* open metafile
*
*for/file 66 plot_resolutions.eps
*meta 66 -113
for/file 66 plot_resolutions.ps
meta 66 -111
*hplot/key 0.2 2.8 22 'R1 occupancy 0.8"Y# '
*
* read data file
*
vector/read layers,dpsig1,dpsig2,dphisig1,dphisig2,dxsig1,dxsig2 'plot_resolutions.dat'
*
vector/print layers
vector/print dpsig1
vector/print dpsig2
vector/print dphisig1
vector/print dphisig2
vector/print dxsig1
vector/print dxsig2
wait
*
*
* fill vector
*
set ndvx 505
set ndvy 505
zone 2 2
*
* delta p / p
*
csize = 0.2
igset chhe 0.2
xmin = 20
xmax = 30
ymin = 0
ymax = 0.02
set pmci 4
hplot/null [xmin] [xmax] [ymin] [ymax]
hplot/symbol layers dpsig1 [npts] 20 [csize]
set pmci 2
hplot/symbol layers dpsig2 [npts] 21 [csize]
hplot/atitle 'Number of Layers' ''
exe window#push
igset chhe 0.3
itx 0.1 0.85 '[D]p/p'
set pmci 4
hplot/key 0.4 0.9 20 '[s]?1!(GeV^-1!)'
set pmci 2
hplot/key 0.4 0.8 21 '[s]?2!'
exe window#pop
*
*
* delta phi (bending angle)
*
csize = 0.2
igset chhe 0.2
xmin = 20
xmax = 30
ymin = 0
ymax = 0.006
set pmci 4
hplot/null [xmin] [xmax] [ymin] [ymax]
hplot/symbol layers dphisig1 [npts] 20 [csize]
set pmci 2
hplot/symbol layers dphisig2 [npts] 21 [csize]
hplot/atitle 'Number of Layers' ''
exe window#push
igset chhe 0.3
itx 0.1 0.85 '[DQ]'
set pmci 4
hplot/key 0.4 0.9 20 '[s]?1!(rad)'
set pmci 2
hplot/key 0.4 0.8 21 '[s]?2!(rad.GeV)'
exe window#pop
*
*
* delta dx
*
csize = 0.2
igset chhe 0.2
xmin = 20
xmax = 30
ymin = 0
ymax = 500
set pmci 4
hplot/null [xmin] [xmax] [ymin] [ymax]
hplot/symbol layers dxsig1 [npts] 20 [csize]
set pmci 2
hplot/symbol layers dxsig2 [npts] 21 [csize]
hplot/atitle 'Number of Layers' ''
exe window#push
igset chhe 0.3
itx 0.1 0.85 '[D]x'
set pmci 4
hplot/key 0.4 0.9 20 '[s]?1!([m]m)'
set pmci 2
hplot/key 0.4 0.8 21 '[s]?2!([m]m.GeV)'
*
itx  1.3 0.6 '[D]p/p=sqrt([s]^2!&?1!p^2!+([s]?2!/[b])^2!)'
itx  1.3 0.4 '[D][Q]=sqrt([s]^2!&?1!+([s]?2!/p[b])^2!)'
itx  1.3 0.2 '[D]x=sqrt([s]^2!&?1!+([s]?2!/p[b])^2!)'
exe window#pop
*
*
* close file
*
*close 1
close 66
*
*
exitm 
return