def R(fig,surfaceA,surfaceB,H,mask):
#produce surface map of YSO distributions
Fig = aplpy.FITSFigure(fig)
levels = (15,45,75,100,150)
cset_H = Fig.show_contour(H, levels=(200,10000),colors='blue')
#cset1 = Fig.show_contour(surfaceA, levels=levels,colors='green')
#cset2 = Fig.show_contour(surfaceB, levels=levels,colors='k')
cset_mask = Fig.show_contour(mask, levels=(0.5,10),colors='k')
#YSO.w40(Fig)
YSO.c2dGBS(Fig)
OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt')
ra = OBmarkers[:,0]
dec = OBmarkers[:,1]
Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=250, alpha=1)
ra = 277.299583333
dec = -2.06388888889
Fig.show_markers(ra,dec, layer='W40a', edgecolor='k', facecolor='yellow', marker='*',s=150, alpha=1)
Fig.show_grayscale(vmax=0.027,vmin=0.0015, stretch='sqrt',invert='TRUE')
#Fig.recenter(277.8995,-2.09200, height=0.6 ,width=0.6)
Fig.recenter(277.4406666667,-1.90194444444, height=0.85 ,width=0.5)
Fig.add_colorbar()
#Fig.colorbar.set_location('top')
Fig.colorbar.set_axis_label_text('850$\mu m$ (Jy/pixel)')
def A(s450,s850,temp,sig450,sig850,surface):
#tri panel map - should tile all maps at once. 1
big = plt.figure(figsize=(12, 17))
#deffine image in APLPY
Fig_850 = aplpy.FITSFigure(s850,figure=big,subplot=[0.1,0.64,0.8,0.26])
Fig_450 = aplpy.FITSFigure(s450,figure=big,subplot=[0.1,0.37,0.8,0.26])
Fig_temp = aplpy.FITSFigure(temp,figure=big,subplot=[0.1,0.1,0.8,0.26])
#Deffine contours
cset1 = Fig_850.show_contour(s850, levels=(sig850,300*sig850), linewidth=3, colors=('black','black'))
cset2 = Fig_450.show_contour(s450, levels=(sig450,300*sig450), linewidth=3, colors=('black','black'))
#cset3 = Fig_temp.show_contour(temp, levels=(12,20,30), linewidth=3, colors=('black'))
#Figure properties
Fig_450.show_grayscale(vmax=10*sig450,vmin=0.0,invert=True, stretch='sqrt')
Fig_450.add_label(0.07, 0.95, r'b) $450\mu m$',size='large', relative=True)
Fig_850.show_grayscale(vmax=10*sig850,vmin=0.0,invert=True, stretch='sqrt')
Fig_850.add_label(0.07, 0.95, r'a) $850\mu m$',size='large', relative=True)
Fig_temp.show_colorscale(vmax=35.0,vmin=10.0, stretch='sqrt',cmap='gist_heat')
Fig_temp.add_label(0.05, 0.95, r'c) $T_{d}$',size='large', relative=True)
#PLOT ysos
YSO.c2dGBS(Fig_850)
#PLOT contours
cset850 = Fig_temp.show_contour(s850, levels=(sig850,100), linewidth=3, colors=('black'))
csetYSO = Fig_temp.show_contour(surface, levels=(20,60,100), linewidth=3, colors=('red'))
Fig_450.tick_labels.set_font(size='medium')
Fig_450.add_colorbar()
#Fig_450.recenter(277.074583333,-3.74686111111,height=0.25 ,width=0.65)
Fig_850.tick_labels.set_font(size='medium')
Fig_850.add_colorbar()
#Fig_850.recenter(277.074583333,-3.74686111111,height=0.25 ,width=0.65)
Fig_temp.tick_labels.set_font(size='medium')
Fig_temp.add_colorbar()
#Fig_temp.recenter(277.074583333,-3.74686111111,height=0.25 ,width=0.65)
Fig_450.hide_xaxis_label()
Fig_850.hide_xaxis_label()
Fig_450.hide_xtick_labels()
Fig_850.hide_xtick_labels()
Fig_850.hide_yaxis_label()
Fig_temp.hide_yaxis_label()
Fig_450.ticks.show()
Fig_450.ticks.set_color('black')
Fig_850.ticks.show()
Fig_850.ticks.set_color('black')
Fig_temp.ticks.show()
Fig_temp.ticks.set_color('black')
big.canvas.draw()
return
def B(s850,temp,sig850,pos):
Fig_temp = aplpy.FITSFigure(temp)
Fig_temp.show_colorscale(vmax=35.0,vmin=10.0, stretch='linear',cmap='gist_heat')
#PLOT ysos
YSO.c2dGBS(Fig_temp)
data = np.loadtxt('OBstars_list.txt')
RA_OB = data[:,0]
DEC_OB = data[:,1]
Fig_temp.show_markers(RA_OB,DEC_OB, edgecolor='k', facecolor='yellow', marker='*',s=100, alpha=1)
#PLOT contours
cset850 = Fig_temp.show_contour(s850, levels=(sig850,1000), linewidth=3, colors=('b'))
Fig_temp.ticks.show()
Fig_temp.ticks.set_color('black')
Fig_temp.tick_labels.set_font(size='medium')
Fig_temp.add_colorbar()
Fig_temp.colorbar.set_axis_label_text('T$_{d}$ (K)')
big = plt.figure(figsize=(25, 6))
Fig = aplpy.FITSFigure(s850,figure=big,subplot=[0.1,0.1,0.375,0.8])
Fig.show_grayscale(vmax=10*sig850,vmin=(1./5.)*sig850,invert=True, stretch='sqrt')
###############
#PLOT-contours#
###############
#OPEN up the region specific clump list
data = fits.open(clumpcat)
DATA = data[1].data
#PLOT ysos
YSO.c2dGBS(Fig)
#OB star position
RA_OB = 328.3708
DEC_OB = 47.2667
Fig.show_markers(RA_OB,DEC_OB, edgecolor='k', facecolor='yellow', marker='*',s=100, alpha=1)
Fig.show_contour(s850, levels=(sig850,300*sig850), linewidth=3, colors=('black','black'))
###########
#PLOT-misc#
###########
Fig.ticks.show()
Fig.ticks.set_color('black')
