Add HalfDome TSZ to pysm

CHANGES.rst

@@ -1,6 +1,7 @@
 Unreleased
 ==========
 - Add NumPy 2 support by replacing deprecated `np.trapz` with `trapezoid` https://github.com/galsci/pysm/pull/232
+- Add HalfDome TSZ maps (new 'tsz2' model) to pysm. Create new documentation and code following the same structural approach used for WebSky.
 
 3.4.3 (2025-10-02)
 ==================

docs/halfdome.rst

@@ -0,0 +1,37 @@
+.. _halfdome:
+
+HalfDome
+******
+
+HalfDome cosmological simulations are a set of simulations tailored for the joint analysis of Stage-IV surveys. For more information, see: [https://halfdomesims.github.io/] (https://halfdomesims.github.io/)
+
+Cosmic Infrared Background
+--------------------------
+
+TODO
+
+Radio Galaxies
+--------------
+
+TODO
+
+Thermal SZ Effect
+-----------------
+
+Provided are eleven realizations of the Compton-y parameter. They were generated with xgpaint bataglia16 profiles at $N_{side}=8192$.
+
+Kinetic SZ Effect
+-----------------
+
+TODO
+
+
+Lensing Convergence
+-------------------
+
+TODO
+
+Primary and lensed CMB
+----------------------
+
+TODO

docs/models.rst

@@ -138,8 +138,11 @@ Sunyaev–Zeldovich emission
 
 - **tsz1**: Thermal SZ emission from WebSky 0.4. Available at $N_{side}=8192$. For more details see :ref:`websky`.
 
+- **tsz2**: Thermal SZ emission from HalfDome 0.1. Generated using xgpaint with Battaglia16 profiles. Eleven realizations available at $N_{side}=8192$. For more details see :ref:`halfdome`.
+
 - **ksz1**: Kinetic SZ emission from WebSky 0.4. Available at $N_{side}=4096$. For more details see :ref:`websky`.
 
+
 Radio galaxies
 ==============
 

src/pysm3/data/presets.cfg

@@ -394,6 +394,13 @@ class = "WebSkySZ"
 version = "0.4"
 sz_type = "thermal"
 max_nside = 4096
+[tsz2]
+class = "HalfDomeSZ"
+version = "0.1"
+sz_type = "thermal"
+map_sz = "y_b16_halo_res1_s{seed}.fits"
+seed = 0
+max_nside = 8192
 [dip1]
 # From NPIPE https://doi.org/10.1051/0004-6361/202038073
 #

src/pysm3/models/__init__.py

@@ -18,5 +18,13 @@
     WebSkySZ,
     WebSkyRadioGalaxies,
 )
+from .halfdome import (
+    HalfDomeSZ,
+#    SPT_CIB_map_scaling,
+#    HalfDomeCMB,
+#    HalfDomeCIB,
+#    HalfDomeSZ,
+#    HalfDomeRadioGalaxies,
+)
 from .catalog import PointSourceCatalog
 from .dipole import CMBDipole