Simulating cross-correlation of planetary emission.
To install from source, run
python3 -m pip install -U pip
python3 -m pip install -U setuptools setuptools_scm pep517
git clone https://github.com/arjunsavel/scope
cd scope
python3 -m pip install -e .
You'll also need to download some data files. Currently, these data files are about 141 MB large. You can download them (to the correct directory, even!) with the following:
cd src/scope
chmod +x download_data.bash
./download_data.bash
This will create a data
directory and plop the relevant files into it.
The bulk of scope
's high-level functionality is contained in scope/run_simulation.py
.
For a detailed tutorial, see the documentation.
To run a large set of models, edit scope/grid.py
to define a parameter grid. This grid is then used to run a set of simulations in scope/run_simulation.py
;
the command python run_simulation.py n
will run a simulation in the defined grid at index n
.
The scope.run_simulation.make_data
function can be used to simulate a single high-resolution dataset. To
simulate detection significances, use scope.run_simulation.calc_log_likelihood
.
Running the script requires an exoplanet spectrum, stellar spectrum, and telluric spectrum. Default parameters are currently correspond to the exoplanet WASP-77Ab.
Once completed, scope.run_simulation.calc_log_likelihood
will output:
simdata
file: the simulated flux cube with PCA performed. That is, the principle components with the largest variance have been removed.nopca_simdata
file: the simulated flux cube, including all spectral components (exoplanet, star, blaze function, tellurics).A_noplanet
file: the simulated flux cube with the lowest-variance principle component removed.lls_
file: the log-likelihood surface for the simulated flux cube, as a Kp--Vsys map.ccfs_
file: the cross-correlation function for the simulated flux cube, as a Kp--Vsys map.