limits_40 = n.log10([Npmin* 9.6 * 10**10. , 5e15]) MPART = n.array([9.63 * 10**7, 1.51 * 10**9, 2.359 * 10**10, 9.6 * 10**10]) zmin = -0.01 zmax = 2.3 #======================= #======================= cos = 'cen' #======================= #======================= # redshift selection zSel = lib.zSelection( data, zmin, zmax ) # mass selection mSel = lib.mSelection(data, qty, limits_04, limits_10, limits_25,limits_40) # minimum number counts selection nSelCen = lib.nSelection(data, NminCount, cos ) # altogether ok = (zSel) & (mSel) & (nSelCen) # selection per box : MD04=(data["boxName"]=='MD_0.4Gpc') MD10=(data["boxName"]=='MD_1Gpc_new_rockS') MD25=(data["boxName"]=='MD_2.5Gpc') MD40=(data["boxName"]=='MD_4Gpc') MD25NW=(data["boxName"]=='MD_2.5GpcNW') MD40NW=(data["boxName"]=='MD_4GpcNW') from hmf import MassFunction from astropy.cosmology import FlatLambdaCDM
cos = 'cen' fw.write("c o s=" + cos + '\n') #======================= #======================= # y coordinates ff = mvir * data["dNdlnM_" + cos] / data["rhom"] / abs(data["dlnsigmaMdlnM"]) ff_c = mvir * data["dNdlnM_" + cos + "_c"] / data["rhom"] / abs( data["dlnsigmaMdlnM"]) log_MF = n.log10(ff) log_MF_c = n.log10(ff_c) # redshift selection zSel = lib.zSelection(data, zmin, zmax) # mass selection mSel = lib.mSelection(data, qty, logNpmin) mSel2_inter = (data["log_mvir"] < 13.2) & (data["redshift"] > 0.) mSel2 = (mSel2_inter == False) # minimum number counts selection nSelCen = lib.nSelection(data, NminCount, cos) # altogether ok1 = (zSel) & (mSel) & (mSel2) & (nSelCen) # selection per box : MD04 = (data["boxName"] == 'MD_0.4Gpc') MD10 = (data["boxName"] == 'MD_1Gpc') MD25 = (data["boxName"] == 'MD_2.5Gpc') MD40 = (data["boxName"] == 'MD_4Gpc') MD25NW = (data["boxName"] == 'MD_2.5GpcNW') MD40NW = (data["boxName"] == 'MD_4GpcNW') DS80 = (data["boxName"] == 'DS_8Gpc')
#======================= #======================= cos = 'cen' #======================= #======================= # y coordinates ff = mvir * data["dNdlnM_"+cos] / data["rhom"] / abs(data["dlnsigmaMdlnM"]) ff_c = mvir * data["dNdlnM_"+cos+"_c"] / data["rhom"] / abs(data["dlnsigmaMdlnM"]) log_MF = n.log10( ff ) log_MF_c = n.log10( ff_c ) # redshift selection zSel = lib.zSelection( data, zmin, zmax ) # mass selection mSel = lib.mSelection(data, qty, logNpmin) mSel2_inter = (data["log_mvir"]<13.2) & (data["redshift"]>0.) mSel2 = (mSel2_inter==False) # minimum number counts selection nSelCen = lib.nSelection(data, NminCount, cos ) # altogether ok = (zSel) & (mSel) & (mSel2) & (nSelCen) # selection per box : MD04=(data["boxName"]=='MD_0.4Gpc') MD10=(data["boxName"]=='MD_1Gpc') MD25=(data["boxName"]=='MD_2.5Gpc') MD40=(data["boxName"]=='MD_4Gpc') MD25NW=(data["boxName"]=='MD_2.5GpcNW') MD40NW=(data["boxName"]=='MD_4GpcNW') x_data = logsig[ok]
limits_40 = n.log10([Npmin * 9.6 * 10**10., 5e15]) MPART = n.array([9.63 * 10**7, 1.51 * 10**9, 2.359 * 10**10, 9.6 * 10**10]) zmin = -0.01 zmax = 2.3 #======================= #======================= cos = 'cen' #======================= #======================= # redshift selection zSel = lib.zSelection(data, zmin, zmax) # mass selection mSel = lib.mSelection(data, qty, limits_04, limits_10, limits_25, limits_40) # minimum number counts selection nSelCen = lib.nSelection(data, NminCount, cos) # altogether ok = (zSel) & (mSel) & (nSelCen) # selection per box : MD04 = (data["boxName"] == 'MD_0.4Gpc') MD10 = (data["boxName"] == 'MD_1Gpc_new_rockS') MD25 = (data["boxName"] == 'MD_2.5Gpc') MD40 = (data["boxName"] == 'MD_4Gpc') MD25NW = (data["boxName"] == 'MD_2.5GpcNW') MD40NW = (data["boxName"] == 'MD_4GpcNW') from hmf import MassFunction from astropy.cosmology import FlatLambdaCDM import astropy.units as u