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-208 1D Redmine <p><em><font color="#505f79"> Created on 2023-02-17 10:12:15 by Vincent Le Brun. % Done: 100</font></em></p> <p>J'ai commencé à travailler sur les spectres incidents EL-COSMOS (/net/CESAM/amazed/vlebrun/EL-COSMOS2020/IncidentSpectra/v3 pour les spectres) et en faisant les tests techniques sur les deux premiers (Id100000_F.fits et Id100001_F.fits) j'ai un comportement étrange :</p> <p>Le linemodel free(répertoire test_LMfree2) marche parfaitement (ce qui est bien le moins vu le niveau des raies) meme si le fit du continu n'est pas bon et trouve les bons redshifts (0.0892 et 0.3467), et donne des mesures de flux correctes. </p> <p>Par contre en tplratio (répertoires test_LMRatio* avec diverses valeurs de la vitesse en emission initiale, de la LSF, du half_window_size) il y a au moins un des deux spectres pour lequel la solution n'est pas bonne, mais pas de loin. Ces spectres sont un peu particuliers dans le sens ou le continu est échantillonné 'à la Euclid' alors que les raies en émission sont bien échantillonnées, mais je ne comprend pas comment un modèle avec aucune des raies en émission peut mieux marcher qu'un tplratio (meme imparfait) au bon redshift , vu que le redshift est presque bon le Khi2 dû au continu ne doit pas varier énormément, et celui dû aux raies en émission ne peut qu'être meilleur ?<br/> Toute idée est bienvenue...</p> REDMINE1D-208

[RM-7882] Pourquoi ?

Task Normal Done Done Redmine-Jira Migtation Redmine-Jira Migtation Wed, 5 Jul 2023 02:44:07 +0000 Mon, 27 Nov 2023 19:11:41 +0000 Wed, 22 Nov 2023 18:55:24 +0000 0 1 <p>Comment by Vincent Le Brun on 2023-03-08 11:56:03:<br/> j'ai fait tourner sur les 700000 spectres et il y a 20% d'erreurs du meme style : raies en émission très fortes mais mismatch à l'arrivée<br/> exemple dans /net/CESAM/amazed/vlebrun/EL-COSMOS2020/IncidentSpectra/outout_0.40_tplratio2<br/> le spectre Id1000320 ou une raie est choisie pour HeI3190<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16875/16875_HeI.png" style="border: 0px solid black" /></span><br/> alors que toutes les raies en émission sont présentes, et la deuxième solution est proche de la bonne (sans doute grec au continu), mais les raies en émission ne sont pas utilisées...<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16876/16876_Sol1.png" style="border: 0px solid black" /></span><br/> Ca voudrait dire que nos tplratio sont tellement loins de ce spectre que ça fit mieux quand on ne prend aucune raie? je n'arrive meme pas à voir comment c'est possible</p> <p>Comment by Didier Vibert on 2023-03-08 13:32:54:<br/> je viens de regarder en direct sous vizu ce spectre, <del>ce qui est bizarre c'est que le modèle n'est pas disponible (pourtant il l'est sur d'autres spectres du même run...) pb vizu</del> ? </p> <p>Et sinon sans doute que le SNR (Fline/Fcontinuum) très élevé des raiess contribue à rendre le fit très sensible aux ratios ? </p> <p>Comment by Didier Vibert on 2023-03-08 13:54:59:<br/> Didier Vibert wrote in #note-4:<br/> &gt; je viens de regarder en direct sous vizu ce spectre, ce qui est bizarre c'est que le modèle n'est pas disponible (pourtant il l'est sur d'autres spectres du même run...) pb vizu ? <br/> &gt; <br/> &gt; Et sinon sans doute que le SNR (Fline/Fcontinuum) très élevé des raiess contribue à rendre le fit très sensible aux ratios ?</p> <p>Didier Vibert wrote in #note-4:<br/> &gt; je viens de regarder en direct sous vizu ce spectre, <del>ce qui est bizarre c'est que le modèle n'est pas disponible (pourtant il l'est sur d'autres spectres du même run...) pb vizu</del> ? <br/> &gt; <br/> &gt; Et sinon sans doute que le SNR (Fline/Fcontinuum) très élevé des raiess contribue à rendre le fit très sensible aux ratios ?</p> <p>y a vraiment un problème: les 2 raies les plus intenses du modèle au best z sont HeI_3190 = 2.92E-18 et MgII = 1.44E-18 donné par la table des flux or on ne voit pas cette dernière ni dans le modèle, ni dans le spectre... et pourtant on voit le double OIII dans le modèle, qui est donné avec un flux de 5.6E-20...</p> <p>Comment by Vincent Le Brun on 2023-03-08 16:52:07:<br/> le SNR est arbitrairement de 50 par pixel, mais en écrivant je réalise que les pixels n'ont pas la meme taille dans le continu et dans les raies, et ça doit être ça qui perturbe le fit. Et ça peut peut-être aussi expliquer les problèmes de flux, si une raie n'est pas à sa place elle peut se retrouver toute entière dans un pixel et avoir une mesure de flux HS</p> <p>Comment by Vincent Le Brun on 2023-03-08 18:43:47:</p> <p>Cela dit il y a aussi des spectres qui ne posent pas de probleme et je ne vois pas vraiment la différence entre les deux<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16877/16877_Id1000323.png" style="border: 0px solid black" /></span></p> <p>Comment by Vincent Le Brun on 2023-03-08 19:01:16:<br/> mais ce genre de fit reste quand meme compliqué à expliquer<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16878/16878_Lyalpha.png" style="border: 0px solid black" /></span></p> <p>Comment by Vincent Le Brun on 2023-03-09 10:47:42:<br/> le même avec le templatefitting. Toujours pas <br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16879/16879_Id1000039_TPL.png" style="border: 0px solid black" /></span></p> <p>et l'autre (Id1000320), ou le templatefiting trouve la bonne solution, mais en 3eme<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16881/16881_Id1000320_Tpl_sol0.png" style="border: 0px solid black" /></span> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16880/16880_Id1000320_Tpl_sol2.png" style="border: 0px solid black" /></span></p> <p>Comment by Vincent Le Brun on 2023-03-09 13:57:58:<br/> le SNR n'est pas le pb, avec SNR de 5 <br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16882/16882_Id1000039_SNR5.png" style="border: 0px solid black" /></span><br/> et pour l'autre la bonne solution est meme passée en 3 alors que la totalité du spectre est compatible avec le modèle (sauf dans les raies)<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16883/16883_Id1000320_SNR5_sol3.png" style="border: 0px solid black" /></span></p> <p>Comment by Mira Sarkis on 2023-03-23 14:29:17:<br/> En échantillonant le spectre Id1000039 avec une résolution constante = 13, on retrouve bien la bonne solution avec linemodel: <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16884/16884_Id1000039_linemodelSolve_Res13.png" style="border: 0px solid black" /></span></p> <p>Ceci n'est toujours pas le cas avec TemplateFitting</p> <p>Comment by Didier Vibert on 2023-03-23 14:39:56:<br/> Mira Sarkis wrote in #note-12:<br/> &gt; En échantillonant le spectre Id1000039 avec une résolution constante = 13, on retrouve bien la bonne solution avec linemodel:<br/> &gt; <br/> le continu est quand même très étonnant...</p> <p>&gt; Ceci n'est toujours pas le cas avec TemplateFitting</p> <p>Comment by Mira Sarkis on 2023-03-23 15:13:15:<br/> Didier Vibert wrote in #note-13:<br/> &gt; Mira Sarkis wrote in #note-12:<br/> &gt; &gt; En échantillonant le spectre Id1000039 avec une résolution constante = 13, on retrouve bien la bonne solution avec linemodel:<br/> &gt; &gt; <br/> &gt; le continu est quand même très étonnant...<br/> &gt; <br/> &gt; &gt; Ceci n'est toujours pas le cas avec TemplateFitting</p> <p>Nous avons eu déjà des cas similaires dans <a href="https://projets.lam.fr/issues/7199" class="external-link" rel="nofollow">https://projets.lam.fr/issues/7199</a>. Selon cette dernière ceci est du au fait que le premier sample du spectre est très petit (de l'ordre de 6*10^-27)<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16885/16885_Id1000039_Flux_Res13.png" style="border: 0px solid black" /></span></p> <p>Comment by Didier Vibert on 2023-03-23 16:00:58:</p> <p>Mira Sarkis wrote in #note-14:<br/> &gt; Nous avons eu déjà des cas similaires dans <a href="https://projets.lam.fr/issues/7199" class="external-link" rel="nofollow">https://projets.lam.fr/issues/7199</a>. Selon cette dernière ceci est du au fait que le premier sample du spectre est très petit (de l'ordre de 6*10^-27)</p> <p>ah ok... mais là ça me semble différent. le spectre est nul sur une plage , à cause de l'igm, normalement les modèles devraient être capable de le faire. Est-ce que l'igm (sur le continu) est bien activé ? car on ne l'active pas sur des données euclid, car les spectres rouge démarrent à 12000 ce qui rend l'igm visible seulement à grand z (7 ou 8 ), là les spectre démarrent à 4000.</p> <p>Comment by Vincent Le Brun on 2023-03-23 16:07:37:<br/> j'ai fait commencer le lambda range à 4000 pour justement éviter la partie concernée par l'IGM, <del>comme ce sont des spectres à z=0</del> c'était impossible à gérer (et en fait la coupure qu'on voit se fait à 912A à l'intérieur de la galaxie, ça correspond aux photons qui ont l'énergie minimale pour ioniser l'hydrogène).</p> <p>Comment by Didier Vibert on 2023-03-23 16:18:16:<br/> Vincent Le Brun wrote in #note-16:<br/> &gt; j'ai fait commencer le lambda range à 4000 pour justement éviter la partie concernée par l'IGM, comme ce sont des spectres à z=0 c'était impossible à gérer (et en fait la coupure qu'on voit se fait à 912A à l'intérieur de la galaxie, ça correspond aux photons qui ont l'énergie minimale pour ioniser l'hydrogène).</p> <p>ok, donc sur le spectre montré par Mira, on ne devrait pas utiliser ce qui est est à gauche de 4000, c'est ça ? pourtant j'ai bien l'impression qu'on voit le spectre à 912... idem sur ton plot #note-10 </p> <p>Comment by Vincent Le Brun on 2023-03-23 16:22:41:<br/> non ça commence bien à 4000 pour contre les spectres ont un redshift, ce ne sont pas des modèles à z=0</p> <p>Comment by Didier Vibert on 2023-03-23 17:46:03:<br/> Vincent Le Brun wrote in #note-18:<br/> &gt; non ça commence bien à 4000 pour contre les spectres ont un redshift, ce ne sont pas des modèles à z=0</p> <p>donc ok ça démarre bien à 4000. Mais le spectre en question est nul à gauche de 5000, et démarre à 700, comment c'est possible ? 5000/4000=1.25 c'est un spectre à z=0.25 ? <br/> bref en tout cas la partie gauche nulle du spectre est gênante, d'autant plus qu'il ne s'agit pas de l'igm...</p> <p>Comment by Vincent Le Brun on 2023-03-23 18:07:40:<br/> pour le spectre Id0000039 c'est bien l'IGM en plus de la coupure&lt;912, la raie à 6700 est Lyman-alpha, donc z=4,5 environ, et à ce redshift là les courbes d'IGM sont toutes à 0en dessous de 912 A donc ça ne devrait pas poser de pb</p> <p>Comment by Mira Sarkis on 2023-03-24 09:22:32:<br/> Didier Vibert wrote in #note-15:</p> <p>&gt;Est-ce que l'igm (sur le continu) est bien activé ? car on ne l'active pas sur des données euclid, car les spectres rouge démarrent à 12000 ce qui rend l'igm visible seulement à grand z (7 ou 8 ), là les spectre démarrent à 4000.</p> <p>L'igm est bien activé sur le continu dans les params.json de Vincent.</p> <p>Voici le plot du modèle de spectre avec IGM désactivé sur le continu, pour un redshift estimé à 7.2010<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16886/16886_Id1000039_linemodelSolve_Res13_noIGM.png" style="border: 0px solid black" /></span></p> <p>Comment by Mira Sarkis on 2023-03-31 09:01:48:<br/> Je vous joins le plot du modèle avec TemplateFitting <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16887/16887_Id1000039_TF_Res13.png" style="border: 0px solid black" /></span></p> <p>Comment by Vincent Le Brun on 2023-03-31 09:16:19:<br/> donc la question de départ 'Pourquoi ?' demeure... est ce que tu peux mettre la PDF du tplfit ?</p> <p>Comment by Mira Sarkis on 2023-03-31 09:21:53:<br/> Vincent Le Brun wrote in #note-23:<br/> &gt; donc la question de départ 'Pourquoi ?' demeure... est ce que tu peux mettre la PDF du tplfit ?<br/> La pdf est aussi bizarre<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16888/16888_Id1000030_TF_Res13_PDF.png" style="border: 0px solid black" /></span></p> <p>Comment by Vincent Le Brun on 2023-03-31 09:38:48:<br/> comme par hasard la PFD plonge autour de 4.3 qui est la valeur attendue...<br/> est ce que le fit d'un spectre qui vaut zéro sur une grande partie par un modèle qui vaut aussi zéro poserait des problèmes (genre si on essaye de calculer une amplitude par une division) ??</p> <p>Comment by Didier Vibert on 2023-03-31 10:49:54:<br/> Vincent Le Brun wrote in #note-25:<br/> &gt; comme par hasard la PFD plonge autour de 4.3 qui est la valeur attendue...<br/> &gt; est ce que le fit d'un spectre qui vaut zéro sur une grande partie par un modèle qui vaut aussi zéro poserait des problèmes (genre si on essaye de calculer une amplitude par une division) ??</p> <p>je ne vois pas trop..., c'est quoi la variance que tu as mise, lorsque le flux est nul ? </p> <p>ça ne serait pas parce que la variance est bcp trop petite, et du coup cette zone du spectre a un poids considérable dans le chi2, or l'extinction igm ne devient pas exactement nulle tout de suite, ce qui l'empêcherait de se caler assez à gauche ? </p> <p>Comment by Didier Vibert on 2023-03-31 10:58:47:<br/> Didier Vibert wrote in #note-26:<br/> &gt; Vincent Le Brun wrote in #note-25:<br/> &gt; &gt; comme par hasard la PFD plonge autour de 4.3 qui est la valeur attendue...<br/> &gt; &gt; est ce que le fit d'un spectre qui vaut zéro sur une grande partie par un modèle qui vaut aussi zéro poserait des problèmes (genre si on essaye de calculer une amplitude par une division) ??<br/> &gt; <br/> &gt; je ne vois pas trop..., c'est quoi la variance que tu as mis, lorsque le flux est nul ? <br/> &gt; <br/> &gt; ça ne serait pas parceque la variance est bcp trop petite, et du coup cette zone du spectre a un poid considérable dans le chi2, hors l'extinction igm ne devient pas exactement nulle tout de suite, ce qui l'empêcherait de se caler assez à gauche ?</p> <p>y a un truc bizarre dans les courbes IGM en dessous de 912:<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16889/16889_clipboard-202303311258-nonkk.png" style="border: 0px solid black" /></span></p> <p>Comment by Vincent Le Brun on 2023-03-31 11:10:23:<br/> bizarre en effet mais pour ce spectre le lambda_rest au bon redshift est d'environ 700 A donc cette zone la ne devrait pas avoir d'impact...<br/> la courbe d'IGM devrait continuer à remonter en dessous de 300A, j'imagine qu'on a coupé en se disant que ça ne servirait jamais...</p> <p>Comment by Didier Vibert on 2023-03-31 12:58:25:<br/> Mira: ça donne quoi si tu mets un lambdarange à 5000 (voire un peu plus) pour supprimer la plage du spectre nulle ? </p> <p>Comment by Mira Sarkis on 2023-03-31 13:15:22:<br/> Didier Vibert wrote in #note-29:<br/> &gt; Mira: ça donne quoi si tu mets un lambdarange à 5000 (voire un peu plus) pour supprimer la plage du spectre nulle ?<br/> J'ai déjà lancé avec lambda_start=6000 et c'est bien mieux comme ca, mais je ne comprend pas comment ca se fait que c'est mieux!</p> <p>Le début du spectre n'est pas vraiment nul, c'est plutot cte = 2*1E-22<br/> <span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/16890/16890_TF_lbdastart6000.png" style="border: 0px solid black" /></span></p> <p>Comment by Didier Vibert on 2023-03-31 14:31:08:<br/> Mira Sarkis wrote in #note-30:<br/> &gt; Didier Vibert wrote in #note-29:<br/> &gt; &gt; Mira: ça donne quoi si tu mets un lambdarange à 5000 (voire un peu plus) pour supprimer la plage du spectre nulle ?<br/> &gt; J'ai déjà lancé avec lambda_start=6000 et c'est bien mieux comme ca, mais je ne comprend pas comment ca se fait que c'est mieux!<br/> &gt; <br/> &gt; Le début du spectre n'est pas vraiment nul, c'est plutot cte = 2*1E-27<br/> ok</p> <p>tu peux essayer avec un lbda range juste après la zone null ie 5000 (6000 c'est déjà pas mal plus loin).... je pense que le pb vient des valeurs de variances dans cette zone (cf #note-26)</p> <p>Comment by Didier Vibert on 2023-03-31 14:35:38:<br/> Mira Sarkis wrote in #note-30:<br/> &gt; J'ai déjà lancé avec lambda_start=6000 et c'est bien mieux comme ca, mais je ne comprend pas comment ca se fait que c'est mieux!</p> <p>y a un truc que je ne comprend pas: entre ce plot lbda_start=6000 #note-30 et le précédent #note-22, la courbe du bruit n'est pas la même ??????</p> <p>Comment by Didier Vibert on 2023-03-31 14:37:41:<br/> Didier Vibert wrote in #note-32:<br/> &gt; Mira Sarkis wrote in #note-30:<br/> &gt; &gt; J'ai déjà lancé avec lambda_start=6000 et c'est bien mieux comme ca, mais je ne comprend pas comment ca se fait que c'est mieux!<br/> &gt; <br/> &gt; y a un truc que je ne comprend pas: entre ce plot lbda_start=6000 #note-30 et le précédent #note-22, la courbe du bruit n'est pas la même ??????</p> <p>quand tu as fait le ré-échantillonnage, tu l'as fait sur le bruit aussi ? </p> <p>Comment by Mira Sarkis on 2023-03-31 14:42:48:<br/> Didier Vibert wrote in #note-32:<br/> &gt; Mira Sarkis wrote in #note-30:<br/> &gt; &gt; J'ai déjà lancé avec lambda_start=6000 et c'est bien mieux comme ca, mais je ne comprend pas comment ca se fait que c'est mieux!<br/> &gt; <br/> &gt; y a un truc que je ne comprend pas: entre ce plot lbda_start=6000 #note-30 et le précédent #note-22, la courbe du bruit n'est pas la même ??????</p> <p>Dans #note-22 il s'agit du spectre dont le flux et la variance sont réechantillonés. <br/> Dans #note-30, c'est le spectre d'origine (avec un zoom sur le flux)</p> <p>Comment by Vincent Le Brun on 2023-03-31 14:43:10:<br/> Didier Vibert wrote in #note-31:</p> <p>&gt; tu peux essayer avec un lbda range juste après la zone null ie 5000 (6000 c'est déjà pas mal plus loin).... je pense que le pb vient des valeurs de variances dans cette zone (cf #note-26)</p> <p>je crois que j'ai compris, c'st à cause du S/N constant : les valeurs de flux sont faibles (1e-27) et donc des variances en 2e-29, mais pas 0, mais du coup ça rend le spectre incompatible avec 0 vu que je n'ai pas mis de bruit de fond constant (qui est toujours présent dans un vrai spectre), et qui permettrait au modele nul de passer par zéro<br/> Je vais refaire une production avec un bruit minimum réaliste et on regarde ce que ça donne</p> <p>Comment by Didier Vibert on 2023-03-31 14:48:32:<br/> Vincent Le Brun wrote in #note-35:<br/> &gt; Didier Vibert wrote in #note-31:<br/> &gt; <br/> &gt; &gt; tu peux essayer avec un lbda range juste après la zone null ie 5000 (6000 c'est déjà pas mal plus loin).... je pense que le pb vient des valeurs de variances dans cette zone (cf #note-26)<br/> &gt; <br/> &gt; je crois que j'ai compris, c'st à cause du S/N constant : les valeurs de flux sont faibles (1e-27) et donc des variances en 2e-29, mais pas 0, mais du coup ça rend le spectre incompatible avec 0 vu que je n'ai pas mis de bruit de fond constant (qui est toujours présent dans un vrai spectre), et qui permettrait au modele nul de passer par zéro<br/> &gt; Je vais refaire une production avec un bruit minimum réaliste et on regarde ce que ça donne</p> <p>+1</p> <p>Comment by Didier Vibert on 2023-03-31 14:49:27:<br/> Mira Sarkis wrote in #note-34:</p> <p>&gt; Dans #note-22 il s'agit du spectre dont le flux et la variance sont réechantillonés. <br/> &gt; Dans #note-30, c'est le spectre d'origine.</p> <p>ça n'explique pas la différence entre les 2 visu du bruit ????</p> <p>Comment by Didier Vibert on 2023-03-31 14:55:33:<br/> Didier Vibert wrote in #note-37:<br/> &gt; Mira Sarkis wrote in #note-34:<br/> &gt; <br/> &gt; &gt; Dans #note-22 il s'agit du spectre dont le flux et la variance sont réechantillonés. <br/> &gt; &gt; Dans #note-30, c'est le spectre d'origine.<br/> &gt; <br/> &gt; ça n'explique pas la différence entre les 2 visu du bruit ????</p> <p>bon si, ça vient des "spikes" dans la variance (notamment celui vers 6700) qui doit disparaitre avec ton rééchantillonnage. </p> <p>Comment by Didier Vibert on 2023-03-31 15:00:11:<br/> Didier Vibert wrote in #note-36:<br/> &gt; Vincent Le Brun wrote in #note-35:<br/> &gt; &gt; je crois que j'ai compris, c'st à cause du S/N constant : les valeurs de flux sont faibles (1e-27) et donc des variances en 2e-29, mais pas 0, mais du coup ça rend le spectre incompatible avec 0 vu que je n'ai pas mis de bruit de fond constant (qui est toujours présent dans un vrai spectre), et qui permettrait au modele nul de passer par zéro<br/> &gt; &gt; Je vais refaire une production avec un bruit minimum réaliste et on regarde ce que ça donne<br/> &gt; <br/> &gt; +1</p> <p>c'est ce que j'essayai de dire #note-26. <br/> Je pense que tu n'as pas forcément besoin de faire une réalisation de bruit de fond, tu peux juste rajouter une constante partout dans la variance, pour lui donner une limite basse. Ensuite c'est pas grave si la partie nulle n'est pas bruitée, elle sera juste un peu trop compatible.</p> <p>Comment by Vincent Le Brun on 2023-04-03 12:29:04:<br/> en effet j'ai rajouté une constante à la variance et pour cet objet là ça marche beaucoup mieux<br/> /net/CESAM/amazed/vlebrun/EL-COSMOS2020/IncidentSpectra/test_noise_LM/<br/> par contre pour le deuxième objet (meme répertoire, Id000320) ça ne marche toujours pas il préfère une solution qui fait tomber toutes les raies au milieu des pixels large, sans doute parce que les line-ratio ne sont pas adapté.<br/> Mira, est ce que tu peux me dire ou est ton script de resampling, je voudrais retenter sur l'échantillon complet (et après on ferme le ticket)</p> <p>Comment by Mira Sarkis on 2023-04-03 13:06:57:<br/> Vincent Le Brun wrote in #note-40:<br/> &gt; en effet j'ai rajouté une constante à la variance et pour cet objet là ça marche beaucoup mieux<br/> &gt; /net/CESAM/amazed/vlebrun/EL-COSMOS2020/IncidentSpectra/test_noise_LM/<br/> &gt; par contre pour le deuxième objet (meme répertoire, Id000320) ça ne marche toujours pas il préfère une solution qui fait tomber toutes les raies au milieu des pixels large, sans doute parce que les line-ratio ne sont pas adapté.</p> <p>&gt; Mira, est ce que tu peux me dire ou est ton script de resampling, je voudrais retenter sur l'échantillon complet (et après on ferme le ticket)<br/> Je viens de rajouter un jupyter notebook sur le cluster: /net/CESAM/amazed/msarkis/resample_spectra_cteRes.ipynb<br/> Attention le resampling se fait avec une resolution = 13A, mais tu peux la changer.</p> <p>Je te joins aussi le spectre Id1000320 resamplé.</p> <p>Comment by Vincent Le Brun on 2023-11-22 17:07:03:<br/> Après échantillonnage des spectres et mise au point des paramètres (quand les raies sont très étroites un échantillonnage fin en z est nécessaire) on arrive à des résultats satisfaisants. On peut fermer ce ticket, pas de code associé</p> Development Rank 0|zzsy8n: