def multipanel_pulsators():
pulsating_allvars = allvars_periods.where(np.in1d(allvars_periods.SOURCEID, pulsating_SOURCEIDs))
pulsating_sdlist = [sdp(photometry_data, x, name='{0}: {1}'.format(y, x)) for x, y in zip(pulsating_allvars.SOURCEID, pulsating_allvars.final_ID)]
# offsets = [0.41, 0.105, 0.1675, -0.3525, -0.305, 0.32, -0.135, -0.08, -0.3025]
plot4.multi_lc_phase_colors(pulsating_sdlist, ['k']*4, pulsating_allvars.best_period, figscale=0.5)#, offsets=offsets)
def multipanel_pulsators():
pulsating_allvars = allvars_periods.where(
np.in1d(allvars_periods.SOURCEID, pulsating_SOURCEIDs))
pulsating_sdlist = [
sdp(photometry_data, x, name='{0}: {1}'.format(y, x))
for x, y in zip(pulsating_allvars.SOURCEID, pulsating_allvars.final_ID)
]
# offsets = [0.41, 0.105, 0.1675, -0.3525, -0.305, 0.32, -0.135, -0.08, -0.3025]
plot4.multi_lc_phase_colors(pulsating_sdlist, ['k'] * 4,
pulsating_allvars.best_period,
figscale=0.5) #, offsets=offsets)
def multipanel_eclipsing_binaries():
EBs_allvars = allvars_periods.where(np.in1d(allvars_periods.SOURCEID, eclipsing_binary_SOURCEIDs))
EB_sdlist = [sdp(photometry_data, x, name='{0}: {1}'.format(y, x)) for x, y in zip(EBs_allvars.SOURCEID, EBs_allvars.final_ID)]
offsets = [0.41, 0.105, 0.1675, -0.3525, -0.305, 0.32, -0.135, -0.08, -0.3025]
fig = plot4.multi_lc_phase_colors(EB_sdlist, ['k']*9, EBs_allvars.best_period, figscale=0.5, offsets=offsets)
return fig
def seven_clean_rotators(cmap=orion_cmap, **kwargs):
rotator_oncvar_ids = [40, 359, 397, 906, 1068, 1091, 1165]
offsets = [0.5, 0.2, 0.1, 0.5, 0, 0.1, 0]
rotator_sourceids = [ukvar_spread['SOURCEID'][ukvar_spread['UKvar_ID'] == oncvar][0]
for oncvar in rotator_oncvar_ids]
rotator_periods = [ukvar_periods[ukvar_spread['UKvar_ID'] == oncvar][0]
for oncvar in rotator_oncvar_ids]
rotator_stardatas = [OrionStarData(variables_photometry, sourceid, name='{}'.format(
oncvar_id)) for sourceid, oncvar_id in zip(rotator_sourceids, rotator_oncvar_ids)]
bands = ['k'] * len(rotator_oncvar_ids)
rotator_fig = multi_lc_phase_colors(
rotator_stardatas, bands, rotator_periods, offsets, cmap=cmap, **kwargs)
for stardata, period, axes_dict in zip(rotator_fig.stardatas, rotator_fig.periods, rotator_fig.axes_dicts):
name = stardata.name
ax_phase = axes_dict['phase']
ax_jhk = axes_dict['jhk']
print "Doing thing with ONCvar {0}!".format(name)
ax_phase.text(0.1, 0.8, "ONCvar {0}".format(
stardata.name), transform=ax_phase.transAxes, fontsize='small')
ax_phase.text(0.6, 0.1, "P = {0:.2f} d".format(
period), transform=ax_phase.transAxes, fontsize='small')
if '359' in stardata.name:
ax_jhk.set_xlim(0.83, 1.04)
ax_jhk.set_ylim(1.623, 1.987)
else:
ax_jhk.set_xlim(0.156, 0.436)
ax_jhk.set_ylim(0.555, 1.055)
# ax_jhk.set_xlim(0.095, 1.05)
# ax_jhk.set_ylim(0.5, 2)
rotator_fig.canvas.draw()
rotator_fig.axes_dicts[0]['phase'].set_title("Period-folded light curve", fontsize='small')
rotator_fig.axes_dicts[0]['lc'].set_title("Light curve (MJD - 54034)", fontsize='small')
rotator_fig.axes_dicts[0]['jhk'].set_title("$J-H$ vs. $H-K$", fontsize='small')
rotator_fig.axes_dicts[0]['khk'].set_title("$K$ vs. $H-K$", fontsize='small')
rotator_fig.canvas.draw()
return rotator_fig
def four_longperiod_variables_bo(cmap=orion_cmap, **kwargs):
fourlong_oncvar_ids = [637, 802, 874, 665]
offsets = [-0.1, 0.3, 0.05, -0.05]
fourlong_sourceids = [ukvar_spread['SOURCEID'][
ukvar_spread['UKvar_ID'] == oncvar][0] for oncvar in fourlong_oncvar_ids]
# try:
# fourlong_periods = [ukvar_periods[ukvar_spread['UKvar_ID'] == oncvar][0] for oncvar in fourlong_oncvar_ids]
fourlong_periods = [42.5, 158.8, 71.5, 44.27]
fourlong_stardatas = [OrionStarData(variables_photometry, sourceid, name='{}'.format(
oncvar_id)) for sourceid, oncvar_id in zip(fourlong_sourceids, fourlong_oncvar_ids)]
bands = ['k'] * 4
four_fig = multi_lc_phase_colors(
fourlong_stardatas, bands, fourlong_periods, offsets, cmap=cmap, **kwargs)
for stardata, period, axes_dict in zip(four_fig.stardatas, four_fig.periods, four_fig.axes_dicts):
name = stardata.name
ax_phase = axes_dict['phase']
ax_jhk = axes_dict['jhk']
print "Doing thing with ONCvar {0}!".format(name)
ax_phase.text(0.1, 0.8, "ONCvar {0}".format(
stardata.name), transform=ax_phase.transAxes, fontsize='small')
ax_phase.text(0.6, 0.8, "P = {0:.2f} d".format(
period), transform=ax_phase.transAxes, fontsize='small')
ax_jhk.set_xlim(0, 2)
ax_jhk.set_ylim(0, 2.5)
four_fig.canvas.draw()
four_fig.axes_dicts[0]['phase'].set_title("Period-folded light curve", fontsize='small')
four_fig.axes_dicts[0]['lc'].set_title("Light curve (MJD - 54034)", fontsize='small')
four_fig.axes_dicts[0]['jhk'].set_title("$J-H$ vs. $H-K$", fontsize='small')
four_fig.axes_dicts[0]['khk'].set_title("$K$ vs. $H-K$", fontsize='small')
four_fig.canvas.draw()
return four_fig
def multipanel_eclipsing_binaries():
EBs_allvars = allvars_periods.where(
np.in1d(allvars_periods.SOURCEID, eclipsing_binary_SOURCEIDs))
EB_sdlist = [
sdp(photometry_data, x, name='{0}: {1}'.format(y, x))
for x, y in zip(EBs_allvars.SOURCEID, EBs_allvars.final_ID)
]
offsets = [
0.41, 0.105, 0.1675, -0.3525, -0.305, 0.32, -0.135, -0.08, -0.3025
]
fig = plot4.multi_lc_phase_colors(EB_sdlist, ['k'] * 9,
EBs_allvars.best_period,
figscale=0.5,
offsets=offsets)
return fig
def aatau_tenpanel(cmap=orion_cmap, **kwargs):
aatau_oncvar_ids = [130, 234, 295, 311, 334, 337, 371, 431, 553, 929]
offsets = [-0.2, -0.15, 0.5, 0.25, 0, -0.25, 0.5, 0.3, 0.1, 0.1]
aatau_sourceids = [ukvar_spread['SOURCEID'][ukvar_spread['UKvar_ID'] == oncvar][0]
for oncvar in aatau_oncvar_ids]
aatau_periods = [ukvar_periods[ukvar_spread['UKvar_ID'] == oncvar][0]
for oncvar in aatau_oncvar_ids]
aatau_stardatas = [OrionStarData(variables_photometry, sourceid, name='{}'.format(
oncvar_id)) for sourceid, oncvar_id in zip(aatau_sourceids, aatau_oncvar_ids)]
bands = ['k'] * 10
aafig = multi_lc_phase_colors(
aatau_stardatas, bands, aatau_periods, offsets, cmap=cmap, **kwargs)
for stardata, period, axes_dict in zip(aafig.stardatas, aafig.periods, aafig.axes_dicts):
name = stardata.name
ax_phase = axes_dict['phase']
ax_jhk = axes_dict['jhk']
print "Doing thing with ONCvar {0}!".format(name)
ax_phase.text(0.1, 0.1, "ONCvar {0}".format(
stardata.name), transform=ax_phase.transAxes, fontsize='small')
ax_phase.text(0.6, 0.1, "P = {0:.2f} d".format(
period), transform=ax_phase.transAxes, fontsize='small')
ax_jhk.set_xlim(0, 1.5)
ax_jhk.set_ylim(0, 2)
aafig.canvas.draw()
aafig.axes_dicts[0]['phase'].set_title("Period-folded light curve", fontsize='small')
aafig.axes_dicts[0]['lc'].set_title("Light curve (MJD - 54034)", fontsize='small')
aafig.axes_dicts[0]['jhk'].set_title("$J-H$ vs. $H-K$", fontsize='small')
aafig.axes_dicts[0]['khk'].set_title("$K$ vs. $H-K$", fontsize='small')
aafig.canvas.draw()
return aafig
def six_old_eclipsing_binaries(cmap=orion_cmap, **kwargs):
oncvar_ids = [122, 427, 518, 1227, 788, 1100]
offsets = [0, 0, -0.15, -0.4, -0.4, -0.15]
sourceids = [ukvar_spread['SOURCEID'][ukvar_spread['UKvar_ID'] == oncvar][0]
for oncvar in oncvar_ids]
periods = [2.635, 20.48, 4.674, 2.65, 9.78, 3.57]
# periods = [ukvar_periods[ukvar_spread['UKvar_ID'] == oncvar][0] for oncvar in oncvar_ids]
stardatas = [OrionStarData(variables_photometry, sourceid, name='{}'.format(
oncvar_id)) for sourceid, oncvar_id in zip(sourceids, oncvar_ids)]
bands = ['k'] * len(oncvar_ids)
fig = multi_lc_phase_colors(stardatas, bands, periods, offsets, cmap=cmap, **kwargs)
for stardata, period, axes_dict in zip(fig.stardatas, fig.periods, fig.axes_dicts):
name = stardata.name
ax_phase = axes_dict['phase']
ax_jhk = axes_dict['jhk']
print "Doing thing with ONCvar {0}!".format(name)
ax_phase.text(0.03, 0.1, "ONCvar {0}".format(
stardata.name), transform=ax_phase.transAxes, fontsize='small')
ax_phase.text(0.7, 0.1, "P = {0:.2f} d".format(
period), transform=ax_phase.transAxes, fontsize='small')
ax_jhk.set_xlim(0, 2)
ax_jhk.set_ylim(0, 2.5)
fig.canvas.draw()
fig.axes_dicts[0]['phase'].set_title("Period-folded light curve", fontsize='small')
fig.axes_dicts[0]['lc'].set_title("Light curve (MJD - 54034)", fontsize='small')
fig.axes_dicts[0]['jhk'].set_title("$J-H$ vs. $H-K$", fontsize='small')
fig.axes_dicts[0]['khk'].set_title("$K$ vs. $H-K$", fontsize='small')
fig.canvas.draw()
return fig
