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