Difference between revisions of "Specially processed maps"

From Planck PLA Wiki
Jump to: navigation, search
(IRAM Maps of the Crab nebula)
 
(31 intermediate revisions by 4 users not shown)
Line 1: Line 1:
<span style="color:Red">NOTE: Text in red to be removed when filling the contents</span>
+
{{DISPLAYTITLE:Additional maps}}
==Product description==
+
==Overview==
<span style="color:Red">A general description of the product, including e.g. figures related to the contents (e.g. maps, tables), and some explanation of its scientific meaning. If there are scientific warnings about the use of the product (User’s caveats), they should also be given here, or at least references to other explanatory documents (papers etc).</span>
 
  
 +
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.
  
 
== Lensing map ==
 
== 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.
  
The fiducial Planck lensing potential map is made from a minimum variance combination of the 143 and 217GHz Planck maps on approximately 70% of the sky, using 857GHz as a dust template. This is the same lens reconstruction on which the Planck lensing likelihood is based. It is distributed in FITS format. The headers of this file are copied at the bottom of this page.
+
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>.
 +
; 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}} '''
 +
|-  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>
 +
|-
 +
|    2 ||    7 || -13.6379  || 15.3409
 +
|-
 +
|    8 ||  20 ||  6.0184  ||  4.8881
 +
|-
 +
|    21 ||  39 ||  -1.0675  ||  3.0940
 +
|-
 +
|    40 ||  65 ||  0.6135  ||  1.8474
 +
|-
 +
|    66 ||  100 ||  1.5030  ||  1.2696
 +
|-
 +
|  101 ||  144 ||  1.3760  ||  0.9950
 +
|-
 +
|  145 ||  198 ||  -1.2289  ||  0.8286
 +
|-
 +
|  199 ||  263 ||  1.1910  ||  0.7001
 +
|-
 +
|  264 ||  338 ||  -0.6567  ||  0.6197
 +
|-
 +
|  339 ||  425 ||  -0.8201  ||  0.5235
 +
|-
 +
|  426 ||  525 ||  -0.7581  ||  0.4850
 +
|-
 +
|  526 ||  637 ||  -0.3201  ||  0.5134
 +
|-
 +
|  638 ||  762 ||  -0.1589  ||  0.4073
 +
|-
 +
|  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 ===
 +
 
 +
The construction PHIBAR, RLPP and NLPP are described in detail in Sec. 2.1 of {{PlanckPapers|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.
 +
 
 +
===Inputs===
 +
 
 +
This product is built from the 143 and 217 GHz Planck [[Frequency Maps|frequency maps]], with 857GHz projected out as a dust template.
 +
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 {{PlanckPapers|planck2013-p12}}.
 +
 
 +
===File names and format===
 +
 
 +
A single file named
 +
*''{{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.
 +
 +
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>.
 +
 
 +
<center>
 +
<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>
  
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 \hat{\phi} with the lensing response function $R_L^{\phi\phi}$. The construction of these quantities are described in detail in Sec. 2.1 of $\fixme{P12_Lensing}$. This map has been band-limited between multipoles $10 \le L \le 2048$. The response function $R_L^{\phi\phi}$ 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. This is discussed further in Appendix A of P12_Lensing.
 
; MASK : This is a NSIDE 2048 HEALPix RING format map, containing the analysis mask used in the lens reconstruction.
 
; RLPP : This column contains the response function $R_L^{\phi\phi}$.
 
; NLPP : This column contains a sky-averaged estimate of the noise power spectrum of PHIBAR. The noise is highly coloured. Note that there is some dependence of the noise power spectrum with the local noise level of the map. Note that the noise power spectrum estimate here is not sufficiently accurate for a power spectrum analysis.
 
  
''[From Duncan Hanson, Feb.2013]''
+
{| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:left"
 +
|+ '''FITS file structure'''
 +
|- bgcolor="ffdead" 
 +
!colspan="4" | 1. EXTNAME = ''LENS-MAP''
 +
|-  bgcolor="ffdead" 
 +
! Column Name || Data Type || Units || Description
 +
|-
 +
|PHIBAR || Real*4 || none || Map of the lensing potential estimate, convolved with RLPP
 +
|-
 +
|MASK || Int || none || Region over which the lensing potential is reconstructed
 +
|- bgcolor="ffdead" 
 +
! Keyword || Data Type || Value || Description
 +
|-
 +
|PIXTYPE  ||  string || HEALPIX || colspan="2"|
 +
|-
 +
|COORDSYS  ||  string || GALACTIC || Coordinate system
 +
|-
 +
|ORDERING  ||  string || NESTED  || Healpix ordering
 +
|-
 +
|NSIDE  || Int*4  ||  2048 ||colspan="2"| Healpix Nside
 +
|-
 +
|FIRSTPIX || Int*4  ||          0 ||
 +
|-
 +
|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
 +
|}
  
==Production process==
+
== IRAM Maps of the Crab nebula ==
<span style="color:Red">Description of the Pipeline used to generate the product. In particular any limitations and approximations used in the data processing should be listed. Avoiding detailed descriptions of methods and referring to other parts of the ES and/or the relevant Planck papers for the details. References however should be quite detailed (i.e. it is not enough to direct the user to a paper, but the relevant section in the paper should be provided).</span>
 
  
==Inputs==
+
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
<span style="color:Red">A list (and brief description to the extent possible) of the input data used to generate this product (down to file names), as well as any external ancillary data sets which were used.</span>
 
==Related products==
 
<span style="color:Red">A description of other products that are related and share some commonalities with the product being described here. E.g. if the description is of a generic product (e.g. frequency maps), all the products falling into that type should be listed and referenced.</span>
 
==Format==
 
<span style="color:Red">A detailed description of the data format of each file, including header keywords for fits files, extension names, column names, formats….</span>
 
===File names===
 
===Extensions===
 
===Keywords===
 
===Numerical Format===
 
  
 +
: [[File:Crab_IRAM_2010.zip]]
  
The file structure is given below
+
See README in the tarball for full details. These data were used in{{BibCite|aumont2010}}
  
XTENSION= 'BINTABLE'          /Written by IDL:  Thu Feb  7 11:29:13 2013     
+
==References==
BITPIX  =                    8 /                                               
+
<References />
NAXIS  =                    2 /Binary table                                   
 
NAXIS1  =            251658240 /Number of bytes per row                       
 
NAXIS2  =                    1 /Number of rows                                 
 
PCOUNT  =                   0 /Random parameter count                         
 
GCOUNT  =                   1 /Group count                                   
 
TFIELDS =                    2 /Number of columns                             
 
COMMENT                                                                       
 
COMMENT  *** End of mandatory fields ***                                       
 
COMMENT                                                                       
 
EXTVER  =                   1 /Extension version                             
 
DATE    = '2013-02-07'        /Creation date                                 
 
COMMENT                                                                       
 
COMMENT  *** Column names ***                                                 
 
COMMENT                                                                       
 
TTYPE1  = 'WFMAP  '          / W-F Lensing map                               
 
TTYPE3  = 'Mask    '          / Mask                                         
 
TTYPE2  = 'MASK    '          /                                               
 
COMMENT                                                                       
 
COMMENT  *** Column formats ***                                               
 
COMMENT                                                                       
 
TFORM1  = '50331648E'          /                                               
 
TFORM2  = '50331648B'          /                                               
 
COMMENT                                                                       
 
COMMENT  *** Column units ***                                                 
 
COMMENT                                                                       
 
TUNIT1  = 'none    '          /                                               
 
TUNIT2  = 'none    '          /                                               
 
EXTNAME = 'LENS-MAP'          / Extension name                               
 
COORSYS = 'GALACTIC'          / Coordinate system                             
 
ORDERING= 'NESTED  '          / Healpix ordering                             
 
NSIDE  =                2048 / Healpix Nside                                 
 
FIRSTPIX=                    0 /                                               
 
LASTPIX =            50331647 /                                               
 
FILENAME= 'COM_Lensing_R1.00.fits' / FITS filename                             
 
CHECKSUM= 'D4iYD4fWD4fWD4fW'  / HDU checksum created 2013-02-07T10:29:14     
 
BAD_DATA=        -1.63750E+30 / bad pixel value                               
 
PROCVER = 'DX9    '          / Product version                               
 
COMMENT                                                                       
 
COMMENT ------------------------------------------------------------------------
 
COMMENT Lensing products:                                                     
 
COMMENT - wfmap (Wiener-filtered phi): map of the C_L^{\phi\phi} \bar{\phi}_{LM}
 
COMMENT - mask: map of the reconstruction mask                                 
 
COMMENT ------------------------------------------------------------------------
 
COMMENT HFI-DMC objects: creation date and group/name                         
 
COMMENT M3 dir: /space/dhanson/share/lensrec_v08_pub/                         
 
COMMENT    ... dx9nom_143_217_857m_g30_simnfg_simbeam_simlens/wfmap/         
 
COMMENT map:  wfmap.fits                                                       
 
COMMENT mask: mask.fits.gz                                                     
 
COMMENT ------------------------------------------------------------------------
 
END                                                                           
 
                                                                               
 
                                                                               
 
 
   
 
   
XTENSION= 'BINTABLE'          /Written by IDL:  Thu Feb  7 11:29:17 2013     
+
 
BITPIX  =                    8 /                                               
+
 
NAXIS  =                    2 /Binary table                                   
+
 
NAXIS1  =                8188 /Number of bytes per row                       
+
 
NAXIS2  =                    1 /Number of rows                                 
+
[[Category:Mission products|010]]
PCOUNT  =                    0 /Random parameter count                         
 
GCOUNT  =                    1 /Group count                                   
 
TFIELDS =                    1 /Number of columns                             
 
COMMENT                                                                       
 
COMMENT  *** End of mandatory fields ***                                       
 
COMMENT                                                                       
 
EXTVER  =                    1 /Extension version                             
 
DATE    = '2013-02-07'        /Creation date                                 
 
COMMENT                                                                       
 
COMMENT  *** Column names ***                                                 
 
COMMENT                                                                       
 
TTYPE1  = 'TLPP    '          / Transfer function                             
 
TUNIT1  = 'NONE    '          /                                               
 
COMMENT                                                                       
 
COMMENT  *** Column formats ***                                               
 
COMMENT                                                                       
 
TFORM1  = '2047E  '          /                                               
 
COMMENT                                                                       
 
COMMENT  *** Planck params ***                                                 
 
COMMENT                                                                       
 
EXTNAME = 'TransFun'          /Transfer Function                             
 
L_MIN  =                    2 /ell min                                       
 
L_MAX  =                2048 /ell max                                       
 
FILENAME= 'COM_Lensing_R1.00.fits' / FITS filename                             
 
CHECKSUM= 'aiAgci7ZaiAfai5Z'  / HDU checksum created 2013-02-07T10:29:17     
 
PROCVER = 'DX9    '          / Product version                               
 
COMMENT                                                                       
 
COMMENT ------------------------------------------------------------------------
 
COMMENT Lensing products:                                                     
 
COMMENT - TLPP: Transfer function T_l between \phi_{LM} and the map           
 
COMMENT ------------------------------------------------------------------------
 
COMMENT HFI-DMC objects: creation date and group/name                         
 
COMMENT M3 dir: /space/dhanson/share/lensrec_v08_pub/                         
 
COMMENT    ... dx9nom_143_217_857m_g30_simnfg_simbeam_simlens/wfmap/         
 
COMMENT tlpp:  tlpp.dat                                                       
 
COMMENT ------------------------------------------------------------------------
 
END
 

Latest revision as of 17:34, 23 July 2014

Overview[edit]

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.

Lensing map[edit]

Description[edit]

Here we present the minimum-variance (MV) lens reconstruction which forms the basis for the main results of Planck-2013-XVII[1]. 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:

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].
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 Planck-2013-XVII[1]. 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 Planck-2013-XVII[1]:

'MV' curl reconstruction bandpowers from Fig A2 of Planck-2013-XVII[1]
[math]L[/math]min [math]L[/math]max [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]
2 7 -13.6379 15.3409
8 20 6.0184 4.8881
21 39 -1.0675 3.0940
40 65 0.6135 1.8474
66 100 1.5030 1.2696
101 144 1.3760 0.9950
145 198 -1.2289 0.8286
199 263 1.1910 0.7001
264 338 -0.6567 0.6197
339 425 -0.8201 0.5235
426 525 -0.7581 0.4850
526 637 -0.3201 0.5134
638 762 -0.1589 0.4073
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[edit]

The construction PHIBAR, RLPP and NLPP are described in detail in Sec. 2.1 of Planck-2013-XVII[1]. 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.

Inputs[edit]

This product is built from the 143 and 217 GHz Planck frequency maps, with 857GHz projected out as a dust template. 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 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, including the ERCSC, SZ and PCCS . The reconstruction was performed using the fiducial beam window functions B(l) from the HFI RIMO . Details of the procedure used to produce a lensing estimate from these inputs are given in Planck-2013-XVII[1].

File names and format[edit]

A single file named

with two BINTABLE extensions containing the items described below.

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].


FITS file structure
1. EXTNAME = LENS-MAP
Column Name Data Type Units Description
PHIBAR Real*4 none Map of the lensing potential estimate, convolved with RLPP
MASK Int none Region over which the lensing potential is reconstructed
Keyword Data Type Value Description
PIXTYPE string HEALPIX
COORDSYS string GALACTIC Coordinate system
ORDERING string NESTED Healpix ordering
NSIDE Int*4 2048 Healpix Nside
FIRSTPIX Int*4 0
LASTPIX Int*4 50331647
2. EXTNAME = TransFun
Column Name Data Type Units Description
RLPP Real*4 none Response function
NLPP Real*4 none Sky-averaged noise power spectrum estimate
Keyword Data Type Value Description
L_MIN Int*4 0 First multipole
L_MAX Int*4 2048 Last multipole

IRAM Maps of the Crab nebula[edit]

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

File:Crab IRAM 2010.zip

See README in the tarball for full details. These data were used in[2]

References[edit]

  1. 1.01.11.21.31.41.5 Planck 2013 results: Gravitational lensing by large-scale structure, Planck Collaboration XVII, A&A, in press, (2014).
  2. Measurement of the Crab nebula polarization at 90 GHz as a calibrator for CMB experiments, J. Aumont, L. Conversi, C. Thum, H. Wiesemeyer, E. Falgarone, J. F. Macías-Pérez, F. Piacentini, E. Pointecouteau, N. Ponthieu, J. L. Puget, C. Rosset, J. A. Tauber, M. Tristram, A&A, 514, A70+, (2010).

(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

Flexible Image Transfer Specification