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

more options for number function

Browse source
This commit is contained in:
P.M. Sutter 2014-05-26 08:09:52 -04:00
parent 792190ccc5
commit e177bddea2
1 changed files with 19 additions and 12 deletions
Download patch file Download diff file Expand all files Collapse all files

31
python_tools/void_python_tools/voidUtil/plotUtil.py
View file

@ -42,12 +42,16 @@ def fill_between(x, y1, y2=0, ax=None, **kwargs):
# ----------------------------------------------------------------------------- # -----------------------------------------------------------------------------
def plotNumberFunction(catalogList, def plotNumberFunction(catalogList,
figDir="./", figDir="./",
plotName="numberfunc"): plotName="numberfunc",
cumulative=True,
binWidth=1):
# plots a cumulative number function # plots a cumulative number function
# catalogList: list of void catalogs to plot # catalogList: list of void catalogs to plot
# figDir: output directory for figures # figDir: output directory for figures
# plotName: name to prefix to all outputs # plotName: name to prefix to all outputs
# cumulative: if True, plots cumulative number function
# binWidth: width of histogram bins in Mpc/h
# returns: # returns:
# numberFuncList: array of len(catalogList), # numberFuncList: array of len(catalogList),
# each element has array of size bins, number, +/- 1 sigma # each element has array of size bins, number, +/- 1 sigma
@ -56,7 +60,11 @@ def plotNumberFunction(catalogList,
plt.clf() plt.clf()
plt.xlabel("$R_{eff}$ [$h^{-1}Mpc$]", fontsize=14) plt.xlabel("$R_{eff}$ [$h^{-1}Mpc$]", fontsize=14)
plt.ylabel(r"log ($n$ (> R) [$h^3$ Gpc$^{-3}$])", fontsize=14)
if cumulative:
plt.ylabel(r"log ($n$ (> R) [$h^3$ Gpc$^{-3}$])", fontsize=14)
else:
plt.ylabel(r"log ($dn/dR$ [$h^3$ Gpc$^{-3}$])", fontsize=14)
numberFuncList = [] numberFuncList = []
@ -70,25 +78,24 @@ def plotNumberFunction(catalogList,
boxVol = vp.getSurveyProps(maskFile, boxVol = vp.getSurveyProps(maskFile,
sample.zBoundary[0], sample.zBoundary[1], sample.zBoundary[0], sample.zBoundary[1],
sample.zRange[0], sample.zRange[1], "all", sample.zRange[0], sample.zRange[1], "all",
selectionFuncFile=None)[0] selectionFuncFile=sample.selFunFile)[0]
#selectionFuncFile=sample.selFunFile)[0]
boxVol *= obsFudgeFactor
else: else:
boxVol = sample.boxLen*sample.boxLen*(sample.zBoundaryMpc[1] - boxVol = sample.boxLen*sample.boxLen*(sample.zBoundaryMpc[1] -
sample.zBoundaryMpc[0]) sample.zBoundaryMpc[0])
boxVol *= 1.e-9 # Mpc->Gpc boxVol *= 1.e-9 # Mpc->Gpc
bins = 100. bins = 100./binWidth
hist, binEdges = np.histogram(data, bins=bins, range=(0., 100.)) hist, binEdges = np.histogram(data, bins=bins, range=(0., 100.))
binCenters = 0.5*(binEdges[1:] + binEdges[:-1]) binCenters = 0.5*(binEdges[1:] + binEdges[:-1])
foundStart = False if cumulative:
for iBin in xrange(len(hist)): foundStart = False
if not foundStart and hist[iBin] == 0: for iBin in xrange(len(hist)):
continue if not foundStart and hist[iBin] == 0:
foundStart = True continue
hist[iBin] = np.sum(hist[iBin:]) foundStart = True
hist[iBin] = np.sum(hist[iBin:])
nvoids = len(data) nvoids = len(data)
var = hist * (1. - hist/nvoids) var = hist * (1. - hist/nvoids)