In [1]:

    

%matplotlib inline

import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
import time

from keputils import koiutils as ku

from starutils.populations import ColormatchMultipleStarPopulation, BGStarPopulation_TRILEGAL


    

In [2]:

    

ra,dec = ku.radec(87) #kepler-22
bgpop = BGStarPopulation_TRILEGAL('kepler22b.h5', ra, dec)


    

In [3]:

    

mags = ku.KICmags(87) #kepler-22
start = time.time()
#pop = ColormatchMultipleStarPopulation(mags, starfield='kepler22b.h5', n=2e4)
pop = ColormatchMultipleStarPopulation(mags, m1=(1,0.1),
                                       age=(9.7,0.1),
                                       feh=(0,0.1), n=2e4)
end = time.time()
print '{} minutes.'.format((end-start)/60)


    

    




16723
20661
0.209011352062 minutes.

In [4]:

    

bgpop.save_hdf('bgpoptest.h5', overwrite=True)
pop.save_hdf('colormatchpoptest.h5', overwrite=True)


    

In [5]:

    

bgpop2 = BGStarPopulation_TRILEGAL().load_hdf('bgpoptest.h5')


    

In [6]:

    

pop2 = ColormatchMultipleStarPopulation().load_hdf('colormatchpoptest.h5')


    

In [ ]: