from lmfit import minimize, Parameters, report_fit from matplotlib.ticker import MaxNLocator import scipy #Inputs folder = "/Users/chrisfuller/Dropbox/phd/herchel/coma/final_outputs/" # input/output folder fname = 'coma_supercluster_cal12.fits' #input name cat = Table(pj(folder,fname)) cat.add_column('all', cat.g) cat.all = 1 cat.add_column('colour', cat.g - cat.r) #cat = cat.where(cat.bptclass ==1) #convert radius vir to log10 cat.RADIUS_VIR = np.log10(cat.RADIUS_VIR) #select detected galaxies #cat = cat.where(np.nan_to_num(cat.K) != 0.0) detected = cat.where((cat.DMASS_TYPE != 0) & (np.nan_to_num(cat.METAL) > 0.0)) undetected = cat.where((cat.DMASS_TYPE == 0) & (np.nan_to_num(cat.METAL) > 0.0)) cats = [detected, undetected] du = ['FIR-detected', 'FIR-undetected'] lisy = ['-', '--'] mark = ['*', '+'] #switchs N_den = 7 #number of density bins N = 60