NonGaussianity 2013

X-ray Cluster Constraints on Non-Gaussianity

Code for evaluating the non-Gaussian mass function

These files contain C code which evaluates the series in equation 2.9 of the paper below for the hierarchical and feeder scalings (equations A.13 and A.14).

The result in sum should be multiplied by the Gaussian mass function to obtain the non-Gaussian one. The variables referenced here should be initialized as follows:

Note that this is C, so the indices of a size N array run from 0 to N-1. Some of the arrays are bigger than strictly necessary to keep the indices from becoming confusing (e.g. M[0], M[1] and M[2] are not used).

Reference:

X-ray cluster constraints on non-Gaussianity.

S. Shandera, A. Mantz, D. Rapetti, and S. W. Allen. JCAP, 8:4, 2013. ADS, arXiv, bibTeX