Sommaire

Comment lancer cosim

Le contenu global du parfile

%I.0 L'identificateur principal
%I Un identificateur secondaire
%I.R Une liste de raw-id
%J Les coordonnées
%M ou %F La magnitude
%V La vitesse
%VLSR La VLSR
%C Le type d'objet principal
%CL Le type d'objet secondaire
%H.p Un parent
%H.c Un enfant
%G ou %GDIM Les dimensions
%P Le mouvement propre
%X La parallaxe
%S Le type spectral
%T Le type morphologique
%N(CLS) Une note
%O.(...) Les data
%COM Commentaire
# Commentaire

L'option .SEARCH rayon
L'option .SEARCH COO ou NOCOO
L'option .NOLIMIT
L'option .SIGMA COO
L'option .SIGMA V
L'option .SIGMA M
Les options .COO .V .M ...
Les options .ACCEPT COO ou .REFUSE COO
L'option .OT+ ou .OT-
L'option .ACRO+ ou .ACRO-
L'option .B+ ou .B-
L'option .HIDE
L'option .IGNORE
L'option .TRUST
L'option .EXPECT
L'option .OK INCORRECT ID
L'option .OK Already Connected
L'option .FORCE FROM ID
Les options .REPLACE NO et .REPLACE EQUAL
L'option .DEBUG

Pour aller plus loin : explication détaillée du fonctionnement de cosim

lancement de cosim

Le programme a besoin d'un fichier d'entrée, le parfile. Ce fichier contient les données extraites de la table.

La ligne de commande à écrire:

cosim parfile

Deux fichiers sont alors créés automatiquement : parfile.out, qui contient entre autres les commandes de mise à jour et parfile.stat, un fichier de statistiques sur les résultats du programme.

Les options ne s'écrivent pas sur la ligne de commande; elles sont écrites dans le parfile.

En cas de récentes modifications qui auraient entraîné un problème, on peut appeler la version précédente de cosim:

cosim-1 parfile

L'outil de visualisation d'un parfile.out qu'Anaïs avait commencé à développer est toujours disponible ici :

http://cacao/RaccordDecoder/WebContent/Decoder.html - Retours bienvenus !

Retour au sommaire

écriture du parfile

le contenu global du parfile

Le fichier parfile doit contenir :

• une en-tête constituée de :
  • une ligne obligatoire avec le bibcode de la table traîtée (si pas de référence à rajouter, remplacer cette ligne par n'importe quelle ligne de commentaire commençant par #)

    .B 2012A&A...544A..81H

    ou

    #no ref

    
    Pour l'instant, on ne peut mettre qu'une seule ligne .B

• • des lignes pour les options (une ligne par option), commençant toujours par un point.

    .OT+1 HII,EmO,PN?,PN
    .SIGMA COO 0,1000

• puis une ligne vide

• puis un paragraphe par objet, avec une ligne pour chaque champ, commençant toujours par un %. Seule la ligne %I.0 est obligatoire.

  %I.0 [SCM2012b] HII001
  %J 0 37 24.12 +40 17 56.2 (Opt) C 2012ApJ...758..133S
  %V v:s -483.6 (Opt) D 2012ApJ...758..133S
  %O.velocities |v| -483.6 | |s,,Opt, | | | |2012ApJ...758..133S|
  %C HII 2012ApJ...758..133S

les champs

l'identificateur principal : %I.0 id

%I.0 IRAS 02074+4541

C'est le seul champ obligatoire. Si l'identificateur de %I.0 est mal écrit, il n'y aura pas d'édition de commandes de mise à jour.

Comme en mise à jour SIMBAD, on peut ajouter un "!" derrière l'identificateur pour chercher le nom exact seulement.

Si il n'y a pas de ligne %I.R, ce sera le contenu de %I.0 qui sera mis en raw-id.

Options qui touchent à ce champ : .OK INCORRECT ID, .ACRO+ ou .ACRO-, (.FORCE FROM ID)

Retour au sommaire

un identificateur secondaire : %I id

%I HIP 10245
%I 2MASS J02114684+4555271
%I NLTT 7250

On peut mettre autant de lignes %I qu'il y a d'identificateurs secondaires à entrer. Attention, les lignes %I ne sont pas reportées dans le raw-id.

Options qui touchent à ce champ : .OK INCORRECT ID, .ACRO+ ou .ACRO-, (.FORCE FROM ID)

la liste des identificateurs à faire apparaître dans le raw-id : %I.R id1;id2

A remplir si différent de %I.0, ou si on veut faire apparaître la liste entière des identificateurs tels qu'ils sont nommés dans la table.

les coordonnées : %J coo (dom) [ellipse d'erreur] Q bibcode

Le format (sexagésimal, degrés décimaux) est automatiquement reconnu.

%J 3 32 38.96 -27 42 43.7 (Opt) C 2013A&A...549A..63K
%J 189.34850810 +62.40285870 C 2011A&A...533A.119E

Pour des coordonnées galactiques, ajouter un G devant les coordonnées :
%J G018.261 -00.297 (IR) D 2012MNRAS.424.2442S

Pour des coordonnées en B1950, ajouter au début du fichier (avec les options) :
.EQ B1950
Attention! Il y a un bug non encore corrigé au 4/3/2021 : COSIM tient bien compte de l'option EQ B1950 pour chercher les objets dans SIMBAD, mais plus pour la mise à jour (il rentre des coordonnées B1950 comme si c'était des J2000).

N.B. : règle pour SIMBAD : les coordonnées doivent être autant que possible des coordonnées pour l'époque 2000. Ceci s'applique à tous les objets dont on connaît l'époque de mesure et le mouvement propre. COSIM ne fait pas de calcul des coordonnées selon l'époque. Il faut donc lui donner directement les coordonnées à l'époque 2000. Pour ce faire, utiliser la colonne Vizier des coordonnées recalculées, ou utiliser astropos.

Options qui touchent à ce champ : .SEARCH rayon, .SEARCH COO ou .SEARCH NOCOO, .SIGMA COO, .COO, .ACCEPT COO ou .REFUSE COO, .IGNORE COO, .REPLACE NO COO et .REPLACE EQUAL COO

une magnitude : %M ou %F bande (AB) mag [erreur] Q bibcode

%M B (AB) 25.8 [0.2] D 2013A&A...549A..63K

%M J 11.214 [0.022] C 2003yCat.2246....0C

Options qui touchent à ce champ : .SIGMA M, .M, .IGNORE M, .TRUST M, .EXPECT M, .REPLACE NO M et .REPLACE EQUAL M

Retour au sommaire

le redshift ou la vitesse : %V [vz]:nat valeur (dom) [erreur] Q bibcode

%V z:a 4.3792 (Opt) [0.0064] D 2013A&A...549A..63K

%V v:s 241.2 (Opt) B [1.5] 2011AJ....142...61C

%V cz 400971 (Opt) C [315] 2011ApJ...728...27O

Options qui touchent à ce champ : .SIGMA V, .V, .IGNORE V, .TRUST V, .EXPECT V, .REPLACE NO V et .REPLACE EQUAL V

la VLSR : %VLSR valeur-centrale minimum maximum [erreur] (dom) bibcode

%VLSR -49.7 -49.8 -49.6 [0.07] (Rad) 2019ApJ...884....4K

%VLSR -49.7 ~ ~ [0.07] (Rad) 2019ApJ...884....4K

N.B.: le champ VLSR a été ajouté à COSIM au printemps 2022; à l'été 2022, il y a des choses à améliorer. (Un score VLSR est bien généré, de -1 à +1, et la commande "a VLSR " ou "c VLSR " est bien générée; les options .IGNORE, .TRUST et .EXPECT sont possibles; mais on ne peut pas encore agir sur les limites; les décisions d'écraser l'ancienne valeur ou pas ne sont pas encore tout à fait au point).

Options qui touchent à ce champ : .IGNORE VLSR, .TRUST VLSR, .EXPECT VLSR, .REPLACE NO VLSR et .REPLACE EQUAL VLSR

le type d'objet principal : %C type bibcode

D'après un tableau de comparaison des types d'objets entre eux, propre à SIMBAD, le type d'objet principal de SIMBAD sera écrasé ou non.

Si le type d'objet principal ne donne pas lieu à la commande c o , il est systématiquement ajouté en type secondaire (commande a o ).

Options qui touchent à ce champ : .OT, .OT+ ou .OT-, .HIDE, .IGNORE OT

un type d'objet secondaire : %CL type bibcode

%CL IR 2012ApJ...761..140C

Ce type secondaire sera systématiquement ajouté (commande a o ).

On peut mettre plusieurs %CL par objet.

un parent : %H.p id_parent!, % [BS] bibcode

%H.p NAME LMC!, 100 B 2011AJ....142...61C

Retour au sommaire

un enfant : %H.c id_enfant!, % [BS] bibcode

les dimensions d'un objet étendu : %G ou %GDIM unit maj-axis min-axis angle lgOnde Q bibcode

le mouvement propre : %P ra dec Q [ellipse d'erreur] bibcode

la parallaxe : %X valeur Q [erreur] bibcode

%X 76.86 C [3.97] 2012ApJ...752...56F

le type spectral : %S type Q bibcode

le type morphologique : %T type Q bibcode

une note attachée à l'objet SIMBAD : %N(CLS)

%NC Ceci est une Note Courte (note publique visible en haut de la page web).

%NL Ceci est une Note Longue (note publique visible en bas de la page web).

%NS Ceci est une Note de Service (note visible seulement au CDS).

une ligne pour les data : %O.[data]

%O.velocities |z| 4.3792 0.0064 | |a,D,Opt, | | | |2013A&A...549A..63K|

%O.X | 76.86 3.97| 2012ApJ...752...56F|

Retour au sommaire

un commentaire interne à cosim : %COM commentaire

Cette ligne est réécrite dans parfile.out. (Mais après la ligne "bye!")

un commentaire interne au parfile : # commentaire

Cette ligne n'est pas réécrite dans parfile.out.

les options

.SEARCH rayon

Pour forcer le rayon de recherche à une valeur donnée :

.SEARCH 5"

.SEARCH 2'

ou pour étendre ou réduire le rayon de recherche d'un facteur :

.SEARCH *2

.SEARCH COO ou .SEARCH NOCOO

Pour forcer dans tous les cas la recherche autour de la position ( COO) ou, au contraire, empêcher toute recherche autour de la position ( NOCOO).

.NOLIMIT

Par défaut la recherche par coordonnées s'arrête à 100 objets trouvés, même si le rayon n'est pas encore atteint. Quand cette limite est atteinte, on voit apparaître dans le parfile.out la ligne :

!... Reach LIMIT 100 objects -> NOLIMIT

En ajoutant l'option :

.NOLIMIT

cette limite est enlevée. La totalité des candidats inclus dans le rayon de recherche sera affichée.

.SIGMA COO

Pour modifier le sigma dans le calcul du score des coordonnées.

Rappel : le sigma est relatif à l'erreur donnée; si il n'y en a pas, il est fixé selon la lettre de qualité :
E → 10 000 mas
D → 3 000 mas
C → 1 000 mas
B → 300 mas
A → 100 mas

Cette option va borner le sigma des coordonnées, toutes lettres confondues, entre une valeur minimale et une valeur maximale, séparées par une virgule et exprimées en mas :

.SIGMA COO 0,1000

ou

.SIGMA COO 1000,10000

Cette option agit aussi bien sur les objets en entrée que sur les objets de SIMBAD.

Attention aux objets étendus dont une taille est donnée : le sigma sera toujours rapporté au minimum au dixième du grand axe.

Option plus fine en cours de réglage : .SIGMA COO lettre par lettre (ne fonctionne pas encore au 17/3/2023)

Retour au sommaire

.SIGMA V

Pour modifier le sigma dans le calcul du score des vitesses.

Rappel : le sigma est égal au dixième de la valeur pour les qualités ABCD et non photométriques, ou à 0.3 x la valeur pour les qualités E ou photométriques, qu'il y ait ou non une erreur donnée. Soit :
σ[ABCD] = a x z
σ[E] = 3 x a x z
a étant par défaut égal à 0.1

Cette option permet de redéfinir la valeur de a :

.SIGMA V 0.01

donnera un sigma de 0.01z pour les qualités ABCD, 0.03z pour la qualité E.

.SIGMA M

Pour modifier le sigma dans le calcul du score des magnitudes.

Rappel : le sigma est une valeur fixe de 0.3 pour les qualités ABCD et de 1 pour la qualité E, qu'il y ait ou non une erreur donnée.

Erreur constatée en décembre 2019 : ce n'est pas comme ça que ça se passe! Apparemment le sigma pour la qualité E est toujours 3 fois le sigma des autres lettres, sauf si on attribue une valeur lettre par lettre comme ci-dessous :

On peut attribuer directement une nouvelle valeur à ce sigma, lettre par lettre :

.SIGMA M D 1
.SIGMA M C 0.5

Par contre, l'indication ci-dessous donnera un sigma de 1 pour les lettres ABCD mais de 3 pour la lettre E :

On peut aussi attribuer une valeur, toutes lettres confondues :

.SIGMA M 1 

Dans les deux cas, l'option agit sur toutes les magnitudes (toutes les longueurs d'onde).

Retour au sommaire

.COO .V .M .OT

Pour modifier les valeurs min et max de chaque score.

Pour chaque score, il existe :
- une valeur min au dessous de laquelle le score est considéré comme BAS (correspond à un MAUVAIS candidat)
- une valeur max au dessus de laquelle le score est considéré comme HAUT (correspond à un BON candidat)
Entre min et max, le score est considéré comme MEDIUM : concrètement dans une même table les scores MEDIUM viendront d'un mélange de BONS candidats et de MAUVAIS candidats.

Par défaut les valeurs min et max de COO sont 0 et 1. Celles des autres paramètres (V et M) sont 0 et 0 (donc pas de score MEDIUM). Cela se note de la façon suivante :

.COO 0,1
.M 0,0
.V 0,0

Pour les scores qui ne prennent que des valeurs entières, comme OT, les valeurs par défaut de min et max sont -1 et +1, 0 étant la valeur MEDIUM. Cela peut aussi bien se noter de cette façon :

.OT -1,1

que de cette façon :

.OT 0.5,0.5

Exemple de modification pour COO:

.COO 0.35,0.55

Avec cette notation, tous les scores COO <= 0.3 seront BAS; tous les scores COO >= 0.6 seront HAUTS. Les scores COO de 0.4 et 0.5 resteront dans les DISPLAY si l'objet a été trouvé par coordonnées, mais n'empêcheront par l'objet de "passer" en UPDATE s'il a été trouvé par identificateur.

.ACCEPT COO ou .REFUSE COO

Pour accepter/refuser les candidats avec des coordonnées en qualité E ou ~ (dont les scores sont tous hauts).

Par défaut un candidat trouvé par coordonnées avec que des scores hauts restera quand même en DISPLAY si ses coordonnées sont en qualité E ou en qualité ~ (=absence de lettre de qualité).

Exemples :

.ACCEPT COO E

pour accepter les candidats trouvés avec des coordonnées en qualité E, comme tout autre candidat aux coordonnées en qualité, A, B, C ou D.

.REFUSE COO ~

pour refuser les candidats trouvés avec des coordonnées sans qualité.

.OT+ ou .OT-

Pour modifier le score OT.

Le résultat de la comparaison des types d'objets principaux donne le score OT. Par défaut il prend les valeurs -1 (types incompatibles) ou 1 (types compatibles). Tant qu'un candidat garde un score OT à -1, il sera considéré comme mauvais et ne pourra passer en mise à jour.

Par exemple, voici deux objets dont le type d'objet (Rad et ERO) est considéré incompatible avec le mien (G) :

!---   2/4: Rad/  G(-1)   0.6"C (  1.0) V( 1.0)       ; from coo
!---  3/10: ERO/  G(-1)   0.2"D (  1.3) V( 1.0) M( 0.5); from coo

Je décide d'accepter ces types d'objets :

.OT+2 Rad,ERO

Cela donnera :

!+++   2/4: Rad/  G(1)   0.6"C (  1.0) V( 1.0)       ; from coo
!+++  3/10: ERO/  G(1)   0.2"D (  1.3) V( 1.0) M( 0.5); from coo

Retour au sommaire

.ACRO+ ou .ACRO-

Pour favoriser/défavoriser certains acronymes.

Exemple 1 :

.ACRO+1 2XMM,XMMM31,XMMU
.ACRO 0,1
.TRUST ACRO

Seuls les candidats contenant au moins un des acronymes suivants : 2XMM, XMMM31, ou XMMU seront retenus (sous réserve que les autres scores OT, COO etc. sont hauts). Tout autre candidat aura un score ACRO de zéro et sera rejeté. !!!Attention : si en entrée il y a aussi un identificateur différent du même acronyme, le candidat trouvé sera rejeté quand même.

Exemple 2 :

.ACRO-1 WARPS,[MEJ2006],[LB2005]
.ACRO -1,0
.TRUST ACRO

Chacun des acronymes suivants : WARPS, [MEJ2006] ou [LB2005], trouvé parmi les identificateurs d'un candidat, fera baisser son score ACRO de 1. En résultat, le candidat sera rejeté. Par contre le score ACRO reste inchangé si l'identificateur était donné en entrée dans la table.

.B+ ou .B-

Pour favoriser/défavoriser certaines références.

!!!!! Attenrion option invalide au 16/12/2021

.HIDE

Pour retirer de l'affichage en DISPLAY certains candidats selon leur type d'objet principal.

Peut être utile lorsque le résultat de la recherche par coordonées est "encombré" d'objets qui ne nous intéressent pas du tout (par exemple des galaxies lorsqu'on cherche des amas de galaxies). Cette option n'a aucune incidence sur la décision prise par cosim; elle n'existe que pour simplifier l'affichage.

Exemple :

.HIDE G,rG,GiG,GiC,BiC,LeG

Retour au sommaire

.IGNORE

Pour ne pas tenir compte d'un score.

Par défaut, pour qu'un candidat soit considéré comme GOOD, tous ses scores doivent être HIGH (éventuellement certains MEDIUM seulement s'il a été trouvé par identificateur). Avec l'option .IGNORE, on décide de passer outre les mauvais résultats d'une comparaison, pour faire la cross-idenfication quand même.

Exemple :

.IGNORE M pour ne pas comparer les magnitudes.

.TRUST

Pour donner plus d'importance à un score.

Par défaut un objet dont un score secondaire (autre que OT ou COO) est LOW mais dont tous les autres scores sont HIGH va rester UNDEFINED, en DISPLAY. Avec l'option .TRUST, on déclare que, ce score étant important, ce même objet doit être considéré comme BAD.

Exemple :

J'entre une table de galaxies avec des redshifts. Toute galaxie trouvée à une distance acceptable mais avec un redshift significativement différent restera par défaut en DISPLAY. Si j'écris :

.TRUST V

c'est que j'attache une importance particulière au redshift. Toute galaxie trouvée à une distance acceptable mais avec un redshift significativement différent sera alors rejetée.

.EXPECT

Pour donner encore plus d'importance à un score.

L'option .EXPECT demande qu'en l'absence d'un paramètre dans le candidat SIMBAD (c'est-à-dire pas de redshift pour .EXPECT V, pas de magnitude comparable pour .EXPECT M...), aucune cross-identification ne soit faite.

Exemple :

J'entre une table de galaxies avec des redshifts. Par défaut toute galaxie trouvée à une distance acceptable, même si elle n'a pas de redshift peut être cross-identifiée. Si j'écris :

.EXPECT V

alors j'attache une grande importance au redshift. Non seulement toute galaxie avec un redshift significativement différent sera rejetée, mais aussi toute galaxie sans redshift.

!! Attention :
- l'option .EXPECT n'a pas d'effet sur les objets trouvés par identificateur.
- l'option .EXPECT n'a pas d'effet dans les cas où le manque de paramètre à comparer est du côté de l'objet de la table d'entrée.

Retour au sommaire

.OK INCORRECT ID

Pour permettre l'édition de commandes de mise à jour même si un identificateur est mal écrit.

.OK ALREADY CONNECTED

Pour passer outre la sécurité "Already Connected".

.FORCE FROM ID

Pour forcer le passage de tous les objets trouvés par identificateur.

A utiliser en dernier recours, quand il n'est plus possible d'augmenter la tolérance sur les scores qui risqueraient d'attirer des cross-identifications non voulues, et après avoir vérifié individuellement chaque DISPLAY restant.

.REPLACE EQUAL et .REPLACE NO

Pour forcer le remplacement d'une donnée à qualités égales ou au contraire empêcher le remplacement d'une donnée même si celle de la table est meilleure.

Exemple 1 :

.REPLACE EQUAL M

donnera :

$c m J 16.443 [0.137] D 2011AcA....61..103G
$!!m J (Vega) 16.442 [0.137] D 2003yCat.2246....0C !!Simbad replaced

Exemple 2 :

.REPLACE NO V

donnera :

$!c V z:s 1.103 (Opt) C 2012ApJ...761..140C
$!!V z:spectroscopic 1.102 (~) D [~] 2012ApJ...761..140C !!Simbad kept

Exemple 3 :

.REPLACE EQUAL COO

donnera :

$c c 150.415665+02.430201 (O) [10 10 90] C 2013yCat.5139....0A
$!!c 10 01 39.802+02 25 49.35 (Xray) [480 480 90] C 2007ApJS..172...29H !!Simbad replaced

.DEBUG

Pour afficher plus de détails sur chaque "décision" de Cosim.

Retour au sommaire

le fonctionnement de cosim

Ce chapitre suit plus ou moins une version PDF qui contient en plus des couleurs et beaucoup de schémas. Pour la voir, cliquer ici.

• Résumé des différentes étapes du programme

Cosim (Comparaison d'Objets pour SIMbad) recherche les objets dans SIMBAD par identificateur et/ou par coordonnées. Puis il évalue chaque candidat trouvé en le comparant à l'objet d'entrée. Pour ce faire, il met en place plusieurs scores, le plus généralement un pour le type d'objet, un pour les coordonnées, un pour la vitesse si il y en a, un pour la magnitude si il y en a. La combinaison de ces scores aboutit à l'évaluation du candidat, pour savoir s'il est bon pour être cross-identifié avec l'objet d'entrée, ou s'il est à rejeter. Dès qu'il y a zéro ou un et un seul bon candidat, cosim édite les commandes pour SIMBAD, de création ou de mise à jour.

• Le vocabulaire utilisé

Chaque score est dit LOW, HIGH ou MEDIUM.
Après combinaison de tous les scores, chaque candidat trouvé est dit BAD, GOOD ou UNDEFINED.
En fonction du nombre de candidats trouvés pour chaque entrée, et de leur statut, le programme va soit créer un nouvel objet (NEW), soit mettre à jour un objet existant (UPDATE), soit rester dans l'indécision, avec affichage des candidats trouvés et de leurs scores (DISPLAY).

• 1ère étape : recherche d'objets dans SIMBAD

Cosim recherche par identificateur. Si aucun objet n'est trouvé par identificateur, il recherche par coordonnées. Le rayon de recherche est fonction de la précision des coordonnées en entrée (comme dans l'ancien raccord).

1° → 3600''
0,1° → 1440''
1' → 300''
0,1' → 36''
1'' → 6''
0,1'' → 2''
0,01'' → 1,2''
1mas → 0,6''
puis borné entre 30'' et 3600''.

Options pouvant agir durant cette étape : .SEARCH rayon, .SEARCH COO ou NOCOO

• 2ème étape : évaluation de chaque candidat trouvé

Chacun des paramètres suivants, s'il est présent à la fois dans l'objet de la table et dans l'objet de SIMBAD, va être comparé et générer un score spécifique : le type d'objet principal, les coordonnées, la vitesse ou le redshift et les magnitudes.

champ %J → score COO
champ %C → score OT
champ %V → score V
champ %M → score M

Dans l'impossibilité de comparer un paramètre, le score correspondant est ignoré.

Il existe deux autres scores, ignorés par défaut mais considérés dès qu'on appelle l'option correspondante : ACRO (pour tester la présence/absence d'un acronyme) et B (pour tester la présence/absence d'un bibcode dans la bibliographie d'un objet).

Options pouvant agir sur le choix des scores : .IGNORE, .ACRO+ ou .ACRO-, .B+ ou .B-

Retour au sommaire

• • Calcul des scores

• • • le score du type d'objet, OT

Le score du type d'objet se base sur un tableau de comparaison des types entre eux utilisé dans SIMBAD.
La comparaison n'est faite qu'entre types d'objets principaux.
Si les deux types sont compatibles, il prend la valeur de 1;
Si les deux types ne sont pas compatibles, il prend la valeur de -1;
[Si les deux types sont neutres, il prend la valeur de 0.] (pour l'instant la comparaison n'aboutit jamais à ce résultat)

Options pouvant agir sur le calcul du score OT : .OT+ ou .OT-

• • • le score des coordonnées, COO

Selon la formule :

1.5-0.5x((valeurSimbad-valeurTable)/√(σSimbad²+σTable²) )

Les limites du score sont définies par la formule: de moins l'infini à +1.5
Le sigma est fonction de l'erreur donnée; sinon il est estimé d'après la lettre de qualité:
E → 10 000 mas
D → 3 000 mas
C → 1 000 mas
B → 300 mas
A → 100 mas
Si une taille est donnée, le sigma est égal au minimum au dixième du grand axe. La taille ira toujours dans le sens d'augmenter un sigma, elle ne peut pas le réduire.

Option pouvant agir sur le calcul du score COO : .SIGMA COO

• • • le score de la vitesse, V

Selon la formule :

1.5-0.5x((valeurSimbad-valeurTable)/√(σSimbad²+σTable²) )

Le sigma est égal au dixième de la valeur pour les qualités ABCD et non photométriques, ou à 0.3 x la valeur pour les qualités E ou photométriques, qu'il y ait ou non une erreur donnée.
Le résultat de la formule est borné entre -1 et +1.

Option pouvant agir sur le calcul du score V : .SIGMA V

Retour au sommaire

• • • le score de la (des) magnitude(s), M

Selon la formule :

1.5-0.5x((valeurSimbad-valeurTable)/√(σSimbad²+σTable²) )

Le sigma est une valeur fixe de 0,3 pour les qualités ABCD et de 1 pour les qualités E, qu'il y ait ou non une erreur donnée entre crochets.
Le résultat de la formule est borné entre -1 et +1.
Quand plusieurs magnitudes sont comparées, un score moyen est calculé à partir des différents scores obtenus.

Option pouvant agir sur le calcul du score M : .SIGMA M

• • • le score de l'acronyme, ACRO

Par défaut, le score ACRO n'existe pas. Il existe dès qu'on utilise l'option .ACRO+, .ACRO- ou .ACRO++.
Il devient alors un score comme les autres qui va entrer en jeu dans la combinaison des scores.
Il ne prend que des valeurs entières.

Options pour utiliser le score ACRO : .ACRO+ ou .ACRO-

• • • le score "bibcode", B

Par défaut, le score B n'existe pas. Il existe dès qu'on utilise l'option .B+ ou .B-.
Il devient alors un score comme les autres qui va entrer en jeu dans la combinaison des scores.
Il ne prend que des valeurs entières.

Options pour utiliser le score B : .B+ ou .B-

• • Classement des scores

A chaque score sont attribués un minimum et un maximum.
Au dessous de min : le score est dit LOW (score a priori pour un candidat BAD).
Au dessus de max : le score est dit HIGH (score a priori pour un candidat GOOD).
Entre les deux, le score est dit MEDIUM.
Valeurs par défaut (notation .SCORE min,max) :
.OT -1,1
.COO 0,1
.V 0,0
.M 0,0
(.ACRO -1,1)
(.B 0,1)
Quand l'objet est trouvé par identificateur, tous les scores MEDIUM deviennent HIGH.

Options pouvant agir sur le classement des scores : .COO, .V, .M ...

Retour au sommaire

• • Combinaison de ces scores

Pour chaque objet trouvé, le programme va décider s'il est bon ou mauvais pour une cross-identification: c'est l'évaluation du candidat.

Tous les scores calculés et non ignorés entrent en compte pour l'évaluation du candidat.
Les scores OT et COO étant considérés comme essentiels, dès que l'un ou l'autre est LOW, le candidat est étiqueté BAD. Ce n'est pas le cas pour les autres scores (V, M...) : si un de ces scores est LOW alors que tous les autres scores sont HIGH, le candidat reste UNDEFINED.
Pour qu'un candidat soit étiqueté GOOD, tous ses scores calculés doivent être HIGH si il a été trouvé par coordonnées; tous ses scores calculés doivent être au moins MEDIUM si il a été trouvé par identificateur.

Options pouvant agir sur l'évaluation du candidat : .IGNORE, .TRUST, .EXPECT, .FORCE FROM ID

Retour au sommaire

• 3ème étape : décision finale

3 grands types de décision sont possibles :

1) les NEW : un nouvel objet sera créé

C'est le cas quand aucun objet n'est trouvé, ou si tous les objets trouvés sont BAD (à condition qu'ils n'aient pas été trouvés par identificateur).

2) les UPDATE : un objet existant sera mis à jour

C'est le cas quand il y a un et un seul candidat GOOD (même si il y en a d'autres UNDEFINED, et à condition qu'il n'y ait pas d'autre candidat, même BAD trouvé par identificateur).

On distingue deux sous-catégories de 'UPDATE' :

Les UPDATE from id quand l'objet a été trouvé par idenficateur,

et les UPDATE from coo quand l'objet a été trouvé par coordonnées.

3) les DISPLAY : le programme ne peut pas décider

Ce sont tous les autres cas. Il n'y a alors pas d'édition de commandes de mise à jour et tout reste affiché (l'objet en entrée, les objets SIMBAD trouvés avec leurs scores).

Option pouvant agir sur l'affichage du DISPLAY : .HIDE

Parmi les DISPLAY, quelques cas se distinguent et sont signalés par l'un des messages suivants :

Possible Merge

...qui est scindé en fait en deux sous-catégories :

Possible Merge (mult. good)

Ce message apparaît dès qu'il y a deux (ou plus) candidats GOOD, quelle que soit la façon dont ils ont été trouvés.

Possible Merge (id+coo)

Ce message apparaît dès qu'il y a un candidat GOOD trouvé par coordonnées et un autre candidat trouvé par identificateur, quel que soit son statut.
Ce message n'est possible qu'avec l'option .SEARCH COO.

Mais la plupart du temps on n'a pas besoin de distinguer entre ces deux sous-catégories. Dans les deux cas c'est qu'il y a du ménage à faire, soit dans SIMBAD, soit dans l'extraction des données de la table.

Conflicting

Ce message apparaît dès que la recherche par identificateur a donné plusieurs objets dans SIMBAD, quelle que soit la valeur de ces candidats.

Already Connected

Ce message apparaît lorsque le programme devrait logiquement faire un UPDATE sur un objet SIMBAD déjà mis à jour plus haut dans le même fichier.

Options pouvant agir durant cette étape : .OK INCORRECT ID, .OK Already Connected

Retour au sommaire

• Commentaires sur les commandes de mise à jour

à propos de la référence bibliographique

La référence bibliographique mise en .B dans le parfile sera ajoutée sur chaque objet mis à jour ou créé. Derrière cette référence sont ajoutées automatiquement trois informations supplémentaires :
La lettre d, qui indique que l'objet apparaît dans une table
Le ou les raw-id, c'est-à-dire la liste des identificateurs tels qu'ils sont écrits dans la table. C'est ce qui est donné dans le parfile par %I.R. Si %I.R n'est pas précisé, le programme va prendre le contenu de %I.0 comme raw-id.
L'occurence +1

à propos du type d'objet

Si le type d'objet principal ne donne pas lieu à la commande 'c o ', il est systématiquement ajouté en type secondaire (commande 'a o ').

données de SIMBAD conservées ou remplacées?

Chaque donnée du parfile donne lieu à une mise à jour éventuelle dans SIMBAD. Dans le cas des UPDATE, le programme va choisir entre garder l'ancienne valeur de SIMBAD ou écraser cette valeur. Pour déterminer si la donnée de la table est plus précise que la donnée de SIMBAD, le programme va suivre une suite de questions. Voici par exemple celui pour le champ vitesse :

J'ai une vitesse en entrée, Vtable.
_Question 1 : Vsimbad existe?_
Réponse Q1 oui → question suivante
Réponse Q1 non → commande "a V ".
Question 2 : Vsimbad a-t-il un bibcode?
Réponse Q2 oui → Question 2.1 : Vtable a-t-il un bibcode?
Réponse Q2.1 oui → question suivante
Réponse Q2.1 non → on ne fait rien.
Réponse Q2 non → Question 2.2 : Vtable a-t-il un bibcode?
Réponse Q2.2 oui → commande "c V ".
Réponse Q2.2 non → question suivante
Question 3 : Comment est la qualité (ABCDE) de Vtable?
Réponse Q3 "moins bonne que Vsimbad" → on ne fait rien.
Réponse Q3 "égale à Vsimbad" → question suivante
Réponse Q3 "meilleure que Vsimbad" → commande "c V ".
Question 4 : Comment est la nature (photo, spectro, vide...) de Vtable?
Réponse Q4 "moins bonne que Vsimbad (p pour s ou s pour sa ou ~ pour p...)" → on ne fait rien.
Réponse Q4 "égale à Vsimbad (~ pour ~ ou s pour s ou p pour p ou sa pour se...)" → question suivante
Réponse Q4 "meilleure que Vsimbad (s pour p ou p pour ~ ou sa pour s...)" → commande "c V ".
Question 5 : Vtable a-t-il une erreur donnée?
Réponse Q5 oui → Question 5.1 : Vsimbad a-t-il une erreur donnée?
Réponse Q5.1 oui → Question 5.2 : erreur de Vtable inférieure à erreur de Vsimbad?
Réponse Q5.2 oui → commande "c V ".
Réponse Q5.2 non → on ne fait rien.
Réponse Q5.1 non → commande "c V ".
Réponse Q5 non → on ne fait rien.

Dans tous les cas, une ligne contenant la commande de mise à jour ( c...) est générée, neutralisée par un '!' si la décision est de ne pas écraser.
Une deuxième ligne au-dessous, en commentaire, rappelle la valeur de SIMBAD. Dans le cas où l'ancienne valeur SIMBAD va être conservée, la ligne contient !!Simbad kept. Dans le cas où l'ancienne valeur SIMBAD va être écrasée, la ligne contient !!Simbad replaced. Par exemple :

!c c 04 43 36.08 -66 43 44.0 (Opt) C [230 230 90] 2011AcA....61..103G
!!c 04 43 36.08-66 43 44.0 (Opt ) [230 230 90] C 2010AcA....60....1P !!Simbad kept

ou

c c 04 40 41.27 -69 15 03.1 (Opt) C [230 230 90] 2011AcA....61..103G
!!c 04 40 41.21-69 15 03.0 (Opt ) [~ ~ ~] D 2011A&A...536A..60S !!Simbad replaced

Pour le type d'objet principal, ces informations sont sur la même ligne que la commande de mise à jour. De plus, en cas de !!Simbad kept, la commande c o n'est pas annulée, mais remplacée par la commande a o. Par exemple :

a o Y*O 2013A&A...549A.130V! # o denseCore ~ !!Simbad kept
ou
c o Y*O 2013A&A...549A.130V! # o Star ~ !!Simbad replaced

Ceci permet d'une part de vérifier rapidement ce que le programme va écraser ou ne pas écraser, d'autre part de récupérer si besoin des commandes de mise à jour qui n'avaient pas été prévues par le programme.

Options pouvant agir durant cette étape : .REPLACE EQUAL et .REPLACE NO

Retour au sommaire

-- CatherineBrunet - 2015-01-23