vec/del *
filecase keep
opt zfl
dir=/wo
field=para
read field
filecase keep
chain -all
*do nr=1,1
ntheta=4
if ( [field] .eq. 'para') then
      vec/create vthlo([ntheta]) r 11.5   12.225  12.775 13.425
      vec/create vthhi([ntheta]) r 12.225 12.775  13.425 15.0
endif
if ( [field] .eq. 'perp') then
      vec/create vthlo([ntheta]) r  11.5    12.4875  13.0875 13.8125
      vec/create vthhi([ntheta]) r  12.4875 13.0875  13.8125 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])
ntu/cut $2 [thlo]<hstheta<=[thhi]
cut $1 $2.and.hsp>.01.and.0.8<hsshtrk/hsp.and.abs(hsdelta)<12.
1dhist 1000 'W ' 40 .7 1.1
ntu/proj 1000 //all/9010.W $1
hist/plot 1000 e
nbin=40
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].ps