In [1]:

    

from astropy.coordinates import SkyCoord, ICRS, CartesianRepresentation, CartesianDifferential, Galactic, Galactocentric
import astropy.units as u
import numpy as np
import healpy as hp


    

In [6]:

    

#low density
c1 = SkyCoord(l = 20.0*u.degree, b = 30*u.degree, frame='galactic')
ra = c1.icrs.ra.value
dec = c1.icrs.dec.value
print("ra in deg:",ra)
print("dec in deg:",dec)
c = SkyCoord(ra=ra*u.degree, dec=dec*u.degree, distance=48.1*u.kpc, frame='icrs',
             pm_ra_cosdec = 1.95*u.mas/u.yr, pm_dec = 0.43*u.mas/u.yr, radial_velocity = 283*u.km/u.s,
             galcen_distance = 8.0*u.kpc, z_sun = 15.0*u.pc,
             galcen_v_sun=CartesianDifferential(d_x=11.1*u.km/u.s, d_y=239.08*u.km/u.s, d_z=7.25*u.km/u.s))
c.galactocentric

hpxlevel = 6
nside = 64
hpx_number = hp.ang2pix(nside = nside, theta = ra, phi = dec, nest=True, lonlat=True)
print("lower and larger source_id value")
print(2**35*np.power(4,12-hpxlevel)*hpx_number)
print(2**35*np.power(4,12-hpxlevel)*(hpx_number+1)-1)


    

    




ra in deg: 250.57034810778302
dec in deg: 3.1817096947567327
lower and larger source_id value
4433793833146253312
4433934570634608639

In [7]:

    

#LMC
c1 = SkyCoord(l = 280.5*u.degree, b = -32.9*u.degree, frame='galactic')
ra = c1.icrs.ra.value
dec = c1.icrs.dec.value
print("ra in deg:",ra)
print("dec in deg:",dec)
c = SkyCoord(ra=ra*u.degree, dec=dec*u.degree, distance=48.1*u.kpc, frame='icrs',
             pm_ra_cosdec = 1.95*u.mas/u.yr, pm_dec = 0.43*u.mas/u.yr, radial_velocity = 283*u.km/u.s,
             galcen_distance = 8.0*u.kpc, z_sun = 15.0*u.pc,
             galcen_v_sun=CartesianDifferential(d_x=11.1*u.km/u.s, d_y=239.08*u.km/u.s, d_z=7.25*u.km/u.s))
c.galactocentric

hpxlevel = 6
nside = 64
hpx_number = hp.ang2pix(nside = nside, theta = ra, phi = dec, nest=True, lonlat=True)
print("lower and larger source_id value")
print(2**35*np.power(4,12-hpxlevel)*hpx_number)
print(2**35*np.power(4,12-hpxlevel)*(hpx_number+1)-1)


    

    




ra in deg: 80.8456130588062
dec in deg: -69.78267074987376
lower and larger source_id value
4657847914607935488
4657988652096290815

In [8]:

    

#Baades Window
c1 = SkyCoord(l = 1.03*u.degree, b = -03.9*u.degree, frame='galactic')
ra = c1.icrs.ra.value
dec = c1.icrs.dec.value
print("ra in deg:",ra)
print("dec in deg:",dec)
c = SkyCoord(ra=ra*u.degree, dec=dec*u.degree, distance=48.1*u.kpc, frame='icrs',
             pm_ra_cosdec = 1.95*u.mas/u.yr, pm_dec = 0.43*u.mas/u.yr, radial_velocity = 283*u.km/u.s,
             galcen_distance = 8.0*u.kpc, z_sun = 15.0*u.pc,
             galcen_v_sun=CartesianDifferential(d_x=11.1*u.km/u.s, d_y=239.08*u.km/u.s, d_z=7.25*u.km/u.s))
c.galactocentric

hpxlevel = 6
nside = 64
hpx_number = hp.ang2pix(nside = nside, theta = ra, phi = dec, nest=True, lonlat=True)
print("lower and larger source_id value")
print(2**35*np.power(4,12-hpxlevel)*hpx_number)
print(2**35*np.power(4,12-hpxlevel)*(hpx_number+1)-1)


    

    




ra in deg: 270.8692734479971
dec in deg: -30.01680074334423
lower and larger source_id value
4050143439889629184
4050284177377984511

In [ ]: