""" 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$")