[REDMINE1D-89] [RM-8117] [linemeas] compute flux uncertainty when using full Gaussian fit Created: 13/Jun/23  Updated: 09/Dec/23

Status: In Progress
Project: 1D Redmine
Component/s: None
Affects Version/s: None
Fix Version/s: None

Type: Task Priority: Normal
Reporter: Redmine-Jira Migtation Assignee: Redmine-Jira Migtation
Resolution: Unresolved Votes: 0
Labels: None
Remaining Estimate: Not Specified
Time Spent: Not Specified
Original Estimate: Not Specified


 Description   

Created on 2023-06-08 14:12:15 by Didier Vibert. % Done: 30

The L-BFGS-B C++ implémentation does not return the estimated Hessian at the fitted parameters which could be used to infer the amplitude/flux uncertainty.

We have to compute it with another way, may be using the residual ? or an analytical estimation following for instance (https://www.osapublishing.org/abstract.cfm?URI=ao-46-22-5374) ?



 Comments   
Comment by Redmine-Jira Migtation [ 09/Dec/23 ]

Comment by Didier Vibert on 2023-06-15 10:12:45:
@vlebrun tu peux me donner ta formule magique

Comment by Redmine-Jira Migtation [ 09/Dec/23 ]

Comment by Vincent Le Brun on 2023-07-18 15:38:54:
j'ai repris la publication que tu donnes et la formule est en gros la meme que celle que j'avais pour la largeur équivalente, et ça donnerait

Var(flux) = 3*sqrt(pi)*error*sigma(pix)

ou error est l'erreur sur le signal (on peut prendre la valeur moyenne sur le support de la raie en 1ere approximation) et sigma(pix) le sigma de la gaussienne exprimé en pixels

Comment by Redmine-Jira Migtation [ 09/Dec/23 ]

Comment by Didier Vibert on 2023-12-01 13:21:39:
Vincent Le Brun wrote in #note-4:
> j'ai repris la publication que tu donnes et la formule est en gros la meme que celle que j'avais pour la largeur équivalente, et ça donnerait
>
> Var(flux) = 3*sqrt(pi)*error*sigma(pix)
>
> ou error est l'erreur sur le signal (on peut prendre la valeur moyenne sur le support de la raie en 1ere approximation) et sigma(pix) le sigma de la gaussienne exprimé en pixels

ok. juste une remarque, la formule du papier correspond au fit d'une Gaussienne seule, c-a-d sans constante ni polynomes de degré >0 ajouté. Or c'est ce qu'on fait et il y a certainement une dépendance de la variance du Flux à ça...

En attendant mieux, je code ça. On verra si on est content avec.

Generated at Sat Feb 10 15:29:12 JST 2024 using Jira 8.3.4#803005-sha1:1f96e09b3c60279a408a2ae47be3c745f571388b.