Retrievals can be run in transmission and emission modes (currently this version only supports transiting planets but directly imaged planets can also be modelled)
To run a retrieval on a spectrum, specify the spectrum path in the parameter file
[Input] spectrum_file = [PATH TO SPECTRUM ASCII FILE] [Output] path = PATH_TO_OUTPUTFOLDER [MultiNest] run = True #alternatively you can run PolyChord (if installed) [PolyChord] run = True
The spectrum file should have 3 to 4 columns:
|wavelength (microns)||(Rp/Rs)2 or Fp/Fs||Error||optional: wavelength bin width (microns)|
The parameter file should contain all necessary attributes to model the spectrum. We have provided examples for 30 planets in transmission in the Tsiaras et al. (2018) paper. The data, parameter files and corresponding retrievals can be found here.
/tests folder provides transmission and emission examples. To run the transmission retrieval case:
python taurex.py -p tests/test_0_transmission/test_retrieval.par --plot
To run the retrieval on multiple cores in parallel, see here.
--plot flag will plot a standard set of output results (posterior distributions, fitted spectrum, opacity contribution breakdown). To plot all available data (e.g. TP-profiles), plots can be
generated after the retrieval has terminated using the created
nest_out.pickle files which contain all results and internal states of the retrieval (see here for plotting and here for
All parameters can be specified in the parameter file or alternatively given as flags in the command line. To print available parameters:
python taurex.py --help