In [1]:
from isochrones import StarModel
from isochrones.dartmouth import Dartmouth_Isochrone
In [2]:
dar = Dartmouth_Isochrone()
In [6]:
dar.radius(1.0, 9.6, 0.0)
Out[6]:
0.97952282264658863
In [7]:
dar.bands
Out[7]:
['g', 'i', 'H', 'K', 'J', 'r', 'W3', 'W2', 'W1', 'Kepler', 'z']
In [8]:
dar.mag['g'](1.0, 9.6, 0.0)
Out[8]:
5.1637306212564864
In [10]:
%%file test_star.ini
maxAV = 1.216
RA = 299.268036
dec = 45.227428
Teff = 4135, 98.0
feh = -0.46, 0.16
logg = 4.711, 0.1
[KIC]
g = 16.5832173, 0.05
r = 15.341168, 0.05
i = 14.917236, 0.05
z = 14.6382752, 0.05
[twomass]
J = 13.513, 0.02
H = 12.845, 0.02
K = 12.693, 0.02
[NIRC2]
resolution = 0.1
separation_1 = 0.6
PA_1 = 100
K_1 = 3.66, 0.05
H_1 = 3.77, 0.03
J_1 = 3.74, 0.05
#separation_2 = 1.2
#PA_2 = 200
#K_2 = 5.1, 0.1
#H_2 = 5.2, 0.1
#J_2 = 5.15, 0.1
[RAO]
resolution=0.15
separation_1=0.65
PA_1 = 105
i_1 = 4.2, 0.02
Overwriting test_star.ini
In [11]:
mod = StarModel.from_ini(Dartmouth_Isochrone, ini_file='test_star.ini')
In [12]:
mod.obs.print_ascii()
root
╚═ KIC g=(16.58, 0.05) @(0.00, 0 [4.00])
╚═ KIC i=(14.92, 0.05) @(0.00, 0 [4.00])
╚═ KIC r=(15.34, 0.05) @(0.00, 0 [4.00])
╚═ KIC z=(14.64, 0.05) @(0.00, 0 [4.00])
╚═ twomass H=(12.85, 0.02) @(0.00, 0 [4.00])
╚═ twomass J=(13.51, 0.02) @(0.00, 0 [4.00])
╚═ twomass K=(12.69, 0.02) @(0.00, 0 [4.00])
╠═ RAO delta-i=(0.00, 0.01) @(0.00, 0 [0.15])
║ ╚═ NIRC2 delta-H=(0.00, 0.01) @(0.00, 0 [0.10])
║ ╚═ NIRC2 delta-J=(0.00, 0.01) @(0.00, 0 [0.10])
║ ╚═ NIRC2 delta-K=(0.00, 0.01) @(0.00, 0 [0.10])
║ ╚═ 0_0, logg=[4.711, 0.1], Teff=[4135.0, 98.0], feh=[-0.46, 0.16]
╚═ RAO delta-i=(4.20, 0.02) @(0.65, 105 [0.15])
╚═ NIRC2 delta-H=(3.77, 0.03) @(0.60, 100 [0.10])
╚═ NIRC2 delta-J=(3.74, 0.05) @(0.60, 100 [0.10])
╚═ NIRC2 delta-K=(3.66, 0.05) @(0.60, 100 [0.10])
╚═ 0_1
In [13]:
mod.n_params
Out[13]:
6
In [14]:
mod.labelstring
Out[14]:
'0_0--0_1'
In [ ]:
# normal parameter order: [mass, log10(age in Gyr), feh, distance, A_V]
# for an associated binary (0_0 + 0_1): [mass_0, mass_1, .. , .., .., ..]
In [15]:
mod.fit()
In [9]:
%matplotlib inline
mod.corner_physical();
In [10]:
mod.corner_observed();
In [16]:
mod2 = StarModel.from_ini(Dartmouth_Isochrone, ini_file='test_star.ini', index=[0,1])
mod2.obs.print_ascii()
root
╚═ KIC g=(16.58, 0.05) @(0.00, 0 [4.00])
╚═ KIC i=(14.92, 0.05) @(0.00, 0 [4.00])
╚═ KIC r=(15.34, 0.05) @(0.00, 0 [4.00])
╚═ KIC z=(14.64, 0.05) @(0.00, 0 [4.00])
╚═ twomass H=(12.85, 0.02) @(0.00, 0 [4.00])
╚═ twomass J=(13.51, 0.02) @(0.00, 0 [4.00])
╚═ twomass K=(12.69, 0.02) @(0.00, 0 [4.00])
╠═ RAO delta-i=(0.00, 0.01) @(0.00, 0 [0.15])
║ ╚═ NIRC2 delta-H=(0.00, 0.01) @(0.00, 0 [0.10])
║ ╚═ NIRC2 delta-J=(0.00, 0.01) @(0.00, 0 [0.10])
║ ╚═ NIRC2 delta-K=(0.00, 0.01) @(0.00, 0 [0.10])
║ ╚═ 0_0, logg=[4.711, 0.1], Teff=[4135.0, 98.0], feh=[-0.46, 0.16]
╚═ RAO delta-i=(4.20, 0.02) @(0.65, 105 [0.15])
╚═ NIRC2 delta-H=(3.77, 0.03) @(0.60, 100 [0.10])
╚═ NIRC2 delta-J=(3.74, 0.05) @(0.60, 100 [0.10])
╚═ NIRC2 delta-K=(3.66, 0.05) @(0.60, 100 [0.10])
╚═ 1_0
In [17]:
mod2.n_params
Out[17]:
10
In [ ]:
mod2.fit()
ERROR: Interrupt received: Terminating
In [3]:
mod3 = StarModel.from_ini(Dartmouth_Isochrone, ini_file='test_star.ini', index=[0,1], N=[2,1])
mod3.obs.print_ascii()
root
╚═ KIC g=(16.58, 0.05) @(0.00, 0 [4.00])
╚═ KIC i=(14.92, 0.05) @(0.00, 0 [4.00])
╚═ KIC r=(15.34, 0.05) @(0.00, 0 [4.00])
╚═ KIC z=(14.64, 0.05) @(0.00, 0 [4.00])
╚═ twomass H=(12.85, 0.02) @(0.00, 0 [4.00])
╚═ twomass J=(13.51, 0.02) @(0.00, 0 [4.00])
╚═ twomass K=(12.69, 0.02) @(0.00, 0 [4.00])
╠═ RAO delta-i=(0.00, 0.01) @(0.00, 0 [0.15])
║ ╚═ NIRC2 delta-H=(0.00, 0.01) @(0.00, 0 [0.10])
║ ╚═ NIRC2 delta-J=(0.00, 0.01) @(0.00, 0 [0.10])
║ ╚═ NIRC2 delta-K=(0.00, 0.01) @(0.00, 0 [0.10])
║ ╠═ 0_0, logg=[4.711, 0.1], Teff=[4135.0, 98.0], feh=[-0.46, 0.16]
║ ╚═ 0_1
╚═ RAO delta-i=(4.20, 0.02) @(0.65, 105 [0.15])
╚═ NIRC2 delta-H=(3.77, 0.03) @(0.60, 100 [0.10])
╚═ NIRC2 delta-J=(3.74, 0.05) @(0.60, 100 [0.10])
╚═ NIRC2 delta-K=(3.66, 0.05) @(0.60, 100 [0.10])
╚═ 1_0
In [4]:
mod3.n_params
Out[4]:
11
In [ ]: