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

Wed, 5 Jul 2023 17:29:21 +0000

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

PFS-JIRA

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