;+
; 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