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 [2]:

    

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)


    

    




0.218014367421 minutes.

In [4]:

    

pop.stars.columns


    

    
Out[4]:





Index([u'B_mag', u'B_mag_A', u'D51_mag', u'D51_mag_A', u'H_mag', u'H_mag_A', u'I_mag', u'I_mag_A', u'J_mag', u'J_mag_A', u'K_mag', u'K_mag_A', u'Kepler_mag', u'Kepler_mag_A', u'R_mag', u'R_mag_A', u'Teff_A', u'U_mag', u'U_mag_A', u'V_mag', u'V_mag_A', u'age_A', u'g_mag', u'g_mag_A', u'i_mag', u'i_mag_A', u'logL_A', u'logg_A', u'mass_A', u'r_mag', u'r_mag_A', u'radius_A', u'z_mag', u'z_mag_A', u'B_mag_B', u'D51_mag_B', u'H_mag_B', u'I_mag_B', u'J_mag_B', u'K_mag_B', u'Kepler_mag_B', u'R_mag_B', u'Teff_B', u'U_mag_B', u'V_mag_B', u'age_B', u'g_mag_B', u'i_mag_B', u'logL_B', u'logg_B', u'mass_B', u'r_mag_B', u'radius_B', u'z_mag_B', u'B_mag_C', u'D51_mag_C', u'H_mag_C', u'I_mag_C', u'J_mag_C', u'K_mag_C', u'Kepler_mag_C', u'R_mag_C', u'Teff_C', u'U_mag_C', u'V_mag_C', u'age_C', u'g_mag_C', u'i_mag_C', u'logL_C', u'logg_C', u'mass_C', u'r_mag_C', u'radius_C', u'z_mag_C', u'distance', u'distmod'], dtype='object')

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 [13]:

    

mags


    

    
Out[13]:





{'H': 10.211,
 'J': 10.523,
 'K': 10.152000000000001,
 'Kepler': 11.664000000000001,
 'Ks': 10.152000000000001,
 'g': 12.042879099999999,
 'h': 10.211,
 'i': 11.4300704,
 'j': 10.523,
 'k': 10.152000000000001,
 'kep': 11.664000000000001,
 'r': 11.596865200000002,
 'z': 11.393061000000001}

In [14]:

    

pop2.ecc_fn


    

    
Out[14]:





<function starutils.utils.draw_eccs>

In [16]:

    

ra, dec


    

    
Out[16]:





(289.21749900000003, 47.884459999999997)

In [ ]: