From 5ab0a34129263286ccc1d631db945bbeb914944c Mon Sep 17 00:00:00 2001 From: "P.M. Sutter" Date: Mon, 10 Feb 2014 03:23:41 -0600 Subject: [PATCH] estimate halo density taken from highest redshiftfor universal minVoidRadius --- python_tools/pipeline_source/prepareCatalogs.in.py | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/python_tools/pipeline_source/prepareCatalogs.in.py b/python_tools/pipeline_source/prepareCatalogs.in.py index 1a805cd..c8dd451 100644 --- a/python_tools/pipeline_source/prepareCatalogs.in.py +++ b/python_tools/pipeline_source/prepareCatalogs.in.py @@ -604,11 +604,12 @@ if (args.script or args.all) and haloFileBase != "": sys.stdout.flush() # estimate number of halos to get density + # take from *highest* redshift to get lowest density if haloFileDummy == '': - dataFile = catalogDir+haloFileBase+fileNums[0] + dataFile = catalogDir+haloFileBase+fileNums[-1] else: dataFile = catalogDir+haloFileBase.replace(haloFileDummy, - fileNums[0]) + fileNums[-1]) numPart = 0 if dataFormat == "sdf": SDFcvt_PATH = "@CMAKE_BINARY_DIR@/external/libsdf/apps/SDFcvt/SDFcvt.x86_64" @@ -865,7 +866,7 @@ if (args.halos or args.all) and haloFileBase != "": exit(-1) numKept = 0 - inFile = open(dataFile, 'r') + inFile = open(tempFile, 'r') outFile = open(outFileName, 'a') for (iLine,line) in enumerate(inFile): if np.random.uniform() > keepFraction: continue