lib.MapMaking.ComponentMapMaking.preset.preset_acquisition module

class lib.MapMaking.ComponentMapMaking.preset.preset_acquisition.PresetAcquisition(seed_noise_qubic, seed_noise_planck, seed_start_pcg, preset_tools, preset_external, preset_qubic, preset_sky, preset_comp, preset_mixing_matrix, preset_gain)

Bases: object

Preset Acquisition.

Instance to initialize the Components Map-Making. It defines the data acquisition variables and methods.

Parameters:

• seed_noise (int) – Seed for random noise generation.

• preset_tools (object) – Class containing tools and simulation parameters.

• preset_external (object) – Class containing external variables and methods.

• preset_qubic (object) – Class containing qubic operator and variables and methods.

• preset_sky (object) – Class containing sky varaibles and methods.

• preset_comp (object) – Class containing component variables and methods.

• preset_mixing_matrix (object) – Class containing mixing-matrix variables and methods.

• preset_gain (object) – Class containing detector gain variables and methods.

seed_noise

Seed for random noise realization.

Type:

int

params_foregrounds

Dictionary containing the paramters associated with foregrounds.

Type:

dict

ites_to_converge

Max iteration number for PCG algorithm.

Type:

int

invN

Inverse noise covariance matrix with input shape \((N_{det}*N_{samples} + N_{pix}*N_{planck}*N_{stokes})\) and output shape \((N_{det}*N_{samples} + N_{pix}*N_{planck}*N_{stokes})\).

Type:

BlockDiagonalOperator

M

Preconditioner for PCG algorithm.

Type:

DiagonalOperator

fwhm_qubic_rec

Resolution of the reconstructed maps.

Type:

array_like

H

Pointing matrix.

Type:

BlockColumnOperator

TOD_qubic

QUBIC simulated TOD.

Type:

array_like

TOD_external

Planck simulated TOD.

Type:

array_like

TOD_obs

Simulated observed TOD.

Type:

array_like

beta_iter

Spectral indices, if d1 \((12 \times Nside_{\beta}^2, N_{comp}-1)\), if not \((N_{comp}-1)\).

Type:

array_like

allbeta

Spectral indeices, if d1 \((iter, 12 \times Nside_{\beta}^2, N_{comp}-1)\), if not \((iter, N_{comp}-1)\).

Type:

array_like

Amm_iter

Mixing matrix.

Type:

array_like

Methods

get_convolution()	Convolutions.
get_noise()	Noise.
get_tod()	TOD.
get_x0()	PCG starting point.

get_convolution()

Convolutions.

Method to define all angular resolutions of the instrument at each frequency.

This method sets the Full Width at Half Maximum (FWHM) for Time-Ordered Data (TOD) and map-making processes. self.fwhm_qubic_tod represents the real angular resolution, and self.fwhm_qubic_mapmaking represents the beams used for reconstruction.

The method checks the convolution_in and convolution_out parameters to determine the appropriate FWHM values. It also calculates the reconstructed FWHM based on these parameters.

Returns:

• fwhm_qubic_tod (array_like) – Resolutions to create input TOD.

• fwhm_qubic_mapmaking (array_like) – Resolutions for map-making.

• fwhm_qubic_rec (array_like) – Resolutions of reconstructed maps.

get_noise()

Noise.

Method to define QUBIC noise, can generate Dual band or Wide Band noise by following:

• Dual Band: n = [Ndet + Npho_150, Ndet + Npho_220]

• Wide Band: n = [Ndet + Npho_150 + Npho_220]

Depending on the instrument type specified in the preset_qubic parameters, this method will instantiate either QubicWideBandNoise or QubicDualBandNoise and return the total noise.

Returns:

noise – The total noise array, flattened.

Return type:

array_like

get_tod()

TOD.

Generate fake observational data from QUBIC and external experiments.

This method simulates observational data, including astrophysical foregrounds contamination using self.beta and systematics using self.g. The data generation follows the formula: \(\vec{d} = H.A.\vec{c} + \vec{n}\).

H

Pointing matrix operator.

Type:

Operator

TOD_qubic

The simulated observational data for QUBIC.

Type:

array_like

TOD_external

The simulated observational data for external experiments.

Type:

array_like

TOD_obs

The combined observational data from QUBIC and external experiments.

Type:

array_like

nsampling_x_ndetectors

The number of samples in self.TOD_qubic.

Type:

int

get_x0()

PCG starting point.

Define starting point of the convergence.

The argument ‘set_comp_to_0’ multiplies the pixel values by a given factor. You can decide to convolve the map by a beam with an FWHM in radians.

This method initializes the beta_iter and Amm_iter attributes based on the foreground model parameters. It also applies convolution and noise to the components based on the preset parameters.

Raises:

TypeError – If an unrecognized component name is encountered.