cl1lt10coords = coordinates.SkyCoord(ysos[0][cl1lt10]['_RAJ2000'], ysos[0][cl1lt10]['_DEJ2000'], frame='fk5').galactic cl1gt10coords = coordinates.SkyCoord(ysos[0][cl1gt10]['_RAJ2000'], ysos[0][cl1gt10]['_DEJ2000'], frame='fk5').galactic fig7 = pl.figure(7) fig7.clf() gray = copy.copy(pl.cm.gray_r) gray.set_bad('white') gray.set_under('white') FMM = FITSFigure(sfmasshdu1e4, figure=fig7, cmap=gray, stretch='log', vmax=99.95) FMM.colorbar._colorbar.set_label("Star-Forming ($n>10^4$ cm$^{-3}$) Gas Surface Density\n[$M_{\odot}$ pc$^{-2}$]", rotation=270, labelpad=50) FMM.save(paths.fpath('StarFormingMassMap.pdf')) FMM.show_markers(cl1lt10coords.l, cl1lt10coords.b, marker='x', edgecolor='r', zorder=1100, layer='kangx') FMM.show_markers(cl1gt10coords.l, cl1gt10coords.b, marker='o', edgecolor='r', zorder=1100, layer='kango') FMM.save(paths.fpath('StarFormingMassMap_Kang2009ClassIYSOs.pdf')) FMM.remove_layer('kangx') FMM.remove_layer('kango') FMM.show_contour(sfmasshdu1e4, levels=[100,1000], zorder=1000, smooth=1, colors=['#AA1111','#FF4411','#BB6666'], linewidths=[2,3], alpha=0.7) FMM.show_contour(sfmasshdu1e5, levels=[100,1000], zorder=1000, smooth=1, colors=['#1111BB','#1144FF','#6699BB'], linewidths=[2,3], alpha=0.9) FMM.save(paths.fpath('StarFormingMassMap_thresholdcontours.pdf')) fig8 = pl.figure(8) fig8.clf() ax8 = fig8.gca() ax8.loglog(totalmassmap[~badsfmass1e4 & ok2cm], sfr2cmd[~badsfmass1e4 & ok2cm], 'r.', alpha=0.1) ax8.loglog(sfmassmap1e4[~badsfmass1e4 & ok2cm], sfr2cmd[~badsfmass1e4 & ok2cm], 'k.', alpha=0.1) ax8.loglog([1,1e4],np.array([1,1e4])*1.2e-8*1e6,'b--', alpha=0.5)