[REDMINE1D-371] [RM-8541] Valeurs snrOII problématiques Created: 14/Dec/23  Updated: 07/Feb/24  Resolved: 07/Feb/24

Status: Done
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: Done Votes: 0
Labels: None
Remaining Estimate: Not Specified
Time Spent: Not Specified
Original Estimate: Not Specified

Attachments: PNG File 00010-00001-1,1-0099546300001228-006-0x6e85dbfdd0e3cf63_spec.png    

 Description   

Created on 2023-12-13 09:00:50 by Vincent Le Brun. % Done: 100

dans les données EDR2 PFS, il y a plusieurs spectres ou la valeur calculée de snrOII est clairement fausse
/net/CESAM/amazed/aallaoui/pfs_weekly/edr/0.46.0/edr2_20231203_bis/
objet 00010-00001-1,1-0099546300001228-006-0x6e85dbfdd0e3cf63
snrOII est donné à 722 alors qu'on a ça



 Comments   
Comment by Redmine-Jira Migtation [ 07/Feb/24 ]

Comment by Didier Vibert on 2023-12-13 18:42:37:
Bon, une explication possible.
Les SNR provenant du fit Gaussien en mode ratio ne sont pas très exploitables. En effet, ils sont basés sur l´incertitude de la mesure d'une amplitude unique (celle du modèle avec toutes les raies fixées), qui est assez contrainte et donc précise. Le SNR est donc trompeur. Si les ratios étaient justes, ce SNR serait bon. Il faudrait avoir des variances sur les ratios pour retrouver quelque chose de plus plausible.
En attendant le SNR directIntegration est correct: 2.2

Comment by Redmine-Jira Migtation [ 07/Feb/24 ]

Comment by Vincent Le Brun on 2023-12-14 07:51:24:
du coup je suggère de ne pas les calculer et encore moins de les afficher. Cela dit certains flux ont l'air très faux aussi (cf #8544)

Comment by Redmine-Jira Migtation [ 07/Feb/24 ]

Comment by Didier Vibert on 2023-12-15 10:57:00:
Après vérif, les SNR (en mode template ratio) de toutes les raies sont bien les mêmes (ils ne sont pas scalés par les amplitudes des ratios, car ce scaling intervient au numérateur et dénominateur du SNR, ie l'amplitude et son incertitude sont scalés par le ratio). Si on run en extended results, les ratios amplitude/amplitude_uncertainties sont bien les mêmes et correspondent au ratio amplitude_globale/amplitude_globale_uncertainty.
D'ailleurs cette dernière quantité n'est pas sortie, je rajoute du code pour ça.
En revanche les snr en flux ne donnent pas la même chose, car la largeur est dépendante de la longueur d'onde, même avec une LSF constante (sigma constant) à cause du terme dispersion de vitesse qui est constant en vitesse et donne donc des largeurs en longueur d'onde dépendantes de la longueur d'onde.

Ensuite à voir si on veut conserver ces quantités individuelles par raies dans le cas des templates-ratio, vu qu'elles n'ont pas grand sens... en revanche le SNR venant du fit global de toutes les raies a un sens.

Comment by Redmine-Jira Migtation [ 07/Feb/24 ]

Comment by Vincent Le Brun on 2023-12-15 11:25:23:
je suis partisan de ne pas calculer ces valeurs individuelles vu qu'elles n'ont pas de sens physique

Comment by Redmine-Jira Migtation [ 07/Feb/24 ]

Comment by Didier Vibert on 2023-12-19 16:36:25:

  • J'ai supprimé des output (mis à NAN) les incertitudes des flux par raie dans le cas tplratio.
  • J'ai ajouté l'incertitude et le SNR de l'amplitude globale du tplratio
  • J'ai ajouté le SNR par raie dans les cas rules ou free (il y avait déjà le flux et l'incertitude), mais ça évite de faire le calcul et permet des filtrages avec vizu
Comment by Redmine-Jira Migtation [ 07/Feb/24 ]

Comment by Didier Vibert on 2023-12-21 12:15:13:
MR: https://gitlab.lam.fr/CPF/cpf-redshift/-/merge_requests/571

Comment by Redmine-Jira Migtation [ 07/Feb/24 ]

Comment by Didier Vibert on 2024-01-26 10:32:41:
some small modifications suggested by @aallaoui

Comment by Redmine-Jira Migtation [ 07/Feb/24 ]

Comment by Pierre-yves Chabaud on 2024-01-31 16:04:11:
Merged into @develop@ (@c84e9cf8@)

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