macro bcal_DRPDE_tradeoff
* 
* Plot parameterizations of Bcal dynamic range
* 05/08/07  ES
*
* set options
*
hi/create/title_global 'Trade off between Dark rate and PDE'
option ndate
option nbox
*set stat 1111111
set stat 111
option nstat
set fit 111
option fit
option grid
*
* plotting options
*
* set *
set xmgl 4.
set ymgl 4.
set asiz 0.5
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.45
set VSIZ 0.5
set TSIZ 0.5
set YHTI 0.9
set HWID 3.0
set BWID 3.0
*
zone
*
* open metafile
*
for/file 66 bcal_DRPDE_tradeoff.ps
meta 66 -111
*for/file 66 bcal_DRPDE_tradeoff.eps
*meta 66 -113
*
npts1 = 101
*
* plot tradeoff between dark rate DR and PDE for fixed resolution.
*
vector/create PDE([npts1]) R [npts1]*0.
vector/create DR1([npts1]) R [npts1]*0.
*
fmin = 0.10
fmax = 0.30
ymin = 10
ymax = 70
* option logx 
* option logy
*
npts = [npts1] - 1
sigma j=array([npts1],0#[npts1]-1)
sigma PDE = [fmin] + ([fmax]-[fmin])*j/[npts]
*
* set contants for resolution GlueX-doc-795 p. 3
*
* b = 0.34
* c = 0.03
* updated 10/29/07
* b = 0.69
* c = 0.12
* updated 4/24/09
b = 0.69
c = 0.085
*
* f=0.12, D=32MHz E=0.06
*
e = 0.06
DeE = 7.5 
sigma DR1 = PDE*[e]*(PDE*[DeE]*[DeE]*[e] - [b]*[b])/([c]*[c])
* convert to MHz
sigma DR1 = DR1*100
* e = 0.06
DeE = 6.5
sigma DR2 = PDE*[e]*(PDE*[DeE]*[DeE]*[e] - [b]*[b])/([c]*[c])
sigma DR2 = DR2*100
* e = 0.06
DeE = 8.5
sigma DR3 = PDE*[e]*(PDE*[DeE]*[DeE]*[e] - [b]*[b])/([c]*[c])
sigma DR3 = DR3*100
*
*
vector/print PDE
wait
vector/print DR1
*
message 'DeE=' [DeE] 'E=' [e]
wait
*
csize = 0.20
igset chhe [csize]
set ndvx 505
set ndvy 505
hplot/null [fmin] [fmax] [ymin] [ymax]
hplot/atitle  'PDE' 'Dark Rate (MHz)'
set pmci 2
set plci 2
hplot/symbol PDE DR1 [npts1] 20 [csize]
set pmci 4
set plci 4
hplot/symbol PDE DR2 [npts1] 20 [csize]
set pmci 3
set plci 3
hplot/symbol PDE DR3 [npts1] 20 [csize]
set plci 1
line 0.21 [ymin] 0.21 [ymax]
line [fmin] 40 [fmax] 40
*
*
set pmci 4
key 0.11 0.9*[ymax] 20 '[D]E/E(stat)=0.075, E=0.06 GeV'
set pmci 2
key 0.11 0.85*[ymax] 20 '[D]E/E(stat)=0.085, E=0.06 GeV (nominal)'
set pmci 3
key 0.11 0.80*[ymax] 20 '[D]E/E(stat)=0.095, E=0.06 GeV'
* exe window#pop
*
*
 close 66
*
exitm 
return