guilhem_lavaux/vide_public
0
0
Fork 0
mirror of https://bitbucket.org/cosmicvoids/vide_public.git synced 2025-07-04 15:21:11 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

user allowed to input number of bins in profile

Browse source
This commit is contained in:
P.M. Sutter 2014-07-18 10:14:21 -04:00
parent 1ae7022f6d
commit b21300a56a
1 changed files with 3 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

5
python_tools/void_python_tools/voidUtil/profileUtil.py
View file

@ -37,12 +37,13 @@ def HSWProfile(r, rs, dc):
return dc * (1 - (r/rs)**alpha) / (1+ (r)**beta) + 1 return dc * (1 - (r/rs)**alpha) / (1+ (r)**beta) + 1
# ----------------------------------------------------------------------------- # -----------------------------------------------------------------------------
def buildProfile(catalog, rMin, rMax): def buildProfile(catalog, rMin, rMax, nBins=10):
# builds a stacked radial density profile from the given catalog # builds a stacked radial density profile from the given catalog
# catalog: void catalog # catalog: void catalog
# rMin: minimum void radius, in Mpc/h # rMin: minimum void radius, in Mpc/h
# rMax: maximum void radius, in Mpc/h # rMax: maximum void radius, in Mpc/h
# nBins: number of bins in profile (detaulf 10)
# #
# returns: # returns:
# binCenters: array of radii in binned profile # binCenters: array of radii in binned profile
@ -71,7 +72,7 @@ def buildProfile(catalog, rMin, rMax):
shiftedPart = shiftPart(localPart, center, periodicLine, catalog.ranges) shiftedPart = shiftPart(localPart, center, periodicLine, catalog.ranges)
dist = np.sqrt(np.sum(shiftedPart[:,:]**2, axis=1)) dist = np.sqrt(np.sum(shiftedPart[:,:]**2, axis=1))
thisProfile, radii = np.histogram(dist, bins=10, range=(0,rMaxProfile)) thisProfile, radii = np.histogram(dist, bins=nBins, range=(0,rMaxProfile))
deltaV = 4*np.pi/3*(radii[1:]**3-radii[0:(radii.size-1)]**3) deltaV = 4*np.pi/3*(radii[1:]**3-radii[0:(radii.size-1)]**3)
thisProfile = np.float32(thisProfile) thisProfile = np.float32(thisProfile)
thisProfile /= deltaV thisProfile /= deltaV