https://pfspipe.ipmu.jp/jira This file is an XML representation of an issue en-us 8.3.4 803005 13-09-2019 https://pfspipe.ipmu.jp/jira/browse/REDMINE1D-89 1D Redmine <p><em><font color="#505f79"> Created on 2023-06-08 14:12:15 by Didier Vibert. % Done: 30</font></em></p> <p>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.</p> <p>We have to compute it with another way, may be using the residual ? or an analytical estimation following for instance (<a href="https://www.osapublishing.org/abstract.cfm?URI=ao-46-22-5374" class="external-link" rel="nofollow">https://www.osapublishing.org/abstract.cfm?URI=ao-46-22-5374</a>) ?</p> REDMINE1D-89

[RM-8117] [linemeas] compute flux uncertainty when using full Gaussian fit

Task Normal In Progress Unresolved Redmine-Jira Migtation Redmine-Jira Migtation Tue, 13 Jun 2023 04:58:25 +0000 Fri, 8 Dec 2023 18:34:17 +0000 0 2 <p>Comment by Didier Vibert on 2023-06-15 10:12:45:<br/> @vlebrun tu peux me donner ta formule magique </p> <p>Comment by Vincent Le Brun on 2023-07-18 15:38:54:<br/> 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</p> <p>Var(flux) = 3*sqrt(pi)*error*sigma(pix)</p> <p>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</p> <p>Comment by Didier Vibert on 2023-12-01 13:21:39:<br/> Vincent Le Brun wrote in #note-4:<br/> &gt; 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<br/> &gt; <br/> &gt; Var(flux) = 3*sqrt(pi)*error*sigma(pix)<br/> &gt; <br/> &gt; 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</p> <p>ok. juste une remarque, la formule du papier correspond au fit d'une Gaussienne seule, c-a-d sans constante ni polynomes de degré &gt;0 ajouté. Or c'est ce qu'on fait et il y a certainement une dépendance de la variance du Flux à ça... </p> <p>En attendant mieux, je code ça. On verra si on est content avec.</p> Development Rank 0|zzsxd3: