def find_20_50_day_periodics_and_save_to_file(): blp = best_long_period(long_periodic_selector(maxvars_spread_per, max_period=100), max_period=100) ukvar_blp = blp.where(np.in1d(blp.SOURCEID, ukvar_spread.SOURCEID)) print "There are {0} stars with nominal periods between 50 and 100.".format(len(ukvar_blp)) for sid, period in zip(ukvar_blp.SOURCEID, ukvar_blp.best_period): sd = OrionStarData(variables_photometry, sid, ukvar_spread.UKvar_ID[ukvar_spread.SOURCEID==sid][0]) fig = lc_and_phase_and_colors(sd, period) fig.ax_j_phase.set_title("ONCvar {}".format(sd.name)) fig.ax_j_lc.set_title(sid) fig.savefig("{0}/{1}/{2}.png".format(path, "50_100", sd.name)) plt.close() blp2 = best_long_period(long_periodic_selector(maxvars_spread_per, min_period=20, max_period=50), min_period=20, max_period=50) ukvar_blp2 = blp2.where(np.in1d(blp2.SOURCEID, ukvar_spread.SOURCEID)) print "There are {0} stars with nominal periods between 20 and 50.".format(len(ukvar_blp2)) for sid, period in zip(ukvar_blp2.SOURCEID, ukvar_blp2.best_period): sd = OrionStarData(variables_photometry, sid, ukvar_spread.UKvar_ID[ukvar_spread.SOURCEID==sid][0]) fig = lc_and_phase_and_colors(sd, period) fig.ax_j_phase.set_title("ONCvar {}".format(sd.name)) fig.ax_j_lc.set_title(sid) fig.savefig("{0}/{1}/{2}.png".format(path, "20_50", sd.name)) plt.close() blp3 = best_long_period(long_periodic_selector(maxvars_spread_per, min_period=10, max_period=20), min_period=10, max_period=20) ukvar_blp3 = blp3.where(np.in1d(blp3.SOURCEID, ukvar_spread.SOURCEID)) print "There are {0} stars with nominal periods between 10 and 20.".format(len(ukvar_blp3)) for sid, period in zip(ukvar_blp3.SOURCEID, ukvar_blp3.best_period): sd = OrionStarData(variables_photometry, sid, ukvar_spread.UKvar_ID[ukvar_spread.SOURCEID==sid][0]) fig = lc_and_phase_and_colors(sd, period) fig.ax_j_phase.set_title("ONCvar {}".format(sd.name)) fig.ax_j_lc.set_title(sid) fig.savefig("{0}/{1}/{2}.png".format(path, "10_20", sd.name)) plt.close()
def find_min_max_day_periodics_and_write(min_period=2, max_period=10, path=path2): blp = best_long_period(long_periodic_selector(maxvars_spread_per, min_period=min_period, max_period=max_period), min_period=min_period, max_period=max_period) ukvar_blp = blp.where(np.in1d(blp.SOURCEID, ukvar_spread.SOURCEID)) print "There are {0} stars with nominal periods between {1} and {2}.".format(len(ukvar_blp), min_period, max_period) for sid, period in zip(ukvar_blp.SOURCEID, ukvar_blp.best_period): sd = OrionStarData(variables_photometry, sid, ukvar_spread.UKvar_ID[ukvar_spread.SOURCEID==sid][0]) fig = lc_and_phase_and_colors(sd, period) fig.ax_j_phase.set_title("ONCvar {} ({})".format(sd.name, megeath_class_column[ukvar_spread.SOURCEID==sid][0])) fig.ax_j_lc.set_title("{0}, Q={1}".format(sid, int((ukvar_spread.autovar+ukvar_spread.strict)[[ukvar_spread.SOURCEID==sid]][0]))) fig.savefig("{0}/{1}_{2}/{3}.png".format(path, min_period, max_period, sd.name)) plt.close()
44508746124658: 15.198723 / 24, 44508746127865: 2.045, 44508746128156: 4.035187, 44508746093223: 4.155218 / 24, # pulsating stars 44508746113922: 13.702932 / 24, 44508746121910: 12.926229 / 24, 44508746105255: 0.213672930845, } # assert that all the overrides are in maxvars! for key in period_override.keys(): assert key in maxvars.SOURCEID periodics = ps.best_long_period(ps.long_periodic_selector( maxvars_spread_per, min_period=2, max_period=75, period_override=period_override), min_period=2, max_period=75, period_override=period_override) periodics.sort('RA') maxvars_periodics = maxvars.where(np.in1d(maxvars.SOURCEID, periodics.SOURCEID)) autovars_true_periodics = autovars_true.where( np.in1d(autovars_true.SOURCEID, periodics.SOURCEID)) autovars_strict_periodics = autovars_strict.where( np.in1d(autovars_strict.SOURCEID, periodics.SOURCEID))
44508746120939 : 3.154922, 44508746124658 : 15.198723 / 24, 44508746127865 : 2.045, 44508746128156 : 4.035187, 44508746093223 : 4.155218 / 24, # pulsating stars 44508746113922 : 13.702932 / 24, 44508746121910 : 12.926229 / 24, 44508746105255 : 0.213672930845, } # assert that all the overrides are in maxvars! for key in period_override.keys(): assert key in maxvars.SOURCEID periodics = ps.best_long_period( ps.long_periodic_selector(maxvars_spread_per, min_period=2, max_period=75, period_override=period_override), min_period=2, max_period=75, period_override=period_override) periodics.sort('RA') maxvars_periodics = maxvars.where( np.in1d(maxvars.SOURCEID, periodics.SOURCEID) ) autovars_true_periodics = autovars_true.where( np.in1d(autovars_true.SOURCEID, periodics.SOURCEID) ) autovars_strict_periodics = autovars_strict.where( np.in1d(autovars_strict.SOURCEID, periodics.SOURCEID) ) c_print( "" ) c_print( "Number of possible variables with detected periods: %d" % len(maxvars_periodics) )