cmateu / galstreams

Milky Way Streams Footprint Library and Toolkit for Python
BSD 3-Clause "New" or "Revised" License
46 stars 17 forks source link

Pal 5 total velocity > 1000km/s in Galactocentric coordintes. #26

Closed appy2806 closed 1 year ago

appy2806 commented 1 year ago

When I use astropy to convert Pal 5 velocites in galactocentric coordaintes using the following code:

import astropy.coordinates as ac
import galstreams
mwsts = galstreams.MWStreams(verbose=False, implement_Off=False)
gsr = ac.Galactocentric()
st_name = 'Pal5-S20'
coords_st = mwsts[st_name].track.transform_to(gsr)

The total velocities for stream particles > 1000 km/s which seems rather high? Even for the highly eccentric orbit it's on?

cmateu commented 1 year ago

Actually these values are meaningless since the Pal5-S20 track does not have any proper motion information, it only has sky coordinates+distance. You can check whether a given track has proper motion info by checking InfoFlag[2], 0=False, 1=True, 2=True_but_caveat(see paper). You can also check the same info in the has_pm column of the summary table (mwsts.summary['has_pm'], same convention). For Pal5 kinematic data better use the PW19 (recommended, see paper) or I21 tracks.

appy2806 commented 1 year ago

Thanks for the prompt response !