vec/del *
ntheta=4
do nth=1,[ntheta]
tar=top
do j=1,2
field=para
do i=1,2
vec/del vw,asymc,acerr,vw2,df,dferr
vec/read vw,asymc,acerr [tar]_[field]h_elas_count_asym_pass3_th[nth].dat
wtype=w_0.7_0.85
*wtype=elas
vec/read vw2,df,dferr  [tar]_[field]h_[wtype]_df_data_data_th[nth].dat
sigma asy=asymc/df
sigma aerr=acerr/asymc*asy
vec/write vw,asy,aerr [tar]_[field]h_[wtype]_asym_df_data_data_th[nth].dat 3(f10.5,1x)
lo=vw(20) 
hi=vw(30)
vec/create tempa(11) r
vec/create tempe(11) r
vec/copy asy(20:30) tempa(1:11)
vec/copy aerr(20:30) tempe(1:11)
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)
message For theta bin [nth] [lo] < W < [hi] [tar] [field] wave=[wt]+-[we]
field=perp
enddo
tar=bot
enddo
enddo