""" Perform a cone search against the main Gaia table and write it to disk """
from gaia.tap import cone_search
# Get all sources within a 0.25 degree radius
sources = cone_search(32.341, -1.4245, 0.25)
print(sources)
# 412
print(sources)
sources.write("cone_search.csv")
# plt.axvline(rc+err,linestyle='--')
# plt.axvline(rc-err,linestyle='--')
plt.title('Clump luminosity: ' + str(rc))
plt.show()
gaia_on_tap = False
if gaia_on_tap:
'''
Kepler FOV is 105 square degrees, so sqrt(115)/2 in each direction from the centre.
RA: 19h 22m 40s / 19.378
Dec: +44 30' 00'' / +44.50
'''
'''
Need the absolute-apparent magnitude conversion
Need to see whether this gaia data is representative
'''
ra = (19. + 22. / 60. + 40. / 3600.) * 15.
dec = 44. + 30. / 60.
r = np.sqrt(105) / 2
sources = gt.cone_search(ra, dec, r, table="gaiadr1.tgas_source")
print(sources)
fig, ax = plt.subplots()
ax.scatter(sources['ra'], sources['dec'], s=1, c="#000000")
ax.set_xlabel(r"$\alpha$")
ax.set_ylabel(r"$\delta$")
plt.show()
""" Perform a cone search against the main Gaia table around M67 """
import astropy.coordinates as coord
import matplotlib.pyplot as plt
from gaia.tap import cone_search
cluster = coord.SkyCoord.from_name("M67")
# Get everything within 1 degree radius of the cluster.
cluster_candidates = cone_search(cluster.ra.deg, cluster.dec.deg, 1.0)
# Plot it.
fig, ax = plt.subplots()
ax.scatter(cluster_candidates["ra"], cluster_candidates["dec"],
s=1, c="#000000")
ax.set_xlabel(r"$\alpha$")
ax.set_ylabel(r"$\delta$")