-
Notifications
You must be signed in to change notification settings - Fork 0
2- and 3-d aperture mass map computation from shear catalogs
License
GPL-2.0, GPL-2.0 licenses found
Licenses found
GPL-2.0
LICENSE
GPL-2.0
COPYING
joergdietrich/mapmap
Folders and files
Name | Name | Last commit message | Last commit date | |
---|---|---|---|---|
Repository files navigation
README for mapmap ----------------- mapmap creates maps of aperture mass (M_ap, Schneider 1996). It takes a FITS_LDAC catalog as input and produces FITS images of aperture mass maps ("mapmaps") and fits cubes of tomographic peak maps (Hennawi & Spergel 2005). Optionally it also produces maps of M_ap significance computed analytically following Eq. 13 of Schirmer et al. (2006, astro-ph/0607022), maps of -M_ap, mapmaps of the cross-component of the shear, and mapmaps for a randomized input catalog. For the latter option the random seed is hardwired to a constant in the program. The input catalog must at least provide columns with (1) right ascencion, (2) declination, (3) 1-component of the shear, and (4) 2-component of the shear. Optionally a fifth column can be used to assign weights to individual shear estimators. Note that although the positions are expected in right ascencion and declination, the program expects shear values measured in a Cartesian frame with the values of the 1-axis increasing from left to right. The table and column names in the input catalog can be configured via command line options or the configuration file. Two filter functions are built in: The polynomial filter of Schneider et al. (1998, option SCHNEIDER) and the truncated pseudo-NFW filter of Schirmer et al. (2006, filter SCHIRMER). Different filter functions can be read from a file (option FILE). The filter radius is set by assuming a fiducial cluster model and specifying the redshift of the model cluster and its radius in kpc. The size of the output grid is determined by the physical scale at the redshift of the fiducial model. The default is that one pixel corresponds to 50 kpc at the model redshift (default 0.3). mapmap is configured via a SExtractor style configuration file or command line options. mapmap -c conf_file specifies the location of the configuration file. If the -c option is not provided mapmap tries to read mapmap_default.conf in the current working directory. If that is not found internal default values are used. mapmap -d prints the internal defaults in the format of the configuration file. Limitations ----------- mapmap computes a very simple transformation from equatorial to pixel coordinates. It will probably fail for very large areas (never tested) and will definitely fail at the 24h -> 0h transition in right ascencion and at the poles. Compilation ----------- Compilation of mapmap requires the cfitsio library, the GNU Scientific Library (gsl) with CBLAS extension. Compilation is done with the usual ./configure make make install Efficient use of parallel processing via OpenMP can be enable during the configure step with the option --enable-openmp. Some machines (notably massive shared ememory Altix) do not maintain cache coherence without further instructions by the user. In this case attaching subprocesses explicitly to a CPU helps. This can be enabled with the Portable Linux Processor Affinity package included in the distribution. It is turned off by default because it decreases performance on most systems but can be enabled with the configuration switch --enable-plpa.
About
2- and 3-d aperture mass map computation from shear catalogs
Resources
License
GPL-2.0, GPL-2.0 licenses found
Licenses found
GPL-2.0
LICENSE
GPL-2.0
COPYING
Stars
Watchers
Forks
Releases
No releases published
Packages 0
No packages published