ecosw = planet_properties['ecosw']
esinw = planet_properties['esinw']
eccentricity = np.sqrt(ecosw**2 + esinw**2)
omega = np.degrees(np.arccos(ecosw/eccentricity))
a, _ = T14b2aRsi(params.per, params.duration, b, params.rp, eccentricity, omega)
params.a = a #semi-major axis (in units of stellar radii)
params.ecc = eccentricity #eccentricity
params.w = omega #longitude of periastron (in degrees)
params.u = map(float, stellar_properties['four_param_limb_darkening'].split(' ')) #limb darkening coefficients
params.limb_dark = "nonlinear" #limb darkening model
return params
transit_params = hat11_params()#get_basic_kepler17_params()
blc = BestLightCurve(best_lc_path, transit_params=transit_params)
blc.plot_whole_lc()
#blc.plot_transits()
#plt.show()
mcmc = MCMCResults(mcmc_paths)
print('Acceptance rate: {0}'.format(str(mcmc.acceptance_rates)))
mcmc.plot_chi2()
mcmc.plot_chains()
#mcmc.plot_star()
#mcmc.plot_corner()
#mcmc.plot_each_spot()
plt.show()
params.rp = 0.13413993 # planet radius (in units of stellar radii)
b = 0.1045441
# inclination = 88.94560#np.arccos(b/params.a)
params.inc = 88.94560 # orbital inclination (in degrees)
inclination = np.radians(params.inc)
params.inclination = params.inc
params.a = b / np.cos(inclination) # semi-major axis (in units of stellar radii)
params.ecc = 0.0 # eccentricity
params.w = 90.0 # longitude of periastron (in degrees)
params.u = [0.1, 0.3] # limb darkening coefficients
params.limb_dark = "quadratic" # limb darkening model
params.duration = (
params.per / np.pi * np.arcsin(np.sqrt((1 + params.rp) ** 2 + b ** 2) / np.sin(inclination) / params.a)
)
return params
transit_params = get_basic_kepler17_params()
blc = BestLightCurve(best_lc_path, transit_params=transit_params)
blc.plot_whole_lc()
blc.plot_transits()
# plt.show()
mcmc = MCMCResults(mcmc_paths)
mcmc.plot_chi2()
mcmc.plot_chains()
# mcmc.plot_star()
# mcmc.plot_corner()
# mcmc.plot_each_spot()
plt.show()
