diff --git a/python_tools/void_python_tools/voidUtil/profileUtil.py b/python_tools/void_python_tools/voidUtil/profileUtil.py index 04d4705..73e61fe 100644 --- a/python_tools/void_python_tools/voidUtil/profileUtil.py +++ b/python_tools/void_python_tools/voidUtil/profileUtil.py @@ -67,7 +67,7 @@ def buildProfile(catalog, rMin, rMax): shiftedPart = shiftPart(localPart, center, periodicLine, catalog.ranges) dist = np.sqrt(np.sum(shiftedPart[:,:]**2, axis=1)) - thisProfile, radii = np.histogram(dist, bins=nBins, range=(0,rMaxProfile)) + thisProfile, radii = np.histogram(dist, bins=10, range=(0,rMaxProfile)) deltaV = 4*np.pi/3*(radii[1:]**3-radii[0:(radii.size-1)]**3) thisProfile = np.float32(thisProfile) thisProfile /= deltaV @@ -75,6 +75,7 @@ def buildProfile(catalog, rMin, rMax): allProfiles.append(thisProfile) binCenters = 0.5*(radii[1:] + radii[:-1]) + nVoid = len(voidsToStack) stackedProfile = np.std(allProfiles, axis=0) / np.sqrt(nVoids) sigmas = np.std(allProfiles, axis=0) / np.sqrt(nVoids)