Difference between revisions of "Subpixel HFI"
(25 intermediate revisions by the same user not shown) | |||
Line 1: | Line 1: | ||
{{DISPLAYTITLE: Very bright regions where the sub-pixel effect in the SRoll CO foreground templates prevent them to be used for cosmology or astrophysics analysis}} | {{DISPLAYTITLE: Very bright regions where the sub-pixel effect in the SRoll CO foreground templates prevent them to be used for cosmology or astrophysics analysis}} | ||
− | |||
− | The bandpass corrections have been optimized for high latitude regions which implied to reduce the noise of the CO and dust bandpass templates to avoid the introduction of significant correlated noise. The effect is negligible for dust but not for CO in bright regions. We simulated the effect to evaluate its level. The difference between the | + | The bandpass corrections have been optimized for high latitude regions which implied to reduce the noise of the CO and dust bandpass templates to avoid the introduction of significant correlated noise. The effect is negligible for dust but not for CO in bright regions. We simulated the effect to evaluate its level. The difference between the downgraded (at N<sub>side</sub>=128) CO simulated map and the input full resolution (at N<sub>side</sub>=2048) CO simulated map gives the simulated bias introduced by this downgrading for each bolometer. For each frequency channel, these bias maps are converted to HPRs, multiplied by the corresponding CO bandpass coefficients computed by SRoll using data, and projected onto a frequency map of the bias for the frequency. Those frequency averaged bias maps are displayed in the table ''(a)''. |
− | For convenience, we also provide ''(b)'' the maps of the logarithm of the ratio of the variance, in each N<sub>side</sub>=128 pixel, | + | For convenience, we also provide ''(b)'' the maps of the logarithm of the ratio of the bias variance, in each N<sub>side</sub>=128 pixel, and the variance of the noise level estimated from the End-to-End simulations. The level of the bias is well determined and specific masks adapted to each scientific analysis can be computed at the appropriate level. The table also gives the sky fraction to be masked for thresholding the bias at 100, 10, and 1 % of the rms of the noise at N<sub>side</sub>=2048 inside the pixel at N<sub>side</sub>=128 referred to the rms of the noise at N<sub>side</sub>=2048 inside the same pixel at N<sub>side</sub>=128. |
− | We also provide ''(c)'' the ratio of the bias maps to the full intensity. | + | We also provide ''(c)'' the ratio of the variance bias maps to the variance of the full intensity signal. |
+ | |||
+ | The effects are largest in intensity and at 100 GHz. It can also be noted that, for polarization, and at all frequencies, only less than 0.04% of the sky is affected by this bias at 10% of the noise level. For intensity, the bias is a factor of 3 lower. If the user wants to reduce the bias at 1% of the noise in polarization, less than 2% of the sky will be affected. Users can also refer to the ratio of the absolute value of the bias with respect to the full intensity maps, which is 2 orders of magnitude lower, to build the relevant mask. For polarization, the two criteria are comparable. | ||
+ | |||
+ | As recommended in {{PlanckPapers|planck2016-l03}}, Galactic science work with the HFI data in these very bright regions should then be done starting from the [[Healpix_Rings_HFI| HFI HPRs]] at the full resolution. | ||
<!-- | <!-- | ||
− | /redtruck/delouis/PROD_RD12_RC4/CHECK/* | + | /redtruck/delouis/PROD_RD12_RC4/CHECK/*_CO_SUBPIX.fits |
− | + | /redtruck/delouis/PROD_RD12_RC4/CHECK/*_CO_SUBPIX_vs_NOISE_LOG10.fits | |
− | + | /redtruck/delouis/PROD_RD12_RC4/CHECK/*_CO_SUBPIX_vs_SIGNAL_LOG10.fits | |
--> | --> | ||
− | + | {| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:center" width=1400px | |
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | {| border="1" cellpadding="3" cellspacing="0" align="center" style="text-align:center" width= | ||
|+ '''CO templates ''' | |+ '''CO templates ''' | ||
|- bgcolor="ffdead" | |- bgcolor="ffdead" | ||
! | ! | ||
!colspan="3" | 100 GHz | !colspan="3" | 100 GHz | ||
− | !colspan=" | + | !colspan="1" | 143 GHz |
!colspan="3" | 217 GHz | !colspan="3" | 217 GHz | ||
!colspan="3" | 353 GHz | !colspan="3" | 353 GHz | ||
Line 35: | Line 30: | ||
|Q | |Q | ||
|U | |U | ||
− | | | + | |. |
− | |||
− | |||
|I | |I | ||
|Q | |Q | ||
Line 45: | Line 38: | ||
|U | |U | ||
|- | |- | ||
− | |''(a)'' bias maps | + | |style="text-align:left;" |''(a)'': frequency averaged bias maps [μK] |
− | |[[File:100- | + | |[[File:100-I_CO_Subpix_Val.png|100px]] |
− | |[[File:100- | + | |[[File:100-Q_CO_Subpix_Val.png|100px]] |
− | |[[File:100- | + | |[[File:100-U_CO_Subpix_Val.png|100px]] |
− | |||
− | |||
|. | |. | ||
− | |[[File:217- | + | |[[File:217-I_CO_Subpix_Val.png|100px]] |
− | |[[File:217- | + | |[[File:217-Q_CO_Subpix_Val.png|100px]] |
− | |[[File:217- | + | |[[File:217-U_CO_Subpix_Val.png|100px]] |
− | |[[File:353- | + | |[[File:353-I_CO_Subpix_Val.png|100px]] |
− | |[[File:353- | + | |[[File:353-Q_CO_Subpix_Val.png|100px]] |
− | |[[File:353- | + | |[[File:353-U_CO_Subpix_Val.png|100px]] |
|- | |- | ||
− | |''(b)'' bias noise | + | |style="text-align:left;" |''(b)'': <math>\log_{10} \left(\frac{<\textrm{bias}>\textrm{map}}{<\textrm{noise}>\textrm{map}}\right)</math> |
|[[File:100-I_CO_Subpix.png|100px]] | |[[File:100-I_CO_Subpix.png|100px]] | ||
|[[File:100-Q_CO_Subpix.png|100px]] | |[[File:100-Q_CO_Subpix.png|100px]] | ||
|[[File:100-U_CO_Subpix.png|100px]] | |[[File:100-U_CO_Subpix.png|100px]] | ||
− | |||
− | |||
|. | |. | ||
|[[File:217-I_CO_Subpix.png|100px]] | |[[File:217-I_CO_Subpix.png|100px]] | ||
Line 73: | Line 62: | ||
|[[File:353-U_CO_Subpix.png|100px]] | |[[File:353-U_CO_Subpix.png|100px]] | ||
|- | |- | ||
− | |''( | + | |style="text-align:left;" |sky fraction above 100% of ''(b)'' |
− | | | + | |<!--0.01 % --> <10<sup>-4</sup> % |
− | | | + | |<!--0.01 % --> <10<sup>-4</sup> % |
− | | | + | |<!--0.07 % --> 0.0005 % |
|. | |. | ||
− | |. | + | |<!--0.00 % --> <10<sup>-4</sup> % |
− | |. | + | |<!--0.01 % --> <10<sup>-4</sup> % |
− | | | + | |<!--0.01 % --> <10<sup>-4</sup> % |
− | | | + | |<!--0.00 % --> <10<sup>-4</sup> % |
− | | | + | |<!--0.01 % --> <10<sup>-4</sup> % |
− | | | + | |<!--0.01 % --> <10<sup>-4</sup> % |
− | |||
− | |||
|- | |- | ||
− | | | + | |style="text-align:left;" |sky fraction above 10% of ''(b)'' |
− | |0. | + | |<!--0.55 % --> 0.0112 % |
− | |0. | + | |<!--0.97 % --> 0.0092 % |
− | + | |<!--1.80 % --> 0.0346 % | |
− | |. | ||
− | |||
|. | |. | ||
− | |0. | + | |<!--0.04 % --> <10<sup>-4</sup> % |
− | |0. | + | |<!--0.29 % --> 0.0061 % |
− | |0. | + | |<!--0.19 % --> 0.0056 % |
− | |0. | + | |<!--0.01 % --> <10<sup>-4</sup> % |
− | |0. | + | |<!--0.38 % --> <10<sup>-4</sup> % |
− | |0. | + | |<!--0.29 % --> 0.0005 % |
|- | |- | ||
− | | | + | |style="text-align:left;" |sky fraction above 1% of ''(b)'' |
− | |0. | + | |<!--7.39 % --> 0.4186 % |
− | |0. | + | |<!--8.63 % --> 0.7370 % |
− | |1. | + | |<!--10.81 % --> 1.4318 % |
− | |||
|. | |. | ||
− | |. | + | |<!--0.90 % --> 0.0270 % |
− | |0. | + | |<!--4.98 % --> 0.1643 % |
− | |0. | + | |<!--2.79 % --> 0.1094 % |
− | |0. | + | |<!--0.46 % --> 0.0071 % |
− | |0. | + | |<!--5.81 % --> 0.0361 % |
− | | | + | |<!--3.38 % --> 0.0478 % |
− | |||
|- | |- | ||
− | | | + | |style="text-align:left;" |''(c)'': <math>\log_{10} \left(\frac{<\textrm{bias}>\textrm{map}}{<\textrm{full intensity}>\textrm{map}}\right)</math> |
− | | | + | |[[File:100-I_CO_Subpix_Val_DIV_Signal.png|100px]] |
− | | | + | |[[File:100-Q_CO_Subpix_Val_DIV_Signal.png|100px]] |
− | | | + | |[[File:100-U_CO_Subpix_Val_DIV_Signal.png|100px]] |
|. | |. | ||
− | |. | + | |[[File:217-I_CO_Subpix_Val_DIV_Signal.png|100px]] |
− | |. | + | |[[File:217-Q_CO_Subpix_Val_DIV_Signal.png|100px]] |
− | | | + | |[[File:217-U_CO_Subpix_Val_DIV_Signal.png|100px]] |
− | | | + | |[[File:353-I_CO_Subpix_Val_DIV_Signal.png|100px]] |
− | | | + | |[[File:353-Q_CO_Subpix_Val_DIV_Signal.png|100px]] |
− | | | + | |[[File:353-U_CO_Subpix_Val_DIV_Signal.png|100px]] |
− | | | ||
− | | | ||
|} | |} | ||
+ | |||
+ | All these files are available on the [[https://www.cosmos.esa.int/web/planck/pla | Planck Legacy Archive]]. | ||
==References== | ==References== | ||
<References /> | <References /> |
Latest revision as of 13:03, 30 May 2018
The bandpass corrections have been optimized for high latitude regions which implied to reduce the noise of the CO and dust bandpass templates to avoid the introduction of significant correlated noise. The effect is negligible for dust but not for CO in bright regions. We simulated the effect to evaluate its level. The difference between the downgraded (at Nside=128) CO simulated map and the input full resolution (at Nside=2048) CO simulated map gives the simulated bias introduced by this downgrading for each bolometer. For each frequency channel, these bias maps are converted to HPRs, multiplied by the corresponding CO bandpass coefficients computed by SRoll using data, and projected onto a frequency map of the bias for the frequency. Those frequency averaged bias maps are displayed in the table (a).
For convenience, we also provide (b) the maps of the logarithm of the ratio of the bias variance, in each Nside=128 pixel, and the variance of the noise level estimated from the End-to-End simulations. The level of the bias is well determined and specific masks adapted to each scientific analysis can be computed at the appropriate level. The table also gives the sky fraction to be masked for thresholding the bias at 100, 10, and 1 % of the rms of the noise at Nside=2048 inside the pixel at Nside=128 referred to the rms of the noise at Nside=2048 inside the same pixel at Nside=128.
We also provide (c) the ratio of the variance bias maps to the variance of the full intensity signal.
The effects are largest in intensity and at 100 GHz. It can also be noted that, for polarization, and at all frequencies, only less than 0.04% of the sky is affected by this bias at 10% of the noise level. For intensity, the bias is a factor of 3 lower. If the user wants to reduce the bias at 1% of the noise in polarization, less than 2% of the sky will be affected. Users can also refer to the ratio of the absolute value of the bias with respect to the full intensity maps, which is 2 orders of magnitude lower, to build the relevant mask. For polarization, the two criteria are comparable.
As recommended in Planck-2020-A3[1], Galactic science work with the HFI data in these very bright regions should then be done starting from the HFI HPRs at the full resolution.
All these files are available on the [| Planck Legacy Archive].
References[edit]
- ↑ Planck 2018 results. III. High Frequency Instrument data processing and frequency maps, Planck Collaboration, 2020, A&A, 641, A3.
(Planck) High Frequency Instrument