cwitools.measurement.luminosity

cwitools.measurement.luminosity(fits_in, redshift=None, mask=None, cosmo=FlatLambdaCDM(name="WMAP9", H0=69.3 km / (Mpc s), Om0=0.286, Tcmb0=2.725 K, Neff=3.04, m_nu=[0. 0. 0.] eV, Ob0=0.0463), var_data=None)

Measure the integrated luminosity from 1D, 2D or 3D data.

Parameters:

• fits_in (astropy HDU or HDUList) – The input data, 2D or 3D. If input is 2D, it is assumed to have units erg/s/cm2/arcsec2. If input is 1D or 3D, units are assumed to be erg/s/cm2/angstrom
• redshift (float) – The redshift of the source.
• mask (numpy.ndarray) – A binary mask of the 3D region to use, where a value of 1 means include and 0 means `exclude.’ If none provided, the entire image or data cube is summed.
• cosmology (astropy FlatLambdaCDM) – An astropy cosmology instance. WMAP9 is used by default.
• var_data (numpy.ndarray) – Array of same dimensions as data and mask, containing variance estimates. Used to propagate error on luminosity.

Returns:

The integrated luminosity of the source in erg/s. float: The error on the luminosity calculation.

Return type:

float