In [1]:
%matplotlib inline

import numpy as np
import matplotlib.pyplot as plt

from isochrones.dartmouth import Dartmouth_Isochrone
from isochrones.starmodel import BinaryStarModel, StarModel
from isochrones.starmodel import TripleStarModel

dar = Dartmouth_Isochrone()

Example: A known EB KIC 6778289, or WISE J192824.57+421508.1 (http://keplerebs.villanova.edu/overview/?k=6778289)


In [2]:
mags = dict( W1=(11.808,0.022), W2=(11.834,0.021), W3=(11.675,0.170), 
            J=(12.119,0.021), H=(11.913,0.017), K=(11.873,0.019))

mod1 = StarModel(dar, maxAV=0.5, **mags)
mod2 = BinaryStarModel(dar, maxAV=0.5, **mags)
mod3 = TripleStarModel(dar, maxAV=0.5, **mags)

In [13]:
#mod1.fit_mcmc()
#mod2.fit_mcmc()
mod3.fit_mcmc(nburn=200, niter=100, nwalkers=300)


Out[13]:
<emcee.ensemble.EnsembleSampler at 0x7f29c6b45d50>

In [4]:
#mod1.triangle_plots();

In [9]:
mod2.triangle_plots();