;+ ; Contains the fond procedure ; ; :Author: ; Philippe Zarka ; ; :History: ; 2006/03/07: Created ; ; 2006/03/07: Last Edit ;- ; ;+ ; Background and sigma of a 1-D distribution of intensities. ; ; :Params: ; tab: in, required, type=sometype ; 1-D array of intensities ; fon: out, required, type=sometype ; background ; sigma: out, required, type=sometype ; fluctuations (1 sigma level) ;- pro FOND, tab, fon,sigma tab=double(tab) if n_elements(tab) gt 1 then sigma=10.*stdev(tab,fon) else begin $ fon=tab(0) & sigma=0 endelse encore: test=where(abs(tab-fon) lt 2.5*sigma) if n_elements(test) eq 1 then return sigma=stdev(tab(test),moy) if moy eq fon then return fon=moy goto,encore end