wale.LoadSimulations module
- wale.LoadSimulations.get_simulation_l1(cosmo_index_to_run, tomobin, edges, centers, snr, R_pixels=30, filter_type='tophat', plot=False)[source]
Load simulation data for a specific cosmology and compute L1 norms and PDFs. Parameters: - cosmo_index_to_run: Index of the cosmology to run (0-9 for 10 different cosmologies). - R_pixels: Physical scale in pixels for smoothing. - filter_type: Type of filter to use for smoothing (‘tophat’ or ‘gaussian’). Returns: - sim_l1_runs: Array of L1 norms for each simulation realization. - sim_pdf_runs: Array of PDF counts for each simulation realization. - avg_sim_l1: Average L1 norm across all realizations. - std_sim_l1: Standard deviation of L1 norms across realizations. - avg_sim_pdf: Average PDF counts across all realizations. - std_sim_pdf: Standard deviation of PDF counts across realizations.
- wale.LoadSimulations.get_smoothed_app_pdf(mass_map, window_radius, binedges, filter_type, **kwargs)[source]
Applies top-hat smoothing in Fourier space at two scales and returns the PDF of the difference map.
The map is filtered with a top-hat window of radius R and 2R, then the difference is computed.
- Parameters:
mass_map – 2D numpy array.
window_radius – The smoothing scale (R) in physical units.
binedges – Bin edges for the histogram.
L – Physical size of the map (default 505 MPC/h).
- Returns:
(bin_edges, pdf_counts, difference_map)
- Return type:
tuple