<jupyter_start><jupyter_text>128311 blind searchWe 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.<jupyter_code>import os

import radvel
import rvsearch
from rvsearch import search, inject, utils, plots, driver<jupyter_output><empty_output><jupyter_text>- Load the data and initialize the search object.- Bin data within 12 hours, and search down to 10 days.<jupyter_code>data = utils.read_from_csv('../../Planets/rvdata/vst128311.csv', binsize=0.5)
searcher = search.Search(data, starname='128311', min_per=60, 
                         workers=8, mcmc=False, verbose=True, mstar=[1., 0.05])<jupyter_output><empty_output><jupyter_text>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.<jupyter_code>searcher.run_search()<jupyter_output>0%|          | 0/769 [00:00<?, ?it/s]<jupyter_text>128311 search with known planets<jupyter_code>data = utils.read_from_csv('../../Planets/rvdata/vst128311.csv', binsize=0.5)
post = radvel.posterior.load('post_final.pkl')
searcher = search.Search(data, post=post, starname='128311', min_per=60, 
                         workers=8, mcmc=False, verbose=True, mstar=[1., 0.05])
searcher.run_search()
args = None
class _args(object):
    def __init__(self):
        self.num_cpus = 8
        self.num_inject = 500
        self.overwrite = True
        self.mstar = 1.
        self.rstar = 1.
        self.teff = 1.
        self.minP = 2.
        self.maxP = 1e5
        self.minK = 0.1
        self.maxK = 1000.0
        self.minE = 0.0
        self.maxE = 1.0
        self.betaE = False

# Run injections.
args = _args()
args.search_dir = '128311' # Directory in which search is saved
args.full_grid = False
args.verbose = False
driver.injections(args)
# Plot injections and completeness.
args.type = ['recovery']
args.fmt = 'png'
args.mstar = 1.
driver.plots(args)<jupyter_output>Creating recovery plot for 128311
Plotting inj_msini vs. inj_au
Recovery plot saved to /Users/lee/Academics/Astronomy/Planets/rvsearch/tutorials/128311/128311_recoveries.png
{'fraction_non_alphanumeric': 0.1045673076923077, 'fraction_numerical': 0.04567307692307692, 'mean_word_length': 3.702448210922787, 'pattern_counts': {'":': 0, '<': 16, '<?xml version=': 0, '>': 15, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '4', 'dir': 'jupyter-structured-clean-dedup', 'id': '148571', 'license': 'MIT', 'n_tokens_mistral': 901, 'n_tokens_neox': 819, 'n_words': 235, 'path': 'tutorials/128311_Search.ipynb', 'repo_name': 'spencerhurt/rvsearch', 'size': '2496'}