Python code computing the probability of planet-planet eclipses (PPE).
This code computes the probability of planet-planet eclipses from a population of objects defined in a CSV file. The format of the CSV file is the same as the one provided by the NASA Exoplanet Archive when downloading the table of planetary systems (e.g. confirmed planets).
The code requires the following packages: pickle, numpy, astropy, pandas, multiprocessing
To run the code, one must first modify the parameter file PARAMS.txt accordingly, and then type python3 main.py PARAMS.txt in a terminal window.
The SELECTION parameter can be set to one of the three following options:
-
LAMBDA_ALL: the code will only consider planet-planet pairs for which both projected obliquities ($\lambda$ ) are known. -
LAMBDA_NOTALL: the code will consider planet-planet pairs for which one of two projected obliquities is known. -
LAMBDA_NONE: the code will only consider planet-planet pairs with both projected obliquities unknown.
Note that the typical computation time for a single planet-planet pair is of several minutes.
To create graphs showing the probability of PPE events as a function of time, type python3 createGraphs.py <inputDir> <outputDir> in a terminal window.
This code has been developed by Maxime Marchand as part of a Master project at the Department of Astronomy of the University of Geneva.