This is the python code for the calculation of:

• the thermal SZ anguler power spectrum (used in Makiya et al. 2018 and Makiya et al. 2020)
• Compton-Y weighted halo bias b_y (used in Chiang et al. 2020)
• redshift dervative of the Compton-Y, dy/dz (used in Chiang et al. 2020)

The tSZ power spectrum calculation is based on the fortran code szfast devoloped by Eiichiro Komatsu, but there are several important updates e.g., treatment of neutrino mass in the calculation of mass function, mass-concentration relation, etc. Please see Makiya et al. 2018 and Makiya et al. 2020 for the details of the updates.

Please cite Komatsu and Seljak 2002 and Makiya et al. 2018 if you use the code in your paper.

Another tSZ power spectrum code, class_sz, is also based on the szfast. I have confirmed that the outputs of pysz and class_sz within the numerical uncertainties.

Any questions and requests are welcome.

Notes

• You can use this script to reproduce the plot of redshift vs Compton-Y weighted halo bias b_y (Figure 10 of Chiang et al. 2020)

0. The code requires the python modules of numpy and classy. You should manually install classy.

1. Move to source directory

cd pysz/pysz/source

2. Compile fortran codes Please specify the fortran compilar you use in the Makefile. Then,

make clean && make

All required fortran modules are at the same directory.

3. then go back to top directory of pysz and install python modules

cd PYSZDIRECTORY
python setup.py install

4. Check the instlattion Run the python interpreter and type

from pysz import pysz

if it doesn't show any errors, installtion is succesfully completed.

See pysz_usage.ipynb for the usage of the code.