Correlate BORG reconstructions with CMB data

Python 100%

Find a file

Deaglan Bartlett 635dd2a862 Add gitattributes file		2026-04-08 10:18:33 +01:00
cmborg	Ensure halo cut-outs are sufficiently wide for beam	2026-04-08 09:58:43 +01:00
notes	Add profile timing notes	2026-03-13 15:54:31 +00:00
scripts	Ensure halo cut-outs are sufficiently wide for beam	2026-04-08 09:58:43 +01:00
tests	Ensure halo cut-outs are sufficiently wide for beam	2026-04-08 09:58:43 +01:00
.gitattributes	Add gitattributes file	2026-04-08 10:18:33 +01:00
.gitignore	Add test suite and update .gitignore	2026-03-13 15:52:00 +00:00
MANIFEST.in	Refactor package structure and clean up codebase	2026-02-20 15:01:03 +00:00
pyproject.toml	Refactor package structure and clean up codebase	2026-02-20 15:01:03 +00:00
README.md	Ensure halo cut-outs are sufficiently wide for beam	2026-04-08 09:58:43 +01:00
requirements-dev.txt	Add code as of Aquila Fall 2025 meeting	2025-12-15 16:32:07 +00:00
requirements.txt	Refactor package structure and clean up codebase	2026-02-20 15:01:03 +00:00

README.md

CMBORG

Cross-correlate BORG density field reconstructions with CMB data (thermal Sunyaev-Zel'dovich effect).

Installation

Clone the repository and install in development mode:

git clone https://git.aquila-consortium.org/Aquila-Consortium/CMBORG
cd CMBORG
pip install -e .

For development (testing, linting):

pip install -e ".[dev]"

Quick start

Run a toy mock inference with scripts/main.py. This generates a synthetic tSZ map from a halo catalog, adds noise, and runs NUTS to recover the profile parameters (P0, xc0, n):

cd scripts
python main.py

Outputs: tsz_map_healpix.png, trace_plot.png, corner_plot.png.

Tests

pytest tests/

Tests cover the interpolation module (1D/2D Catmull-Rom splines, correction factor accuracy vs scipy), tSZ profile models (physical sanity, interpolated vs analytic agreement, JAX differentiability), and cosmological distances.

License

MIT