hk_empty, hk_filled, filter_color_slopes)
from tablemate_comparisons import (mated_ukvar, ukvar_spread,
ukvar_periods, source_period_digger)
from tablemate_script import (Megeath2012, Megeath_P, Megeath_D,
Megeath_Full, Megeath_Allgoodsources,
XMM_north, Rice_UKvars)
from tablemate_core import index_secondary_by_primary, tablemater
# from montage_script import conf_subj_periodics, conf_subj_nonpers
dropbox_bo = os.path.expanduser("~/Dropbox/Bo_Tom/")
output_directory = dropbox_bo+"paper/publication_tables/"
# Let's grab IRAC colors from Megeath.
megeath2012_by_ukvar = index_secondary_by_primary(mated_ukvar,
Megeath2012)
megeath2012_full_by_ukvar = index_secondary_by_primary(mated_ukvar,
Megeath_Full)
megeath2012_all_by_ukvar = index_secondary_by_primary(mated_ukvar,
Megeath_Allgoodsources)
# XMM catalog gets matched too.
XMM_north_by_ukvar = index_secondary_by_primary(
tablemater(Rice_UKvars, [XMM_north]), XMM_north)
# And let's make some color slope references that we like.
jhk_slope_reference = filter_color_slopes(autovars_strict, 'jhk',
slope_confidence=0.5)
jh_slope_reference = filter_color_slopes(autovars_true, 'jh',
slope_confidence=0.5)
hk_slope_reference = filter_color_slopes(autovars_true, 'hk',
from color_slope_filtering import (jhk_empty, jhk_filled, jh_empty, jh_filled,
hk_empty, hk_filled)
from tablemate_comparisons import (mated_ukvar, ukvar_spread,
ukvar_periods, source_period_digger)
from tablemate_script import (Megeath2012, Megeath_P, Megeath_D)
from tablemate_core import index_secondary_by_primary
from variables_data_filterer import filter_by_tile, variables_photometry
from table_maker import make_megeath_class_column
from montage_script import conf_subj_periodics, conf_subj_nonpers
from plot2 import plot_trajectory_vanilla
from helpers3 import band_cut
import robust as rb
# Let's grab IRAC colors from Megeath.
megeath2012_by_ukvar = index_secondary_by_primary(mated_ukvar, Megeath2012)
color_dict = {}
color_dict['disk'] = '#e41a1c' # red
color_dict['nondisk'] = '#377eb8' # blue
color_dict['protostar'] = '#4daf4a' # green
def f_hist_periods():
"""
A figure showing a histogram of periods reported for periodic variables.
"""
fig = plt.figure()
# First group
plt.hist(np.concatenate((conf_subj_periodics.best_period,
