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-213 1D Redmine <p><em><font color="#505f79"> Created on 2022-11-23 10:42:15 by Mira Sarkis. % Done: 50</font></em></p> <p>Following the addition of qp8+linemodel variation with igmfit activated (#7647): </p> <p>The activation of igmfit for spectrum qp8 11095-58439-0081 degenerates the redshift estimation:</p> <ul> <li>z_ref = 5.198454</li> </ul> <ul> <li>For IgmFit <b>deactivated</b> : z_estimated = 5.198518<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/17115/17115_igmfit_false_qp8.png" style="border: 0px solid black" /></span></li> <li>For IgmFit <b>activated</b> : z_estimated = 5.1737, with Meiksin_Var_curves_5.0.txt used and igmIndex = 3.<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/17116/17116_igmfit_true_qp8.png" style="border: 0px solid black" /></span></li> </ul> <p>Remark from Vincent: <a href="#note-7]">[https://projets.lam.fr/issues/7647#note-7]</a>: </p> <p><em>Alors il faut vérifier si l'affichage est OK dans le cas igmfit=true, ou si le probleme vient de l'utilisation de la courbe, pour cette courbe, il doit rester 20% du flux juste en dessous de 1215 et là il semble rester à peine 10% (sauf si indice = 3 veut dire qu'on prend la 3eme courbe la plus absorbée??)</em></p> REDMINE1D-213

[RM-7688] unlink continuum igm and Lya igm

Task Normal In Progress Unresolved Redmine-Jira Migtation Redmine-Jira Migtation Wed, 5 Jul 2023 05:33:32 +0000 Mon, 5 Feb 2024 18:29:19 +0000 0 1 <p>Comment by Vincent Le Brun on 2023-04-26 11:57:39:<br/> same remark, check the Ly-alpha reference wavelength...</p> <p>Comment by Vincent Le Brun on 2023-04-27 10:10:45:<br/> après la discussion en réunion, il faut tester des paramètres d'IGM différents pour la raie et le continu</p> <p>Comment by Didier Vibert on 2024-01-11 15:11:12:<br/> Vincent Le Brun wrote in #note-4:<br/> &gt; après la discussion en réunion, il faut tester des paramètres d'IGM différents pour la raie et le continu</p> <p>@vlebrun tu pourrais nous ré-expliquer pourquoi l'extinction de l'igm ne serait pas la même pour le continu et Lya ? </p> <p>je laisse en planif car à mon avis pas trop de contrainte avec le multi-obs, en plus je viens de vérifier: le code est buggé ! et l'index igm est refité systématiquement sur Lya, même si le continu est non nul et inclu l'igm... il faudrait s'assurer que ce bug était déjà là avec l'exemple donné ! </p> <p>cf <a href="https://gitlab.lam.fr/CPF/cpf-redshift/-/blob/0.46-RC4/RedshiftLibrary/src/lib/linemodel/abstractfitter.cpp?ref_type=tags#L194" class="external-link" rel="nofollow">https://gitlab.lam.fr/CPF/cpf-redshift/-/blob/0.46-RC4/RedshiftLibrary/src/lib/linemodel/abstractfitter.cpp?ref_type=tags#L194</a> <br/> @isSymIgm()@ devrait être remplacé par @! isSymIgmFit()@ </p> <p>Comment by Didier Vibert on 2024-01-12 15:56:24:<br/> après vérif: </p> <ul> <li>Le bug mentionné ci-dessus #note-7 n'était pas présent lors de la découverte du pb par Mira. Donc à l'époque (au cours du dev de la 0.40) l'igm Lya était bien égal à celui du continu si disponible.</li> </ul> <ul> <li>Il est probable que ce soit moi, lors d'un refacto qui ait introduit le bug qui détruit donc involontairement le lien.</li> </ul> <p>Malheureusement, celà n'améliore pas vraiment la situation:<br/> <img src="https://pfspipe.ipmu.jp/jira/images/icons/attach/noimage.png" imagetext="Capture%20d%E2%80%99%C3%A9cran%20du%202024-01-12%2014-45-18.png" align="absmiddle" border="0" /></p> <p>comparé à la situation précédente où l'indice igm du continu était 3 (absorption moyenne) qui donnait une extinction sur le continu et la raie Lya de 10% à z=5 &amp; lambda=1216 angströms on a maintenant, avec le fit de l'igm libre sur Lya, l'extinction la plus faible possible, en choisissant la première courbe (indice 1), mais qui correspond quand même à une extinction de 13,6% à z=5 &amp; lambda=1216 angströms (donc 13,6% du flux de la raie et 10% du continu).</p> <p>Bizarrement, le seul fit qui marchait bien, montré par Mira, c'est une absorption du continu seul et pas d'absorption du tout de la raie, ie un profil Gaussien symétrique posé sur un continu avec un break...</p> <p>Bref, en attendant, il faut quand même corriger le bug et ajouter un param booleen pour laisser le fit de l'igm libre sur Lya même avec un igm fitté sur un continu significatif.<br/> Ensuite je ne sais pas ce qu'il faut faire pour faire rentrer ce cas dans le moule ! </p> <p>Comment by Didier Vibert on 2024-01-12 15:59:33:<br/> ça ne serait pas tout simplement la LSF qui n'est pas bonne, ie trop étroite, le spectre qp8 c'est un spectre qui vient d'où ? SDSS ? </p> <p>Comment by Vincent Le Brun on 2024-01-12 16:09:36:<br/> vu que c'est un spectre de quasar, je ne pense pas que la LSF soit le probleme. Est ce que la vitesse est au max possible ?<br/> Cela dit je n'avais pas réalisé que c'était un QSO, la précision attendue n'est pas la même et le profils peuvent varier. Et il peut y avoir des absorptions indépendantes qui modifient le profil et peuvent donner l'impression que le fit n'est pas parfait. On s'en sort bien vu que ça a permis de détecter un bug mais au final je pense qu'on peut fermer l'issue et laisser ce cas à ses imperfections...</p> <p>Comment by Didier Vibert on 2024-01-25 09:22:25:<br/> le bug identifié #note-7 est corrigé dans #8475 : ie l'igm sur la raie Lya est celui du continu, sans fit sur la raie, si le continu est significatif et l'option igmfit est true.</p> <p>reste à ajouter un param pour defaire le lien volontairement, ie fitter un coeff igm différent sur le profil Lya, même lorsque le continu permet d'estimer indépendamment l'igm.<br/> On pourra dans la foulée modifier les param Lya:</p> <ul> <li>créer une section dédiée @lya@</li> <li>déplacer et renommer la section @lyafit@ en @lya.asymProfile@</li> <li>déplacer et renommer @lyaforcefit@ en @lya.asymProfile.alwaysFit@ (ie lorsque le profil est aymfixed on le force en asymfit)</li> <li>déplacer et renommer @lyaforcedisablefit@ en @lya.asymProfile.neverFit@ (ie lorsque le profil est asymfit on le force en asymfixed)</li> <li>renommer @igmfit@ en @lya.igmProfile@ (true: switch from asymmetric Skew-Gaussian profile to Gaussian profile with igm extinction)</li> <li>ajouter @lya.igmFromContinuum@ (true/false, if false igm is fitted independently on continuum and Lya profile)</li> </ul> <p>pour le check des param, même situation que #8435, ie @asymProfile@ (old @asmyfit@) section necessary only when @igmProfile@ is false, and add the following condition @igmFromContinuum@ required when @igmProfile@ is true (useless otherwise).</p> <p>rmq: une autre possibilité est d'introduire @lya.profile (= asym | igm)@ et @lya.igmProfile@ devient une section qui contient juste @lya.igmProfile.igmFromContinuum@. </p> <p>Comment by Fanny Dufresne on 2024-02-05 14:03:04:<br/> Je propose de faire la partie renommage dans le ticket #8435, et de ne garder dans ce ticket que la partie ajout d'un paramètre pour fitter un coeff igm différent sur le profil Lya.</p> <p>Je propose également comme renommage pour:</p> <ul class="alternate" type="square"> <li>@lyaforcefit@ : @lya.asymProfile.switchFixedToFit@ plutôt que @lya.asymProfile.alwaysFit@</li> <li>@lyadisablefit@: @lya.asymProfile.switchFitToFixed@ plutôt que @lya.asymProfile.neverFit@</li> </ul> <p>Comment by Didier Vibert on 2024-02-05 15:35:18:<br/> Fanny Dufresne wrote in #note-14:<br/> &gt; Je propose de faire la partie renommage dans le ticket #8435, et de ne garder dans ce ticket que la partie ajout d'un paramètre pour fitter un coeff igm différent sur le profil Lya.<br/> &gt; <br/> &gt; Je propose également comme renommage pour:<br/> &gt; - @lyaforcefit@ : @lya.asymProfile.switchFixedToFit@ plutôt que @lya.asymProfile.alwaysFit@<br/> &gt; - @lyadisablefit@: @lya.asymProfile.switchFitToFixed@ plutôt que @lya.asymProfile.neverFit@</p> <p>parfait ça me va ! </p> <p>on reste aussi sur la proposition @lya.profile (= asym | igm)@ et les deux sections @lya.igmProfile@ @lya.asymProfile@</p> Development Rank 0|zzsy9r: