Completeness map

[1]:

import matplotlib.pyplot as plt
import numpy as np

[2]:

from exogaia.data import EpochAstrometry
from exogaia.limits import CompletenessMap

The archive is unstable and may perform below expectations. If launching multiple, consecutive, heavy queries through Python, please space them out (e.g., using sleep(1)) to avoid overloading the system. Please contact the Gaia helpdesk in case of questions (https://www.cosmos.esa.int/web/gaia/gaia-helpdesk). Workaround solutions for the issues following the December 2025 infrastructure upgrade: https://www.cosmos.esa.int/web/gaia/news#WorkaroundArchive

[3]:

np.random.seed(42)

[4]:

source_id = 6199395656645840384  # HD 135344 A

[5]:

primary_mass = (2.2, 0.1)  # (Msun)

[6]:

cmap = CompletenessMap(source_id=source_id, primary_mass=primary_mass, gaia_release='DR5')

Gaia ID = 6199395656645840384
Gaia release = DR5

Primary mass = 2.20 +/- 0.10
INFO: Query finished. [astroquery.utils.tap.core]

[7]:

fig = cmap.calc_completeness(
    n_sigma=3,
    n_samples=20,
    mass_points=np.linspace(0.0001, 0.02, 30),
    filter_sigma=1.0,
    plot_file='cmap.png',
)


----------------
Completeness map
----------------
../_images/tutorials_completeness_7_2.png 
[8]:

fig.axes[0].plot(20., 0.01, 'X', ms=8., mec='white', mfc='tomato')
fig.axes[0].text(20., 0.008, 'HD 135344 Ab', fontsize=10., color='tomato', va='top', ha='center')
fig.savefig('cmap.png')
fig

[8]:
../_images/tutorials_completeness_8_0.png