'*', '*', '*', '*', '*', '*', '*', '*', '*', '*', '*', '*' ] data_J_b = [ '*', '*', '*', 4, '*', '*', '*', '*', 9, 10, '*', '*', '*', '*', 15, '*', '*', '*', 19, '*', '*', 22, '*', 24, 25, 26, '*' ] #import 850um data ID = [] RA = [] DEC = [] flux_850 = [] for line in cat_850.readlines()[2:]: tmp = line.split() ID.append(tmp[0]) RA.append(co.convHMS(tmp[1])) DEC.append(co.convDMS(tmp[2])) flux_850.append(float(tmp[3])) #import MIPS IDs ID_MIPS = [] RA_MIPS = [] DEC_MIPS = [] for line in cat_MIPS.readlines()[2:]: tmp = line.split() if tmp[1][0:2] == '17': ID_MIPS.append(tmp[0]) RA_MIPS.append(co.convHMS(tmp[1])) DEC_MIPS.append(co.convDMS(tmp[2])) else: pass
sub2_1700 = data_1700[2] #but don't use these or you'll double count Rband_1700 = data_1700[ 3] #undetected in 24um have estimates of their 24um flux from the R-band mag all_1700 = data_1700[4] #everything detected at 24um width_1700 = np.logspace(0.05 + 0.15, 3.5 + 0.15, len( data_1700[0])) - bins_1700 #850um data ID_850 = [] RA_850 = [] DEC_850 = [] cat_850 = open('../h1549/4sig.gaia', 'r') for line in cat_850.readlines()[2:]: tmp = line.split() ID_850.append(tmp[0]) RA_850.append(co.convHMS(tmp[1])) DEC_850.append(co.convDMS(tmp[2])) ID_850 = np.array(ID_850) RA_850 = np.array(RA_850) DEC_850 = np.array(DEC_850) cat_850.close() #proto-z LBG data for HS1549 ID = [] RA = [] DEC = [] RA_1549 = co.convHMS('15:51:52.604') DEC_1549 = co.convDMS('19:11:22.90') flux = [] z = [] cat = open('../h1549/MIPS/h1549_24um_proto_cat.gaia', 'r')
S) - 25452.43730 * erf(.6324555320 * np.sqrt(S)) + 25870 # From Geach et al. 2017 S2-CLS paper def boost(snr): return 1.0 + 0.2 * (snr / 5.0)**(-2.3) def S2CLS_core(S, a, b, c): return a * (S / b)**(c) * np.exp(-S / b) ############# # HS1549 DATA ############# RA_map_1549 = co.convHMS('15:51:53.582') DEC_map_1549 = co.convDMS('19:10:50.30') data_1549 = open('../h1549/flux_4sig.gaia', 'r') ID_1549 = [] RA_1549 = [] DEC_1549 = [] flux_1549 = [] flux_1549_deboosted = [] snr_1549 = [] for line in data_1549.readlines()[2:]: tmp = line.split() if tmp[0] == '1': #SMA resolved source ID_1549.append('1.1') ID_1549.append('1.2') ID_1549.append('1.3')
#sources without multi-lambda data data_24_b = ['*','*','*','*',5,6,7,'*','*','*','*','*','*',14,'*','*','*',18,19,20,'*','*','*','*','*','*','*','*','*']#29??? data_ch4_b = ['*','*','*','*',5,6,7,'*','*','*','*','*','*',14,'*','*','*',18,19,20,'*','*','*','*','*','*','*','*',29] data_ch2_b = ['*','*','*','*',5,6,7,'*','*','*','*','*','*',14,'*','*','*',18,19,20,'*','*','*','*','*','*','*','*',29] data_J_b = ['*','*','*','*',5,6,7,8,'*','*','*',12,'*',14,'*',16,17,18,19,20,'*','*',23,'*',25,'*','*',28,29] #import 850um data ID = [] RA = [] DEC = [] flux_850 = [] for line in cat_850.readlines()[2:]: tmp = line.split() ID.append(tmp[0]) RA.append(co.convHMS(tmp[1])) DEC.append(co.convDMS(tmp[2])) flux_850.append(float(tmp[3])) #import IRAC IDs ID_IRAC = [] RA_IRAC = [] DEC_IRAC = [] for line in cat_IRAC.readlines()[2:]: tmp = line.split() if tmp[1][0:2] == '15': ID_IRAC.append(tmp[0]) RA_IRAC.append(co.convHMS(tmp[1])) DEC_IRAC.append(co.convDMS(tmp[2])) else: pass
import xlrd import matplotlib.patches as patches from PIL import Image # Define a Gaussian function def Gaussian(x, A, mu, sigma): return A*np.exp(-0.5*((x-mu)/sigma)**2) # Physical properties from Ned Wright cosmo calc. DL = 101.6*1e6 #pc scale = 0.471 #kpc/" # Coma Cluster properties from Zwicky+1942 z_coma = 0.023100 v_coma = 6925.0 #km/s sigma_coma = 1000.0 #km/s RA_coma=co.convHMS('12:59:48.7') #deg DEC_coma=co.convDMS('+27:58:50') #deg ################ # Importing Data ################ loc = ('coma.xlsx') wb=xlrd.open_workbook(loc) sheet=wb.sheet_by_index(0) # Sorting data and converting to useful units ID = [] RA = [] DEC = [] type = []
import coords as co mom0_map = '../moment_maps/spt0348_C+_dirty_contsub_briggs_robust05.image.mom0.fits' mom1_map = '../moment_maps/spt0348_C+_dirty_contsub_briggs_robust05.image.mom1.fits' mom2_map = '../moment_maps/spt0348_C+_dirty_contsub_briggs_robust05.image.mom2.fits' ############ #MOMENT MAPS ############ #set up figure fig = pl.figure() pl.rc('font', size=30) pl.rc('mathtext', default='regular') RA_c = co.convHMS('03:48:42.271') DEC_c = co.convDMS('-62:20:50.85') #moment-0 map mom0 = ap.FITSFigure(mom0_map, figure=fig, subplot=[0, 0, 0.75, 1]) mom0.tick_labels.set_xformat('hh:mm:ss') mom0.tick_labels.set_yformat('dd:mm:ss') mom0.hide_xaxis_label() mom0.hide_xtick_labels() mom0.hide_yaxis_label() mom0.hide_ytick_labels() mom0.recenter(RA_c, DEC_c, width=5. / 3600, height=5. / 3600) mom0.show_colorscale(cmap='jet', vmax=12, vmin=0.0, interpolation='nearest') mom0.add_colorbar(axis_label_text=r'S$_{[CII]}$ (Jy beam$^{-1}$ km s$^{-1}$)', axis_label_rotation=-90, axis_label_pad=50)