[REDMINE1D-298] [RM-8259] SNR négatif under lines Created: 07/Sep/23  Updated: 07/Nov/23  Resolved: 07/Nov/23

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


 Description   

Created on 2023-09-06 12:35:46 by Vincent Le Brun. % Done: 100

dans l'échantillon /net/CESAM/amazed/validation_tests/release_0.44/output_pfc_0.44-RC5 dans un des objets classé quasar ( 00232106 ) le paramètre StrongEmissionLinesSNR est négatif (ce qui ne devrait pas pouvoir exister), je ne sais pas si c'est un problème librairie ou vizu.

MR: https://gitlab.lam.fr/CPF/cpf-redshift/-/merge_requests/559



 Comments   
Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Didier Vibert on 2023-10-17 08:26:49:
a fix to add: the @StrongEmissionLinesSNR@ is duplicated by the output (in same HDF5 dataset ie candidate/model_parameters) @StrongELSNR@. The latter should be removed from result_specification.csv

Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Didier Vibert on 2023-10-17 09:56:24:
pour info: cette quantité @StrongEmissionLinesSNR@ est une somme de SNR des raies en émission avec un SNR > 8 (hrad-codé) ou dont le SNR est > 80% de la somme de tous les SNR des raies en émission. Le SNR de chaque raie étant calculé comme: Flux intégré du spectre moins le continu / sqrt(somme des pixels variance).

Je trouve ça assez bizarre de faire des sommes de SNR, ça n'a pas de sens à mon avis...

si on tient à conserver ce type d'info, je propose de calculer non pas la somme des SNR, mais le SNR de la somme des raies, c-a-d somme_des_flux / sqrt(somme_des_pixels_variance), en conservant à peu près la même sélection des raies: SNR individuel > 8, ou SNR individuel > 80% du SNR total (de toutes les raies en émission visibles)

@vlebrun, ça te va ?
on peut aussi virer tout ça !

Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Vincent Le Brun on 2023-10-17 10:01:33:
clairement la somme des SNR n'a aucun sens, et ça aurait un sens de garder le SNR calculé sur les raies strong du linecatalog (ce qui éviterait d'avoir un critère sur le SNR de chaque raie qui a une origine assez trouble pour moi)

Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Didier Vibert on 2023-10-17 10:13:32:
Vincent Le Brun wrote in #note-7:
> clairement la somme des SNR n'a aucun sens, et ça aurait un sens de garder le SNR calculé sur les raies strong du linecatalog (ce qui éviterait d'avoir un critère sur le SNR de chaque raie qui a une origine assez trouble pour moi)

ok du coup je propose 2 sorties:

  • StrongEmissionLinesSNR: SNR calculé sur les raies strong du linecatalog
  • EmissionLinesSNR: SNR calculé sur toutes les raies en émission du linecatalog

ça peut être utile, si on trouve des solutions avec un SNR total plus for que le SNR strong...

Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Vincent Le Brun on 2023-10-17 10:14:25:
en effet, c'est bon pour moi

Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Didier Vibert on 2023-10-17 10:17:13:
@aallaoui ça va sans doute conflicter un peu avec ton issue #8312

Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Didier Vibert on 2023-10-17 10:35:42:
"StrongEmissionLinesSNR" , strict ou pas ?
dans le code actuel, pour chaque SNR individuel on agrège les raies qui se chevauchent, et par exemple pour Euclid, on comptera NII dans le SNR Ha.

Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Vincent Le Brun on 2023-10-17 10:36:59:
vu que c'est rien d'officiel, on va au plus simple, si les raies se chevauchent ça n'aurait pas forcément de sens de séparer..

Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Ali Allaoui on 2023-10-17 12:18:43:
Didier Vibert wrote in #note-10:
> @aallaoui ça va sans doute conflicter un peu avec ton issue #8312

c'est urgent ? parce que je peux aussi modifier légèrement 8312 pour ne rien casser au monoobs et la valider maintenant
J'allais justement commencer cette partie, @dvibert tu peux vérifier que 8312 en l'état ne conflicte pas trop avec cette issue ? et si on ne merge pas que je travaille sur des éléments qui ne conflictent pas

Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Didier Vibert on 2023-10-17 13:17:32:
Ali Allaoui wrote in #note-13:
> Didier Vibert wrote in #note-10:
> > @aallaoui ça va sans doute conflicter un peu avec ton issue #8312
>
> c'est urgent ? parce que je peux aussi modifier légèrement 8312 pour ne rien casser au monoobs et la valider maintenant
> J'allais justement commencer cette partie, @dvibert tu peux vérifier que 8312 en l'état ne conflicte pas trop avec cette issue ? et si on ne merge pas que je travaille sur des éléments qui ne conflictent pas

rien de très urgent... je la passerai après toi, et donc je gèrerai les conflits, si tu ne tardes pas trop sur #8312. héhé... t'as pas un peu la pression là ? 😀

Comment by Redmine-Jira Migtation [ 07/Nov/23 ]

Comment by Pierre-yves Chabaud on 2023-10-27 12:22:33:
Merged into @develop@ (@05b54892@)

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