COSMOMC (v2010 Jan) 2dF support?

[Michael Drinkwater]

Is analysis of the 2dF mpk data still supported in CosmoMC ver Jan 2010?

I have compiled with EXTDATA = LRG and it works fine just using the LRG, but gives errors when I try to add the 2dF data.

I put this in params.ini; the rest is unchanged from default settings

#filenames for matter power spectrum datasets, incl twodf
mpk_numdatasets = 2
mpk_dataset1 = data/lrgDR7kmax02kmin02newmaxLv2ALL_MAGCOVv3.dataset
#mpk_dataset1 = data/sdss_lrgDR4.dataset
mpk_dataset2 = data/2df_2005.dataset
use_mpk = T

The log file shows it reading LRG and 2dF data files OK, but then gives this error and appears to hang with no further output to any files:

k/h out of bounds in MatterPowerAt ( 0.7870000 )

Cheers, Michael.

[Antony Lewis]

I think it should work if you don't compile with that option.

[Michael Drinkwater]

Yes that works; thanks for the quick reply!

I compiled without the LRG option in the Makefile and it runs OK reading the 2dF data.

I was trying to read both LRG and 2dF at the same time to investigate my options for adding our new WiggleZ mpk data into the fits. Ideally we would have the option of fitting both at the same time since ours are at significantly higher redshift. Is fitting two mpk sets at once a sensible (or possible) approach with the code?

We want to write the necessary code to feed the WiggleZ mpk into cosmomc. I'm hoping we can modify Beth's code for the LRGs, but do you have any general advice before we start, especially so what we do is compatible with the way the LRG data are currently used?

Cheers, Michael