exosim.tasks.foregrounds.estimateZodi

Classes

EstimateZodi

It estimate the zodiacal radiance in the target direction for a specific wavelength range

Module Contents

class EstimateZodi

Bases: exosim.tasks.task.Task

It estimate the zodiacal radiance in the target direction for a specific wavelength range

Returns:

zodiacal radiance

Return type:

Radiance

Examples

>>> estimateZodi = EstimateZodi()
>>> wavelength = np.logspace(np.log10(0.45), np.log10(2.2), 6000) * u.um
>>> zodi = estimateZodi(wavelength=wavelength, zodiacal_factor=1)

or, given the pointing direction

>>> zodi = self.estimateZodi(wavelength=wavelength,  coordinates=(90 * u.deg, -66 * u.deg))

execute()

Class execution. It runs on call and executes all the task actions returning the outputs. It requires the input with correct keywords

model(a, wl)

The used zodiacal model is based on the zodiacal light model presented in Glasse et al. 2010

\[I_{zodi}(\lambda) = a \left( 3.5 \cdot 10^{-14} BB(\lambda, 5500 \, K) + 3.52 \cdot 10^{-8} BB(\lambda, 270 \, K) \right)\]

where \(BB(\lambda, T)\) is the Planck black body law and \(a\) is the fitted coefficient.

Parameters:

• a (float) – zodiacal multiplicative factor. Default is 0.

• wl (ndarray or Quantity) – wavelength grid. If no units are attached is considered as expressed in um

Returns:

zodiacal radiance

Return type:

Radiance

zodiacal_fit_direction(coordinates, zodi_map=None)

In this case the \(A\) coefficient is selected by a precompiled grid. The grid has been estimated by fitting our model with Kelsall et al. (1998) data.

A custom map can be provided, to replace the default one, as long as it matches the format.

Parameters:

• coordinates ((float, float)) – pointing coordinates as (ra, dec)

• zodi_map (str (optional)) – map file containing the zodiacal coefficients per sky positions. A default map is inclluded in ExoSim, that contains the coefficients fitted over Kelsall et al. 1998 model.

Returns:

zodiacal factor for Glasse et al. 2010 zodiacal model

Return type:

float