Esempio n. 1
0
from __future__ import (absolute_import, division, print_function,
                        unicode_literals)

# Import dev version of friedrich:
import sys
sys.path.insert(0, '../')
from friedrich.analysis import MCMCResults
from friedrich.lightcurve import hat11_params_morris
from glob import glob
import matplotlib.pyplot as plt
import numpy as np

archive_paths = sorted(glob('/local/tmp/friedrich/hat11/chains???.hdf5'))

#archive_paths = ['/local/tmp/friedrich/hat11/chains043.hdf5']

n_significant_spots = []
indices = []
for archive_path in archive_paths:
    m = MCMCResults(archive_path, hat11_params_morris())
    delta_bics = m.get_spots_delta_chi2(plots=False)
    indices.append(m.index)
    n_significant_spots.append(np.count_nonzero(delta_bics > 10))
    #plt.savefig('spots/{0}.png'.format(m.index))
    #plt.close()
indices = [int(i[6:]) for i in indices]
np.savetxt('nspots.txt', np.vstack([indices, n_significant_spots]).T)
plt.show()
Esempio n. 2
0
elif any(map(lambda x: 'hyak' in x, os.uname())):
    mt_wilson_paths = glob('/gscratch/stf/bmmorris/Mt_Wilson_Tilt/*/sspot??.dat')
    archive_paths = sorted(glob('/gscratch/stf/bmmorris/hat11_flip_lambda/chains???.hdf5'))
else:
    mt_wilson_paths = glob('/astro/store/scratch/tmp/bmmorris/friedrich/Mt_Wilson_Tilt/*/sspot??.dat')
    archive_paths = sorted(glob('/astro/store/scratch/tmp/bmmorris/friedrich/hat11_flip_lambda/chains???.hdf5'))

dropbox_path = '/astro/users/bmmorris/Dropbox/Apps/ShareLaTeX/STSP_HAT-P-11/figures/'

transit_params = hat11_params_morris_experiment()
labels = "a1, mean_latitude, lnv1, lnv2, new_i_s".split(', ')

transits = []
all_times = []
for archive_path in archive_paths:
    m = MCMCResults(archive_path, transit_params)
    all_times.extend(m.lc.times.jd)
    spots = m.get_spots()
    transits.append(Transit(spots))

from friedrich.orientation import times_to_occulted_lat_lon
from friedrich.lightcurve import hat11_params_morris_flip_lambda

all_lats = []
all_lons = []
all_amps = []
all_lats_errors = []
all_spot_times = []
all_BICs = []
for transit in transits:
    for spot in transit.spots:
Esempio n. 3
0
import matplotlib.pyplot as plt

#archive_path = ('/gscratch/stf/bmmorris/friedrich/chains{0:03d}.hdf5'
#                .format(int(sys.argv[1])))
# archive_path = ('/media/PASSPORT/friedrich/chains{0:03d}.hdf5'
#                 .format(int(sys.argv[1])))
#
# m = MCMCResults(archive_path)
# # m.plot_corner()
# # m.plot_lnprob()
# # m.plot_max_lnp_lc(hat11_params_morris())
# # m.plot_lat_lon()
# # m.plot_star()
# m.plot_star_projected()
# plt.savefig('tmp/{0:03d}.png'.format(int(sys.argv[1])))
# #plt.show()

archive_paths = sorted(glob('/local/tmp/friedrich/k17/chains???.hdf5'))
#archive_paths = ['chains017.hdf5']
for archive_path in archive_paths:
    m = MCMCResults(archive_path, k17_params_morris())
    m.plot_star_projected()
    #m.plot_corner()
    #m.plot_lnprob()
    # m.plot_max_lnp_lc()
    transit_number = m.index.split('chains')[1]
    plt.savefig('tmp/{0:03d}.png'.format(int(transit_number)))
    # plt.show()
    plt.close()
#plt.show()
Esempio n. 4
0
import matplotlib.pyplot as plt

#archive_path = ('/gscratch/stf/bmmorris/friedrich/chains{0:03d}.hdf5'
#                .format(int(sys.argv[1])))
# archive_path = ('/media/PASSPORT/friedrich/chains{0:03d}.hdf5'
#                 .format(int(sys.argv[1])))
#
# m = MCMCResults(archive_path)
# # m.plot_corner()
# # m.plot_lnprob()
# # m.plot_max_lnp_lc(hat11_params_morris())
# # m.plot_lat_lon()
# # m.plot_star()
# m.plot_star_projected()
# plt.savefig('tmp/{0:03d}.png'.format(int(sys.argv[1])))
# #plt.show()

archive_paths = sorted(glob('/local/tmp/friedrich/hat11/chains???.hdf5'))
#archive_paths = ['/astro/users/bmmorris/Desktop/chains000.hdf5']
for archive_path in archive_paths:
    m = MCMCResults(archive_path, hat11_params_morris())
    #m.max_lnp_theta_phi()
    m.plot_lnprob()
    m.plot_corner()
    #m.plot_star_projected()
    transit_number = m.index.split('chains')[1]
    #plt.savefig('tmp/{0:03d}.png'.format(int(transit_number)))
    plt.show()
    #plt.close()
#plt.show()