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

masking of mock catalogs now works

Browse source
This commit is contained in:
P.M. Sutter 2012-12-02 13:08:20 -06:00
parent 7484bebbfa
commit f3f96b0dcd
1 changed files with 12 additions and 11 deletions
Download patch file Download diff file Expand all files Collapse all files

23
python_tools/pipeline_source/prepareCatalogs.in.py
View file

@ -417,7 +417,7 @@ RESOLUTION {numPartPerSide}
BOX_SIZE {boxSize} BOX_SIZE {boxSize}
% output % output
root_filename hod root_filename {workDir}/hod
""" """
if (args.script or args.all) and haloFileBase != "": if (args.script or args.all) and haloFileBase != "":
@ -449,17 +449,18 @@ if (args.hod or args.all) and haloFileBase != "":
haloFile=haloFile, haloFile=haloFile,
haloFileFormat=dataFormat, haloFileFormat=dataFormat,
numPartPerSide=numPart**(1/3.), numPartPerSide=numPart**(1/3.),
boxSize=lbox)) boxSize=lbox,
workDir=catalogDir))
parFile.close() parFile.close()
os.system(hodPath+" "+parFileName+"") os.system(hodPath+" "+parFileName+">& /dev/null")
#os.system(hodPath+" "+parFileName+">& /dev/null")
sampleName = getSampleName(prefix+"hod_dr72dim2", redshift, False) sampleName = getSampleName(prefix+"hod_dr72dim2", redshift, False)
outFileName = catalogDir+"/"+sampleName+".dat" outFileName = catalogDir+"/"+sampleName+".dat"
os.system("mv hod.mock" + " " + outFileName) os.system("mv %s/hod.mock %s" % (catalogDir, outFileName))
os.system("rm ./hod.*") os.system("rm %s/hod.*" % catalogDir)
os.system("rm %s/hod-*" % catalogDir)
# ----------------------------------------------------------------------------- # -----------------------------------------------------------------------------
# now the BOSS HOD # now the BOSS HOD
@ -493,15 +494,15 @@ if (args.hod or args.all) and haloFileBase != "":
haloFile=haloFile, haloFile=haloFile,
haloFileFormat=dataFormat, haloFileFormat=dataFormat,
numPartPerSide=numPart**(1/3.), numPartPerSide=numPart**(1/3.),
boxSize=lbox)) boxSize=lbox,
workDir=catalogDir))
parFile.close() parFile.close()
os.system(hodPath+" "+parFileName+">& /dev/null") os.system(hodPath+" "+parFileName+">& /dev/null")
sampleName = getSampleName(prefix+"hod_dr9mid", redshift, False) sampleName = getSampleName(prefix+"hod_dr9mid", redshift, False)
outFileName = catalogDir+"/"+sampleName+".dat" outFileName = catalogDir+"/"+sampleName+".dat"
os.system("mv hod.mock" + " " + outFileName) os.system("mv %s/hod.mock %s" % (catalogDir, outFileName))
os.system("rm ./hod.*")
os.system("rm %s/hod.*" % catalogDir)
os.system("rm %s/hod-*" % catalogDir)