ReuV11mar2014 < Documentaliste/VizierDocumentaliste

Réunion VizieR du 11/03/2014

Présents : Marianne B., Sylvain G., Gilles L., François O., Pierre O., Emmanuelle P., Patricia V.

Question Gilles :

Cat. J/ApJ/752/42 : graphe SQL fonctionnait, pourquoi être revenu sur un .graph classique ?

Problème avec la période qui ne modifie pas la phase.
=> Gilles regarde

Sinon, graphe peut-être un peu grand en taille de pixels ?
Possibilité de changer les "+" en gros ronds "o" ou en petits ronds "." dans le .graph_sql avant la couleur. :
Exemple du .graph_sql (fonctionne !) :

Set asu="asu -source=J/ApJ/752/42/table2 -out=B-HJD,Bmag,V-HJD,Vmag"
B light curve:.b:SELECT $1 as 'HJD-2454117'  , $2 as 'Bmag' FROM "${asu}"
V light curve:.g:SELECT $3 as 'HJD-2454116'  , $4 as 'Vmag' FROM "${asu}"

Commande du .status associée (fonctionne ; normalement !) :
\vizExplain{ table2 }{}{+ \vizContent{timeSerie} (\vPlotTS{@{@cat}}{.graph_sql}{Light curves}{-P=0.6593}{Light crves}) }

Demande (si possible) d'Emmanuelle :

Afficher les objets "not found" clairement dans la sortie asu lorsqu'on interroge VizieR avec une liste d'identificateurs (Cat. J/ApJ/754/121)

Dans le répertoire du cat. J/ApJ/754/121, le fichier "asu.tsv" est le résultat d'une interrogation du cat. J/ApJS/199/18 par liste de Trigger Name (format YYMMDDNNN) [pour avoir la correspondance avec les noms GRB (YYYYaa)]. (Dans VizieR, on a interrogé J/ApJS/199/18/GBM, en cochant "list of constraints" pour envoyer la liste de noms (trigger names) venant de J/ApJ/754/121, options "verbose" et "add your input" cochées, avec en sortie : "| - separated values" demandé).

En plus, dans cette sortie, j'ai commenté avec un "#" les lignes qui commencent par underscore ou par "---------" et j'ai décommenté, les lignes qui commencent par un chiffre pour avoir tous mes identificateurs.
J'obtiens le fichier "corr_grb" (via un grep) dans lequel je dois enlever les lignes qui ne comprennent qu'un trigger name sans autre info seulement dans le cas où cette ligne est suivie d'une ligne avec le même trigger name et des infos. Sinon je veux qu'il reste afin d'obtenir le même nombre d'objet que dans ma liste initiale aux mêmes endroits...
La ligne de commande : "perl -ane'$x{$F[0]}=$_;END{print foreach sort values%x}' in corr_grb > grb" me donne ce que je veux.

Ce serait plus simple, si on pouvait distinguer dès le départ, dans la sortie asu, les objets "not found" pour qu'ils puissent apparaître simplement dans la liste de résultat même si le résultat est vide.

=> Gilles voit si on peut ajouter une ligne après ces objets là...

Réponse Gilles (va arriver bientôt sur VizieR) : "no DATA" ajouté dans le paragraphe des "No object found". Par exemple :

#RESOURCE=
#Name:
#Title:
#INFO recno=23
#
#Table I_239_hip_main:I_239_hip_main
#Name: I/239/hip_main
#Title:
HIP|RAhms|DEdms|Vmag|RA(ICRS)|DE(ICRS)|Plx|pmRA|pmDE|e_Plx|B-V|Notes
| | |mag|deg|deg|mas|mas/yr|mas/yr|mas|mag|
------|-----------|-----------|-----|------------|------------|-------|--------|--------|------|------|-
23|00 00 17.86|+13 18 44.0| 7.57| 0.07442930| 13.31221083| 12.21| 54.15| 9.65| 0.95| 0.456|

#RESOURCE=
#Name:
#Title:
#INFO recno=0
#
#Table I_239_hip_main:I_239_hip_main
#Name: I/239/hip_main
#Title:
#INFO result=0: no DATA part

Ce qui permet d'utiliser la commande : egrep 'no DATA|^ *[0-9].*' asu.tsv pour retrouver toutes les lignes interrogées.

Sinon, dans le même genre d'idée ; en ligne de commande, on peut faire comme pour le Cat. J/ApJ/748/49, où on interroge un cat. par liste de coordonnées :
asu -f '%v%u|' -source=J/ApJ/715/429/table1 -out.add=_q,_1 -c.rm=60 -c='<<'coo
#-out.add : ajoute un num seq a la requete (_q) et l'input 1ere col. (_1) de la requete
#-c'<<'coo : prend la contrainte du fichier "coo" (qui contient liste de positions en JHHMM+DDMM par exemple)
#-f '%*v%u|' : | pour separer les valeurs et %u pour ajouter les unites.
#Résultats : si "No object" trouvé => input seulement (mais sans numéro de requête... Pas forcément pratique...) et reste de la ligne vide ; si plusieurs résultats trouvés => même num. de requête et input repétés.
Il suffit ensuite d'enlever les lignes de "---", de faire attention aux numéros de requêtes doubles et au décalage pour les "not found" (qui n'ont pas de numéros de requêtes)...

Réponse François :

Avec asu -f '%v%u|' -source=J/ApJ/715/429/table1 -out.add=_1 -c.rm=60 -c='<<'coo, on obtient une table utilisable comme je veux avec une ligne par objet...

------------------------ - ------------ ------------------------------ --------- --------- ------ ------ ----- ---- --- -------- ---- ----- -----
J0042+2320 | |J0041.9+2318|PKS 0039+230 |010.51896deg|+23.33367deg| 0.043deg| 0.201deg| 0.98|FSRQ|
------------------------ - ------------ ------------------------------ --------- --------- ------ ------ ----- ---- --- -------- ---- ----- -----
J0044-8422 | |J0048.0-8412|PKS 0044-84 |011.11192deg|-84.37781deg| 0.191deg| 0.229deg| 0.96|FSRQ|
------------------------ - ------------ ------------------------------ --------- --------- ------ ------ ----- ---- --- -------- ---- ----- -----
J0048+2235
------------------------ - ------------ ------------------------------ --------- --------- ------ ------ ----- ---- --- -------- ---- ----- -----

J0058+3311 | |J0058.0+3314|CRATES J0058+3311 |014.63363deg|+33.18811deg| 0.104deg| 0.110deg| 0.95|BLL |

Questions Patricia :

Ref. de Mar/14 : le setUCD donne le message d'erreur "#+++Column ${recno} could not be found!"; est-ce normal ?

=> Oui ! Dans le .status, il y a :

\vizLink{ table1 }{ EW -Name }{EW}{ \vMore{-source=${$cat}/table2\&-out.max=301\&-out=\htarg{#1,#2,#3,#4,#${*recno}+4}}{${}}}{Display the equivalent width for that star }

qui permet de faire les liens vers les EW correspondants aux bons objets (col. recno+4) mais le recno n'existe pas tant que la ref. n'a pas été chargée dans VizieR...

RAVE-DR4 (Cat. III/272) : quels paramètres afficher par défaut ?

De manière générale, il faut préférer les paramètres mesurés pour ce catalogue (pas les info venant de xid et donc d'autres papiers), en l'occurrence : vitesses radiales et abondances (métallicité ou Fe/H et éléments alpha tels que le Mg en particulier) et flag donnant la qualité du spectre.

En plus :

utiliser la dernière commande \originalcolumnnames (voir réunion du 18 février 2014) pour afficher les labels originaux des colonnes.
Enlever les XXXXXXXXXXXX lorsqu'il n'y a pas l'info en utilisant, soit :
- \vizSQL{ ravedr4 }{set TYCHO2=null where TYCHO2=X% ; \
  set UCAC2=null where UCAC2=X% } etc. pour les autres colonnes. Le "%" remplace une suite de caractères.
- \vizConvert{ ravedr4 }{ trcol -t180-204 'X' ' ' } par exemple (en prenant les colonnes concernées).

Question Sylvain :

Cat. J/AJ/145/15 : une magnitude apparente négative - est-ce normal ?

=> Oui ! Il s'agit d'une magnitude Spitzer/MIPS à 70um. En infrarouge, la comparaison avec la magnitude de l'étoile de référence (Vega - qui n'est pas forcément très brillante en IR) peut être négative.

Voir, par exemple, le tableau des magnitudes apparentes sur la page de Wikipédia où Rigel peut aller jusqu'à -38 (en mag V apparente)...

Précision François :

Attention pour Rigel, le -38 c'est la magnitude qu'aurait Rigel si elle était à ma place du Soleil (à 1 unité astronomique). Ce n'est donc pas la magnitude apparente -- ni sa magnitude absolue (la magnitude absolue est celle qu'aurait l'étoile si elle était à 10pc (environ 2millions d'unités astronomiques)).

Proposition de Gilles :

Suite à plusieurs questions des utilisateurs sur le compte question pour savoir pourquoi il n'y a pas de RAJ2000 dans le ReadMe (par exemple) : pourquoi ne pas récupérer les données de VizieR pour faire un "ReadMe VizieR" ?

=> Bonne idée ! Cela permettrait de faire la différence entre le ReadMe qui concerne les tables FTP et la description des colonnes telles qu'elles sont dans VizieR qui peuvent être différentes (labels différents, colonnes additionnelles, parfois noms de tables différents, explications modifiées, etc.)

Peut-être en reprenant les données via le METAcola et en choisissant les colonnes lisibles et en créant un 2è onglet ReadMe _VizieR en plus du ReadMe _FTP et/ou en faisant un deuxième lien vers ce ReadMe recréé... Une sortie "parsable" de ce ReadMe sera nécessaire pour les utilisateurs...