Difference between revisions of "Specially processed maps"

From Planck Legacy Archive Wiki
Jump to: navigation, search
(→‎Description)
 
(70 intermediate revisions by 10 users not shown)
Line 1: Line 1:
−
{{DISPLAYTITLE:Additional maps}}
+
{{DISPLAYTITLE: Additional maps}}
−
==Overview==
+
==Introduction==
  
−
This section describes products that require special processing. Only one such product is available at this time; this section will be expanded with time as more products are added.
+
This section describes the map-based products that required special processing.
  
−
== Lensing map ==
+
== 2015 lensing map ==
−
=== Description ===
 
  
−
Here we present the minimum-variance (MV) lens reconstruction which forms the basis for the main results of {{PlanckPapers|planck2013-p12}}. This map is produced using a combination of the 143 and 217 GHz Planck maps on approximately 70% of the sky, and is the same map on which the Planck lensing likelihood is based.
+
We distribute the minimum-variance (MV) lensing potential estimate presented in {{PlanckPapers|planck2014-a17}} as part of the 2014 data release. This map represents an estimate of the CMB lensing potential on approximately 70% of the sky, and also forms the basis for the Planck 2014 lensing likelihood. It is produced using filtered temperature and polarization data from the SMICA DX11 CMB map; its construction is discussed in detail in {{PlanckPapers|planck2014-a11}}.
  
−
We distribute:
 
  
−
; PHIBAR : A (transfer-function convolved) map of the lensing potential, in NSIDE 2048 HEALPix RING format. It is obtained by convolving the lensing potential estimate <math>\hat{\phi}</math> with the lensing response function <math>R_L^{\phi\phi}</math>. This map has been band-limited between multipoles <math>10 \le L \le 2048</math>.
+
The estimate is contained in a single gzipped tarball named  ''{{PLASingleFile|fileType=map|name=COM_CompMap_Lensing_2048_R2.00.tgz|link=COM_CompMap_Lensing_2048_R2.00.tgz}}''. Its contents are described below. The convergence map "dat_klm.fits" that can be found in the tarball, has been categorized as COM_Lensing-Convergence-dat-klm_2048_R2.00.fits in the Lensing Products section of the archive.
−
; MASK : This is a NSIDE = 2048 HEALPix map, containing the analysis mask used in the lens reconstruction. ''Note'': the lensing map PHIBAR may take small but non-zero values inside the masked regions because it has been bandlimited.
 
−
; RLPP : This column contains the response function <math>R_L^{\phi\phi}</math>.
 
−
; NLPP : This column contains a sky-averaged estimate of the noise power spectrum of PHIBAR, <math>N_L^{\phi\phi}</math>. The noise is highly coloured. There is a dependence of the noise power spectrum on the local noise level of the map, discussed in Appendix A of {{PlanckPapers|planck2013-p12}}. Note that the noise power spectrum estimate here is not sufficiently accurate for a power spectrum analysis.
 
  
−
Also, the table below gives the lensing curl-mode power spectrum data used to produce Figure A2 of {{PlanckPapers|planck2013-p12}}:
 
  
−
{| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:center"
+
 
−
|+ ''''MV' curl reconstruction bandpowers from Fig A2 of {{PlanckPapers|planck2013-p12}} '''
+
{| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:left"
 +
|+ ''' Contents of Lensing package '''
 
|-  bgcolor="ffdead"   
 
|-  bgcolor="ffdead"   
−
! <math>L</math><sub>min</sub> ||  <math>L</math><sub>max</sub> ||  <math>{\hat{C}_L}^{\psi\psi} (L(L+1))^2/(2\pi)</math> || <math> std({\hat{C}_L}^{\psi\psi})  (L(L+1))^2 /(2\pi)</math>
+
! Filename ||  Format ||  Description
 
|-
 
|-
−
|     2 ||    7 || -13.6379 || 15.3409
+
| dat_klm.fits || HEALPix FITS format alm, with <math>L_{\rm max}=2048</math> || Contains the estimated lensing convergence <math> \hat{\kappa}_{LM} = \frac{1}{2} L(L+1)\hat{\phi}_{LM} </math>.
 
|-
 
|-
−
|     8 ||  20 ||   6.0184 || 4.8881
+
| mask.fits.gz || HEALPix FITS format map, with <math>N_{\rm side}=2048</math> || Contains the lens reconstruction analysis mask.
 
|-
 
|-
−
|   21 ||   39 ||  -1.0675 ||  3.0940
+
| nlkk.dat || ASCII text file, with columns = (<math>L</math>, <math>N_L </math>, <math>C_L+N_L</math>) || The approximate noise <math>N_L</math> (and signal+noise, <math>C_L+N_L</math>) power spectrum of <math> \hat{\kappa}_{LM} </math>, for the fiducial cosmology used in {{PlanckPapers|planck2014-a15}}.
 +
|}
 +
 
 +
== 2015 Compton <i>y</i> parameter map ==
 +
 
 +
We distribute here the Planck full mission Compton parameter maps (<i>y</i>-maps hereafter) obtained using the NILC and MILCA component-separation algorithms as described in {{PlanckPapers|planck2014-a28}}. We also provide the ILC weights per scale and per frequency that were used to produce these <i>y</i>-maps. IDL routines are also provided to allow the user to apply those weights. Compton parameters produced by keeping either the first or the second half of stable pointing periods are also provided; we call these the FIRST and LAST <i>y</i>-maps. Additionally we construct noise estimates of full mission Planck <i>y</i>-maps from the half difference of the FIRST and LAST <i>y</i>-maps. These estimates are used to construct standard deviation maps of the noise in the full mission Planck <i>y</i>-maps, which are also provided. To complement this we also provide the power spectra of the noise estimate maps after correcting for inhomogeneities using the standard deviation maps. We also deliver foreground masks including point-source and Galactic masks.
 +
 
 +
The full data set is contained in a single gzipped tarball named ''COM_CompMap_YSZ_R2.00.fits.tgz''. Its contents are described below.
 +
''Temporarily the tarball file can be found in the Supplement Data area of the PLA (go to pla.esac.esa.int/pla, go to Subsection Supplementary Data and type "YSZ" in the search form)''
 +
 
 +
 
 +
{| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:left"
 +
|+ ''' Contents of COM_CompMap_YSZ_R2.00.fits.tgz} '''
 +
|- bgcolor="ffdead" 
 +
! Filename ||  Format ||  Description
 
|-
 
|-
−
|   40 ||   65 ||  0.6135 || 1.8474
+
| nilc_ymaps.fits || HEALPix FITS format map in Galactic coordinates with <math>N_{\rm side}=2048 </math>|| Contains the NILC full mission, FIRST and LAST <i>y</i>-maps.
 
|-
 
|-
−
|   66 ||  100 ||   1.5030 || 1.2696
+
| milca_ymaps.fits || HEALPix FITS format map in Galactic coordinates with <math> N_{\rm side} = 2048 </math> || Contains the MILCA full mission, FIRST and LAST <i>y</i>-maps.
 
|-
 
|-
−
|   101 ||  144 ||   1.3760  ||  0.9950
+
| nilc_weights_BAND.fits || HEALPix FITS format map in Galactic coordinates with <math> N_{\rm side} = 128 </math>|| Contains the NILC ILC weights for the full mission <i>y</i>-map for band BAND 0 to 9. For each band we provide a weight map per frequency.
 
|-
 
|-
−
|   145 ||  198 || -1.2289  ||  0.8286
+
| milca_FREQ_Csz.fits || HEALPix FITS format map in Galactic coordinates with <math> N_{\rm side} = 2048 </math> || Contains the MILCA ILC weights for the full mission <i>y</i>-map for frequency FREQ (100, 143, 217, 353, 545, 857). For each frequency we provide a weight map per filter band.
 
|-
 
|-
−
|   199 ||  263 ||   1.1910 || 0.7001
+
| nilc_stddev.fits || HEALPix FITS format map in Galactic coordinates with <math> N_{\rm side} = 2048 </math>|| Contains the stddev map for the NILC full mission <i>y</i>-map.
 
|-
 
|-
−
|   264 ||  338 || -0.6567  ||  0.6197
+
| milca_stddev.fits || HEALPix FITS format map in Galactic coordinates with <math> N_{\rm side} = 2048 </math> || Contains the stddev maps for the MILCA full mission <i>y</i>-map.
 
|-
 
|-
−
|   339 || 425 || -0.8201  ||  0.5235
+
| nilc_homnoise_spect.fits || ASCII table FITS format || Contains the angular power spectrum of the homogeneous noise in the NILC full mission <i>y</i>-map.
 
|-
 
|-
−
|   426 || 525 || -0.7581  ||  0.4850
+
| milca_homnoise_spect.fits || ASCII table FITS format || Contains the angular power spectrum of the homogeneous noise in the MILCA full mission <i>y</i>-map.
 
|-
 
|-
−
|   526 ||  637 || -0.3201  ||  0.5134
+
| masks.fits || HEALPix FITS format map, with <math> N_{\rm side} = 2048 </math> || Contains foreground masks.
 
|-
 
|-
−
|   638 ||  762 ||  -0.1589  ||  0.4073
+
| nilc_bands.fits || ASCII table FITS format || Contains NILC wavelet bands in multipole space
−
|-
 
−
|   763 || 901 ||  -0.6451  ||  0.4044
 
−
|-
 
−
|  902 || 1054 ||  0.4910  ||  0.3718
 
−
|-
 
−
|  1055 || 1221 ||  -0.2186  ||  0.3702
 
−
|-
 
−
|  1222 || 1404 ||  -0.3295  ||  0.4146
 
−
|-
 
−
|  1405 || 1602 ||  -0.3647  ||  0.4703
 
−
|-
 
−
|  1603 || 1816 ||  -0.1060  ||  0.5904
 
−
|-
 
−
|  1817 || 2020 ||  -0.7887  ||  0.8507
 
 
|}
 
|}
  
−
=== Production process ===
+
== 2015 lensing-induced <i>B</i>-mode map ==
 +
 
 +
We distribute the lensing-induced <i>B</i>-mode map presented in {{PlanckPapers|planck2015-XLI}}. The lensing <i>B</i>-mode Stokes parameter maps are produced by combining the lensing potential reconstruction from the SMICA CMB temperature map with <i>E</i>-mode data from the SMICA CMB polarization maps. The SMICA temperature and polarization products are described in {{PlanckPapers|planck2014-a11}}. The lensing-induced <i>B</i>-mode polarization maps are used in cross-correlation with the SMICA CMB polarization maps to obtain a lensing <i>B</i>-mode power spectrum measurement from approximately 70% of the sky, as described in {{PlanckPapers|planck2015-XLI}}.
 +
 
 +
We provide both raw products, which can be utilized to generate products adapted to one's specific needs in term of mask, filtering, etc., and "ready-to-use" products for cross-correlation study purposes.
 +
 
 +
=== Raw products ===
 +
 
 +
We deliver the non-normalized lensing-induced Stokes parameter maps, labelled <math> \bar{Q}^{\rm{lens}} </math> and <math> \bar{U}^{\rm{lens}} </math>, which form the basis of the final lensing <i>B</i>-mode estimator defined in equation (6) of {{PlanckPapers|planck2015-XLI}}. They are defined as
 +
 
 +
<math> \begin{eqnarray}
 +
\bar Q^{\rm{lens}}({\bf n}) &=& \widetilde Q^{E}({\bf n}) \cdot \nabla \widetilde \phi({\bf n}), \\
 +
\bar U^{\rm{lens}}({\bf n}) &=& \widetilde U^{E}({\bf n}) \cdot \nabla \widetilde \phi({\bf n}),
 +
\end{eqnarray} </math>
 +
 
 +
where <math> \widetilde Q^{E} </math> and  <math> \widetilde U^{E} </math> are the filtered pure <i>E</i>-mode polarization maps given in equation (5) of {{PlanckPapers|planck2015-XLI}}, and <math> \widetilde \phi</math> is the filtered lensing potential estimate. 
 +
 
 +
We also provide the normalization transfer function <math> \mathcal{B}_\ell </math> defined in equation (11) of {{PlanckPapers|planck2015-XLI}}, as well as the "L70" mask <math> M({\bf n}) </math> that retains 69% of the sky before
 +
apodization, and its apodized version <math> \tilde{M}({\bf n}) </math>, which has an effective sky fraction <math> f_{\rm{sky}}^{\rm{eff}} = 65\% </math>.
 +
 
 +
As an example of the utilization of these products, the lensing <i>B</i>-mode maps that are shown in figure 4 of {{PlanckPapers|planck2015-XLI}} are generated from
 +
 +
<math> Q^{\rm{lens}} \pm i U^{\rm{lens}} = \sum_{\ell m}  \left( G_\ell \mathcal{B}_\ell^{-1}  \int d{\bf n} {\, }_{\pm 2}Y_{\ell m}^*({\bf n}) \left(\bar{Q}^{\rm{lens}} \pm i \bar{U}^{\rm{lens}} \right)  \right)  {\, }_{\pm 2}Y_{\ell m}({\bf n}) </math>,
 +
 
 +
where <math>G_\ell</math> is a Gaussian of 60 arcmin FWHM (introduced for highlighting large angular scales, although it can be removed or replaced by any other filter). This can be practically done by ingesting <math>\bar{Q}^{\rm{lens}} </math> and <math> \bar{U}^{\rm{lens}} </math> in the HEALPix "smoothing" routine, and using the product <math> G_\ell\mathcal{B}_\ell^{-1} </math> as an input filtering function.
  
−
The construction PHIBAR, RLPP and NLPP are described in detail in Sec. 2.1 of {{BibCite|planck2013-p12}}. The response function <math>R_L^{\phi\phi}</math> here is analogous to the the beam transfer function in a CMB temperature or polarization map. We have chosen to distribute this transfer-function convolved map rather than the normalized lens reconstruction as it is a significantly more localized function of the CMB temperature map from which it is derived, and therefore more useful for cross-correlation studies.
+
=== Specific products ===
  
−
===Inputs===
+
We provide the lensing <i>B</i>-mode spherical harmonic coefficient estimate <math> B_{\ell m}^{\rm{lens}} </math> over approximately 70% of the sky.
  
−
This product is built from the 143 and 217 GHz Planck [[Frequency Maps|frequency maps]], with 857GHz projected out as a dust template.
+
It can also be constructed using the raw products described above from  
−
The analysis mask is constructed from a combination of thresholding in the 857GHz map (to remove the regions which are most contaminated by Galactic dust) and the [[CMB_and_astrophysical_component_maps#CO_emission_maps | Type2 CO map]] (to reduce contamination from CO lines at 217GHz). This is joined with a compact object mask synthesized from several Planck source [[Catalogues | catalogues]], including the [[ Catalogues#ERCSC | ERCSC]], [[ Catalogues#SZ | SZ ]] and [[ Catalogues#The Catalogue of Compact Sources | PCCS ]]. The reconstruction was performed using the fiducial beam window functions B(l) from the [[ The RIMO | HFI RIMO ]]. Details of the procedure used to produce a lensing estimate from these inputs are given in {{BibCite|planck2013-p12}}.
 
  
−
===File names and format===
+
<math> B_{\ell m}^{\rm{lens}} = f_{10 \rightarrow 2000} \,  \mathcal{B}_\ell^{-1} \, \, {\, }_{\pm 2}\mathcal{Y} \left[ \tilde{M}({\bf n}) \left( \bar{Q}^{\rm{lens}}({\bf n}) \pm i \bar{U}^{\rm{lens}}({\bf n}) \right) \right] </math>,
  
−
A single file named
+
where <math> f_{10 \rightarrow 2000} </math> is a function for producing band-powers over the range <math> 10 \le \ell \le 2000 </math>, and <math> {\, }_{\pm2}\mathcal{Y} </math> is a short-hand notation for transforming a map into spin-weighted spherical harmonic coefficients <math> {\, }_{+2}a_{\ell m}</math>, <math>{\, }_{-2}a_{\ell m} </math> and forming <math>1/(2i)\left({\, }_{+2}a_{\ell m} - {\, }_{-2}a_{\ell m}\right)</math>. This can be done using, e.g., the HEALPix "anafast" tool.
−
*''{{PLASingleFile|fileType=map|name=COM_CompMap_Lensing_2048_R1.10.fits|link=COM_CompMap_Lensing_2048_R1.10.fits}}''
+
 
−
with two BINTABLE extensions containing the items described below.
+
The lensing <i>B</i>-mode power spectrum estimate <math> \hat{C}_\ell^{BB^{\rm{lens}}} </math> discussed in {{PlanckPapers|planck2015-XLI}} is obtained by forming the cross-correlation power spectrum of <math> B_{\ell m}^{\rm{lens}} </math> and the <i>B</i>-mode data from the SMICA polarization maps <math> B_{\ell m} </math>:
 
   
 
   
−
For illustration, we show in the figures below the maps of the Wiener-filtered CMB lensing potential in Galactic coordinates using orthographic projection. The reconstruction was bandpass filtered to <math>L \in [10, 2048]</math>. Note that the lensing reconstruction, while highly statistically significant, is still noise dominated for every individual mode, and is at best <math>S/N \simeq 0.7</math> around <math>L = 30</math>.
+
<math> \hat{C}_\ell^{BB^{\rm{lens}}} = \frac{\left(f_{\rm{sky}}^{\rm{eff}}\right)^{-1}}{2 \ell +1} G_\ell^{-2} \sum_m B_{\ell m}^* B_{\ell m}^{\rm{lens}}</math>,
  
−
<center>
+
where <math> G_\ell </math> is the 5 arcmin Gaussian beam that convolves the SMICA CMB maps.  
−
<gallery perrow=3 widths=260px heights=170px>  
 
−
File: analysis_lens_pub_map_orth_north_dat_p.png | Galactic north
 
−
File: analysis_lens_pub_map_orth_south_dat_p.png | Galactic south
 
−
</gallery></center>
 
  
 +
 +
The products are contained in a single gzipped tarball named  ''{{PLASingleFile|fileType=map|name=COM_Lensing-Bmode_R2.00.tgz|link=COM_Lensing-Bmode_R2.00.tgz}}''. Its contents are described below.
  
  
 
{| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:left"
 
{| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:left"
−
|+ '''FITS file structure'''
+
|+ ''' Contents of Lensing B-mode package '''
−
|- bgcolor="ffdead" 
 
−
!colspan="4" | 1. EXTNAME = ''LENS-MAP''
 
 
|-  bgcolor="ffdead"   
 
|-  bgcolor="ffdead"   
−
! Column Name || Data Type || Units || Description
+
! Filename || Format || Description
 
|-
 
|-
−
|PHIBAR || Real*4 || none || Map of the lensing potential estimate, convolved with RLPP
+
| bar_q_lens_map.fits || HEALPix FITS format map in Galactic coordinates with <math> N_{\rm  side} = 2048 </math>  || Contains the non-normalized lensing-induced <i>Q</i> Stokes parameter map <math> \bar Q^{\rm{lens}}({\bf n}) </math>.
 
|-
 
|-
−
|MASK || Int || none || Region over which the lensing potential is reconstructed
+
| bar_u_lens_map.fits || HEALPix FITS format map in Galactic coordinates with <math> N_{\rm  side} = 2048 </math>  || Contains the non-normalized lensing-induced <i>U</i> Stokes parameter map <math> \bar U^{\rm{lens}}({\bf n}) </math>.
−
|- bgcolor="ffdead" 
 
−
! Keyword || Data Type || Value || Description
 
 
|-
 
|-
−
|PIXTYPE  ||  string || HEALPIX || colspan="2"|
+
| mask.fits || HEALPix FITS format map in Galactic coordinates with <math> N_{\rm side} = 2048 </math> || The L70 mask.
 
|-
 
|-
−
|COORDSYS  ||  string || GALACTIC || Coordinate system
+
| mask_noapo.fits || HEALPix FITS format map in Galactic coordinates with <math> N_{\rm side} = 2048 </math> || The L70 mask without apodization.
 
|-
 
|-
−
|ORDERING  || string || NESTED  || Healpix ordering
+
| transfer_function_b_l.dat || ASCII text file, with columns = (<math>\ell</math>, <math> \mathcal{B}_\ell </math>) || The transfer function.
 
|-
 
|-
−
|NSIDE  || Int*4 || 2048 ||colspan="2"| Healpix Nside
+
| lensing_bmode_b_lm.fits || HEALPix FITS format alm, with <math> \ell_{\rm max} = 2000 </math> || Contains the lensing <i>B</i>-mode harmonic coefficients <math> B_{\ell m}^{\rm{lens}} </math>.
 
|-
 
|-
−
|FIRSTPIX || Int*4 ||          0 ||
+
| lensing_bmode_bandpowers.dat || ASCII text file, with columns = (<math>\ell_{\rm min}</math>, <math>\ell_{\rm b} </math>, <math>\ell_{\rm max} </math>, <math> \hat{C}_{\ell_{\rm b}}^{BB^{\rm{lens}}} </math>, <math> \Delta \hat{C}_{\ell_{\rm b}}^{BB^{\rm{lens}}} </math> ) || The lensing <i>B</i>-mode bandpower estimate on approximativily 70% of the sky and over the multipole range from 10 to 2000 shown in figure 9 of {{PlanckPapers|planck2015-XLI}} (for plotting purposes only).
−
|-
 
−
|LASTPIX || Int*4 || 50331647 ||
 
−
|- bgcolor="ffdead"  
 
−
!colspan="4"| 2. EXTNAME = ''TransFun''
 
−
|- bgcolor="ffdead"  
 
−
! Column Name || Data Type || Units || Description
 
−
|-
 
−
|RLPP || Real*4 || none || Response function
 
−
|-
 
−
|NLPP || Real*4 || none || Sky-averaged noise power spectrum estimate
 
−
|- bgcolor="ffdead" 
 
−
! Keyword || Data Type || Value || Description
 
−
|-
 
−
|L_MIN || Int*4  ||          0 || First multipole
 
−
|-
 
−
|L_MAX || Int*4  ||  2048 || Last multipole
 
 
|}
 
|}
 +
 +
== 2015 Integrated Sachs-Wolfe effect map ==
 +
 +
We distribute estimates of the integrated Sachs-Wolfe (ISW) maps presented in {{PlanckPapers|planck2014-a26}} as part of the 2015 data release. These map represents an estimate of the ISW anisotropies using different data sets:
 +
 +
* SEVEM DX11 CMB map, together with all the large-scale structure tracers considered in the ISW paper, namely: NVSS, SDSS, WISE, and the Planck lensing map
 +
* Using only the large-scale structure tracers mentioned above
 +
* SEVEM DX11 CMB map, together with NVSS and the Planck lensing maps (since these two tracers capture most of the information, as compared to SDSS and WISE)
  
  
−
== IRAM Maps of the Crab nebula ==
+
For all the three cases, the reconstruction is provided on approximately 85% of the sky, and they are produced using the LCB filter described in the Planck ISW paper (Section 5), described in detail in [http://cdsads.u-strasbg.fr/abs/2008ISTSP...2..747B| Barreiro et al. 2008] and [http://cdsads.u-strasbg.fr/doi/10.1093/mnras/stw415| Bonavera et al. 2016].
  
−
Maps of the Crab nebula at  89.189 GHz (HCO+(1-0) transition) in both temperature and polarization, prodouced from observations performed at the IRAM 30m telescope from January 9th to January 12th 2009, are delivered as a tarball of 416 KB in the file
+
These ISW maps, together with their corresponding uncertainties maps and masks, are given in a file named  ''{{PLASingleFile|fileType=map|name=COM_CompMap_ISW_0064_R2.00.fits|link=COM_CompMap_ISW_0064_R2.00.fits}}''. Its contents are described below.
  
−
: [[File:Crab_IRAM_2010.zip]]
 
  
−
See README in the tarball for full details. These data were used in {{BibCite|aumont2010}}
+
{| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:left"
 +
|+ ''' Contents of the ISW maps file: COM_CompMap_ISW_0064_R2.00.fits '''
 +
|-  bgcolor="ffdead"
 +
! Extension ||  Format ||  Description || Used data sets
 +
|-
 +
| 0 || HEALPix FITS format map with three components, <math>N_{\rm  side}=64</math>, Ordering='Nest' || Contains three components: i) ISW map [Kelvin], ii) Error map [Kelvin], iii) Mask map ||  SEVEM DX11 CMB + NVSS + SDSS + WISE + Planck lensing.
 +
|-
 +
| 1 || HEALPix FITS format map with three components, <math>N_{\rm  side}=64</math>, Ordering='Nest' || Contains three components: i) ISW map [Kelvin], ii) Error map [Kelvin], iii) Mask map ||  NVSS + SDSS + WISE + Planck lensing.
 +
|-
 +
| 2 || HEALPix FITS format map with three components, <math>N_{\rm  side}=64</math>, Ordering='Nest' || Contains three components: i) ISW map [Kelvin], ii) Error map [Kelvin], iii) Mask map ||  SEVEM DX11 CMB + NVSS + Planck lensing.
 +
|}
  
 
==References==
 
==References==

Latest revision as of 06:56, 26 April 2016

Introduction[edit]

This section describes the map-based products that required special processing.

2015 lensing map[edit]

We distribute the minimum-variance (MV) lensing potential estimate presented in Planck-2015-A15[1] as part of the 2014 data release. This map represents an estimate of the CMB lensing potential on approximately 70% of the sky, and also forms the basis for the Planck 2014 lensing likelihood. It is produced using filtered temperature and polarization data from the SMICA DX11 CMB map; its construction is discussed in detail in Planck-2015-A09[2].


The estimate is contained in a single gzipped tarball named COM_CompMap_Lensing_2048_R2.00.tgz. Its contents are described below. The convergence map "dat_klm.fits" that can be found in the tarball, has been categorized as COM_Lensing-Convergence-dat-klm_2048_R2.00.fits in the Lensing Products section of the archive.


Contents of Lensing package
Filename Format Description
dat_klm.fits HEALPix FITS format alm, with [math]L_{\rm max}=2048[/math] Contains the estimated lensing convergence [math] \hat{\kappa}_{LM} = \frac{1}{2} L(L+1)\hat{\phi}_{LM} [/math].
mask.fits.gz HEALPix FITS format map, with [math]N_{\rm side}=2048[/math] Contains the lens reconstruction analysis mask.
nlkk.dat ASCII text file, with columns = ([math]L[/math], [math]N_L [/math], [math]C_L+N_L[/math]) The approximate noise [math]N_L[/math] (and signal+noise, [math]C_L+N_L[/math]) power spectrum of [math] \hat{\kappa}_{LM} [/math], for the fiducial cosmology used in Planck-2015-A13[3].

2015 Compton y parameter map[edit]

We distribute here the Planck full mission Compton parameter maps (y-maps hereafter) obtained using the NILC and MILCA component-separation algorithms as described in Planck-2015-A22[4]. We also provide the ILC weights per scale and per frequency that were used to produce these y-maps. IDL routines are also provided to allow the user to apply those weights. Compton parameters produced by keeping either the first or the second half of stable pointing periods are also provided; we call these the FIRST and LAST y-maps. Additionally we construct noise estimates of full mission Planck y-maps from the half difference of the FIRST and LAST y-maps. These estimates are used to construct standard deviation maps of the noise in the full mission Planck y-maps, which are also provided. To complement this we also provide the power spectra of the noise estimate maps after correcting for inhomogeneities using the standard deviation maps. We also deliver foreground masks including point-source and Galactic masks.

The full data set is contained in a single gzipped tarball named COM_CompMap_YSZ_R2.00.fits.tgz. Its contents are described below. Temporarily the tarball file can be found in the Supplement Data area of the PLA (go to pla.esac.esa.int/pla, go to Subsection Supplementary Data and type "YSZ" in the search form)


Contents of COM_CompMap_YSZ_R2.00.fits.tgz}
Filename Format Description
nilc_ymaps.fits HEALPix FITS format map in Galactic coordinates with [math]N_{\rm side}=2048 [/math] Contains the NILC full mission, FIRST and LAST y-maps.
milca_ymaps.fits HEALPix FITS format map in Galactic coordinates with [math] N_{\rm side} = 2048 [/math] Contains the MILCA full mission, FIRST and LAST y-maps.
nilc_weights_BAND.fits HEALPix FITS format map in Galactic coordinates with [math] N_{\rm side} = 128 [/math] Contains the NILC ILC weights for the full mission y-map for band BAND 0 to 9. For each band we provide a weight map per frequency.
milca_FREQ_Csz.fits HEALPix FITS format map in Galactic coordinates with [math] N_{\rm side} = 2048 [/math] Contains the MILCA ILC weights for the full mission y-map for frequency FREQ (100, 143, 217, 353, 545, 857). For each frequency we provide a weight map per filter band.
nilc_stddev.fits HEALPix FITS format map in Galactic coordinates with [math] N_{\rm side} = 2048 [/math] Contains the stddev map for the NILC full mission y-map.
milca_stddev.fits HEALPix FITS format map in Galactic coordinates with [math] N_{\rm side} = 2048 [/math] Contains the stddev maps for the MILCA full mission y-map.
nilc_homnoise_spect.fits ASCII table FITS format Contains the angular power spectrum of the homogeneous noise in the NILC full mission y-map.
milca_homnoise_spect.fits ASCII table FITS format Contains the angular power spectrum of the homogeneous noise in the MILCA full mission y-map.
masks.fits HEALPix FITS format map, with [math] N_{\rm side} = 2048 [/math] Contains foreground masks.
nilc_bands.fits ASCII table FITS format Contains NILC wavelet bands in multipole space

2015 lensing-induced B-mode map[edit]

We distribute the lensing-induced B-mode map presented in Planck-2015-XLI[5]. The lensing B-mode Stokes parameter maps are produced by combining the lensing potential reconstruction from the SMICA CMB temperature map with E-mode data from the SMICA CMB polarization maps. The SMICA temperature and polarization products are described in Planck-2015-A09[2]. The lensing-induced B-mode polarization maps are used in cross-correlation with the SMICA CMB polarization maps to obtain a lensing B-mode power spectrum measurement from approximately 70% of the sky, as described in Planck-2015-XLI[5].

We provide both raw products, which can be utilized to generate products adapted to one's specific needs in term of mask, filtering, etc., and "ready-to-use" products for cross-correlation study purposes.

Raw products[edit]

We deliver the non-normalized lensing-induced Stokes parameter maps, labelled [math] \bar{Q}^{\rm{lens}} [/math] and [math] \bar{U}^{\rm{lens}} [/math], which form the basis of the final lensing B-mode estimator defined in equation (6) of Planck-2015-XLI[5]. They are defined as

[math] \begin{eqnarray} \bar Q^{\rm{lens}}({\bf n}) &=& \widetilde Q^{E}({\bf n}) \cdot \nabla \widetilde \phi({\bf n}), \\ \bar U^{\rm{lens}}({\bf n}) &=& \widetilde U^{E}({\bf n}) \cdot \nabla \widetilde \phi({\bf n}), \end{eqnarray} [/math]

where [math] \widetilde Q^{E} [/math] and [math] \widetilde U^{E} [/math] are the filtered pure E-mode polarization maps given in equation (5) of Planck-2015-XLI[5], and [math] \widetilde \phi[/math] is the filtered lensing potential estimate.

We also provide the normalization transfer function [math] \mathcal{B}_\ell [/math] defined in equation (11) of Planck-2015-XLI[5], as well as the "L70" mask [math] M({\bf n}) [/math] that retains 69% of the sky before apodization, and its apodized version [math] \tilde{M}({\bf n}) [/math], which has an effective sky fraction [math] f_{\rm{sky}}^{\rm{eff}} = 65\% [/math].

As an example of the utilization of these products, the lensing B-mode maps that are shown in figure 4 of Planck-2015-XLI[5] are generated from

[math] Q^{\rm{lens}} \pm i U^{\rm{lens}} = \sum_{\ell m} \left( G_\ell \mathcal{B}_\ell^{-1} \int d{\bf n} {\, }_{\pm 2}Y_{\ell m}^*({\bf n}) \left(\bar{Q}^{\rm{lens}} \pm i \bar{U}^{\rm{lens}} \right) \right) {\, }_{\pm 2}Y_{\ell m}({\bf n}) [/math],

where [math]G_\ell[/math] is a Gaussian of 60 arcmin FWHM (introduced for highlighting large angular scales, although it can be removed or replaced by any other filter). This can be practically done by ingesting [math]\bar{Q}^{\rm{lens}} [/math] and [math] \bar{U}^{\rm{lens}} [/math] in the HEALPix "smoothing" routine, and using the product [math] G_\ell\mathcal{B}_\ell^{-1} [/math] as an input filtering function.

Specific products[edit]

We provide the lensing B-mode spherical harmonic coefficient estimate [math] B_{\ell m}^{\rm{lens}} [/math] over approximately 70% of the sky.

It can also be constructed using the raw products described above from

[math] B_{\ell m}^{\rm{lens}} = f_{10 \rightarrow 2000} \, \mathcal{B}_\ell^{-1} \, \, {\, }_{\pm 2}\mathcal{Y} \left[ \tilde{M}({\bf n}) \left( \bar{Q}^{\rm{lens}}({\bf n}) \pm i \bar{U}^{\rm{lens}}({\bf n}) \right) \right] [/math],

where [math] f_{10 \rightarrow 2000} [/math] is a function for producing band-powers over the range [math] 10 \le \ell \le 2000 [/math], and [math] {\, }_{\pm2}\mathcal{Y} [/math] is a short-hand notation for transforming a map into spin-weighted spherical harmonic coefficients [math] {\, }_{+2}a_{\ell m}[/math], [math]{\, }_{-2}a_{\ell m} [/math] and forming [math]1/(2i)\left({\, }_{+2}a_{\ell m} - {\, }_{-2}a_{\ell m}\right)[/math]. This can be done using, e.g., the HEALPix "anafast" tool.

The lensing B-mode power spectrum estimate [math] \hat{C}_\ell^{BB^{\rm{lens}}} [/math] discussed in Planck-2015-XLI[5] is obtained by forming the cross-correlation power spectrum of [math] B_{\ell m}^{\rm{lens}} [/math] and the B-mode data from the SMICA polarization maps [math] B_{\ell m} [/math]:

[math] \hat{C}_\ell^{BB^{\rm{lens}}} = \frac{\left(f_{\rm{sky}}^{\rm{eff}}\right)^{-1}}{2 \ell +1} G_\ell^{-2} \sum_m B_{\ell m}^* B_{\ell m}^{\rm{lens}}[/math],

where [math] G_\ell [/math] is the 5 arcmin Gaussian beam that convolves the SMICA CMB maps.


The products are contained in a single gzipped tarball named COM_Lensing-Bmode_R2.00.tgz. Its contents are described below.


Contents of Lensing B-mode package
Filename Format Description
bar_q_lens_map.fits HEALPix FITS format map in Galactic coordinates with [math] N_{\rm side} = 2048 [/math] Contains the non-normalized lensing-induced Q Stokes parameter map [math] \bar Q^{\rm{lens}}({\bf n}) [/math].
bar_u_lens_map.fits HEALPix FITS format map in Galactic coordinates with [math] N_{\rm side} = 2048 [/math] Contains the non-normalized lensing-induced U Stokes parameter map [math] \bar U^{\rm{lens}}({\bf n}) [/math].
mask.fits HEALPix FITS format map in Galactic coordinates with [math] N_{\rm side} = 2048 [/math] The L70 mask.
mask_noapo.fits HEALPix FITS format map in Galactic coordinates with [math] N_{\rm side} = 2048 [/math] The L70 mask without apodization.
transfer_function_b_l.dat ASCII text file, with columns = ([math]\ell[/math], [math] \mathcal{B}_\ell [/math]) The transfer function.
lensing_bmode_b_lm.fits HEALPix FITS format alm, with [math] \ell_{\rm max} = 2000 [/math] Contains the lensing B-mode harmonic coefficients [math] B_{\ell m}^{\rm{lens}} [/math].
lensing_bmode_bandpowers.dat ASCII text file, with columns = ([math]\ell_{\rm min}[/math], [math]\ell_{\rm b} [/math], [math]\ell_{\rm max} [/math], [math] \hat{C}_{\ell_{\rm b}}^{BB^{\rm{lens}}} [/math], [math] \Delta \hat{C}_{\ell_{\rm b}}^{BB^{\rm{lens}}} [/math] ) The lensing B-mode bandpower estimate on approximativily 70% of the sky and over the multipole range from 10 to 2000 shown in figure 9 of Planck-2015-XLI[5] (for plotting purposes only).

2015 Integrated Sachs-Wolfe effect map[edit]

We distribute estimates of the integrated Sachs-Wolfe (ISW) maps presented in Planck-2015-A21[6] as part of the 2015 data release. These map represents an estimate of the ISW anisotropies using different data sets:

  • SEVEM DX11 CMB map, together with all the large-scale structure tracers considered in the ISW paper, namely: NVSS, SDSS, WISE, and the Planck lensing map
  • Using only the large-scale structure tracers mentioned above
  • SEVEM DX11 CMB map, together with NVSS and the Planck lensing maps (since these two tracers capture most of the information, as compared to SDSS and WISE)


For all the three cases, the reconstruction is provided on approximately 85% of the sky, and they are produced using the LCB filter described in the Planck ISW paper (Section 5), described in detail in Barreiro et al. 2008 and Bonavera et al. 2016.

These ISW maps, together with their corresponding uncertainties maps and masks, are given in a file named COM_CompMap_ISW_0064_R2.00.fits. Its contents are described below.


Contents of the ISW maps file: COM_CompMap_ISW_0064_R2.00.fits
Extension Format Description Used data sets
0 HEALPix FITS format map with three components, [math]N_{\rm side}=64[/math], Ordering='Nest' Contains three components: i) ISW map [Kelvin], ii) Error map [Kelvin], iii) Mask map SEVEM DX11 CMB + NVSS + SDSS + WISE + Planck lensing.
1 HEALPix FITS format map with three components, [math]N_{\rm side}=64[/math], Ordering='Nest' Contains three components: i) ISW map [Kelvin], ii) Error map [Kelvin], iii) Mask map NVSS + SDSS + WISE + Planck lensing.
2 HEALPix FITS format map with three components, [math]N_{\rm side}=64[/math], Ordering='Nest' Contains three components: i) ISW map [Kelvin], ii) Error map [Kelvin], iii) Mask map SEVEM DX11 CMB + NVSS + Planck lensing.

References[edit]

  1. ↑ Planck 2015 results. XV. Gravitational Lensing, Planck Collaboration, 2016, A&A, 594, A15.
  2. ↑ 2.02.1 Planck 2015 results. XI. Diffuse component separation: CMB maps, Planck Collaboration, 2016, A&A, 594, A9.
  3. ↑ Planck 2015 results. XIII. Cosmological parameters, Planck Collaboration, 2016, A&A, 594, A13.
  4. ↑ Planck 2015 results. XXII. A map of the thermal Sunyaev-Zeldovich effect, Planck Collaboration, 2016, A&A, 594, A22.
  5. ↑ 5.05.15.25.35.45.55.65.7 Planck intermediate results. XLI. A map of lensing-induced B-modes, Planck Collaboration Int. XLI A&A, 596, A102, (2016).
  6. ↑ Planck 2015 results. XXI. The integrated Sachs-Wolfe effect, Planck Collaboration, 2016, A&A, 594, A21.

Cosmic Microwave background

(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).

Flexible Image Transfer Specification

Planck Legacy Archive

Full-Width-at-Half-Maximum

Retrieved from ‘https://wiki.cosmos.esa.int/planck-legacy-archive/index.php?title=Specially_processed_maps&oldid=12226’