128311 Blind search¶
We demonstrate how to conduct a blind search in pre- and post-uprade-separate HIRES data. NOTE: If notebook is not working, can just open iPython, import rvsearch.search and rvsearch.utils, and run the search with the 3 commands below.
[3]:
import os
import rvsearch
from rvsearch import search, utils
Load the data and initialize the search object.
Bin data within 12 hours, and search down to 10 days.
[4]:
data = utils.read_from_csv(os.path.join(rvsearch.DATADIR, 'HD128311.csv'), binsize=0.5)
searcher = search.Search(data, starname='128311', min_per=10,
workers=1, mcmc=True, verbose=True)
Run search on a single core. In general, you’ll want to allocate more CPUs for the search, but multi-threading within the Jupyter notebook environment is buggy. Run MCMC after search is done (default). With our current configuration, outputs will be saved in in the current working directory. Will save periodogram for each successive search, orbit plots, and corner plot.
[ ]:
searcher.run_search()