pro offband2,fno
;use only for 309 nm filter
!p.multi=[0,1,2]
!p.charsize=0.8
j=fno

!p.psym=10
!x.style=1
!y.range=[0,2000]
!x.range=[3000,7000]
!x.title='Wavelength (A)'
!y.title='Brightness (Rayleighs/A)'
!p.title='Hartley 2 FOS composite'

file1='bq02t'
groups_np,file1,w1,ray1
plot,w1,ray1
file2='b303t'
groups_np,file2,w2,ray2
oplot,w2,ray2*2.
file3='b502t'
groups_np,file3,w3,ray3
oplot,w3,ray3

openr,1,'~/contour/filters/Filter_20Pred.dat'
dummy=''
for i=0,1 do readf,1,dummy
data=fltarr(10,900)
readf,1,data
close,1
name=[' ','309 nm','344.8 nm','387 nm','445 nm']
wave=10.*data(0,*)

index1=where(w1 ge 3000)
len1=n_elements(index1)
r1=ray1(index1)
ww1=w1(index1)
q1=0.2
t1=interpol(data(j,*),wave,ww1)
for i=0,len1-1 do r1(i)=r1(i)*q1*t1(i)
diff1=0.1*(ww1(10)-ww1(0))
sum1=total(r1)*diff1
index1a=where(ww1 le 3200)
sum1a=total(r1(index1a))*diff1
print,diff1,sum1,sum1a

index2=where(w2 ge 3300 and w2 le 4600)
len2=n_elements(index2)
r2=ray2(index2)*2.
ww2=w2(index2)
q2=ww2/1.e4
t2=interpol(data(j,*),wave,ww2)
for i=0,len2-1 do r2(i)=r2(i)*t2(i)*q2(i)
diff2=-0.1*(ww2(10)-ww2(0))
sum2=total(r2)*diff2
index2a=where(ww2 ge 3380 and ww2 le 3520)
sum2a=total(r2(index2a))*diff2
index2b=where(ww2 ge 3800 and ww2 le 3950)
sum2b=total(r2(index2b))*diff2
print,diff2,sum2,sum2a,sum2b

index3=where(w3 ge 4600)
len3=n_elements(index3)
r3=ray3(index3)
ww3=w3(index3)
q3=ww3/1.e4
index6=where(ww3 gt 6000)
q3(index6)=0.6
t3=interpol(data(j,*),wave,ww3)
for i=0,len3-1 do r3(i)=r3(i)*t3(i)*q3(i)
diff3=-0.1*(ww3(10)-ww3(0))
sum3=total(r3)*diff3
print,diff3,sum3

if (j eq 1) then rmax=max(r1)
if (j ne 1) then rmax=max(r2)
!y.range=[1.e-6,1.]
!y.title='Relative signal'
!y.style=1
!p.title='CONTOUR Barr estimate '+name(j)
plot_io,ww1,r1/rmax
;plot,ww1,r1,/ylog
oplot,ww2,r2/rmax
oplot,ww3,r3/rmax
if (j eq 1) then oob=100.*(sum1-sum1a+sum2+2*sum3)/sum1a
if (j eq 2) then oob=100.*(sum1-sum2a+sum2+1.5*sum3)/sum2a
if (j eq 3) then oob=100.*(sum1-sum2b+sum2+2*sum3)/sum2b
soob=string(format='(f7.4)',oob)
xyouts,5000.,0.1,'Out of band = '+soob+'%',size=1.0
datetag
end