Package can be downloaded here (size 34.7 MB)

NEW: Mathematica notebook including code for Fisher matrix calculation

If you use the Fisher matrix code you can additionally cite arXiv:1303.7212

*** Note for arbitrary primordial power spectrum ***

Arbitrary input power spectrum can be specified by replacing the dQdz with an integral over the power spectrum. See the "frun[]" function in the new mathematica notebook with Fisher matrix calculation (idistort_fisher_beta.nb). Replace the integrand in frun with your ( power spectrum * Exp[-2 k^2 /(mparsec*kd[z])^2] ) and copy the next few lines until dQdz[z_]=... into the main code replacing the dQdz for silk damping. The definition of power spectrum in the integrand is such that the standard power spectrum is (k/k0)^ns, with k0 the pivot point and ns=0.96 the Planck value. This routine numerically integrates the modes dissipating in different redshift bins to calculate the power injected into the CMB at those redshifts.

Typical runtime for a spectrum calculation is < 30 seconds