vec/del *
filecase keep
opt zfl
dir=/wo
field=para
read field
filecase keep
chain -all
*do nr=1,1
if ( [field] .eq. 'para') then
nrun=43775
endif
if ( [field] .eq. 'perp') then
nrun=43407
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
do i=1,[nchain]
nrun=runlist([i])
nseg=1
dir=/work/hallc/e93026/rss/pass3
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 1 2
txt = [fname]  Fit to background
title_gl [txt]
cut $1 hsp>.01.and.0.8<hsshtrk/hsp<1.4.and.abs(hsdelta)<8.
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
sum=$hinfo(1000,'sum')
vec/create par(6) r 5000. 1.1 .05 0. 0. 0.
set fcol 2.
hist/fit 1000 g+p2 SQ 6 par
vec/write par ' ' 6e15.5
set fcol 1.
do i=1,6
p[i]=par([i])
enddo
sigma y=[p1]*exp(-0.5*(wv-[p2])*(wv-[p2])/[p3]/[p3])+[p4]+[p5]*wv+[p6]*wv*wv
SET PLCI 4 
pline 40 wv y
SET PLCI 1
sigma df=1-y/tot
sigma edf=y/tot/tot*etot
vec/create evw([nbin]) r [nbin]*.001
hplot/errors wv df evw edf 200 24 .1 W
pict/print fit_bg_g+p2_data_[nrun].ps