Abstract
Broadband receiver data need color corrections applying to correct for the different source spectra across their wide bandwidths. The full integration over a receiver bandpass may be computationally expensive and redundant when repeated many times. Color corrections can be applied, however, using a simple quadratic fit based on the full integration instead. Here we describe fastcc and interpcc, quick Python and IDL codes that return, respectively, color correction coefficients for different power-law spectral indices and modified blackbodies for various Cosmic Microwave Background related experiments. The codes are publicly available, and can be easily extended to support additional telescopes.
Export citation and abstract BibTeX RIS
![](https://content.cld.iop.org/journals/2515-5172/6/12/252/revision1/rnaasaca6eblicense.gif)
Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
1. Background
Cosmic Microwave Background (CMB) observations require telescopes using high sensitivity receivers with wide bandwidths that are a large fraction of the observing frequency. For example, the 28.4 GHz channel of the Planck satellite was sensitive to photons at frequencies between ∼24 and 34 GHz (see Figure 1), a fractional bandwidth (i.e., Δν/ν) of ∼35%. With such high fractional bandwidths, the measured flux density of an astronomical source depends on the convolution of the observing spectral bandpass with the source spectrum. For example, a bandpass that is most sensitive at low frequencies will result in an increase in the measured flux density of a steep spectrum radio source compared to a flatter spectrum source. Color corrections avoid these effects by correcting the flux density to that which a monochromatic receiver operating at a reference frequency ν0 would have measured.
Figure 1. Left: example bandpasses from Planck LFI, showing how these can be complex structures that extend significantly beyond the reference frequency for broadband detectors (dashed lines). Right: corresponding color correction fits for various spectral indexes.
Download figure:
Standard image High-resolution imageThe color correction thus depends on the spectral index of the source, and if the value for this changes (e.g., while fitting a model with multiple spectral components, or when adding new observations), then the color correction also changes. For this reason, raw flux densities are published without color corrections, but color corrections are applied for any scientific analysis and for the inclusion of data in data visualization and subsequent analysis.
The fastcc and interpcc codes provide fast color corrections for CMB data sets, particularly in MCMC fitting techniques.
2. Formalism
We follow and expand on the formalism from Planck Collaboration (2016). When the source spectrum is described by a power-law, S ∝ να , the color correction is the ratio of integrals over frequency ν, with the numerator integrating over the bandpass g(ν) (if measured in temperature units) and the spectrum the data is calibrated to, δ, and the denominator integrating over the combination of the bandpass and the source spectrum with spectral index α, as
![Equation (1)](https://content.cld.iop.org/journals/2515-5172/6/12/252/revision1/rnaasaca6ebeqn1.gif)
where ν0 is the reference frequency. If the instrument bandpass has been measured in intensity units, τ(ν), this becomes:
![Equation (2)](https://content.cld.iop.org/journals/2515-5172/6/12/252/revision1/rnaasaca6ebeqn2.gif)
For instruments calibrated in CMB thermodynamic units, this becomes:
![Equation (3)](https://content.cld.iop.org/journals/2515-5172/6/12/252/revision1/rnaasaca6ebeqn3.gif)
where , with x = h
ν/(kB
TCMB) to adjust for the CMB spectrum in Rayleigh–Jeans temperature units, with Planck constant h, Boltzmann constant kB, and CMB temperature TCMB = 2.7255 K. However, the ratio of ηΔT
is usually small, only becoming important when ν0 is far from the center of the bandpass.
The resulting color corrections are very smooth (see Figure 1, right), and can be fitted with a quadratic function:
![Equation (4)](https://content.cld.iop.org/journals/2515-5172/6/12/252/revision1/rnaasaca6ebeqn4.gif)
The source flux density is then corrected as:
![Equation (5)](https://content.cld.iop.org/journals/2515-5172/6/12/252/revision1/rnaasaca6ebeqn5.gif)
The color corrections for thermal dust spectra can be calculated using the same equations as for power laws, but changing the source spectrum for a graybody, with the dust temperature Tdust and spectral index βdust as input parameters. For equations, see the fastcc explanatory supplement. A quadratic fit becomes increasingly different from the color corrections at extreme values. As such, we also include interpcc, which uses interpolation over a set of cached color correction values.
There are two equivalent approaches to color correction: either change the flux density of the source at a fixed reference frequency, or change the reference frequency for the observations. When characterizing multiple sources within a survey, the analysis is normally simpler if all sources are at the same frequency with corrected flux density. Alternatively, effective frequencies νeff, for which C(νeff, α) = 1, can be calculated from the color corrections (in thermodynamic units: convert to intensity at νeff) via:
![Equation (6)](https://content.cld.iop.org/journals/2515-5172/6/12/252/revision1/rnaasaca6ebeqn6.gif)
3. Code
The code takes input α for fastcc, or Tdust, and βdust for interpcc, and returns the factor that the measured flux density must be multiplied to obtain the corrected flux density. We include both Python and IDL versions of fastcc, and interpcc in Python. Test/example scripts are also included for both fastcc and interpcc. The code is available at https://github.com/mpeel/fastcc and 10.5281/zenodo.7376510 , along with a more detailed explanatory supplement. An IDL version for Planck LFI only was released in the Planck Legacy Archive. 5 While no dependencies are required for fastcc, interpcc and the color correction calculations use Astropy (Astropy Collaboration et al. 2018), Numpy (Harris et al. 2020), and scipy (Virtanen et al. 2020).
The code now includes color corrections for Planck LFI (Planck Collaboration 2014a, 2014b, 2016) and HFI (Planck Collaboration 2014b), WMAP (Bennett et al. 2013), IRAS (IRAS 1988), DIRBE, 6 QUIJOTE MFI (Rubiño-Martín et al. 2022, MNRAS in press; Genova-Santos et al. 2022, in preparation), and C-BASS (King et al. 2014; Taylor et al. 2022, in preparation). Additional experiments can be added on request and provision of bandpasses.
Footnotes
- 5
- 6