f=occ_dbase('occ_discsp*', '*', '*')
f = f(where( strpos(strlowcase(f),'occ_discsp') ne -1))
break_file, f,disk,dir,fnam,ext
date=strmid(fnam,11,6)
sdate=date
det=strmid(ext,5,1)
det=fix(det)
for i=0,n_elements(date)-1 do date(i)= arr2str(string(reform(byte(date(i)),2,3)),'/')
udates=date(uniq(date,sort(date)))
usdate= sdate(uniq(sdate,sort(sdate)))
get_occ_onday,udates(0),occout=occout,det_id=det_id 
for j=0,n_elements(occout)-1 do begin
        for m=0,2 do begin
                plot_occs, occout=occout(j),step=step,tdur=60.+30*m
                occout(j).counts(*,*,m) = step
                endfor
        endfor


inradec = [83.5160, 22.0197]
gro_point, [occout(0).zradec,occout(0).xradec],inradec, cos_lad,cos_spec 
print,cos_spec
;    -0.338114     0.211976    -0.444681     0.105408    -0.105408     0.444681    -0.211976     0.338114
edges=discsp_edges(occout(0).(0),det=occout(0).det_id)
print,edges
;      25.8513      52.0713
;      52.0713      14809.3
;      14809.3      40000.0
;      40000.0      100000.
apar=[.0025,2.12,fltarr(2)]
resp2=fltarr(2)
for i=0,1 do begin 
	discsp_resp, edges(*,i) <1000. >3.0, $
		!radeg*acos(cos_spec(occout(0).det_id)) <80., apar,resp 
	resp2(i)=resp 
endfor
print,resp2
;      8.52851      12.6790
y=fltarr(2,n_elements(occout))
for i=0,n_elements(occout)-1 do $
  for j=0,1 do begin
	w =  (reverse(sort(occout(i).counts(1,j,*))))(0)
	y(j,i) = occout(i).counts(0,j,w)
  endfor
obs=fltarr(2)
cut = 2 ;2 standard deviations
for j=0,1 do begin
	resistant_mean, (y(j,*))(*), cut, mean
	obs(j) = mean
endfor
end