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

cleaner handling of HOD errors

Browse source
This commit is contained in:
Paul Matthew Sutter 2014-01-14 11:36:19 +01:00
parent 89082b2f1e
commit 65cfa8a073
1 changed files with 17 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

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

@ -812,6 +812,8 @@ if (args.hod or args.all) and haloFileBase != "":
print " ", thisHod['name'] print " ", thisHod['name']
sys.stdout.flush() sys.stdout.flush()
sampleName = getSampleName(prefix+"hod_"+thisHod['name'], redshift, False)
parFileName = "./hod.par" parFileName = "./hod.par"
parFile = open(parFileName, 'w') parFile = open(parFileName, 'w')
parFile.write(parFileText.format(omegaM=omegaM, parFile.write(parFileText.format(omegaM=omegaM,
@ -831,19 +833,29 @@ if (args.hod or args.all) and haloFileBase != "":
sampleName=sampleName)) sampleName=sampleName))
parFile.close() parFile.close()
sampleName = getSampleName(prefix+"hod_"+thisHod['name'], redshift, False)
tempFile = "./hod.out_"+sampleName tempFile = "./hod.out_"+sampleName
output = open(tempFile, 'w') output = open(tempFile, 'w')
HOD_PATH = "/home/psutter2/projects/Voids/vide/c_tools/hod/hod" HOD_PATH = "@CMAKE_BINARY_DIR@/c_tools/hod/hod"
#os.system(HOD_PATH+" "+parFileName+">& " + tempFile) #os.system(HOD_PATH+" "+parFileName+">& " + tempFile)
subprocess.call(HOD_PATH+" "+parFileName, stdout=output, stderr=output, shell=True) subprocess.call(HOD_PATH+" "+parFileName, stdout=output, stderr=output,
shell=True)
output.close() output.close()
hodWorked = False
for line in open(tempFile): for line in open(tempFile):
if "MLO" in line: if "MLO" in line:
print " (minimum halo mass = ", line.split()[1], ")" print " (minimum halo mass = ", line.split()[1], ")"
hodWorked = True
break break
os.unlink(tempFile)
if hodWorked:
os.unlink(tempFile)
else:
print "HOD Failed! Log follows:"
for line in open(tempFile):
print line
exit(-1)
outFileName = catalogDir+"/"+sampleName+".dat" outFileName = catalogDir+"/"+sampleName+".dat"
os.system("mv %s/hod_%s.mock %s" % (catalogDir, sampleName, outFileName)) os.system("mv %s/hod_%s.mock %s" % (catalogDir, sampleName, outFileName))