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

added ability to do many HOD runs easily

Browse source
This commit is contained in:
P.M. Sutter 2013-02-11 14:17:39 -06:00
parent 85e05aa696
commit ee87db1517
1 changed files with 35 additions and 85 deletions
Download patch file Download diff file Expand all files Collapse all files

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

@ -419,7 +419,7 @@ if (args.halos or args.all) and haloFileBase != "":
outFile.close() outFile.close()
# ----------------------------------------------------------------------------- # -----------------------------------------------------------------------------
# now the SDSS HOD # now the HOD
parFileText = """ parFileText = """
% cosmology % cosmology
OMEGA_M {omegaM} OMEGA_M {omegaM}
@ -461,103 +461,53 @@ BOX_SIZE {boxSize}
root_filename {workDir}/hod root_filename {workDir}/hod
""" """
## TEMP: TURN OFF DR7
if 1 == 0:
#if (args.script or args.all) and haloFileBase != "":
print " Doing DR7 HOD scripts"
setName = prefix+"hod_dr72dim2"
writeScript(setName, prefix+"hod_dr72dim2_z", "multidark",
scriptDir, catalogDir, fileNums, redshifts,
numSubvolumes, numSlices, False, lbox, 5, omegaM)
writeScript(setName, prefix+"hod_dr72dim2_z", "multidark",
scriptDir, catalogDir, fileNums, redshifts,
numSubvolumes, numSlices, True, lbox, 5, omegaM)
if 1 == 0:
#if (args.hod or args.all) and haloFileBase != "":
print " Doing DR7 HOD"
for (iRedshift, redshift) in enumerate(redshifts):
print " z = ", redshift
parFileName = "./hod.par"
parFile = open(parFileName, 'w')
if haloFileDummy == '':
haloFile = catalogDir+haloFileBase+fileNums[iRedshift]
else:
haloFile = catalogDir+haloFileBase.replace(haloFileDummy,
fileNums[iRedshift])
parFile.write(parFileText.format(omegaM=omegaM,
hubble=hubble,
redshift=redshift,
Mmin=1.99526e12,
M1=3.80189e13,
sigma_logM=0.21,
alpha=1.12,
Mcut=6.91831e11,
galden=0.02,
haloFile=haloFile,
haloFileFormat=dataFormat,
numPartPerSide=numPart**(1/3.),
boxSize=lbox,
workDir=catalogDir))
parFile.close()
os.system(hodPath+" "+parFileName+">& /dev/null")
sampleName = getSampleName(prefix+"hod_dr72dim2", redshift, False)
outFileName = catalogDir+"/"+sampleName+".dat"
os.system("mv %s/hod.mock %s" % (catalogDir, outFileName))
os.system("rm %s/hod.*" % catalogDir)
# -----------------------------------------------------------------------------
# now the BOSS HOD
if (args.script or args.all) and haloFileBase != "": if (args.script or args.all) and haloFileBase != "":
print " Doing DR9 HOD scripts" print " Doing HOD scripts"
setName = prefix+"hod_dr9mid" for thisHod in hodParmList:
writeScript(setName, prefix+"hod_dr9mid_z", "multidark", print " ", thisHod['name']
scriptDir, catalogDir, fileNums, redshifts, setName = prefix+"hod_"+thisHod['name']
numSubvolumes, numSlices, False, lbox, 15, omegaM) writeScript(setName, prefix+"hod_"+thisHod['name']+"_z", "multidark",
writeScript(setName, prefix+"hod_dr9mid_z", "multidark", scriptDir, catalogDir, fileNums, redshifts,
scriptDir, catalogDir, fileNums, redshifts, numSubvolumes, numSlices, False, lbox, 15, omegaM)
numSubvolumes, numSlices, True, lbox, 15, omegaM) writeScript(setName, prefix+"hod_"+thisHod['name']+"_z", "multidark",
scriptDir, catalogDir, fileNums, redshifts,
numSubvolumes, numSlices, True, lbox, 15, omegaM)
if (args.hod or args.all) and haloFileBase != "": if (args.hod or args.all) and haloFileBase != "":
print " Doing DR9 HOD" print " Doing HOD"
for (iRedshift, redshift) in enumerate(redshifts): for thisHod in hodParmList:
print " z = ", redshift print " ", thisHod['name']
for (iRedshift, redshift) in enumerate(redshifts):
print " z = ", redshift
# these parameters come from Manera et al 2012, eq. 26 parFileName = "./hod.par"
# modified to match the observed galaxy density parFile = open(parFileName, 'w')
parFileName = "./hod.par" if haloFileDummy == '':
parFile = open(parFileName, 'w') haloFile = catalogDir+haloFileBase+fileNums[iRedshift]
if haloFileDummy == '': else:
haloFile = catalogDir+haloFileBase+fileNums[iRedshift] haloFile = catalogDir+haloFileBase.replace(haloFileDummy,
else:
haloFile = catalogDir+haloFileBase.replace(haloFileDummy,
fileNums[iRedshift]) fileNums[iRedshift])
parFile.write(parFileText.format(omegaM=omegaM, parFile.write(parFileText.format(omegaM=omegaM,
hubble=hubble, hubble=hubble,
redshift=redshift, redshift=redshift,
Mmin=0.0, Mmin=thisHod['Mmin'],
#Mmin=1.23e13, #Mmin=1.23e13,
M1=1.e14, M1=thisHod['M1'],
sigma_logM=0.596, sigma_logM=thisHod['sigma_logM'],
alpha=1.0127, alpha=thisHod['alpha'],
Mcut=1.19399e13, Mcut=thisHod['Mcut'],
galden=galDens, galden=thisHod['galDens'],
#galden=0.000225, #galden=0.000225,
haloFile=haloFile, haloFile=haloFile,
haloFileFormat=dataFormat, haloFileFormat=dataFormat,
numPartPerSide=numPart**(1/3.), numPartPerSide=numPart**(1/3.),
boxSize=lbox, boxSize=lbox,
workDir=catalogDir)) 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_"+thisHod['name'], redshift, False)
outFileName = catalogDir+"/"+sampleName+".dat" outFileName = catalogDir+"/"+sampleName+".dat"
os.system("mv %s/hod.mock %s" % (catalogDir, outFileName)) os.system("mv %s/hod.mock %s" % (catalogDir, outFileName))
os.system("rm %s/hod.*" % catalogDir)
os.system("rm %s/hod.*" % catalogDir)