Cosmomc SDSS Lyman-Alpha patch.

A completely rewritten likelihood evaluator, now all in Fortran 90.

Please not that this new code has not been extensively tested and any problems should be reported to mclovin@berkeley.edu.

DISCLAIMER: This code works by calculating the amplitude, slope and curvature of power spectrum at the k=0.009s/km and then uses table to interpolate to a corresponding chi2. This table is not intended for use without WMAP, or for models with power spectra qualitatively different in shape from LCDM.

Download this patch. Apply it as follows:

patch -p0 <LyA070329.patch

and compile as usual. The patch might alter your make file so you might need to reapply some changes (WMAP directory, etc.). To use this LyAlpha constraints specify LyAlphaSDSS as an mpk dataset.

This is the dataset that was used in Seljak, Slosar and McDonald 2006.

Please cite

@Article{McDonald:2004eu,
     author    = "McDonald, Patrick and others",
     title     = "The Lyman-alpha Forest Power Spectrum from the Sloan
                  Digital Sky Survey",
     journal   = "Astrophys. J. Suppl.",
     volume    = "163",
     year      = "2006",
     pages     = "80-109",
     eprint    = "astro-ph/0405013",
     SLACcitation  = "%%CITATION = ASTRO-PH 0405013;%%"
} 

and

@Article{McDonald:2004xn,
     author    = "McDonald, Patrick and others",
     title     = "The Linear Theory Power Spectrum from the Lyman-alpha
                  Forest in the Sloan Digital Sky Survey",
     journal   = "Astrophys. J.",
     volume    = "635",
     year      = "2005",
     pages     = "761-783",
     eprint    = "astro-ph/0407377",
     SLACcitation  = "%%CITATION = ASTRO-PH 0407377;%%"}

Acknowledgement of me (A Slosar) is appreciated but not necessary.

A Slosar, 29/03/07