Examples

kima comes with a few examples ready to run, showing some of the common use cases.

51 Peg

Here we analyse RV data from the first exoplanet host.
The getting started guide goes into more detail on this example.

BL2009

In the /BL2009 example we analyse the simulated datasets created by Balan & Lahav (2009).
The folder contains two datasets with 1 and 2 simulated planets, which were used to test the ExoFit package.

The kima_setup.cpp file sets the main options for the model: we use a standard sum-of-Keplerians model (no GP) without an offset due to the HARPS change of optical fibers, and no linear trend. The number of planets in the model is fixed to 1. Inside the RVmodel() constructor, we define the priors for the model parameters, the same as those used by Balan & Lahav (2009).

To compile and run, type

make
./run 

CoRoT7

The /CoRoT7 example uses HARPS radial velocity measurements of the active star CoRoT-7.
This example reproduces the analysis done by Faria et al. (2016), where we successfully recovered both the orbits of CoRoT-7b and CoRot-7c, and the activity-induced signal.

In the kima_setup.ccp file, we set a GP model with hyperpriors, no linear trend, and no offset because of the HARPS fibers. The number of planets is free, with a uniform prior between 0 and 5.

Note: in Faria et al. (2016) we considered a uniform prior for Np between 0 and 10.

To compile and run, type

make
./run 

This example takes a considerably longer time to run because of the GP model.

Many Planets

The /many_planetsexample uses the radial velocity measurements of the star HD10180 where seven possible planetary signals were found by Lovis et al. (2011).

Here we also use the sum-of-Keplerians model, with no HARPS offset and no linear trend, but this time we consider hyperpriors for the orbital periods and semi-amplitudes. The number of planets is a free parameter in the MCMC, with a uniform prior between 0 and 10.

To compile and run, type

make
./run