function hist_stats, vals, binsize=binsize, minval=minval, maxval=maxval

if n_elements(binsize) eq 0 then binsize = .01
if n_elements(minval) eq 0 then minval = 0.0
if n_elements(maxval) eq 0 then maxval = 1.0

hist = histogram(vals,min=minval,max=maxval,binsize=binsize,reverse_indices=ri,locations=loc)
nbins = n_elements(hist)
res = fltarr(3,nbins)
res[0,*] = loc + binsize/2.0

for i=0,nbins-1 do begin
	if ri[i+1]-ri[i] lt 3 then continue
	inds = ri[ri[i]:ri[i+1]-1]
	res[1:2,i] = [mean(vals[inds]),stddev(vals[inds])]
endfor

return,res

end