planetplanet is a general photodynamical code for modeling exoplanet transits, secondary eclipses, phase curves, and exomoons, as well as eclipsing binaries, circumbinary planets, and more. The code was originally developed to model planet-planet occultation (PPO) light curves for the TRAPPIST-1 system. During a PPO, a planet occults (transits) the disk of another planet in the same planetary system, blocking its thermal (and reflected) light, which can be measured photometrically by a distant observer.

planetplanet is coded in C and wrapped in a user-friendly Python interface. Once installed, generating light curves is as easy as

import planetplanet as pp
import numpy as np
import matplotlib.pyplot as pl

star = pp.Star('A', m = 0.1, r = 0.1, limbdark = [0.4, 0.26])
planet = pp.Planet('b', m = 1, r = 1, t0 = 0., per = 3.)
system = pp.System(star, planet)
time = np.arange(-1., 1., 0.0001)
system.compute(time)
system.plot_occultation('A', time = 0)
pl.show()

Please check out the documentation or read the paper for more information.

The planetplanet code is now pip-installable:

pip install planetplanet

Alternatively, to install from source:

git clone git@github.com:rodluger/planetplanet.git
cd planetplanet
git submodule init && git submodule update
python setup.py develop

Note that you may need to install the GNU Scientific Library. On a Mac, it's as simple as

brew install gsl

Here's a an example of a planet-planet occultation [code]:

And here's a wacky example of a transit of a circumbinary exomoon [code]: