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
