I've created reconstructions that are similar to Alexia's configuration. They are in the following directories:
workingDir = 'run201077_1424'
workingDir = 'run2010712_1255'
And the following log files have the code I used to generate these:
../logs/100714log.py
../logs/100714_2log.py
Specs of the runs two runs
Run 1:workingDir = 'run201077_1424'
#Angular Correlation
oversample = 5.
corrBins = 40.0 #should be one more than you want it to be
mincorr = 6.0
maxcorr = 12.0
convo = 180./pi
corrspace = (maxcorr - mincorr)/corrBins
corrarray = makecorrarray(mincorr, maxcorr, corrspace)
#Auto Correlation
nbins = 40
boxside = 1000.0 #in Mpc/h
xicorrmin = 5.0
xicorrmax = 50.0
nxicorr = 40
#make redshift bins
#boxside = 1000.0 #in Mpc/h
rbinmin = 0.0 #in Mpc/h
rbinspace = boxside/10 #in Mpc/h
rbinmax = boxside/2.0 + 1.0 #in Mpc/h
rbins = makeRbins(rbinmin,rbinmax,rbinspace)
nrbins = rbins.size - 1
#Catalog sizes
psubsize=100000
ssubsize=100000
Run 2 (larger):workingDir = 'run2010712_1255'
#Angular Correlation
oversample = 5.
corrBins = 40.0 #should be one more than you want it to be
mincorr = 6.0
maxcorr = 12.0
convo = 180./pi
#Auto Correlation
boxside = 1000.0 #in Mpc/h
xicorrmin = 5.0
xicorrmax = 50.0
nxicorr = 40
#make redshift bins
#boxside = 1000.0 #in Mpc/h
rbinmin = 0.0 #in Mpc/h
rbinspace = boxside/20 #in Mpc/h
rbinmax = boxside/2.0 + 1.0 #in Mpc/h
rbins = makeRbins(rbinmin,rbinmax,rbinspace)
nrbins = rbins.size - 1
#Catalog sizes
psubsize=1000000
ssubsize=1000000
I've uploaded the run directories onto the repository for Alexia to look at, but I'm getting the same problems as before.
I've also updated the files for running and re-running the reconstruction. The primary change is that I now save all the run constants into the run directory, including the masks, so that everything can be recovered when you go back to an old run without looking at the initial run directory.
The code to run a full reconstruction in a new directory is here:
../pythonjess/jessreconstructionFull.py
The code to read in all the run constants from an old run (where the correlation functions have already been computed) and perform the reconstruction is here:
../pythonjess/jessreconstruction.py
Simply change first line of the code, to the working directory of the old run:
workingDir = 'run2010712_1255'
You can also change line 35 to skip redshift bins with too noisy correlation functions:
skip = 0
If all you want to do is rebin the reconstruction, use the following code:
../pythonjess/rebinjessrecontruction.py
And change the workingDir, npbins, rmin and rmax to be the appropriate run directory, number of reconstruction bins, min redshift bin, and max redshift bin.
Or you can fun the following function from correlationData:
reconstructionrebin(workingDir, npbins, rmin, rmax)
Alexia is moving and is still looking to see if she can figure out where the bug is, but it seems to me I am at the point where I need to possibly start from her working reconstruction, and step by step implement my changes to see where things break down.
In the mean time, I need to write the likelihood paper, so I'm going to work on that for now until Alexia is back online. This sucks.
