Difference between revisions of "Map-making LFI"

From Planck Legacy Archive Wiki
Jump to: navigation, search
(Low-resolution maps and Noise Covariance Matrices)
(Low-resolution maps and noise covariance matrices)
 
(35 intermediate revisions by 2 users not shown)
Line 1: Line 1:
==Map-making==
+
==Mapmaking==
The input of the map-making step consists of the calibrated timelines, along with the corresponding pointing information.
+
The inputs to the mapmaking procedure consist of the calibrated timelines, along with the corresponding pointing information.
 
The main output consists of temperature and polarization maps.  
 
The main output consists of temperature and polarization maps.  
An important part of the map-making step is the removal of correlated 1/<math>f</math> noise.
+
An important part of the mapmaking step is the removal of correlated 1/<i>f</i> noise.
  
LFI maps were constructed with the Madam map-making code, version 3.7.4.
+
LFI maps were constructed with the MADAM mapmaking code, version 3.7.4.
The code is based on generalized destriping technique, where the correlated noise component is modeled as a sequence of constant offset, "baselines".  
+
The code is based on a generalized destriping technique, where the correlated noise component is modelled as a sequence of constant offset "baselines".  
 
The baseline solution is constrained by a noise filter.
 
The baseline solution is constrained by a noise filter.
As auxiliary information the code produces a hit count map and a white noise covariance matrix.
+
As auxiliary information, the code produces a hit-count map and a white noise covariance matrix.
No beam information is used, but the signal is assigned to the pixel where the center of the beam falls.
+
No beam information is used, with the signal being simply assigned to the pixel where the centre of the beam falls.
  
The chosen baseline length was one second for the 44GHz and 70GHz maps, 0.25sec for the 30GHz map. This gives a good noise removal,
+
The chosen baseline length was 1s for the 44GHz and 70GHz maps, 0.25s for the 30GHz map. This gives good noise removal,
without being computationally heavy.
+
without being computationally burdensome.
The noise filter was built according to the noise parameters listed in Table 1.
+
The noise filter was built according to the noise parameters (see noise section).
 +
Flagged samples were excluded from the analysis. The Galaxy region was masked out in the destriping phase, to reduce errors arising from strong signal gradients. Radiometers were combined according to the horn-uniform weighting scheme to minimize systematics. The polarization component was included in the analysis and is part of this release.
 +
A detailed description of the mapmaking procedure is given in {{PlanckPapers|planck2013-p02}}, {{PlanckPapers|planck2014-a03}}, {{PlanckPapers|planck2014-a07}} no relevant changes has been applyed in the 2018 release. See also the section on [[Frequency Maps]].
  
Flagged samples were excluded from the analysis. The galaxy region was masked out in the destriping phase, to reduce error arising from strong signal gradients. Radiometers were combined according to the horn-uniform weighting scheme to minimize systematics. The polarization component was included in the analysis, although only the temperature maps are released.
+
The maps are in HEALPix format, at resolution <i>N</i><sub>side</sub>=1024 for all frequencies with an additional map at <i>N</i><sub>side</sub>=2048 for the LFI 70GHz channel, in the nested pixelization scheme.
A detailed description of the map-making procedure is given in {{PlanckPapers|planck2013-p02}} {{PlanckPapers|planck2014-a03||Planck-2015-A03}} {{PlanckPapers|planck2014-a07||Planck-2015-A07}}. See also section [[Frequency Maps]].
 
 
 
The maps are in Healpix format, at resolution <math>N_{\rm side}=</math>1024 for all frequencies with an additional map at <math>N_{\rm side}</math>=2048 for the LFI 70 GHz channel, in nested pixeling scheme.
 
 
Unobserved pixels are marked by a special value.
 
Unobserved pixels are marked by a special value.
  
The released maps are in galactic coordinates.
+
The released maps are in Galactic coordinates.
The conversion from ecliptic to galactic coordinate system is described by rotation matrix
+
The conversion between ecliptic and Galactic coordinates is described by the rotation matrix
 
<math>
 
<math>
 
   \newcommand{\Re}{\mathrm{Re}\,}
 
   \newcommand{\Re}{\mathrm{Re}\,}
Line 39: Line 38:
 
The conversion was applied to the input pointing data, prior to the construction of the map.
 
The conversion was applied to the input pointing data, prior to the construction of the map.
  
[[File:lfi_noise.jpg|650px]]
+
==Low-resolution maps and noise covariance matrices==
:'''Table 1: Knee frequencies and slopes for the LFI radiometers.'''
 
<!--
 
{| class="wikitable"  align="center" style="text-align:center" border="1" cellpadding="5" cellspacing="0"
 
|+ '''Table 1: Noise parameters used in the construction of the noise filter, and the radiometer weights.'''
 
|-bgcolor="ffdead"
 
!Radiometer || White noise sigma [mK<math>_{\rm CMB}</math>] || Slope ||Knee frequency [mHz] || <math>f_{\rm min}</math> [Hz] || Horn-uniform weight
 
|-
 
|'''LFI18M''' || 4.542 ||  -1.07 || 15.3 || 1.15e-5 || 5.297E+04
 
|-
 
|'''LFI18S'''  || 4.139 ||  -1.20 || 18.3 || 1.15e-5 || 5.297E+04
 
|-
 
|'''LFI19M''' || 5.146  ||    -1.22 || 11.9 ||  1.15e-5 || 3.948E+04
 
|-
 
|'''LFI19S'''  || 4.918 ||  -1.12 || 14.6 || 1.15e-5 || 3.948E+04
 
|-
 
|'''LFI20M''' || 5.214 ||  -1.31 || 8.4 || 1.15e-5 || 3.482E+04
 
|-
 
|'''LFI20S'''  || 5.500 ||  -1.34 || 6.0 || 1.15e-5 || 3.482E+04
 
|-
 
|'''LFI21M''' || 4.008 ||  -1.26 || 39.3 || 1.15e-5 || 4.897E+04
 
|-
 
|'''LFI21S'''  || 4.978 ||  -1.24 || 14.0 || 1.15e-5 || 4.897E+04
 
|-
 
|'''LFI22M''' || 4.348 ||  -1.53 || 10.1 || 1.15e-5 || 4.864E+04
 
|-
 
|'''LFI22S'''  || 4.713 ||  -1.20 || 15.9 || 1.15e-5 || 4.864E+04
 
|-
 
|'''LFI23M''' || 4.468 ||  -1.07 || 30.2 || 1.15e-5 || 4.670E+04
 
|-
 
|'''LFI23S'''  || 4.782 ||  -1.21    || 58.8 || 1.15e-5 || 4.670E+04
 
|-
 
|'''LFI24M''' || 3.145 ||  -0.94 || 26.9  || 1.15e-5 || 1.158E+05
 
|-
 
|'''LFI24S'''  || 2.717 ||  -0.91 || 73.0  || 1.15e-5 || 1.158E+05
 
|-
 
|'''LFI25M''' || 2.821 ||  -0.85 || 20.1 ||  1.15e-5 || 1.319E+05
 
|-
 
|'''LFI25S'''  || 2.684 || -0.90 ||  46.1 ||  1.15e-5 || 1.319E+05
 
|-
 
|'''LFI26M''' || 3.280 ||  -0.92 ||  64.4 ||  1.15e-5 || 1.056E+05
 
|-
 
|'''LFI26S'''  || 2.859 || -0.88 ||  43.8 ||  1.15e-5 || 1.056E+05
 
|-
 
|'''LFI27M''' || 1.609 ||  -0.93 || 175.1 ||  1.15e-5 || 3.567E+05
 
|-
 
|'''LFI27S'''  || 1.737 ||  -0.91 ||  109.6 ||  1.15e-5 || 3.567E+05
 
|-
 
|'''LFI28M''' || 1.814 ||  -0.93 ||  127.9 ||  1.15e-5 || 3.353E+05
 
|-
 
|'''LFI28S'''  || 1.635 ||  -0.91 || 43.9 ||  1.15e-5 || 3.353E+05
 
|}
 
-->
 
  
==Low-resolution maps and Noise Covariance Matrices==
+
To fully exploit the information contained in the large-scale structure of the microwave sky, pixel-pixel covariances are needed in the maximum likelihood estimation of the CMB power spectrum. However, full covariance matrices are impossible to employ at the native map resolution due to resource limitations. A low-resolution data set is therefore required for the low-&#8467; analysis. This data set has been packed into three different files, one per frequency, called "LFI_NoiseCovMat_0??_0016_R3.00.tgz", that can be downloaded from the Cosmology section of the Planck Legacy Archive.  
To fully exploit the information contained in the large scale structure of the microwave sky, pixel-pixel covariances are needed in the maximum likelihood estimation of the CMB power spectrum. However, full covariance matrices are impossible to employ at the native map resolution due to resource limitations. A low-resolution dataset is therefore required for the low-<math>l</math> analysis. This dataset '''has been packed into three different files, one per frequency, called "LFI_NoiseCovMat_0??_0016_R2.00.tgz", that can be downloaded from the Cosmology section of the Planck Legacy Archive.  
+
They consist of low-resolution maps, and descriptions of residual noise present in those maps given by pixel-pixel noise covariance matrices (NCVMs).
They consist of low-resolution maps, and descriptions of residual noise present in those maps given by pixel-pixel noise covariance matrices (NCVMs).'''
+
Note the in the 2018 release the Low-resolution maps full mission coverage, excluding Surveys 2 and 4 has not been used. We release them for crosschecking purposes with respect 2015 release.
  
 
{| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:left" width=800px
 
{| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:left" width=800px
|+ Contents of '''LFI_NoiseCovMat_0??_0016_R2.00.tgz'''
+
|+ Contents of '''LFI_NoiseCovMat_0??_0016_R3.20.tgz'''
 
|- bgcolor="ffdead"
 
|- bgcolor="ffdead"
 
! Filenames || Comment
 
! Filenames || Comment
 
|-   
 
|-   
| LFI_SkyMap_0??_0016_nobs_DX11D_full_regnoise.fits || Low resolution maps. Full mission coverage
+
| LFI_SkyMap_0??_0016_coswin_DX12_full_regnoise.fits || Low-resolution maps. Full mission coverage.
 +
|- 
 +
| LFI_SkyMap_0??_0016_coswin_DX12_Corrected_full_regnoise.fits || Low-resolution maps BandPass Corrected. Full mission coverage. Corrected for the BandPass.
 +
|-
 +
| LFI_SkyMap_0??_0016_coswin_DX12_s1-s3-s5-s6-s7-s8_regnoise.fits || Low-resolution maps. Full mission coverage, excluding Surveys 2 and 4.
 +
|-
 +
| LFI_SkyMap_0??_0016_coswin_DX12_Corrected_s1-s3-s5-s6-s7-s8_regnoise.fits || Low-resolution maps. Full mission coverage, excluding Surveys 2 and 4. Corrected for the BandPass.
 
|-  
 
|-  
| LFI_SkyMap_0??_0016_nobs_DX11D_s1-s3-s5-s6-s7-s8_regnoise.fits || Low resolution maps. Full mission coverage, excluding Survey 2 and 4
+
| offset_covmat_toast_nside64T16_coswin_025sec_??GHz_DX12_full_bin.dat || Low-resolution noise covariance matrices. Regularization noise added. Format: C unformatted.
 
|-  
 
|-  
| offset_covmat_toast_nside64T16_nobs_025sec_??GHz_DX11D_full_bin.dat || Low resolution noise covariance matrices. Regularisation noise added. Format: C unformatted.
+
| offset_covmat_toast_nside64T16_coswin_025sec_??GHz_DX12_s1-s3-s5-s6-s7-s8_bin.dat || Low-resolution noise covariance matrices. Regularization noise added. Format: C unformatted.
 
|-  
 
|-  
| offset_covmat_toast_nside64T16_nobs_025sec_??GHz_DX11D_s1-s3-s5-s6-s7-s8_bin.dat || Low resolution noise covariance matrices. Regularisation noise added. Format: C unformatted.
+
| covmat_toast_nside64T16_coswin_025sec_??GHz_DX12_full_bin.dat || Low-resolution noise covariance matrices. No regularization noise added. Format: C unformatted.
 
|-  
 
|-  
| covmat_toast_nside64T16_nobs_025sec_??GHz_DX11D_full_bin.dat || Low resolution noise covariance matrices. No regularisation noise added. Format: C unformatted.
+
| covmat_toast_nside64T16_coswin_025sec_??GHz_DX12_s1-s3-s5-s6-s7-s8_bin.dat || Low-resolution noise covariance matrices. No regularization noise added. Format: C unformatted.
 
|-  
 
|-  
| covmat_toast_nside64T16_nobs_025sec_??GHz_DX11D_s1-s3-s5-s6-s7-s8_bin.dat || Low resolution noise covariance matrices. No regularisation noise added. Format: C unformatted.
+
| mask_comm17Tv2_pol_R2.2x_ns16_nest.fits || Low-resolution Mask 2.2.
 
|-  
 
|-  
| DX11D_0??_016_full_conv060_VarianceBPcorrection_v2.fits ||
+
| mask_comm17Tv2_pol_R1.8x_ns16_nest.fits || Low-resolution Mask 1.8.
 
|-  
 
|-  
| DX11D_0??_016_full_conv060_BPcorrection_v2.fits ||
+
| offset_RESCALED_base_bestfit_cov_IQU.da || Rescaled noise ONLY for 70 GHz. C unformatted.
 
|-  
 
|-  
| QUmask_DX11_16_0??GHz_R1.50.fits ||
 
|-
 
 
|}
 
|}
  
  
  
The low-resolution dataset can currently be utilized efficiently only at resolution <math>N_{\rm side} = 16</math>, or lower. All the low-resolution data products are produced at this target resolution.
+
The low-resolution data set can currently be utilized efficiently only at resolution <i>N</i><sub>side</sub> = 16, or lower. All the low-resolution data products are produced at this target resolution.
===Low Resolution Maps===
+
===Low-resolution maps===
A number of different schemes to obtain the low resolution maps are discussed in {{BibCite|keskitalo2013}}. We chose to downgrade the maps using the inverse noise weighting. See {{PlanckPapers|planck2013-p02}} {{PlanckPapers|planck2014-a07||Planck-2015-A07}} for discussion.
+
A number of different schemes to obtain the low-resolution maps are discussed in {{BibCite|keskitalo2013}}. We chose to downgrade the maps using the inverse noise weighting, no changes on the procedure has been applyed to the 2018 release. This is discussed further in {{PlanckPapers|planck2013-p02}} {{PlanckPapers|planck2014-a07}}.
 
====Inputs====
 
====Inputs====
We took the high resolution maps described in [[Map-making LFI#Map-making|Map-making]] and [[Frequency Maps]], and the corresponding <math>3 \times 3</math> matrices as an input for this analysis step.
+
We took the high-resolution maps described in [[Map-making LFI#Map-making|Map-making]] and [[Frequency Maps]], and the corresponding 3&times;3 matrices as an input for this analysis step.
  
 
====Production====
 
====Production====
The high resolution maps were downgraded to <math>N_\mathrm{side} = 16</math> using inverse noise weights (given by the <math>3 \times 3</math> matrices), and subsequently the temperature part was smoothed with a symmetric Gaussian beam with <math>\mbox{FWHM} = 440\,\mbox{arcmin}</math>.
+
The high-resolution maps were downgraded to <i>N</i><sub>side</sub> = 16 using inverse noise weights (given by the 3&times;3 matrices), and subsequently the temperature part was smoothed with a symmetric Gaussian beam with FWHM = 440arcmin.
  
===Noise Covariance Matrices===
+
===Noise covariance matrices===
  
 
The statistical description of the residual noise in the maps is given in the form of a pixel-to-pixel noise covariance matrix (NCVM), as described in {{BibCite|keskitalo2013}}.  
 
The statistical description of the residual noise in the maps is given in the form of a pixel-to-pixel noise covariance matrix (NCVM), as described in {{BibCite|keskitalo2013}}.  
Line 139: Line 90:
 
====Inputs====
 
====Inputs====
  
The noise model was determined by three noise parameters: white noise level <math>\sigma</math>, slope, and knee frequency <math>f_{\rm knee}</math>. We actually used three sets of noise parameters one for the entire mission (noise parameters are listed in Table 1), and one for each sky survey (SS1 and SS2).
+
The noise model was determined by three noise parameters: the white noise level &sigma;; slope; and knee frequency <i>f</i><sub>knee</sub>. We actually used three sets of noise parameters, one for the entire mission (noise parameters are listed in Table 1), and one for each sky survey (SS1 and SS2).
  
 
We used the same pointing as in the noise Monte Carlo simulations. See the description in [[Map-making LFI#Noise Monte Carlo Simulation#Inputs|Noise Monte Carlo Simulation Inputs]].
 
We used the same pointing as in the noise Monte Carlo simulations. See the description in [[Map-making LFI#Noise Monte Carlo Simulation#Inputs|Noise Monte Carlo Simulation Inputs]].
  
We used the gap files produced during the making of the flight maps to leave out samples that were flagged bad for various reasons.
+
We used the gap files produced during the making of the flight maps to leave out samples that were flagged as bad for various reasons.
  
 
====Production====
 
====Production====
  
The output of the NCVM module of Madam map-maker are inverse NCVMs.  Since the inverse matrices are additive, we divided the computations into a number of small chunks to save computational resources. We first calculated one inverse NCVM per radiometer per survey at resolution <math>N_\mathrm{side}=32</math>, and then combined these individual inverse matrices to form the actual inverse matrices. The map-making parameters were almost identical to the standard map-making runs. The differing parameter values are listed below:
+
The output of the NCVM module of MADAM mapmaker are inverse NCVMs.  Since the inverse matrices are additive, we divided the computations into a number of small chunks to save computational resources. We first calculated one inverse NCVM per radiometer per survey at resolution <i>N</i><sub>side</sub>=32, and then combined these individual inverse matrices to form the actual inverse matrices. The mapmaking parameters were almost identical to the standard mapmaking runs. The differing parameter values are listed below:
* Baseline lengths were 0.25 s (8 samples), 0.50 s (24 samples), and 0.50 s (39 samples) for 30 GHz, 44 GHz, and 70GHz, respectively.
+
* baseline lengths were 0.25s for 30GHzand 1.0s  for 44GHz, and 70GHz;
* The calculations were performed at resolution <math>N_\mathrm{side} = 32</math>.
+
* the calculations were performed at resolution <i>N</i><sub>side</sub> = 64;
* No destriping mask was applied.
+
* no destriping mask was applied;
* The horns were weighted optimally.
+
* the horns were weighted optimally.
  
 +
To obtain the noise covariance from its inverse, the matrices are inverted using the eigen decomposition of a matrix.  The monopole of the temperature map cannot be resolved by the mapmaker, and thus the matrix becomes singular. This ill-determined mode is left out of the analysis.
  
To get the noise covariance from its inverse, the matrices are inverted using the eigen decomposition of a matrix.  The monopole of the temperature map cannot be resolved by the map-maker, and thus the matrix becomes singular. This ill-determined mode is left out of the analysis.
+
Having calculated the eigen decomposition in the previous step, we can apply the same linear operators to modify the eigenvectors as were applied to the high-resolution maps while downgrading them. The eigenvectors are downgraded to <i>N</i><sub>side</sub> = 16 using inverse noise weights, and subsequently the temperature part is smoothed with a symmetric Gaussian beam with FWHM = 440arcmin.
 
 
Having calculated the eigen decomposition in the previous step, we can apply the same linear operators to modify the eigenvectors as were applied to the high resolution maps while downgrading them. The eigenvectors are downgraded to <math>N_\mathrm{side} = 16</math> using inverse noise weights, and subsequently the temperature part is smoothed with a symmetric Gaussian beam with <math>\mbox{FWHM} = 440\,\mbox{arcmin}</math>.
 
  
 
The final matrices are then recomposed from the original eigenvalues and modified eigenvectors.
 
The final matrices are then recomposed from the original eigenvalues and modified eigenvectors.
  
The low resolution noise covariance matrices
+
The low-resolution noise covariance matrices:
* are in C binary format files.
+
* are C binary format files;
 
* are organized in block form,
 
* are organized in block form,
 
<math>
 
<math>
Line 174: Line 124:
 
             QI & QQ & QU \\
 
             QI & QQ & QU \\
 
             UI & UQ & UU
 
             UI & UQ & UU
             \end{matrix} \right).
+
             \end{matrix} \right);
 
\end{align}   
 
\end{align}   
 
</math>
 
</math>
  
* are in HEALPix nested pixelisation scheme. Resolution is <math>N_\mathrm{side} = 16</math>, and thus there are <math>N_\mathrm{pix} = 3072</math> pixels.
+
* are in the HEALPix nested pixelisation scheme (with resolution is <i>N</i><sub>side</sub> = 16, and thus there are <i>N</i><sub>pix</sub> = 3072 pixels);
* are in Galactic coordinates.
+
* are in Galactic coordinates;
* have units <math>\mathrm{K}_{\mathrm{CMB}}</math>.
+
* have K<sub>CMB</sub> units.
  
==Half-ring jackknife Noise Maps==
+
==Half-ring jackknife noise maps==
  
 
===Overview===
 
===Overview===
 +
In the 2018 release we follow the same procedure as in {{PlanckPapers|planck2014-a03||Planck-2015-A03}} in order to estimate the noise directly at the map level and in the angular power spectra.
  
In order to estimate the noise directly at the map level and in the angular power spectra, we divided the time-ordered data into two halves and produced half-ring jackknife maps as described in Planck First Results: II. The Low Frequency Instrument performance and data processing.
+
Briefly, instead of using the full time-ordered data as described above, we produced two sets of maps using either only the first half of each pointing period (map named <b>j</b><sub>1</sub> below) or only the second half of each pointing period (map named <b>j</b><sub>2</sub>). At each pixel <i>p</i>, these half-ring jackknife maps <b>j</b><sub>1</sub> and <b>j</b><sub>2</sub> contain the same sky signal, since they result from the same scanning pattern on the sky.  However, because of instrumental noise, the maps <b>j</b><sub>1</sub> and <b>j</b><sub>2</sub> are not identical.
 
 
Briefly: Instead of using the full time ordered data as described above, we produced two sets of maps using either only the first half of each pointing period (map named <math>\mathbf{j_1}</math> below) or only the second half of each pointing period (map named <math> \mathbf{j_2}</math>). At each pixel <math>p</math>, these half-ring jackknife maps <math>\mathbf{j_1}</math> and <math>\mathbf{j_2}</math> contain the same sky signal, since they result from the same scanning pattern on the sky.  However, because of instrumental noise, the maps <math>\mathbf{j_1}</math> and <math>\mathbf{j_2}</math> are not identical.
 
  
We estimated the noise level in each map <math>\mathbf{m}</math> made using the full(ring) data, by constructing a half-ring difference map
+
We estimated the noise level in each map <b>m</b> made using the full (ring) data, by constructing a half-ring difference map
  
 
<math>\mathbf{n_{m}}(p) =  [ \mathbf{j_1}(p) - \mathbf{j_2}(p)] \ / \ \mathbf{w_{\rm hit}}(p)\,,</math>
 
<math>\mathbf{n_{m}}(p) =  [ \mathbf{j_1}(p) - \mathbf{j_2}(p)] \ / \ \mathbf{w_{\rm hit}}(p)\,,</math>
Line 199: Line 148:
 
       \frac{1}{\mathbf{hit_2}(p)} \right]}\,</math>.
 
       \frac{1}{\mathbf{hit_2}(p)} \right]}\,</math>.
  
Here <math>\mathbf{hit_{full}}(p) = \mathbf{hit_1}(p) + \mathbf{hit_2}(p)</math>
+
Here <b>hit</b><sub>full</sub>(<i>p</i>) = <b>hit</b><sub>1</sub>(<i>p</i>) + <b>hit</b><sub>2</sub>(<i>p</i>)
is the hit count at pixel <math>p</math> in the full map <math>\mathbf{m}</math>, while <math>\mathbf{hit_1}</math> and <math>\mathbf{hit_2}</math> are the hit counts of <math>\mathbf{j_1}</math> and <math>\mathbf{j_2}</math>, respectively.  The weight factor <math>\mathbf{w_{hit}}(p)</math> is equal to <math>2</math> only in those pixels where <math>\mathbf{hit_1}(p) = \mathbf{hit_2}(p)</math>.  In a typical pixel, <math>\mathbf{hit_1}(p)</math> will differ slightly from
+
is the hit count at pixel <i>p</i> in the full map <b>m</b>, while <b>hit</b><sub>1</sub> and <b>hit</b><sub>2</sub> are the hit counts of <b>j</b><sub>1</sub> and <b>j</b><sub>2</sub>, respectively.  The weight factor <b>w</b><sub>hit</sub>(<i>p</i>) is equal to 2 only in those pixels where <b>hit</b><sub>1</sub>(<i>p</i>) = <b>hit</b><sub>2</sub>(<i>p</i>) .  In a typical pixel, <b>hit</b><sub>1</sub>(<i>p</i>) will differ slightly from
<math>\mathbf{hit_2}(p)</math> and hence the weight factor is <math>\mathbf{w_{hit}}(p) > 2</math>.
+
<b>hit</b><sub>2</sub>(<i>p</i>)  and hence the weight factor is <b>w</b><sub>hit</sub>(<i>p</i>)>2.
 
 
The half-ring difference maps <math>\mathbf{n_{m}}</math> are the most direct measure of the noise in the actual maps. The other noise estimates (NCVM and noise Monte Carlo) rely on specific modelling of the noise and this modelling can be validated by comparing to the half-ring difference maps. However, the half-ring difference maps can only capture the noise that varies faster than half of the duration of the pointing period, i.e., the noise whose frequency is <math>f \gtrsim 1/20\,\mbox{min} = 0.85\,</math>mHz.
 
 
 
We calculated the noise maps <math>\mathbf{n_{m}}</math>, from half-ring jackknife maps for temperature (I) and polarization (Q and U) and as a first quality check of the maps (and as one of the tests of the whole data processing pipeline up to the maps) tested both numerically and visually that these noise maps divided pixel-by-pixel by square root of the white noise covariance maps were approximately Gaussian with variance near to unityTemperature noise maps for the nominal survey and for the first and second sky surveys are shown in the next subsection. Further we calculated from the noise maps the temperature and polarization (E and B mode) auto-correlation and cross-correlation noise angular power spectra by anafast and compared to these the results from the white noise covariance matrices and from the noise Monte Carlo simulations. A similar comparison was made between downgraded half-ring noise maps, downgraded noise Monte Carlo maps and the low resolution noise covariance maps.  Detailed results are presented in the Systematic Effects paper.
 
 
 
===Examples of Half-ring Difference Maps and Noise Angular Power Spectra===
 
 
 
====Hit Count Weighted Half-ring Difference Maps====
 
 
 
Here we show some hitcount weighted half-ring difference maps, i.e., noise maps for temperature. The columns are for different LFI frequencies: 30, 44, and 70 GHz. The rows are the nominal survey (a bit more than 1 year of observations), survey 1 (the first sky survey, approx first half a year of observations) and survey 2 (the second sky survey, approx the second half a year of observations). Some features are visible in particular in the galactic plane. These are due to "gradient leakage". (In regions where the gradient in the sky signal is very large even a tiny difference in the pointing of the first and second half of each pointing period causes the signal to "leak" to the half-ring difference map. In practice this is not a problem for noise estimation, since these regions - the galaxy, orion, crab nebula, etc - will be masked in the cosmology analysis.
 
 
 
<B>Half-ring difference maps</B> calculated at the native <math>N_{\rm side}</math> 1024 resolution (note for the LFI 70 GHz channel maps are also generated at <math>N_{\rm side}</math>=2048). Columns: frequency (30, 44, 70 GHz), rows: sky survey (nominal, survey_1, survey_2).
 
[[File:LFI_4_5_5_3_DX9_Delta_TemperatureNoiseMap_1s_frequency_and_SScomparison_NoSmoothing.png|thumb|900px|center]]
 
 
 
Half-ring difference maps (the same as above, but) <B>smoothed</B> with 60 arcmin fwhm Gaussian.  Columns: frequency (30, 44, 70 GHz), rows: sky survey (nominal, survey_1, survey_2).
 
[[File:LFI_4_5_5_3_DX9_Delta_TemperatureNoiseMap_1s_frequency_and_SScomparison_fwhm60arcmin.png|thumb|900px|center]]
 
 
 
====Hit Count Weighted Half-ring Difference Maps Normalized by sqrt of white noise variance at each pixel====
 
 
 
Now we show the same as above, but divided by the square root of the estimate of white noise variance in each pixel. These normalized noise maps should be approximately Gaussian with a unit variance (at the native resolution), apart from some stripes that are due to correlated (non-Gaussian) <math>1/f</math> noise. The large-scale <math>1/f</math> noise is more apparent in the smoothed version of the figure that follows after the native resolution version.
 
 
 
<B>Normalized Half-ring difference maps</B> calculated at the native <math>N_{\rm side}</math> 1024 resolution (note for the LFI 70 GHz channel maps are also generated at <math>N_{\rm side}</math>=2048).  Columns: frequency (30, 44, 70 GHz), rows: sky survey (nominal, survey_1, survey_2).
 
[[File:LFI_4_5_5_3_DX9_Delta_NormalizedTemperatureNoiseMap_1s_frequency_and_SScomparison_NoSmoothing.png|thumb|900px|center]]
 
 
 
Normalized Half-ring difference maps (the same as above, but) <B>smoothed</B> with 60 arcmin fwhm Gaussian.  Columns: frequency (30, 44, 70 GHz), rows: sky survey (nominal, survey_1, survey_2).
 
[[File:LFI_4_5_5_3_DX9_Delta_NormalizedTemperatureNoiseMap_1s_frequency_and_SScomparison_fwhm60arcmin.png|thumb|900px|center]]
 
 
 
====Noise Angular Power Spectra from Half-ring Difference Maps====
 
 
 
The noise angular power spectra calculated by anafast from half-ring difference temperature maps and normalized by the sky coverage to estimate the noise level if there was a full sky coverage in order to make a comparison of different surveys easier. Columns: frequency (30, 44, 70 GHz). Colors: black = nominal, red = survey_1, blue = survey_2.
 
[[File:LFI_4_5_5_3_DX9_Delta_SSnoiseComparison_loglin_TT.png|thumb|900px|center]]
 
 
 
===Comparison of Noise Calculated from Half-ring Difference and from Other Noise Estimates===
 
 
 
Here we compare noise angular power spectra estimated from half-ring difference maps (RED) to the estimate from white noise covariance maps (BLUE) and the full noise Monte Carlo simulations (BLACK, top curve 16% quantile, middle curve 50% quantile, i.e., median, and bottom curve 84% quantile) - see the next section for the details of noise Monte Carlo.  
 
  
 +
The half-ring difference maps <b>n</b><sub>m</sub> are the most direct measure of the noise in the actual maps. The other noise estimates (NCVM and noise Monte Carlo) rely on specific modelling of the noise and this modelling can be validated by comparing to the half-ring difference maps. However, the half-ring difference maps can only capture the noise that varies faster than half of the duration of the pointing period, i.e., the noise whose frequency is approximately <i>f</i> > 1/20min = 0.85mHz.
  
[[File:DX9_Delta_noise_comparison_nominal_nomask.png|thumb|900px|center|Nominal Survey]]
+
We calculated the noise maps <b>n</b><sub>m</sub>, from half-ring "jackknife" maps for temperature (<i>I</i>) and polarization (<i>Q</i> and <i>U</i>) and as a first quality check of the maps (and as one of the tests of the whole data processing pipeline up to the map level) tested both numerically and visually that these noise maps divided pixel-by-pixel by the square root of the white noise covariance maps were approximately Gaussian with variance near to unity. Temperature noise maps for the nominal survey and for the first and second sky surveys are shown in the next subsection. Furthermore we calculated from the noise maps the temperature and polarization (E and B mode) auto-correlation and cross-correlation noise angular power spectra using HEALPix anafast and compared to these the results from the white noise covariance matrices and from the noise Monte Carlo simulations. A similar comparison was made between downgraded half-ring noise maps, downgraded noise Monte Carlo maps, and the low-resolution noise covariance maps.  Detailed results are presented in {{PlanckPapers|planck2016-l02}}.
  
[[File:DX9_Delta_noise_comparison_survey1_nomask.png|thumb|900px|center|Survey_1]]
+
===Comparison of noise estimates using Half-Ring===
  
[[File:DX9_Delta_noise_comparison_survey2_nomask.png|thumb|900px|center|Survey_2]]
+
Here we compare noise angular power spectra estimated from half-ring difference maps (red), white noise covariance maps (black dash-dotted lines), and 100 full noise Monte Carlo simulations (grey band showing range for 16th  and 84th quantiles of noise simulations, and the black solid lines giving the median, i.e., 50th quantile, of distributions). See the next section for details of noise Monte Carlo simulations. From top to botto we show ''TT'', ''EE'' and ''BB'' power spectra for 30 GHz (left), 44 GHz (centre), and 70 GHz (right). Half-ring spectra are binned with <math>\mathbf{\Delta}l = 25 </math> for <math> l \mathbf{\geq} 75 </math>.
  
====High-ell average (<math>\ell =1150..1800</math>) noise relative to the white noise estimate====
+
[[File:test_ffp10.png|thumb|800px|center|<b>Consistency Check</b>]]
  
The same as previous figures, but here the noise comparison is made from the high ell tails of the angular power spectra, where the white noise dominates. We have taken the average of <math>C_\ell</math> from the multipoles between 1150 and 1800, and subtracted the white noise estimate. Half-ring difference noise estimate is RED, the full noise Monte Carlo estimate is BLACK, and the white noise estimate from noise Monte Carlo is GREEN. The white noise estimate from the white noise covariance map (WNC), BLUE, has been subtracted from all the results.  
+
Below the null-tests comparing power spectra from survey differences to those from teh half-ring maps are showed. Difefrences are: left, Survey 1 - Survey 2; Middle , Survey 1- Survey 3; and right, Survey 1 - Survey4. these are for 30 GHz (top), 44 GHz (middle), and 70 GHz (bottom), for both ''TT'' and ''EE'' power spectra. There is a significant improvement in Surve1 - Survey 2 and Survey 1 - Survey 4 at 30 GHz, especially in ''EE''. See {{PlanckPapers|planck2016-l02}} for further details.
  
 +
[[File:nullhr1517.png|thumb|800px|center|<b>Null Test</b>]]
  
[[File:Comparison_Cls_NMC_WNC_JN_WNCmean_removed_nominal_nomask.png|thumb|700px|center|Nominal survey]]
 
  
 +
====High-ell average noise relative to white noise estimate====
  
[[File:Comparison_Cls_NMC_WNC_JN_WNCmean_removed_survey1_nomask.png|thumb|700px|center|Survey_1]]
+
The figure below is the same as the previous figures, but here the noise comparison is made from the high &#8467; tails of the angular power spectra, where the white noise dominates. We have taken the average of <i>C</i><sub>&#8467;</sub> from multipoles between 1150 and 1800 for both temperature and polarization and tehn comparing with the WNCVM. As already shown in previous releases, there is still an excess of 1/f noise, meaning tha both the real data and the noise MCs predict slightly larger noise than the WNCVM. It is important to note that such noise excess is reduced considerably with respect to the 2015 release.
  
 +
[[File:high.png|thumb|800px|center|<b>Ratio at high multipoles</b>]]
  
[[File:Comparison_Cls_NMC_WNC_JN_WNCmean_removed_survey2_nomask.png|thumb|700px|center|Survey_2]]
+
==Noise Monte Carlo simulations==
 
 
==Noise Monte Carlo Simulation==
 
  
 
===Overview===
 
===Overview===
 
Calculating and handling full pixel-to-pixel noise covariance matrices for Planck maps if feasible only at low resolution.
 
Calculating and handling full pixel-to-pixel noise covariance matrices for Planck maps if feasible only at low resolution.
To support the analysis of high-resolution maps, a Monte Carlo set of noise maps were produced.  These maps were produced from noise timelines using the same map-making procedure as for the flight data.  In the noise Monte Carlo it was possible to follow exactly the map-making procedure used for the flight maps, whereas for the calculation of noise covariance matrices some approximations had to be made.
+
To support the analysis of high-resolution maps, a Monte Carlo set of noise maps were produced.  These maps were produced from noise timelines using the same map-making procedure as for the flight data.  In the noise Monte Carlo it was possible to follow exactly the mapmaking procedure used for the flight maps, whereas for the calculation of noise covariance matrices some approximations had to be made.
Such noise Monte Carlos were produced at two levels of the analysis:  1) LFI Monte Carlo (MC) as part of the LFI data processing, and 2) Full Focal Plane (FFP) Monte Carlo as part of the joint HFI/LFI data processing. This page describes the LFI noise MC.  For the FFP MC, see [[HL-sims]] and [[Simulation data]].
+
Such noise Monte Carlos were produced at two levels of the analysis:  (1) LFI Monte Carlo (MC) as part of the LFI data processing procedure; and (2) Full Focal Plane (FFP) Monte Carlos as part of the joint HFI/LFI data processing. This page describes the LFI noise MCs.  For the FFP MC, see [[HL-sims]] and [[Simulation data]].
  
 
===Inputs===
 
===Inputs===
The noise MC uses a three-parameter (white noise level (<math>\sigma</math>), slope, and knee frequency (<math>f_\mathrm{knee}</math>)) noise model, where the noise consists of white noise and correlated (<math>1/f</math>) noise and the latter has a power spectrum
+
The noise MC uses a three-parameter noise model, consisting of white noise level (&sigma;), slope, and knee frequency (<i>f</i><sub>knee</sub>)). Here the noise consists of white noise and correlated 1/<i>f</i> noise, with a power spectrum
  
:<math> P(f) = \frac{2\sigma^2}{f_\mathrm{sample}}\left(\frac{f}{f_\mathrm{knee}}\right)^\mathrm{slope} </math>.
+
:<math> P(f) = \frac{2\sigma^2}{f_\mathrm{sample}}\left(\frac{f}{f_\mathrm{knee}}\right)^\mathrm{slope} </math>,
  
Here <math>f_\mathrm{sample}</math> is the sampling frequency of the instrument.  The noise parameters were determined separately for each radiometer as described in the section [[TOI processing LFI#Noise| Noise]] above, assuming they stayed constant over the mission.   
+
where <i>f</i><sub>sample</sub> is the sampling frequency of the instrument.  The noise parameters were determined separately for each radiometer, as described in the section [[TOI processing LFI#Noise| Noise]] above, assuming they stayed constant over the mission.   
  
The detector pointing was reconstructed from satellite pointing, focal-plane geometry, pointing correction (tilt angle), and sample timing, using Level-S simulation software.  The same pointing solution (two focal planes) was used as for the LFI flight maps.  Due to numerical accuracy, the detector pointing in the noise MC was not exactly the same as for the flight maps, so some data samples (of the order of one in a thousand) whose pointing was near the pixel boundary ended up assigned to the neighboring pixel. During the map-making from the flight data, a gap file was produced to represent the samples that were omitted from map-making due to various flags.  This gap file was used in the noise MC instead of the full set of flags.  The flight map-making used destriping masks to exclude regions of strong signal gradients from contributing to the noise baseline solution.  These same destriping masks (one for each frequency channel) were used for the noise MC.
+
The detector pointing was reconstructed from satellite pointing information, focal-plane geometry, pointing correction (tilt angle), and sample timing, using Level-S simulation software.  The same pointing solution (two focal planes) was used as for the LFI flight maps.  Due to numerical accuracy, the detector pointing in the noise MC was not exactly the same as for the flight maps, so some data samples (of the order of one in a thousand) whose pointing was near the pixel boundary ended up assigned to the neighbouring pixel. During mapmaking from the flight data, a "gap file" was produced to represent the samples that were omitted from mapmaking due to various flags.  This gap file was used in the noise MC instead of the full set of flags.  The flight mapmaking procedure used destriping masks to prevent regions of strong signal gradients from contributing to the noise baseline solution.  These same destriping masks (one for each frequency channel) were used for the noise MC.
  
 
===Production===
 
===Production===
The noise was generated internally in the Madam map-making code using a Stochastic Differential Equation (SDE) method, to avoid time-consuming writing and reading noise timelines to/from disk.  Noise for each pointing period was generated separately, using a double-precision random number seed constructed from the realization number, radiometer number, and the pointing period number; to allow regeneration of the same noise realization when needed.  White noise and <math>1/f</math> noise were generated separately.   
+
The noise was generated internally in the Madam mapmaking code using a stochastic differential equation (SDE) method, to avoid time-consuming writing and reading of noise timelines to and from disk.  Noise for each pointing period was generated separately, using a double-precision random number seed constructed from the realization number, radiometer number, and the pointing period number; this allowed for regeneration of the same noise realization when needed.  White noise and 1/<i>f</i> noise were generated separately.   
  
The same map-making code (Madam) with the same parameter settings was used for the noise MC as for the flight maps.
+
The same mapmaking code (Madam) with the same parameter settings was used for the noise MC as for the flight maps.
In addition to the destriped maps from the full noise (output maps), also binned maps from just the white noise (binned white noise maps) were produced; they represent the white noise part of the output maps.  The difference between these two maps represent the residual correlated noise in the output map.  The maps were made at Healpix resolution <math>N_\mathrm{side} = 1024</math> for all LFI frequency channels and also at Healpix resolution <math>N_\mathrm{side} = 2048</math> for the 70 GHz channel.
+
In addition to the destriped maps from the full noise (output maps), binned maps from just the white noise (binned white noise maps) were produced; they represent the white noise part of the output maps.  The difference between these two maps represents the residual correlated noise in the output map.  The maps were made at HEALPix resolution <i>N</i><sub>side</sub> = 1024 for all LFI frequency channels and also at HEALPix resolution <i>N</i><sub>side</sub> = 2048 for the 70 GHz channel.
For low-resolution analysis, these maps were downgraded (and the temperature part was smoothed) to <math>N_\mathrm{side} = 32</math> and <math>N_\mathrm{side} = 16</math>.
+
For low-resolution analysis, these maps were downgraded (and the temperature part was smoothed) to <i>N</i><sub>side</sub> = 32 and <i>N</i><sub>side</sub> = 16.
  
In addition to frequency maps for the nominal survey, also single-survey and 70 GHz horn-pair maps were produced in the noise MC.  For each case 102--1026 realizations were produced.
+
In addition to frequency maps for the nominal survey, also single-survey and 70 GHz horn-pair maps were produced in the noise MC.  For each case 102-1026 realizations were produced.
  
 
===Usage===
 
===Usage===
  
These noise Monte Carlo maps were used for a number of things in LFI data analysis.  They were compared to the low-resolution noise covariance matrices, generated for the same noise model, in order to reveal the impact of the approximations in the noise covariance matrix calculation.  They were compared to the half-ring noise maps to see how well the noise model matches the noise in the flight maps (noting, however, that the half-ring noise maps misrepresent the lowest noise frequencies in the flight maps, and contain some effects from the sky signal).  They were also used in power spectrum estimation and non-Gaussianity estimation.
+
These noise Monte Carlo maps were used for a number of purposes in LFI data analysis.  They were compared to the low-resolution noise covariance matrices, generated for the same noise model, in order to reveal the impact of the approximations in the noise covariance matrix calculation.  They were compared to the half-ring noise maps to see how well the noise model matches the noise in the flight maps (noting, however, that the half-ring noise maps misrepresent the lowest noise frequencies in the flight maps, and contain some effects from the sky signal).  They were also used in power spectrum estimation and non-Gaussianity estimation.
  
 
===Examples===
 
===Examples===
  
As an example, we show below images of the first realization of the 70 GHz frequency map noise for the nominal survey.  The images are in order: destriped full noise, binned white noise, residual correlated noise.  Note that it is difficult to see any difference between the first two images, since the residual correlated noise is more than an order of magnitude below the white noise level. The units are CMB K.  
+
As an example, we show below images of the first realization of the 70GHz frequency map noise for the nominal survey.  The images are in the order: destriped full noise; binned white noise; and residual correlated noise.  Note that it is difficult to see any difference between the first two images, since the residual correlated noise is more than an order of magnitude below the white noise level. The units here are K<sub>CMB</sub>.  
  
  
  
[[File:LFI_4_5_5_4_madam_mask_1sec_DB10_1_noise_70GHz_all_DX9delta_nom_1024outmap.00000.gif|thumb|800px|center|destriped full noise]]
+
[[File:LFI_4_5_5_4_madam_mask_1sec_DB10_1_noise_70GHz_all_DX9delta_nom_1024outmap.00000.gif|thumb|800px|center|<b>Destriped full noise.</b>]]
[[File:LFI_4_5_5_4_madam_mask_1sec_DB10_1_white_70GHz_all_DX9delta_nom_1024binmap.00000.gif|thumb|800px|center|binned white noise]]
+
[[File:LFI_4_5_5_4_madam_mask_1sec_DB10_1_white_70GHz_all_DX9delta_nom_1024binmap.00000.gif|thumb|800px|center|<b>Binned white noise.</b>]]
[[File:LFI_4_5_5_4_madam_mask_1sec_DB10_1_rcnoise_70GHz_all_DX9delta_nom_1024map.00000.gif|thumb|800px|center|residual correlated noise]]
+
[[File:LFI_4_5_5_4_madam_mask_1sec_DB10_1_rcnoise_70GHz_all_DX9delta_nom_1024map.00000.gif|thumb|800px|center|<b>Residual correlated noise.</b>]]
  
  
The following two images show the statistics of the angular power spectra of 101 realizations of the 70 GHz frequency map noise for the nominal survey. The thick black line shows the median <math>C_\ell</math>, and the green line the mean <math>C_\ell</math>. Thin black lines show the minimum, 16% quantile, 84% quantile, and the maximum <math>C_\ell</math>. The red line is the 102<math>^{nd}</math> realization.  The first plot is for the full noise in the output map, the second plot is for the residual correlated noise.
+
The following two images show the statistics of the angular power spectra of 101 realizations of the 70 GHz frequency map noise for the nominal survey. The thick black line shows the median <i>C</i><sub>&#8467;</sub>, while the green line the mean <i>C</i><sub>&#8467;</sub>. Thin black lines show the minimum, 16th percentile, 84th percentile, and the maximum <i>C</i><sub>&#8467;</sub>. The red line is the 102nd realization.  The first plot is for the full noise in the output map, while the second plot is for the residual correlated noise.
  
[[File:LFI_4_5_5_4_cl_TT_stat_noisemc101_1sec_70GHz_all_DX9delta_nom.png|thumb|800px|center|full noise]]
+
[[File:LFI_4_5_5_4_cl_TT_stat_noisemc101_1sec_70GHz_all_DX9delta_nom.png|thumb|800px|center|<b>Full noise.</b>]]
[[File:LFI_4_5_5_4_cl_TT_stat_rcnoisemc101_1sec_70GHz_all_DX9delta_nom.png|thumb|800px|center|residual correlated noise]]
+
[[File:LFI_4_5_5_4_cl_TT_stat_rcnoisemc101_1sec_70GHz_all_DX9delta_nom.png|thumb|800px|center|<b>Residual correlated noise.</b>]]
  
 
==References==
 
==References==
Line 307: Line 222:
  
 
<References />
 
<References />
 +
  
  
  
 
[[Category:LFI data processing|004]]
 
[[Category:LFI data processing|004]]

Latest revision as of 16:04, 6 July 2018

Mapmaking[edit]

The inputs to the mapmaking procedure consist of the calibrated timelines, along with the corresponding pointing information. The main output consists of temperature and polarization maps. An important part of the mapmaking step is the removal of correlated 1/f noise.

LFI maps were constructed with the MADAM mapmaking code, version 3.7.4. The code is based on a generalized destriping technique, where the correlated noise component is modelled as a sequence of constant offset "baselines". The baseline solution is constrained by a noise filter. As auxiliary information, the code produces a hit-count map and a white noise covariance matrix. No beam information is used, with the signal being simply assigned to the pixel where the centre of the beam falls.

The chosen baseline length was 1s for the 44GHz and 70GHz maps, 0.25s for the 30GHz map. This gives good noise removal, without being computationally burdensome. The noise filter was built according to the noise parameters (see noise section). Flagged samples were excluded from the analysis. The Galaxy region was masked out in the destriping phase, to reduce errors arising from strong signal gradients. Radiometers were combined according to the horn-uniform weighting scheme to minimize systematics. The polarization component was included in the analysis and is part of this release. A detailed description of the mapmaking procedure is given in Planck-2013-II[1], Planck-2015-A02[2], Planck-2015-A06[3] no relevant changes has been applyed in the 2018 release. See also the section on Frequency Maps.

The maps are in HEALPix format, at resolution Nside=1024 for all frequencies with an additional map at Nside=2048 for the LFI 70GHz channel, in the nested pixelization scheme. Unobserved pixels are marked by a special value.

The released maps are in Galactic coordinates. The conversion between ecliptic and Galactic coordinates is described by the rotation matrix [math] \newcommand{\Re}{\mathrm{Re}\,} \newcommand{\pFq}[5]{{}_{#1}\mathrm{F}_{#2} \left( \genfrac{}{}{0pt}{}{#3}{#4} \bigg| {#5} \right)} [/math]

[math] \begin{align} \label{def:Rot_matrix} \left (\begin{matrix} -0.054882486 & 0.494116468 & -0.867661702\\ -0.993821033 & -0.110993846 & -0.000346354\\ -0.096476249 & 0.86228144 & 0.497154957 \end{matrix} \right). \end{align} [/math]

The conversion was applied to the input pointing data, prior to the construction of the map.

Low-resolution maps and noise covariance matrices[edit]

To fully exploit the information contained in the large-scale structure of the microwave sky, pixel-pixel covariances are needed in the maximum likelihood estimation of the CMB power spectrum. However, full covariance matrices are impossible to employ at the native map resolution due to resource limitations. A low-resolution data set is therefore required for the low-ℓ analysis. This data set has been packed into three different files, one per frequency, called "LFI_NoiseCovMat_0??_0016_R3.00.tgz", that can be downloaded from the Cosmology section of the Planck Legacy Archive. They consist of low-resolution maps, and descriptions of residual noise present in those maps given by pixel-pixel noise covariance matrices (NCVMs). Note the in the 2018 release the Low-resolution maps full mission coverage, excluding Surveys 2 and 4 has not been used. We release them for crosschecking purposes with respect 2015 release.

Contents of LFI_NoiseCovMat_0??_0016_R3.20.tgz
Filenames Comment
LFI_SkyMap_0??_0016_coswin_DX12_full_regnoise.fits Low-resolution maps. Full mission coverage.
LFI_SkyMap_0??_0016_coswin_DX12_Corrected_full_regnoise.fits Low-resolution maps BandPass Corrected. Full mission coverage. Corrected for the BandPass.
LFI_SkyMap_0??_0016_coswin_DX12_s1-s3-s5-s6-s7-s8_regnoise.fits Low-resolution maps. Full mission coverage, excluding Surveys 2 and 4.
LFI_SkyMap_0??_0016_coswin_DX12_Corrected_s1-s3-s5-s6-s7-s8_regnoise.fits Low-resolution maps. Full mission coverage, excluding Surveys 2 and 4. Corrected for the BandPass.
offset_covmat_toast_nside64T16_coswin_025sec_??GHz_DX12_full_bin.dat Low-resolution noise covariance matrices. Regularization noise added. Format: C unformatted.
offset_covmat_toast_nside64T16_coswin_025sec_??GHz_DX12_s1-s3-s5-s6-s7-s8_bin.dat Low-resolution noise covariance matrices. Regularization noise added. Format: C unformatted.
covmat_toast_nside64T16_coswin_025sec_??GHz_DX12_full_bin.dat Low-resolution noise covariance matrices. No regularization noise added. Format: C unformatted.
covmat_toast_nside64T16_coswin_025sec_??GHz_DX12_s1-s3-s5-s6-s7-s8_bin.dat Low-resolution noise covariance matrices. No regularization noise added. Format: C unformatted.
mask_comm17Tv2_pol_R2.2x_ns16_nest.fits Low-resolution Mask 2.2.
mask_comm17Tv2_pol_R1.8x_ns16_nest.fits Low-resolution Mask 1.8.
offset_RESCALED_base_bestfit_cov_IQU.da Rescaled noise ONLY for 70 GHz. C unformatted.


The low-resolution data set can currently be utilized efficiently only at resolution Nside = 16, or lower. All the low-resolution data products are produced at this target resolution.

Low-resolution maps[edit]

A number of different schemes to obtain the low-resolution maps are discussed in [4]. We chose to downgrade the maps using the inverse noise weighting, no changes on the procedure has been applyed to the 2018 release. This is discussed further in Planck-2013-II[1]Planck-2015-A06[3].

Inputs[edit]

We took the high-resolution maps described in Map-making and Frequency Maps, and the corresponding 3×3 matrices as an input for this analysis step.

Production[edit]

The high-resolution maps were downgraded to Nside = 16 using inverse noise weights (given by the 3×3 matrices), and subsequently the temperature part was smoothed with a symmetric Gaussian beam with FWHM = 440arcmin.

Noise covariance matrices[edit]

The statistical description of the residual noise in the maps is given in the form of a pixel-to-pixel noise covariance matrix (NCVM), as described in [4].

Inputs[edit]

The noise model was determined by three noise parameters: the white noise level σ; slope; and knee frequency fknee. We actually used three sets of noise parameters, one for the entire mission (noise parameters are listed in Table 1), and one for each sky survey (SS1 and SS2).

We used the same pointing as in the noise Monte Carlo simulations. See the description in Noise Monte Carlo Simulation Inputs.

We used the gap files produced during the making of the flight maps to leave out samples that were flagged as bad for various reasons.

Production[edit]

The output of the NCVM module of MADAM mapmaker are inverse NCVMs. Since the inverse matrices are additive, we divided the computations into a number of small chunks to save computational resources. We first calculated one inverse NCVM per radiometer per survey at resolution Nside=32, and then combined these individual inverse matrices to form the actual inverse matrices. The mapmaking parameters were almost identical to the standard mapmaking runs. The differing parameter values are listed below:

  • baseline lengths were 0.25s for 30GHzand 1.0s for 44GHz, and 70GHz;
  • the calculations were performed at resolution Nside = 64;
  • no destriping mask was applied;
  • the horns were weighted optimally.

To obtain the noise covariance from its inverse, the matrices are inverted using the eigen decomposition of a matrix. The monopole of the temperature map cannot be resolved by the mapmaker, and thus the matrix becomes singular. This ill-determined mode is left out of the analysis.

Having calculated the eigen decomposition in the previous step, we can apply the same linear operators to modify the eigenvectors as were applied to the high-resolution maps while downgrading them. The eigenvectors are downgraded to Nside = 16 using inverse noise weights, and subsequently the temperature part is smoothed with a symmetric Gaussian beam with FWHM = 440arcmin.

The final matrices are then recomposed from the original eigenvalues and modified eigenvectors.

The low-resolution noise covariance matrices:

  • are C binary format files;
  • are organized in block form,

[math] \newcommand{\Re}{\mathrm{Re}\,} \newcommand{\pFq}[5]{{}_{#1}\mathrm{F}_{#2} \left( \genfrac{}{}{0pt}{}{#3}{#4} \bigg| {#5} \right)} [/math] [math] \begin{align} \label{def:Block_form} \left (\begin{matrix} II & IQ & IU \\ QI & QQ & QU \\ UI & UQ & UU \end{matrix} \right); \end{align} [/math]

  • are in the HEALPix nested pixelisation scheme (with resolution is Nside = 16, and thus there are Npix = 3072 pixels);
  • are in Galactic coordinates;
  • have KCMB units.

Half-ring jackknife noise maps[edit]

Overview[edit]

In the 2018 release we follow the same procedure as in Planck-2015-A03[2] in order to estimate the noise directly at the map level and in the angular power spectra.

Briefly, instead of using the full time-ordered data as described above, we produced two sets of maps using either only the first half of each pointing period (map named j1 below) or only the second half of each pointing period (map named j2). At each pixel p, these half-ring jackknife maps j1 and j2 contain the same sky signal, since they result from the same scanning pattern on the sky. However, because of instrumental noise, the maps j1 and j2 are not identical.

We estimated the noise level in each map m made using the full (ring) data, by constructing a half-ring difference map

[math]\mathbf{n_{m}}(p) = [ \mathbf{j_1}(p) - \mathbf{j_2}(p)] \ / \ \mathbf{w_{\rm hit}}(p)\,,[/math]

with weights

[math]\mathbf{w_{hit}}(p) = \sqrt{ \mathbf{hit_{full}}(p) \left[ \frac{1}{\mathbf{hit_1}(p)} + \frac{1}{\mathbf{hit_2}(p)} \right]}\,[/math].

Here hitfull(p) = hit1(p) + hit2(p) is the hit count at pixel p in the full map m, while hit1 and hit2 are the hit counts of j1 and j2, respectively. The weight factor whit(p) is equal to 2 only in those pixels where hit1(p) = hit2(p) . In a typical pixel, hit1(p) will differ slightly from hit2(p) and hence the weight factor is whit(p)>2.

The half-ring difference maps nm are the most direct measure of the noise in the actual maps. The other noise estimates (NCVM and noise Monte Carlo) rely on specific modelling of the noise and this modelling can be validated by comparing to the half-ring difference maps. However, the half-ring difference maps can only capture the noise that varies faster than half of the duration of the pointing period, i.e., the noise whose frequency is approximately f > 1/20min = 0.85mHz.

We calculated the noise maps nm, from half-ring "jackknife" maps for temperature (I) and polarization (Q and U) and as a first quality check of the maps (and as one of the tests of the whole data processing pipeline up to the map level) tested both numerically and visually that these noise maps divided pixel-by-pixel by the square root of the white noise covariance maps were approximately Gaussian with variance near to unity. Temperature noise maps for the nominal survey and for the first and second sky surveys are shown in the next subsection. Furthermore we calculated from the noise maps the temperature and polarization (E and B mode) auto-correlation and cross-correlation noise angular power spectra using HEALPix anafast and compared to these the results from the white noise covariance matrices and from the noise Monte Carlo simulations. A similar comparison was made between downgraded half-ring noise maps, downgraded noise Monte Carlo maps, and the low-resolution noise covariance maps. Detailed results are presented in Planck-2020-A2[5].

Comparison of noise estimates using Half-Ring[edit]

Here we compare noise angular power spectra estimated from half-ring difference maps (red), white noise covariance maps (black dash-dotted lines), and 100 full noise Monte Carlo simulations (grey band showing range for 16th and 84th quantiles of noise simulations, and the black solid lines giving the median, i.e., 50th quantile, of distributions). See the next section for details of noise Monte Carlo simulations. From top to botto we show TT, EE and BB power spectra for 30 GHz (left), 44 GHz (centre), and 70 GHz (right). Half-ring spectra are binned with [math]\mathbf{\Delta}l = 25 [/math] for [math] l \mathbf{\geq} 75 [/math].

Consistency Check

Below the null-tests comparing power spectra from survey differences to those from teh half-ring maps are showed. Difefrences are: left, Survey 1 - Survey 2; Middle , Survey 1- Survey 3; and right, Survey 1 - Survey4. these are for 30 GHz (top), 44 GHz (middle), and 70 GHz (bottom), for both TT and EE power spectra. There is a significant improvement in Surve1 - Survey 2 and Survey 1 - Survey 4 at 30 GHz, especially in EE. See Planck-2020-A2[5] for further details.

Null Test


High-ell average noise relative to white noise estimate[edit]

The figure below is the same as the previous figures, but here the noise comparison is made from the high ℓ tails of the angular power spectra, where the white noise dominates. We have taken the average of C from multipoles between 1150 and 1800 for both temperature and polarization and tehn comparing with the WNCVM. As already shown in previous releases, there is still an excess of 1/f noise, meaning tha both the real data and the noise MCs predict slightly larger noise than the WNCVM. It is important to note that such noise excess is reduced considerably with respect to the 2015 release.

Ratio at high multipoles

Noise Monte Carlo simulations[edit]

Overview[edit]

Calculating and handling full pixel-to-pixel noise covariance matrices for Planck maps if feasible only at low resolution. To support the analysis of high-resolution maps, a Monte Carlo set of noise maps were produced. These maps were produced from noise timelines using the same map-making procedure as for the flight data. In the noise Monte Carlo it was possible to follow exactly the mapmaking procedure used for the flight maps, whereas for the calculation of noise covariance matrices some approximations had to be made. Such noise Monte Carlos were produced at two levels of the analysis: (1) LFI Monte Carlo (MC) as part of the LFI data processing procedure; and (2) Full Focal Plane (FFP) Monte Carlos as part of the joint HFI/LFI data processing. This page describes the LFI noise MCs. For the FFP MC, see HL-sims and Simulation data.

Inputs[edit]

The noise MC uses a three-parameter noise model, consisting of white noise level (σ), slope, and knee frequency (fknee)). Here the noise consists of white noise and correlated 1/f noise, with a power spectrum

[math] P(f) = \frac{2\sigma^2}{f_\mathrm{sample}}\left(\frac{f}{f_\mathrm{knee}}\right)^\mathrm{slope} [/math],

where fsample is the sampling frequency of the instrument. The noise parameters were determined separately for each radiometer, as described in the section Noise above, assuming they stayed constant over the mission.

The detector pointing was reconstructed from satellite pointing information, focal-plane geometry, pointing correction (tilt angle), and sample timing, using Level-S simulation software. The same pointing solution (two focal planes) was used as for the LFI flight maps. Due to numerical accuracy, the detector pointing in the noise MC was not exactly the same as for the flight maps, so some data samples (of the order of one in a thousand) whose pointing was near the pixel boundary ended up assigned to the neighbouring pixel. During mapmaking from the flight data, a "gap file" was produced to represent the samples that were omitted from mapmaking due to various flags. This gap file was used in the noise MC instead of the full set of flags. The flight mapmaking procedure used destriping masks to prevent regions of strong signal gradients from contributing to the noise baseline solution. These same destriping masks (one for each frequency channel) were used for the noise MC.

Production[edit]

The noise was generated internally in the Madam mapmaking code using a stochastic differential equation (SDE) method, to avoid time-consuming writing and reading of noise timelines to and from disk. Noise for each pointing period was generated separately, using a double-precision random number seed constructed from the realization number, radiometer number, and the pointing period number; this allowed for regeneration of the same noise realization when needed. White noise and 1/f noise were generated separately.

The same mapmaking code (Madam) with the same parameter settings was used for the noise MC as for the flight maps. In addition to the destriped maps from the full noise (output maps), binned maps from just the white noise (binned white noise maps) were produced; they represent the white noise part of the output maps. The difference between these two maps represents the residual correlated noise in the output map. The maps were made at HEALPix resolution Nside = 1024 for all LFI frequency channels and also at HEALPix resolution Nside = 2048 for the 70 GHz channel. For low-resolution analysis, these maps were downgraded (and the temperature part was smoothed) to Nside = 32 and Nside = 16.

In addition to frequency maps for the nominal survey, also single-survey and 70 GHz horn-pair maps were produced in the noise MC. For each case 102-1026 realizations were produced.

Usage[edit]

These noise Monte Carlo maps were used for a number of purposes in LFI data analysis. They were compared to the low-resolution noise covariance matrices, generated for the same noise model, in order to reveal the impact of the approximations in the noise covariance matrix calculation. They were compared to the half-ring noise maps to see how well the noise model matches the noise in the flight maps (noting, however, that the half-ring noise maps misrepresent the lowest noise frequencies in the flight maps, and contain some effects from the sky signal). They were also used in power spectrum estimation and non-Gaussianity estimation.

Examples[edit]

As an example, we show below images of the first realization of the 70GHz frequency map noise for the nominal survey. The images are in the order: destriped full noise; binned white noise; and residual correlated noise. Note that it is difficult to see any difference between the first two images, since the residual correlated noise is more than an order of magnitude below the white noise level. The units here are KCMB.


Destriped full noise.
Binned white noise.
Residual correlated noise.


The following two images show the statistics of the angular power spectra of 101 realizations of the 70 GHz frequency map noise for the nominal survey. The thick black line shows the median C, while the green line the mean C. Thin black lines show the minimum, 16th percentile, 84th percentile, and the maximum C. The red line is the 102nd realization. The first plot is for the full noise in the output map, while the second plot is for the residual correlated noise.

Full noise.
Residual correlated noise.

References[edit]

  1. 1.01.1 Planck 2013 results. II. Low Frequency Instrument data processing, Planck Collaboration, 2014, A&A, 571, A2.
  2. 2.02.1 Planck 2015 results. II. LFI processing, Planck Collaboration, 2016, A&A, 594, A2.
  3. 3.03.1 Planck 2015 results. VI. LFI mapmaking, Planck Collaboration, 2016, A&A, 594, A6.
  4. 4.04.1 Residual noise covariance for Planck low-resolution data analysis, R. Keskitalo, M. A. J. Ashdown, P. Cabella, T. Kisner, T. Poutanen, R. Stompor, ArXiv e-prints, (2013).
  5. 5.05.1 Planck 2018 results. II. Low Frequency Instrument data processing, Planck Collaboration, 2020, A&A, 641, A2.

(Planck) Low Frequency Instrument

(Hierarchical Equal Area isoLatitude Pixelation of a sphere, <ref name="Template:Gorski2005">HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere, K. M. Górski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke, M. Bartelmann, ApJ, 622, 759-771, (2005).

Cosmic Microwave background

Full-Width-at-Half-Maximum

(Planck) High Frequency Instrument