%matplotlib inline
import numpy as np
from isochrones.dartmouth import Dartmouth_Isochrone
from isochrones.starmodel import BinaryStarModel, StarModel
from isochrones.starmodel import TripleStarModel
dar = Dartmouth_Isochrone()
mass = np.ones(100)
age, feh = (9.5,0.0)
AV = np.ones(100)*0.1
dist = np.arange(100)+1
dar(mass,age, feh, AV=AV, distance=dist)
B_mag | D51_mag | H_mag | I_mag | J_mag | K_mag | Kepler_mag | R_mag | Teff | U_mag | ... | W3_mag | age | g_mag | i_mag | logL | logg | mass | r_mag | radius | z_mag | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
0 | 0.762614 | 0.205111 | -1.549342 | -0.741485 | -1.216058 | -1.592913 | -0.232721 | -0.378845 | 5751.548729 | 1.063375 | ... | -1.620139 | 9.5 | 0.359039 | -0.316348 | -0.048279 | 4.478884 | 1 | -0.182274 | 0.954492 | -0.345894 |
1 | 2.267764 | 1.710261 | -0.044192 | 0.763665 | 0.289092 | -0.087763 | 1.272429 | 1.126305 | 5751.548729 | 2.568525 | ... | -0.114989 | 9.5 | 1.864189 | 1.188802 | -0.048279 | 4.478884 | 1 | 1.322876 | 0.954492 | 1.159256 |
2 | 3.148220 | 2.590717 | 0.836264 | 1.644121 | 1.169548 | 0.792693 | 2.152885 | 2.006761 | 5751.548729 | 3.448981 | ... | 0.765467 | 9.5 | 2.744645 | 2.069258 | -0.048279 | 4.478884 | 1 | 2.203332 | 0.954492 | 2.039713 |
3 | 3.772914 | 3.215411 | 1.460958 | 2.268815 | 1.794242 | 1.417387 | 2.777579 | 2.631455 | 5751.548729 | 4.073675 | ... | 1.390161 | 9.5 | 3.369339 | 2.693952 | -0.048279 | 4.478884 | 1 | 2.828026 | 0.954492 | 2.664406 |
4 | 4.257464 | 3.699961 | 1.945508 | 2.753365 | 2.278792 | 1.901937 | 3.262129 | 3.116005 | 5751.548729 | 4.558225 | ... | 1.874711 | 9.5 | 3.853889 | 3.178502 | -0.048279 | 4.478884 | 1 | 3.312576 | 0.954492 | 3.148956 |
5 | 4.653370 | 4.095867 | 2.341414 | 3.149271 | 2.674698 | 2.297843 | 3.658035 | 3.511911 | 5751.548729 | 4.954131 | ... | 2.270617 | 9.5 | 4.249795 | 3.574408 | -0.048279 | 4.478884 | 1 | 3.708482 | 0.954492 | 3.544863 |
6 | 4.988104 | 4.430601 | 2.676148 | 3.484005 | 3.009432 | 2.632577 | 3.992769 | 3.846645 | 5751.548729 | 5.288865 | ... | 2.605351 | 9.5 | 4.584529 | 3.909142 | -0.048279 | 4.478884 | 1 | 4.043216 | 0.954492 | 3.879597 |
7 | 5.278064 | 4.720561 | 2.966108 | 3.773965 | 3.299392 | 2.922537 | 4.282729 | 4.136605 | 5751.548729 | 5.578825 | ... | 2.895311 | 9.5 | 4.874489 | 4.199102 | -0.048279 | 4.478884 | 1 | 4.333176 | 0.954492 | 4.169556 |
8 | 5.533827 | 4.976324 | 3.221871 | 4.029727 | 3.555154 | 3.178299 | 4.538491 | 4.392367 | 5751.548729 | 5.834588 | ... | 3.151073 | 9.5 | 5.130251 | 4.454864 | -0.048279 | 4.478884 | 1 | 4.588938 | 0.954492 | 4.425319 |
9 | 5.762614 | 5.205111 | 3.450658 | 4.258515 | 3.783942 | 3.407087 | 4.767279 | 4.621155 | 5751.548729 | 6.063375 | ... | 3.379861 | 9.5 | 5.359039 | 4.683652 | -0.048279 | 4.478884 | 1 | 4.817726 | 0.954492 | 4.654106 |
10 | 5.969578 | 5.412075 | 3.657621 | 4.465478 | 3.990905 | 3.614050 | 4.974242 | 4.828118 | 5751.548729 | 6.270338 | ... | 3.586824 | 9.5 | 5.566002 | 4.890615 | -0.048279 | 4.478884 | 1 | 5.024689 | 0.954492 | 4.861070 |
11 | 6.158520 | 5.601017 | 3.846564 | 4.654421 | 4.179848 | 3.802993 | 5.163185 | 5.017061 | 5751.548729 | 6.459281 | ... | 3.775767 | 9.5 | 5.754945 | 5.079558 | -0.048279 | 4.478884 | 1 | 5.213632 | 0.954492 | 5.050013 |
12 | 6.332331 | 5.774828 | 4.020375 | 4.828232 | 4.353659 | 3.976804 | 5.336996 | 5.190871 | 5751.548729 | 6.633092 | ... | 3.949577 | 9.5 | 5.928755 | 5.253369 | -0.048279 | 4.478884 | 1 | 5.387443 | 0.954492 | 5.223823 |
13 | 6.493254 | 5.935751 | 4.181298 | 4.989155 | 4.514582 | 4.137727 | 5.497919 | 5.351795 | 5751.548729 | 6.794015 | ... | 4.110501 | 9.5 | 6.089679 | 5.414292 | -0.048279 | 4.478884 | 1 | 5.548366 | 0.954492 | 5.384747 |
14 | 6.643070 | 6.085567 | 4.331114 | 5.138971 | 4.664398 | 4.287543 | 5.647735 | 5.501611 | 5751.548729 | 6.943831 | ... | 4.260317 | 9.5 | 6.239495 | 5.564108 | -0.048279 | 4.478884 | 1 | 5.698182 | 0.954492 | 5.534563 |
15 | 6.783214 | 6.225711 | 4.471258 | 5.279115 | 4.804542 | 4.427687 | 5.787879 | 5.641755 | 5751.548729 | 7.083975 | ... | 4.400461 | 9.5 | 6.379639 | 5.704252 | -0.048279 | 4.478884 | 1 | 5.838326 | 0.954492 | 5.674706 |
16 | 6.914859 | 6.357356 | 4.602903 | 5.410759 | 4.936187 | 4.559331 | 5.919523 | 5.773399 | 5751.548729 | 7.215620 | ... | 4.532105 | 9.5 | 6.511283 | 5.835896 | -0.048279 | 4.478884 | 1 | 5.969971 | 0.954492 | 5.806351 |
17 | 7.038977 | 6.481474 | 4.727021 | 5.534877 | 5.060304 | 4.683449 | 6.043641 | 5.897517 | 5751.548729 | 7.339738 | ... | 4.656223 | 9.5 | 6.635401 | 5.960014 | -0.048279 | 4.478884 | 1 | 6.094088 | 0.954492 | 5.930469 |
18 | 7.156382 | 6.598879 | 4.844426 | 5.652283 | 5.177710 | 4.800855 | 6.161047 | 6.014923 | 5751.548729 | 7.457143 | ... | 4.773629 | 9.5 | 6.752807 | 6.077420 | -0.048279 | 4.478884 | 1 | 6.211494 | 0.954492 | 6.047874 |
19 | 7.267764 | 6.710261 | 4.955808 | 5.763665 | 5.289092 | 4.912237 | 6.272429 | 6.126305 | 5751.548729 | 7.568525 | ... | 4.885011 | 9.5 | 6.864189 | 6.188802 | -0.048279 | 4.478884 | 1 | 6.322876 | 0.954492 | 6.159256 |
20 | 7.373711 | 6.816208 | 5.061755 | 5.869611 | 5.395038 | 5.018183 | 6.378375 | 6.232251 | 5751.548729 | 7.674472 | ... | 4.990957 | 9.5 | 6.970135 | 6.294748 | -0.048279 | 4.478884 | 1 | 6.428822 | 0.954492 | 6.265203 |
21 | 7.474728 | 6.917225 | 5.162771 | 5.970628 | 5.496055 | 5.119200 | 6.479392 | 6.333268 | 5751.548729 | 7.775488 | ... | 5.091974 | 9.5 | 7.071152 | 6.395765 | -0.048279 | 4.478884 | 1 | 6.529839 | 0.954492 | 6.366220 |
22 | 7.571253 | 7.013750 | 5.259297 | 6.067154 | 5.592581 | 5.215726 | 6.575918 | 6.429794 | 5751.548729 | 7.872014 | ... | 5.188500 | 9.5 | 7.167678 | 6.492291 | -0.048279 | 4.478884 | 1 | 6.626365 | 0.954492 | 6.462746 |
23 | 7.663670 | 7.106167 | 5.351714 | 6.159571 | 5.684998 | 5.308143 | 6.668335 | 6.522211 | 5751.548729 | 7.964431 | ... | 5.280917 | 9.5 | 7.260095 | 6.584708 | -0.048279 | 4.478884 | 1 | 6.718782 | 0.954492 | 6.555163 |
24 | 7.752314 | 7.194811 | 5.440358 | 6.248215 | 5.773642 | 5.396787 | 6.756979 | 6.610855 | 5751.548729 | 8.053075 | ... | 5.369561 | 9.5 | 7.348739 | 6.673352 | -0.048279 | 4.478884 | 1 | 6.807426 | 0.954492 | 6.643807 |
25 | 7.837481 | 7.279978 | 5.525525 | 6.333382 | 5.858809 | 5.481954 | 6.842146 | 6.696021 | 5751.548729 | 8.138242 | ... | 5.454727 | 9.5 | 7.433905 | 6.758519 | -0.048279 | 4.478884 | 1 | 6.892593 | 0.954492 | 6.728973 |
26 | 7.919433 | 7.361930 | 5.607477 | 6.415334 | 5.940761 | 5.563906 | 6.924098 | 6.777973 | 5751.548729 | 8.220194 | ... | 5.536679 | 9.5 | 7.515858 | 6.840471 | -0.048279 | 4.478884 | 1 | 6.974545 | 0.954492 | 6.810925 |
27 | 7.998404 | 7.440901 | 5.686448 | 6.494305 | 6.019732 | 5.642877 | 7.003069 | 6.856945 | 5751.548729 | 8.299165 | ... | 5.615651 | 9.5 | 7.594829 | 6.919442 | -0.048279 | 4.478884 | 1 | 7.053516 | 0.954492 | 6.889897 |
28 | 8.074604 | 7.517101 | 5.762648 | 6.570505 | 6.095932 | 5.719077 | 7.079269 | 6.933145 | 5751.548729 | 8.375365 | ... | 5.691851 | 9.5 | 7.671029 | 6.995642 | -0.048279 | 4.478884 | 1 | 7.129716 | 0.954492 | 6.966096 |
29 | 8.148220 | 7.590717 | 5.836264 | 6.644121 | 6.169548 | 5.792693 | 7.152885 | 7.006761 | 5751.548729 | 8.448981 | ... | 5.765467 | 9.5 | 7.744645 | 7.069258 | -0.048279 | 4.478884 | 1 | 7.203332 | 0.954492 | 7.039713 |
... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
70 | 10.018906 | 9.461403 | 7.706950 | 8.514807 | 8.040234 | 7.663379 | 9.023571 | 8.877446 | 5751.548729 | 10.319667 | ... | 7.636152 | 9.5 | 9.615330 | 8.939944 | -0.048279 | 4.478884 | 1 | 9.074018 | 0.954492 | 8.910398 |
71 | 10.049277 | 9.491774 | 7.737321 | 8.545177 | 8.070604 | 7.693749 | 9.053941 | 8.907817 | 5751.548729 | 10.350038 | ... | 7.666523 | 9.5 | 9.645701 | 8.970314 | -0.048279 | 4.478884 | 1 | 9.104388 | 0.954492 | 8.940769 |
72 | 10.079228 | 9.521725 | 7.767272 | 8.575129 | 8.100556 | 7.723701 | 9.083893 | 8.937769 | 5751.548729 | 10.379989 | ... | 7.696475 | 9.5 | 9.675653 | 9.000266 | -0.048279 | 4.478884 | 1 | 9.134340 | 0.954492 | 8.970721 |
73 | 10.108773 | 9.551270 | 7.796817 | 8.604673 | 8.130100 | 7.753245 | 9.113437 | 8.967313 | 5751.548729 | 10.409534 | ... | 7.726019 | 9.5 | 9.705197 | 9.029810 | -0.048279 | 4.478884 | 1 | 9.163885 | 0.954492 | 9.000265 |
74 | 10.137920 | 9.580417 | 7.825964 | 8.633821 | 8.159248 | 7.782393 | 9.142585 | 8.996461 | 5751.548729 | 10.438681 | ... | 7.755167 | 9.5 | 9.734345 | 9.058958 | -0.048279 | 4.478884 | 1 | 9.193032 | 0.954492 | 9.029413 |
75 | 10.166682 | 9.609179 | 7.854726 | 8.662583 | 8.188010 | 7.811155 | 9.171347 | 9.025223 | 5751.548729 | 10.467443 | ... | 7.783929 | 9.5 | 9.763107 | 9.087720 | -0.048279 | 4.478884 | 1 | 9.221794 | 0.954492 | 9.058174 |
76 | 10.195068 | 9.637565 | 7.883112 | 8.690968 | 8.216396 | 7.839540 | 9.199732 | 9.053608 | 5751.548729 | 10.495829 | ... | 7.812314 | 9.5 | 9.791492 | 9.116105 | -0.048279 | 4.478884 | 1 | 9.250180 | 0.954492 | 9.086560 |
77 | 10.223087 | 9.665584 | 7.911131 | 8.718988 | 8.244415 | 7.867560 | 9.227752 | 9.081628 | 5751.548729 | 10.523848 | ... | 7.840334 | 9.5 | 9.819512 | 9.144125 | -0.048279 | 4.478884 | 1 | 9.278199 | 0.954492 | 9.114579 |
78 | 10.250750 | 9.693247 | 7.938794 | 8.746650 | 8.272077 | 7.895222 | 9.255414 | 9.109290 | 5751.548729 | 10.551510 | ... | 7.867996 | 9.5 | 9.847174 | 9.171787 | -0.048279 | 4.478884 | 1 | 9.305861 | 0.954492 | 9.142242 |
79 | 10.278064 | 9.720561 | 7.966108 | 8.773965 | 8.299392 | 7.922537 | 9.282729 | 9.136605 | 5751.548729 | 10.578825 | ... | 7.895311 | 9.5 | 9.874489 | 9.199102 | -0.048279 | 4.478884 | 1 | 9.333176 | 0.954492 | 9.169556 |
80 | 10.305039 | 9.747536 | 7.993083 | 8.800940 | 8.326367 | 7.949512 | 9.309704 | 9.163580 | 5751.548729 | 10.605800 | ... | 7.922286 | 9.5 | 9.901464 | 9.226077 | -0.048279 | 4.478884 | 1 | 9.360151 | 0.954492 | 9.196532 |
81 | 10.331683 | 9.774180 | 8.019727 | 8.827584 | 8.353011 | 7.976156 | 9.336348 | 9.190224 | 5751.548729 | 10.632444 | ... | 7.948930 | 9.5 | 9.928108 | 9.252721 | -0.048279 | 4.478884 | 1 | 9.386795 | 0.954492 | 9.223176 |
82 | 10.358005 | 9.800502 | 8.046049 | 8.853905 | 8.379332 | 8.002477 | 9.362669 | 9.216545 | 5751.548729 | 10.658765 | ... | 7.975251 | 9.5 | 9.954429 | 9.279042 | -0.048279 | 4.478884 | 1 | 9.413116 | 0.954492 | 9.249497 |
83 | 10.384011 | 9.826508 | 8.072054 | 8.879911 | 8.405338 | 8.028483 | 9.388675 | 9.242551 | 5751.548729 | 10.684771 | ... | 8.001257 | 9.5 | 9.980435 | 9.305048 | -0.048279 | 4.478884 | 1 | 9.439122 | 0.954492 | 9.275503 |
84 | 10.409709 | 9.852206 | 8.097753 | 8.905609 | 8.431037 | 8.054181 | 9.414373 | 9.268249 | 5751.548729 | 10.710470 | ... | 8.026955 | 9.5 | 10.006133 | 9.330746 | -0.048279 | 4.478884 | 1 | 9.464821 | 0.954492 | 9.301201 |
85 | 10.435106 | 9.877603 | 8.123150 | 8.931007 | 8.456434 | 8.079579 | 9.439771 | 9.293647 | 5751.548729 | 10.735867 | ... | 8.052353 | 9.5 | 10.031531 | 9.356144 | -0.048279 | 4.478884 | 1 | 9.490218 | 0.954492 | 9.326599 |
86 | 10.460210 | 9.902707 | 8.148254 | 8.956111 | 8.481538 | 8.104683 | 9.464875 | 9.318751 | 5751.548729 | 10.760971 | ... | 8.077457 | 9.5 | 10.056635 | 9.381248 | -0.048279 | 4.478884 | 1 | 9.515322 | 0.954492 | 9.351703 |
87 | 10.485028 | 9.927524 | 8.173071 | 8.980928 | 8.506355 | 8.129500 | 9.489692 | 9.343568 | 5751.548729 | 10.785788 | ... | 8.102274 | 9.5 | 10.081452 | 9.406065 | -0.048279 | 4.478884 | 1 | 9.540139 | 0.954492 | 9.376520 |
88 | 10.509564 | 9.952061 | 8.197608 | 9.005465 | 8.530892 | 8.154037 | 9.514229 | 9.368105 | 5751.548729 | 10.810325 | ... | 8.126811 | 9.5 | 10.105989 | 9.430602 | -0.048279 | 4.478884 | 1 | 9.564676 | 0.954492 | 9.401057 |
89 | 10.533827 | 9.976324 | 8.221871 | 9.029727 | 8.555154 | 8.178299 | 9.538491 | 9.392367 | 5751.548729 | 10.834588 | ... | 8.151073 | 9.5 | 10.130251 | 9.454864 | -0.048279 | 4.478884 | 1 | 9.588938 | 0.954492 | 9.425319 |
90 | 10.557821 | 10.000318 | 8.245865 | 9.053722 | 8.579149 | 8.202294 | 9.562486 | 9.416362 | 5751.548729 | 10.858582 | ... | 8.175068 | 9.5 | 10.154246 | 9.478859 | -0.048279 | 4.478884 | 1 | 9.612933 | 0.954492 | 9.449313 |
91 | 10.581553 | 10.024050 | 8.269597 | 9.077454 | 8.602881 | 8.226026 | 9.586218 | 9.440094 | 5751.548729 | 10.882314 | ... | 8.198800 | 9.5 | 10.177978 | 9.502591 | -0.048279 | 4.478884 | 1 | 9.636665 | 0.954492 | 9.473046 |
92 | 10.605029 | 10.047526 | 8.293073 | 9.100930 | 8.626357 | 8.249502 | 9.609694 | 9.463569 | 5751.548729 | 10.905790 | ... | 8.222275 | 9.5 | 10.201453 | 9.526067 | -0.048279 | 4.478884 | 1 | 9.660141 | 0.954492 | 9.496521 |
93 | 10.628253 | 10.070750 | 8.316297 | 9.124154 | 8.649581 | 8.272726 | 9.632918 | 9.486794 | 5751.548729 | 10.929014 | ... | 8.245500 | 9.5 | 10.224678 | 9.549291 | -0.048279 | 4.478884 | 1 | 9.683365 | 0.954492 | 9.519746 |
94 | 10.651232 | 10.093729 | 8.339276 | 9.147133 | 8.672560 | 8.295705 | 9.655897 | 9.509773 | 5751.548729 | 10.951993 | ... | 8.268479 | 9.5 | 10.247657 | 9.572270 | -0.048279 | 4.478884 | 1 | 9.706344 | 0.954492 | 9.542725 |
95 | 10.673970 | 10.116467 | 8.362014 | 9.169871 | 8.695298 | 8.318443 | 9.678635 | 9.532511 | 5751.548729 | 10.974731 | ... | 8.291217 | 9.5 | 10.270395 | 9.595008 | -0.048279 | 4.478884 | 1 | 9.729082 | 0.954492 | 9.565463 |
96 | 10.696473 | 10.138970 | 8.384517 | 9.192373 | 8.717801 | 8.340945 | 9.701137 | 9.555013 | 5751.548729 | 10.997234 | ... | 8.313719 | 9.5 | 10.292897 | 9.617511 | -0.048279 | 4.478884 | 1 | 9.751585 | 0.954492 | 9.587965 |
97 | 10.718745 | 10.161241 | 8.406788 | 9.214645 | 8.740072 | 8.363217 | 9.723409 | 9.577285 | 5751.548729 | 11.019505 | ... | 8.335991 | 9.5 | 10.315169 | 9.639782 | -0.048279 | 4.478884 | 1 | 9.773856 | 0.954492 | 9.610237 |
98 | 10.740790 | 10.183287 | 8.428834 | 9.236691 | 8.762118 | 8.385263 | 9.745455 | 9.599331 | 5751.548729 | 11.041551 | ... | 8.358037 | 9.5 | 10.337215 | 9.661828 | -0.048279 | 4.478884 | 1 | 9.795902 | 0.954492 | 9.632282 |
99 | 10.762614 | 10.205111 | 8.450658 | 9.258515 | 8.783942 | 8.407087 | 9.767279 | 9.621155 | 5751.548729 | 11.063375 | ... | 8.379861 | 9.5 | 10.359039 | 9.683652 | -0.048279 | 4.478884 | 1 | 9.817726 | 0.954492 | 9.654106 |
100 rows × 23 columns
Example: A known EB KIC 6778289, or WISE J192824.57+421508.1 (http://keplerebs.villanova.edu/overview/?k=6778289)
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)
#mod1.fit_mcmc()
mod2.fit_mcmc()
#mod3.fit_mcmc()
<emcee.ensemble.EnsembleSampler at 0x1022a74d0>
mod2._make_samples(lnprob_thresh=0.01)
mod2.triangle_plots(extent=0.99);
mod1.triangle_plots();
mod3.triangle_plots();