def sgr(): """ Top-down plot of Sgr particles, with selected stars and then re-observed """ cached_data_file = os.path.join(plot_path, 'sgr_kde.pickle') fig,axes = plt.subplots(1, 2, figsize=(14,6.5), sharex=True, sharey=True) # read in all particles as a table pdata = particle_table(N=0, expr="(Pcol<7) & (abs(Lmflag)==1)") extent = {'x':(-90,55), 'y':(-58,68)} if not os.path.exists(cached_data_file): Z = sgr_kde(pdata, extent=extent) with open(cached_data_file, 'w') as f: pickle.dump(Z, f) with open(cached_data_file, 'r') as f: Z = pickle.load(f) axes[1].imshow(Z**0.5, interpolation="nearest", extent=extent['x']+extent['y'], cmap=cm.Blues, aspect=1) np.random.seed(42) particles = particles_today(N=100, expr="(Pcol<7) & (Pcol>0) & (abs(Lmflag)==1)") # read in Stripe 82 RR Lyrae s82 = ascii.read(stripe82_catalog) stripe82_xyz = ra_dec_dist_to_xyz(s82['ra']*u.deg, s82['dec']*u.deg, s82['dist']*u.kpc) catalina = ascii.read(os.path.join(project_root, "data/spitzer_sgr_sample.txt")) catalina = catalina[catalina['dist'] < 35.] catalina_xyz = ra_dec_dist_to_xyz(catalina['ra']*u.deg, catalina['dec']*u.deg, catalina['dist']*u.kpc) ''' quest = ascii.read(quest_catalog) quest = quest[quest['dist'] < 35] quest = quest[quest['dist'] > 10] quest_xyz = ra_dec_dist_to_xyz(quest['ra']*u.deg, quest['dec']*u.deg, quest['dist']*u.kpc) ''' rcparams = {'lines.linestyle' : 'none', 'lines.color' : 'k', 'lines.marker' : 'o'} with rc_context(rc=rcparams): axes[0].plot(pdata['x'], pdata['z'], marker='.', alpha=0.1, ms=5) axes[1].plot(stripe82_xyz[:,0], stripe82_xyz[:,2], marker='.', color='#111111', alpha=0.85, markersize=12, label="Stripe 82", markeredgewidth=0) axes[1].plot(catalina_xyz[:,0], catalina_xyz[:,2], marker='^', color='#111111', alpha=0.85, markersize=7, label="Catalina", markeredgewidth=0) #axes[1].plot(quest_xyz[:,0], quest_xyz[:,2], marker='s', # color='#111111', alpha=0.85, markersize=6, # label="QUEST", markeredgewidth=0) # add solar symbol axes[0].text(-8., 0., s=r"$\odot$") axes[1].text(-8., 0., s=r"$\odot$") axes[0].set_xlabel("$X_{GC}$ [kpc]") axes[1].set_xlabel("$X_{GC}$ [kpc]") axes[0].set_ylabel("$Z_{GC}$ [kpc]") axes[1].legend(loc='upper left') plt.tight_layout() axes[0].set_xlim(extent['x']) axes[0].set_ylim(extent['y']) # turn off right, left ticks respectively axes[0].yaxis.tick_left() axes[1].yaxis.tick_right() # turn off top ticks axes[0].xaxis.tick_bottom() axes[1].xaxis.tick_bottom() fig.subplots_adjust(wspace=0.) fig.savefig(os.path.join(plot_path, "lm10.png"))