https://github.com/sflender/szmap
Making a Sunyaev Zel'dovich map from a halo catalog
https://github.com/sflender/szmap
Last synced: 8 months ago
JSON representation
Making a Sunyaev Zel'dovich map from a halo catalog
- Host: GitHub
- URL: https://github.com/sflender/szmap
- Owner: sflender
- Created: 2017-03-21T17:37:11.000Z (over 8 years ago)
- Default Branch: master
- Last Pushed: 2018-03-22T21:28:00.000Z (over 7 years ago)
- Last Synced: 2025-01-21T12:34:00.852Z (10 months ago)
- Language: C++
- Size: 14.6 KB
- Stars: 0
- Watchers: 3
- Forks: 0
- Open Issues: 0
-
Metadata Files:
- Readme: README
Awesome Lists containing this project
README
szmap - a fast C++ code for generating a mock tSZ map from a halo catalog
by Samuel Flender
Added beta profile by Vinu Vikraman
currently the code is using the pressure profile from Battaglia+12 (arxiv 1109.3711) & Birkinshaw 1999
You need to install Healpix, cfitsio, MPI, and gsl. To compile, type "make". To run, type "./szmap szmap.cfg".
The config parameters are:
profiletype = 1/2 for Battaglia/Beta
order = order of the Healpix map
order_subpix = if this is larger than order, then each pixel is subpixelized into pixels of the higher order. This can reduce some shotnoise at very high ell.
convert_map_to_ring = Set this to 1 if you want to have the output map in RING ordering (otherwise NESTED)
write_output_file = 1 for writing output file
output_tsz_name = name of the output file
cluster_catalog = name of the cluster catalog
nu = frequency for converting y to T
Cosmology parameters:
h, Omega_M, Omega_b, Omega_k, w_DE
max_int_radius = radius out to which to integrate the pressure profile.
trunc_radius = where to truncate the 2D SZ signal. This should be smaller than max_int_radius.
beta = Beta value in the profile
Mmin = lower mass cut
Mmax = upper mass cut
zmin = lower z cut
zmax = upper z cut