exoplanet-light-curve-analysis's Introduction

This code is no longer being maintained and has been merged into EXOTIC.

The api syntax has remained the same but requires a new import

from elca.tools import ... -> from exotic.api.elca import ...

Exoplanet Light Curve Analyzer

A python 3 package for modeling exoplanet light curves. The transit function is based on the analytic expressions of Mandel and Agol 2002. Estimate parameters and uncertainties using nested sampling.

The data points are color coded to the likelihood value and the contours represent different sigma levels. One sigma uncertainties are reported above each histogram.

Dependencies

Numpy, Matplotlib, Cython, Ultranest

Installation

pip install git+https://github.com/pearsonkyle/Exoplanet-Light-Curve-Analysis.git

Examples

Quickstart

import numpy as np
import matplotlib.pyplot as plt

from elca.tools import transit, lc_fitter

if __name__ == "__main__":

    prior = { 
        'rprs':0.05,        # Rp/Rs
        'ars':14.5,         # a/Rs
        'per':3.336817,     # Period [day]
        'inc':86.5,         # Inclination [deg]
        'u1': 0.3, 'u2': 0.1, # limb darkening
        'ecc':0,            # Eccentricity
        'omega':0,          # Arg of periastron
        'tmid':0.75         # time of mid transit [day]
    } 

    # GENERATE NOISY DATA
    time = np.linspace(0.65,0.85,500) # [day]
    data = transit(time, prior) + np.random.normal(0, 2e-4, len(time))
    dataerr = np.random.normal(300e-6, 50e-6, len(time))

    mybounds = {
        'rprs':[0,2*prior['rprs']],
        'tmid':[min(time),max(time)],
        'ars':[13,16], 
        #'inc':[85,87]
    }

    myfit = lc_fitter(time, data, dataerr, prior, mybounds)
    
    for k in myfit.bounds.keys():
        print("{:.6f} +- {}".format( myfit.parameters[k], myfit.errors[k]))

    fig,axs = myfit.plot_bestfit()
    plt.show()

    myfit.plot_triangle()
    plt.show()