Exemplos de mSelection em Python

Exemplo n.º 1

Arquivo: mvir-orders-of-magnitude.py Projeto: JohanComparat/pySU

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

Exemplo n.º 2

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')

Exemplo n.º 3

Arquivo: mvir-3-joint_fit.py Projeto: JohanComparat/pySU

#=======================
#=======================
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]

Exemplo n.º 4

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