[RM-5547] [linemeas] estim deltaFlux from residual

Created on 2020-01-20 10:51:25 by Didier Vibert. % Done: 0

Pour la mesure du Flux des raies ou du continu, modifier l'estimation de l'erreur de mesure du Flux: remplacer le bruit (variance associée au spectre en entrée) par la variance mesurée sur le résidu du fit

ie si l'amplitude est mesurée par a = Dt N-1 T / Tt N-1 T (D vecteur de données et T profil de la raie)
alors l'estimation de var(a) utilise actuellement: var(a) = 1/ Tt N-1 T = 1/sum(square(T/sigma))
sachant que F = a * Sum(T) et donc var(F) = Square(Sum(T)) * var(a)

à remplacer par var(a) = 1 / Tt var(D - aT)^-1 T ou var(D-aT) est un estimateur de la variance du résidu
(i.e. on remplace le bruit N donné en entrée par le bruit mesuré dans le résidu du fit)

• pour le line model, considérer la variance comme uniforme sur le support de la raie, l'estimateur de la variance est donc la moyenne des carrés du résidu

• pour le template ratio, ou le continu, utiliser un estimateur local de la variance, e.g. moyenne glissante ou par morceau des carrés du résidu

=> à modifier aussi pour comparer les flux dans le linemodel rules