def rwa_guys(): """ Creates JJH plots for RWA sources. """ print "Now for RWA sources:" for s, n in zip(rwa_sources, rwa_names): fig = plot3.jjh(rwa_data, s, name=str(n), color_slope=True, date_offset=54579, outfile=path2+'RWA_sources/'+str(n), png_too=True) if fig == None: print "dude %s failed to plot right" % str(s)
def wise_guys(): """ Creates JJH plots for WISE and Aspin sources. """ # wise disks for s, n in zip(wise_disks[rwd], wise_disks_names[rwd]): fig = plot3.jjh(watso, s, name=str(n), color_slope=True, date_offset=54579, outfile=path2+'Wise/'+str(n), png_too=True) if fig == None: print "dude %s failed to plot right" % str(s) # wise special case for s, n in zip(wise_v, wise_v_names): fig = plot3.jjh(data, s, name=str(n), color_slope=True, date_offset=54579, outfile=path2+'Wise/'+str(n), png_too=True) # aspins for s, n in zip(aspin_sources[ras], aspin_sources_names[ras]): fig = plot3.jjh(watso, s, name=str(n), color_slope=True, date_offset=54579, outfile=path2+'Aspin/'+str(n), png_too=True) if fig == None: print "dude %s failed to plot right" % str(s) # wise trans for s, n in zip(wise_trans[rwt], wise_trans_names[rwt]): fig = plot3.jjh(watso, s, name=str(n), color_slope=True, date_offset=54579, outfile=path2+'transition/'+str(n), png_too=True) if fig == None: print "dude %s failed to plot right" % str(s) # wise extras for s, n in zip(wise_extras[rwe], wise_extras_names[rwe]): fig = plot3.jjh(data, s, name=str(n), color_slope=True, date_offset=54579, outfile=path2+'Wise_extras/'+str(n), png_too=True) if fig == None: fig2 = plot3.jjh(welo, s, name=str(n), color_slope=True, date_offset=54579, outfile=path2+'Wise_extras/'+str(n), png_too=True) if fig2 == None: print "dude %s failed to plot right" % str(s) print "did that."
def gen_ukvar_all_jjh(start=0, stop=len(ukvar)): """ Creates all the UKvar JJH plots in one shot. Uses the UKvar ID as the primary identifier; lists info in the filename distinguishing auto/strict/subj, plus per/nonper I'm thinking: 'a' for autovar, 't' for strict, 'j' for subj, 'p' for periodic, 'n' for nonperiodic. Each dude gets two letters. The glued plots only consist of (a) colored-by-time and (b) colored-by-phase. No periodograms here, go see the lightcurves for that. if periodic: do the gluedvars stuff else: just a normal lightcurve! """ for s, id, i in zip(ukvar.SOURCEID, ukvar.UKvar_ID, range(len(ukvar)))[start:stop]: # Each plot gets a suffix: ('a'|'t'|'j')+('p'|'n') suffix = suffix_generator(ukvar, i) # Periodics first # if ukvar.periodic[i] == 1: if False: # uncomment the above line when jjh_phase is functional. # Dig up the best period for this dude! 3 main cases. if ukvar.autovar[i] == 1: if s in autovars_true_periods.SOURCEID: t = autovars_true_periods elif s in autovars_true_periods_s1.SOURCEID: t = autovars_true_periods_s1 else: raise Exception("Something weird happened!") else: t = conf_subj_periodics best_period = t.best_period[t.SOURCEID == s][0] # Let's make 2 plots. LC and folded. Save em into a place. # print out the names as ID_fs_lc.png plot3.jjh(data, s, color_slope=True, date_offset=54034, # timecolor=True, name = "%s: UKvar %s (%s)" % (str(s), str(id), suffix), outfile=jjh_path+"%s_%s_lc.png" % (str(id), suffix)) # ID_fs_phase.png plot3.jjh_phase(data, s, timecolor='time', color_slope=True, date_offset=54034, period=best_period, name = "%s: UKvar %s (%s)" % (str(s), str(id), suffix), outfile=jjh_path+"%s_%s_phase.png" % (str(id), suffix)) # now glue em together! try: call(["montage","-mode", "concatenate", "-tile", "2x", jjh_path+"%s_%s*.png" % (str(id), suffix), jjh_path_g+"%s_%s.png" % (str(id), suffix) ]) except Exception, e: print "Why did montage fail?" raise e else: # Just make the lightcurve. plot3.jjh(data, s, color_slope=True, date_offset=54034, # timecolor=True, name = "%s: UKvar %s (%s)" % (str(s), str(id), suffix), outfile=jjh_path_g+"%s_%s.png" % (str(id), suffix)) print "Completed jjh: UKvar %s" % str(id)