Difference between revisions of "Lensing map reconstruction from 02.00 sims w/ and w/o foreground+inhomogeneous noise"

From CMB-S4 wiki
Jump to: navigation, search
Line 31: Line 31:
  
 
(Note that what I actually use is the mirror of the above data at Odyssey.)  
 
(Note that what I actually use is the mirror of the above data at Odyssey.)  
 +
 +
 +
Fig.1 shows the angular power spectra of the CMB alms.
 +
 +
[[File:fig_cmbcl_TT.jpg|400px]]
 +
[[File:fig_cmbcl_EE.jpg|400px]]
 +
[[File:fig_cmbcl_BB.jpg|400px]]
 +
 +
Fgi.1: The CMB temperature, EE, and BB power spectra from the signal only or signal+noise+foreground sims.
 +
  
 
----
 
----
Line 40: Line 50:
 
I first compute the diagonal filtering to be applied to the CMB harmonic coefficients (see text around Eq.(17) of BKVIII). The T, E and B at 500<L<3000 are used for the reconstruction. The unnormalized kappa estimator is then computed (e.g., for EB, fullsky counterpart of Eq.(19) of BKVIII). The analytic fullsky normalization is computed and multiplied to the unnormalized estimator. The reconstructed kappa map is finally cross-correlated with the input kappa map.  
 
I first compute the diagonal filtering to be applied to the CMB harmonic coefficients (see text around Eq.(17) of BKVIII). The T, E and B at 500<L<3000 are used for the reconstruction. The unnormalized kappa estimator is then computed (e.g., for EB, fullsky counterpart of Eq.(19) of BKVIII). The analytic fullsky normalization is computed and multiplied to the unnormalized estimator. The reconstructed kappa map is finally cross-correlated with the input kappa map.  
  
 
[[File:fig_cmbcl_TT.jpg|400px|The CMB temperature power spectrum from the signal only or signal+noise+foreground sims]]
 
[[File:fig_cmbcl_EE.jpg|400px|The CMB EE power spectrum from the signal only or signal+noise+foreground sims]]
 
[[File:fig_cmbcl_BB.jpg|400px|The CMB BB power spectrum from the signal only or signal+noise+foreground sims]]
 
  
 
----
 
----
Line 49: Line 55:
 
'''Results'''
 
'''Results'''
  
[[File:fig_kmap_input.jpg|400px|Example of the input kappa map after applying the mask. The kappa multipole at 50<L<200 is included in the fluctuations.]]
+
[[File:fig_kmap_input.jpg|400px]]
[[File:fig_kmap_recon.jpg|400px|Same as the above but for the reconstructed kappa map using EB estimator from signal+noise+foreground simulation.]]
+
[[File:fig_kmap_recon.jpg|400px]]
 +
 
 +
Fig.2: Example of the input kappa map after applying the mask. The kappa multipole at 50<L<200 is included in the fluctuations.
 +
The right panel is the same as the left panel but for the reconstructed kappa map using EB estimator from signal+noise+foreground simulation.
 +
 
 +
 
 +
[[File:fig_kk_TT.png|400px]]
 +
[[File:fig_kk_EB.png|400px]]
 +
 
 +
 
 +
Fig.3: The cross spectrum between input and reconstructed kappa (TT or EB estimator).
  
[[File:fig_kk_TT.png|400px|The cross spectrum between input and reconstructed kappa (TT estimator).]]
+
[[File:fig_vkk_TT.png|400px]]
[[File:fig_kk_EB.png|400px|The cross spectrum between input and reconstructed kappa (EB estimator).]]
+
[[File:fig_vkk_EB.png|400px]]
  
[[File:fig_vkk_TT.png|400px|The 1sigma error of the reconstructed kappa auto spectrum (TT estimator).]]
+
Fig.4: The 1 sigma error of the reconstructed kappa auto spectrum (TT or EB estimator).
[[File:fig_vkk_EB.png|400px|The 1sigma error of the reconstructed kappa auto spectrum (EB estimator).]]
 
  
  

Revision as of 20:44, 26 September 2018

Septemper 26, 2018 (Toshiya Namikawa posted)



Summary

In this posting, I show some results of the reconstruction of the CMB lensing kappa map from the signal only or signal+noise+foreground simulations. I used two quadratic estimators, TT and EB, and show that the reconstructed kappa map looks agree with the input kappa map even after including the noise and foreground. However, the reconstruction noise becomes larger at smaller scales of the reconstructed kappa. This could affect the delensing efficiency.



Data

Following [1], I use the following simulation data set:

- The high res masked lensed-LCDM single frequency map

 /project/projectdirs/cmbs4/data_xx.yy/02.00/cmbs4_02_llcdm_f145_b04_ellmin30_map_2048_mc_0000.fits

- The high res masked lensed-LCDM + noise + foreground single frequency map

 /project/projectdirs/cmbs4/data_xx.yy/02.00/cmbs4_02.00_comb_f145_b04_ellmin30_map_2048_mc_0000.fits

- Relative hits mask

 /project/projectdirs/cmbs4/expt_xx/02/rhits/n2048.fits.

- The true input phi map

 /project/projectdirs/cmb/data/generic/cmb/ffp10/mc/scalar/ffp10_unlensed_scl_cmb_000_tebplm_mc_0000.fits (4th binary table)

- The input lensed-LCDM angular power spectra (TT,TE,EE,BB)

 /project/projectdirs/cmb/data/generic/cmb/ffp10/cls/ffp10_lensedCls.dat

(Note that what I actually use is the mirror of the above data at Odyssey.)


Fig.1 shows the angular power spectra of the CMB alms.

Fig cmbcl TT.jpg Fig cmbcl EE.jpg Fig cmbcl BB.jpg

Fgi.1: The CMB temperature, EE, and BB power spectra from the signal only or signal+noise+foreground sims.



Method

The method of the lensing reconstruction is very similar to that of BICEP2/Keck Array Results Paper VIII [[2]] (hereafter BKVIII), except the curved sky analysis.

I first compute the diagonal filtering to be applied to the CMB harmonic coefficients (see text around Eq.(17) of BKVIII). The T, E and B at 500<L<3000 are used for the reconstruction. The unnormalized kappa estimator is then computed (e.g., for EB, fullsky counterpart of Eq.(19) of BKVIII). The analytic fullsky normalization is computed and multiplied to the unnormalized estimator. The reconstructed kappa map is finally cross-correlated with the input kappa map.



Results

Fig kmap input.jpg Fig kmap recon.jpg

Fig.2: Example of the input kappa map after applying the mask. The kappa multipole at 50<L<200 is included in the fluctuations. The right panel is the same as the left panel but for the reconstructed kappa map using EB estimator from signal+noise+foreground simulation.


Fig kk TT.png Fig kk EB.png


Fig.3: The cross spectrum between input and reconstructed kappa (TT or EB estimator).

Fig vkk TT.png Fig vkk EB.png

Fig.4: The 1 sigma error of the reconstructed kappa auto spectrum (TT or EB estimator).



Output/Usage

The input/reconstructed kappa map is currently located at

 /n/holylfs02/LABS/kovac_lab/cmbs4/nersc_mirror/cmbs4dat/reanalysis/phi_recons/02.00_namikawa_180920/kmap_inp/
 /n/holylfs02/LABS/kovac_lab/cmbs4/nersc_mirror/cmbs4dat/reanalysis/phi_recons/02.00_namikawa_180920/kmap_rec/

The file format is the fortran binary file. An example of reading the files (and showing the kappa map) is written in showkmap.py.