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()