vec/del *
filecase keep
opt zfl
dir=/wo
field=para
read field
filecase keep
chain -all
*do nr=1,1
ntheta=2
if ( [field] .eq. 'para') then
      vec/create vthlo([ntheta]) r 11.5    13.0875
      vec/create vthhi([ntheta]) r 13.0875 15.0
endif
if ( [field] .eq. 'perp') then
      vec/create vthlo([ntheta]) r  11.5    12.7625
      vec/create vthhi([ntheta]) r  12.7625 16.0
endif
if ( [field] .eq. 'para') then
nrun=43775
nchain=3
vec/create runlist([nchain]) i 43775 43776 43799
endif
if ( [field] .eq. 'perp') then
nrun=43407
nchain=6
vec/create runlist([nchain]) i 43227 43230 43285 43324 43395 43396
endif
dir=/work/hallc/e93026/rss/pass3
do i=1,[nchain]
nrun=runlist([i])
nseg=1
fname=[dir]/ntup/hms[nrun].[nseg].hbook
message check [fname]
while ($FEXIST([fname])) do
message add [fname] to chain
chain all [fname]
nseg=[nseg]+1
fname=[dir]/ntup/hms[nrun].[nseg].hbook
endwhile
enddo
zone 2 2
txt = [fname]  
title_gl [txt]
do nth=1,[ntheta]
thlo=vthlo([nth])
thhi=vthhi([nth])
if ([field] .eq. 'perp') then
ntu/cut $11 abs(hsyptar)<.04.and.abs(hsxptar+.05)<.12
endif
if ([field] .eq. 'para') then
ntu/cut $11 abs(hsyptar)<.04.and.abs(hsxptar)<.12
endif
ntu/cut $2 [thlo]<hstheta<=[thhi].and.$11
cut $1 $2.and.hsp>.01.and.0.8<hsshtrk/hsp.and.abs(hsdelta)<12.
nbin=50
1dhist 1000 'W ' [nbin] .6 1.1
ntu/proj 1000 //all/9010.W $1
hist/plot 1000 e
vec/create wv([nbin]) r
vec/create tot([nbin]) r
vec/create etot([nbin]) r
hist/get_vect/contents 1000 tot
hist/get_vect/error 1000 etot
hist/get_vect/abscissa 1000 wv
vec/write wv,tot,etot [field]_bg_[nrun]_th[nth].dat 3(f15.5,1x)
enddo
pict/print extract_bg_data_[nrun]_th_[nth].ps