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

streamlined merger tree processing

Browse source
This commit is contained in:
P.M. Sutter 2012-12-20 13:38:59 -06:00
parent 3c8d7b2fb2
commit d4febc6942
2 changed files with 12 additions and 12 deletions
Download patch file Download diff file Expand all files Collapse all files

6
crossCompare/analysis/datasetsToAnalyze.py
View file

@ -2,14 +2,14 @@
workDir = "/home/psutter2/workspace/Voids/"
figDir = "./figs"
dataDir = "./data"
dataDir = "/home/psutter2/workspace/Voids/multidark/crossCompare/mergerTree/"
CTOOLS_PATH = "../../c_tools/"
baseSampleDir = "multidark/md_ss0.000175/sample_md_ss0.000175_z0.56_d00/"
sampleDirList = [
"multidark/md_ss0.000175/sample_md_ss0.000175_z0.56_d00/",
"multidark/md_ss0.0004/sample_md_ss0.0004_z0.56_d00/",
]

18
crossCompare/analysis/mergerTree.py
View file

@ -30,30 +30,30 @@ if not os.access(filename, os.F_OK):
parms = imp.load_source("name", filename)
globals().update(vars(parms))
if not os.access(figDir, os.F_OK):
os.makedirs(figDir)
if not os.access(dataDir, os.F_OK):
os.makedirs(dataDir)
outFileName = dataDir + "/" + dataNameBase #+ ".dat"
with open(workDir+baseSampleDir+"/sample_info.dat", 'rb') as input:
baseSample = pickle.load(input)
for (iSample, sampleDir) in enumerate(sampleDirList):
if iSample == 0: continue
with open(workDir+sampleDir+"/sample_info.dat", 'rb') as input:
sample = pickle.load(input)
with open(workDir+sampleDirList[0]+"/sample_info.dat", 'rb') as input:
baseSample = pickle.load(input)
print " Working with", sample.fullName, "..."
print " Working with", sample.fullName, "...",
sampleName = sample.fullName
binPath = CTOOLS_PATH+"/analysis/voidOverlap"
logFile = os.getcwd()+"/mergerTree.log"
stepOutputFileName = outFileName + "_" + sampleName + "_"
stepOutputFileName = outFileName + "_" + baseSample.fullName + "_" + \
sampleName + "_"
#stepOutputFileName = os.getcwd()+"/data/overlap_"
launchVoidOverlap(baseSample, sample, workDir+sampleDirList[0],
launchVoidOverlap(baseSample, sample, workDir+baseSampleDir,
workDir+sampleDir, binPath,
thisDataPortion="central", logFile=logFile,
continueRun=False, workDir=workDir,