[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: |
| Comment by Redmine-Jira Migtation [ 09/Dec/23 ] |
|
Comment by Vincent Le Brun on 2023-07-18 15:38:54: 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: 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. |