macro comp_data
field=perp
read field
ttype=h
pmom=4.723
tar=top
read tar
fname=[tar]_[field][ttype]_[pmom].out
vec/del  vrun,vtar,vtpol,vbpol,vcasy,ptime,ntime,pcharge,ncharge 
vec/read vrun,vtar,vtpol,vbpol,vcasy,ptime,ntime,pcharge,ncharge [fname]
nrtot=$vlen(vrun,1)
filecase keep
chain -all
nst=3
nrtot=8
tot_charge=0
tot_eff=0
do nr=[nst],[nrtot]
nrun=vrun([nr])
sh tclsh get_info_run.tcl [nrun] > t.out
sh cat t.out
vec/del ch,t,f
vec/read ch,t,f t.out
charge=ch(1)
time=t(1)
current=[charge]/[time]*1000.
fideff=f(1)/100.
sh rm t.out
message Total charge = [charge] mC
message Total time = [time] sec
message ave current = [current] nA
message Tracking eff = [fideff]
tot_charge=[tot_charge]+[charge]
tot_eff=[tot_eff]+[fideff]
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
zone 1 2
enddo
nr=[nrtot]-[nst]+1
ave_fideff=[tot_eff]/[nr]
cf=1./[tot_charge]/[ave_fideff]
message Total charge = [tot_charge] and averarge fideff= [ave_fideff]
* background run
if ( [field] .eq. 'para') then
nrun=43775
* C+He during parah running
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
chain -back
btot_charge=0
btot_eff=0
do nr=1,[nchain]
nrun=runlist([nr])
sh tclsh get_info_run.tcl [nrun] > t.out
sh cat t.out
vec/del ch,t,f
vec/read ch,t,f t.out
bcharge=ch(1)
btime=t(1)
bcurrent=[charge]/[time]*1000.
bfideff=f(1)/100.
btot_charge=[btot_charge]+[bcharge]
btot_eff=[btot_eff]+[bfideff]
sh rm t.out
message Background Total charge = [bcharge] mC
message Background Total time = [btime] sec
message Background ave current = [bcurrent] nA
message Background Tracking eff = [bfideff]

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 back [fname]
nseg=[nseg]+1
fname=[dir]/ntup/hms[nrun].[nseg].hbook
endwhile
zone 1 2
enddo
bavefideff=[btot_eff]/[nchain]
bcf=1./[btot_charge]/[bavefideff]
message Background Total charge = [btot_charge] and averarge fideff= [bavefideff]
*
zone 1 1
cut $1 hsp>.01.and.0.8<hsshtrk/hsp.and.abs(hsdelta)<8.
nbin=50
1dhist 1000 'W ' [nbin] .6 1.1
ntu/proj 1000 //all/9010.W $1
hrin 1000
vec/create wv([nbin]) r
vec/create wverr([nbin]) r [nbin]*0.0001
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
sigma tot=tot*[cf]
sigma etot=etot*[cf]
1dhist 1000 'W ' [nbin] .6 1.1
ntu/proj 1000 //back/9010.W $1
hrin 1000 
vec/create btot([nbin]) r
vec/create ebtot([nbin]) r
hist/get_vect/contents 1000 btot
hist/get_vect/error 1000 ebtot
sigma btot=btot*[bcf]
sigma ebtot=ebtot*[bcf]
sigma rat=btot/tot
sigma errrat=rat*sqrt(ebtot*ebtot/btot/btot+etot*etot/tot/tot)
hplot/errors wv tot wverr etot [nbin] 24 .1 W
set pcol 2.
hplot/errors wv btot wverr ebtot [nbin] 20 .1 
set pcol 1.
wait
hplot/errors wv rat wverr errrat [nbin] 24 .1 W
vec/write wv,rat,errrat [tar]_[field]_data_Che_NH3_comp.dat 3(f15.5,1x)
*
* get weigthed ave 0.6 < W < .85
vec/create tempa(25) r
vec/create tempe(25) r
vec/copy rat(1:25) tempa(1:25)
vec/copy errrat(1:25) tempe(1:25)
sigma wserr=vsum(1/tempe/tempe)
sigma werr=sqrt(1./wserr)
sigma wave=vsum(tempa/tempe/tempe)
wt=wave(1)/wserr(1)
we=werr(1)
wlo=wv(1)
whi=wv(25)
message  [wlo] < W < [whi] [field] wave=[wt]+-[we]