|
17 | 17 | from cluster_toolkit import _dcast, _lib |
18 | 18 | import numpy as np |
19 | 19 | from scipy.integrate import quad |
20 | | -import scipy.special as spec |
21 | 20 |
|
22 | 21 |
|
23 | 22 | __BBPS_params_P_0 = (18.1, 0.154, -0.758) |
@@ -269,42 +268,6 @@ def projected_y_BBPS(r, M, z, omega_b, omega_m, |
269 | 268 | epsrel=epsrel) |
270 | 269 |
|
271 | 270 |
|
272 | | -def Cxl(M, z, omega_b, omega_m, da, ell, |
273 | | - epsrel=1e-3): |
274 | | - ''' |
275 | | - The Fourier-transform of a cluster's Compton y parameter on the sky. Assumes |
276 | | - the flat-sky approximation. |
277 | | -
|
278 | | - Args: |
279 | | - M (float): Cluster :math:`M_{\\Delta}`, in Msun. |
280 | | - z (float): Cluster redshift. |
281 | | - omega_b (float): Baryon fraction. |
282 | | - omega_m (float): Matter fraction. |
283 | | - d_a (float): Angular diameter distance at redshift `z`. Should be \ |
284 | | - in Mpc. |
285 | | - ell (float): The Fourier-transform wavenumber. |
286 | | -
|
287 | | - Returns: |
288 | | - float: :math:`C_{x, l}`. Unitless. |
289 | | - ''' |
290 | | - return quad(lambda theta: 2 * np.pi * theta * spec.j0(ell * theta) |
291 | | - * projected_y_BBPS(theta * da, M, z, omega_b, omega_m, |
292 | | - epsrel=epsrel), |
293 | | - 0, 2 * np.pi, |
294 | | - epsrel=epsrel)[0] |
295 | | - |
296 | | - |
297 | | -def smoothed_xi(theta, M, z, omega_b, omega_m, da, |
298 | | - epsrel=1e-3, maxl=10000): |
299 | | - # Convert from arcmin to radians |
300 | | - theta = theta * 60 * np.pi / 180 |
301 | | - return quad(lambda ell: 1 / (2 * np.pi) * ell * spec.j0(ell * theta) |
302 | | - * Cxl(M, z, omega_b, omega_m, da, ell, |
303 | | - epsrel=epsrel), |
304 | | - 0, maxl, |
305 | | - epsrel=epsrel)[0] |
306 | | - |
307 | | - |
308 | 271 | ################################################## |
309 | 272 | # The following functions are for testing only!! # |
310 | 273 | ################################################## |
|
0 commit comments