Difference between revisions of "Effective Beams"

Revision as of 16:45, 17 March 2013

Product description[edit]

The effective beam is the average of all scanning beams pointing at a certain direction within a given pixel of the sky map for a given scan strategy. It takes into account the coupling between azimuthal asymmetry of the beam and the uneven distribution of scanning angles across the sky. It captures the complete information about the difference between the true and observed image of the sky. They are, by definition, the objects whose convolution with the true CMB sky produce the observed sky map.

Details of the beam processing are given in the respective pages for HFI and LFI.

The full algebra involving the effective beams for temperature and polarisation was presented in [Mitra, Rocha, Gorski et al.] #mitra2010, and a discussion of its application to Planck data is given in the appropriate LFI #planck2013-p02d and HFI #planck2013-p03c papers. Relevant details of the processing steps are given in the Effective Beams section of this document.

Comparison of the images of compact sources observed by Planck with FEBeCoP products[edit]

We show here a comparison of the FEBeCoP derived effective beams, and associated point spread functions, PSF (the transpose of the beam matrix), to the actual images of a few compact sources observed by Planck, for all LFI and HFI frequency channels, as an example. We show below a few panels of source images organized as follows:

Row #1- DX9 images of four ERCSC objects with their galactic (l,b) coordinates shown under the color bar
Row #2- linear scale FEBeCoP PSFs computed using input scanning beams, Grasp Beams, GB, for LFI and B-Spline beams,BS, Mars12 apodized for the CMB channels and the BS Mars12 for the sub-mm channels, for HFI (see section Inputs below).
Row #3- log scale of #2; PSF iso-contours shown in solid line, elliptical Gaussian fit iso-contours shown in broken line

Comparison images of compact sources and effective beams, PSFs
30GHz
44GHz
70GHz
100GHz
143GHz
217GHz
353GHz
545GHz
857GHz

Histograms of the effective beam parameters[edit]

Here we present histograms of the three fit parameters - beam FWHM, ellipticity, and orientation with respect to the local meridian and of the beam solid angle. The shy is sampled (pretty sparsely) at 3072 directions which were chosen as HEALpix nside=16 pixel centers for HFI and at 768 directions which were chosen as HEALpix nside=8 pixel centers for LFI to uniformly sample the sky.

Where beam solid angle is estimated according to the definition: [math] 4 \pi \sum[/math](effbeam)/max(effbeam) i.e.,

Histograms for LFI effective beam parameters

Histograms for HFI effective beam parameters

Sky variation of effective beams solid angle and ellipticity of the best-fit Gaussian[edit]

The discontinuities at the Healpix domain edges in the maps are a visual artifact due to the interplay of the discretized effective beam and the Healpix pixel grid.

Sky variation of effective beams ellipticity of the best-fit Gaussian
ellipticity - 30GHz
ellipticity - 44GHz
ellipticity - 70GHz
ellipticity - 100GHz
ellipticity - 143GHz
ellipticity - 217GHz
ellipticity - 353GHz
ellipticity - 545GHz
ellipticity - 857GHz

Sky (relative) variation of effective beams solid angle of the best-fit Gaussian
beam solid angle (relative) variations wrt scanning beam - 30GHz
beam solid angle (relative) variations wrt scanning beam - 44GHz
beam solid angle (relative) variations wrt scanning beam - 70GHz
beam solid angle (relative) variations wrt scanning beam - 100GHz
beam solid angle (relative) variations wrt scanning beam - 143GHz
beam solid angle (relative) variations wrt scanning beam - 217GHz
beam solid angle (relative) variations wrt scanning beam - 353GHz
beam solid angle (relative) variations wrt scanning beam - 545GHz
beam solid angle (relative) variations wrt scanning beam - 857GHz

Sky (relative) variation of effective beams fwhm of the best-fit Gaussian
fwhm (relative) variations wrt scanning beam - 30GHz
fwhm (relative) variations wrt scanning beam - 44GHz
fwhm (relative) variations wrt scanning beam - 70GHz
fwhm (relative) variations wrt scanning beam - 100GHz
fwhm (relative) variations wrt scanning beam - 143GHz
fwhm (relative) variations wrt scanning beam - 217GHz
fwhm (relative) variations wrt scanning beam - 353GHz
fwhm (relative) variations wrt scanning beam - 545GHz
fwhm (relative) variations wrt scanning beam - 857GHz

Sky variation of effective beams $\psi$ angle of the best-fit Gaussian
$\psi$ - 30GHz
$\psi$ - 44GHz
$\psi$ - 70GHz
$\psi$ - 100GHz
$\psi$ - 143GHz
$\psi$ - 217GHz
$\psi$ - 353GHz
$\psi$ - 545GHz
$\psi$ - 857GHz

Statistics of the effective beams computed using FEBeCoP[edit]

We tabulate the simple statistics of FWHM, ellipticity (e), orientation ([math] \psi[/math]) and beam solid angle, ([math] \Omega [/math]), for a sample of 3072 and 768 directions on the sky for HFI and LFI data respectively. Statistics shown in the Table are derived from the histograms shown above.

The derived beam parameters are representative of the DPC NSIDE 1024 and 2048 healpix maps (they include the pixel window function).
The reported FWHM_eff are derived from the beam solid angles, under a Gaussian approximation. These are best used for flux determination while the the Gaussian fits to the effective beam maps are more suited for source identification.

**Statistics of the FEBeCoP Effective Beams Computed with the BS Mars12 apodized for the CMB channels and oversampled**
frequency	mean(fwhm) [arcmin]	sd(fwhm) [arcmin]	mean(e)	sd(e)	mean([math] \psi[/math]) [degree]	sd([math] \psi[/math]) [degree]	mean([math] \Omega [/math]) [arcmin[math]^{2}[/math]]	sd([math] \Omega [/math]) [arcmin[math]^{2}[/math]]	FWHM_eff [arcmin]
030	32.239	0.013	1.320	0.031	-0.304	55.349	1189.513	0.842	32.34
044	27.005	0.552	1.034	0.033	0.059	53.767	832.946	31.774	27.12
070	13.252	0.033	1.223	0.026	0.587	55.066	200.742	1.027	13.31
100	9.651	0.014	1.186	0.023	-0.024	55.400	105.778	0.311	9.66
143	7.248	0.015	1.036	0.009	0.383	54.130	59.954	0.246	7.27
217	4.990	0.025	1.177	0.030	0.836	54.999	28.447	0.271	5.01
353	4.818	0.024	1.147	0.028	0.655	54.745	26.714	0.250	4.86
545	4.682	0.044	1.161	0.036	0.544	54.876	26.535	0.339	4.84
857	4.325	0.055	1.393	0.076	0.876	54.779	24.244	0.193	4.63

Beam solid angles for the PCCS[edit]

- [math]\Omega_{eff}[/math] - is the mean beam solid angle of the effective beam, where beam solid angle is estimated according to the definition: 4pi*sum(effbeam)/max(effbeam), i.e. as an integral over the full extent of the effective beam, i.e. .

- from [math]\Omega_{eff}[/math] we estimate the [math]fwhm_{eff}[/math], under a Gaussian approximation - these are tabulated above
- [math]\Omega^{(1)}_{eff}[/math] is the beam solid angle estimated up to a radius equal to one [math]fwhm_{eff}[/math] and [math]\Omega^{(2)}_{eff}[/math] up to a radius equal to twice the [math]fwhm_{eff}[/math].
  - These were estimated according to the procedure followed in the aperture photometry code for the PCCS: if the pixel centre does not lie within the given radius it is not included (so inclusive=0 in query disc).

**Band averaged beam solid angles**
Band	[math]\Omega_{eff}[/math][arcmin[math]^{2}[/math]]	spatial variation [arcmin[math]^{2}[/math]]	[math]\Omega^{(1)}_{eff}[/math] [arcmin[math]^{2}[/math]]	spatial variation-1 [arcmin[math]^{2}[/math]]	[math]\Omega^{(2)}_{eff}[/math] [arcmin[math]^{2}[/math]]	spatial variation-2 [arcmin[math]^{2}[/math]]
30	1189.513	0.842	1116.494	2.274	1188.945	0.847
44	832.946	31.774	758.684	29.701	832.168	31.811
70	200.742	1.027	186.260	2.300	200.591	1.027
100	105.778	0.311	100.830	0.410	105.777	0.311
143	59.954	0.246	56.811	0.419	59.952	0.246
217	28.447	0.271	26.442	0.537	28.426	0.271
353	26.714	0.250	24.827	0.435	26.653	0.250
545	26.535	0.339	24.287	0.455	26.302	0.337
857	24.244	0.193	22.646	0.263	23.985	0.191

Production process[edit]

FEBeCoP, or Fast Effective Beam Convolution in Pixel space [Mitra, Rocha, Gorski et al.], is an approach to representing and computing effective beams (including both intrinsic beam shapes and the effects of scanning) that comprises the following steps:

identify the individual detectors' instantaneous optical response function (presently we use elliptical Gaussian fits of Planck beams from observations of planets; eventually, an arbitrary mathematical representation of the beam can be used on input)
follow exactly the Planck scanning, and project the intrinsic beam on the sky at each actual sampling position
project instantaneous beams onto the pixelized map over a small region (typically <2.5 FWHM diameter)
add up all beams that cross the same pixel and its vicinity over the observing period of interest
create a data object of all beams pointed at all N'_pix_' directions of pixels in the map at a resolution at which this precomputation was executed (dimension N'_pix_' x a few hundred)
use the resulting beam object for very fast convolution of all sky signals with the effective optical response of the observing mission

Computation of the effective beams at each pixel for every detector is a challenging task for high resolution experiments. FEBeCoP is an efficient algorithm and implementation which enabled us to compute the pixel based effective beams using moderate computational resources. The algorithm used different mathematical and computational techniques to bring down the computation cost to a practical level, whereby several estimations of the effective beams were possible for all Planck detectors for different scanbeam models and different lengths of datasets.

Pixel Ordered Detector Angles (PODA)[edit]

The main challenge in computing the effective beams is to go through the trillion samples, which gets severely limited by I/O. In the first stage, for a given dataset, ordered lists of pointing angles for each pixels---the Pixel Ordered Detector Angles (PODA) are made. This is an one-time process for each dataset. We used computers with large memory and used tedious memory management bookkeeping to make this step efficient.

effBeam[edit]

The effBeam part makes use of the precomputed PODA and unsynchronized reading from the disk to compute the beam. Here we tried to made sure that no repetition occurs in evaluating a trigonometric quantity.

One important reason for separating the two steps is that they use different schemes of parallel computing. The PODA part requires parallelisation over time-order-data samples, while the effBeam part requires distribution of pixels among different computers.

Computational Cost[edit]

The computation of the effective beams has been performed at the NERSC Supercomputing Center. The table below shows the computation cost for FEBeCoP processing of the nominal mission.

Computational cost for PODA, Effective Beam and single map convolution. Wall-clock time is given as a guide, as found on the NERSC supercomputers.
Channel	030	044	070	100	143	217	353	545	857
PODA/Detector Computation time (CPU hrs)	85	100	250	500	500	500	500	500	500
PODA/Detector Computation time (wall clock hrs)	7	10	20	20	20	20	20	20	20
Beam/Channel Computation time (CPU hrs)	900	2000	2300	2800	3800	3200	3000	900	1100
Beam/Channel Computation time (wall clock hrs)	0.5	0.8	1	1.5	2	1.2	1	0.5	0.5
Convolution Computation time (CPU hr)	1	1.2	1.3	3.6	4.8	4.0	4.1	4.1	3.7
Convolution Computation time (wall clock sec)	1	1	1	4	4	4	4	4	4
Effective Beam Size (GB)	173	123	28	187	182	146	132	139	124

The computation cost, especially for PODA and Convolution, is heavily limited by the I/O capacity of the disc and so it depends on the overall usage of the cluster done by other users.

Inputs[edit]

In order to fix the convention of presentation of the scanning and effective beams, we show the classic view of the Planck focal plane as seen by the incoming CMB photon. The scan direction is marked, and the toward the center of the focal plane is at the 85 deg angle w.r.t spin axis pointing upward in the picture.

Planck Focal Plane

The Focal Plane DataBase (FPDB)[edit]

The FPDB contains information on each detector, e.g., the orientation of the polarisation axis, different weight factors, ... (see the instrument RIMOs):

HFI - LFI_RIMO_DX9_PTCOR6 - The HFI RIMO
LFI - HFI-RIMO-3_16_detilt_t2_ptcor6.fits - The LFI RIMO

Template:PLADoc

The scanning strategy[edit]

The scanning strategy, the three pointing angle for each detector for each sample: Detector pointings for the nominal mission covers about 15 months of observation from Operational Day (OD) 91 to OD 563 covering 3 surveys and half.

The scanbeam[edit]

The scanbeam modeled for each detector through the observation of planets. Which was assumed to be constant over the whole mission, though FEBeCoP could be used for a few sets of scanbeams too.

LFI: GRASP scanning beam - the scanning beams used are based on Radio Frequency Tuned Model (RFTM) smeared to simulate the in-flight optical response.
HFI: B-Spline, BS based on 2 observations of Mars.

(see the instrument RIMOs)

HFI - LFI_RIMO_DX9_PTCOR6 - The HFI RIMO
LFI - HFI-RIMO-3_16_detilt_t2_ptcor6.fits - The LFI RIMO

LFI effective beams

Beam cutoff radii[edit]

N times geometric mean of FWHM of all detectors in a channel, where N

**Beam cut off radius**
channel	Cutoff Radii in units of fwhm	fwhm of full beam extent
30 - 44 - 70	2.5
100	2.25	23.703699
143	3	21.057402
217-353	4	18.782754
sub-mm	4	18.327635(545GHz) ; 17.093706(857GHz)

Map resolution for the derived beam data object[edit]

[math]N_{side} = 1024 [/math] for LFI frequency channels
[math]N_{side} = 2048 [/math] for HFI frequency channels

Related products[edit]

Monte Carlo simulations[edit]

FEBeCoP software enables fast, full-sky convolutions of the sky signals with the Effective beams in pixel domain. Hence, a large number of Monte Carlo simulations of the sky signal maps map convolved with realistically rendered, spatially varying, asymmetric Planck beams can be easily generated. We performed the following steps:

generate the effective beams with FEBeCoP for all frequencies for dDX9 data and Nominal Mission
generate 100 realizations of maps from a fiducial CMB power spectrum
convolve each one of these maps with the effective beams using FEBeCoP
estimate the average of the Power Spectrum of each convolved realization, C'_\ell_'^out^'}, and 1 sigma errors

As FEBeCoP enables fast convolutions of the input signal sky with the effective beam, thousands of simulations are generated. These Monte Carlo simulations of the signal (might it be CMB or a foreground (e.g. dust)) sky along with LevelS+Madam noise simulations were used widely for the analysis of Planck data. A suite of simulations were rendered during the mission tagged as Full Focalplane simulations, FFP#. For example FFP6

Beam Window Functions[edit]

The Transfer Function or the Beam Window Function [math] W_l [/math] relates the true angular power spectra [math]C_l [/math] with the observed angular power spectra [math]\widetilde{C}_l [/math]:

Note that, the window function can contain a pixel window function (depending on the definition) and it is {\em not the angular power spectra of the scanbeams}, though, in principle, one may be able to connect them though fairly complicated algebra.

The window functions are estimated by performing Monte-Carlo simulations. We generate several random realisations of the CMB sky starting from a given fiducial [math] C_l [/math], convolve the maps with the pre-computed effective beams, compute the convolved power spectra [math] C^\text{conv}_l [/math], divide by the power spectra of the unconvolved map [math]C^\text{in}_l [/math] and average over their ratio. Thus, the estimated window function

For subtle reasons, we perform a more rigorous estimation of the window function by comparing C^{conv}_l with convolved power spectra of the input maps convolved with a symmetric Gaussian beam of comparable (but need not be exact) size and then scaling the estimated window function accordingly.

Beam window functions are provided in the RIMO.

Beam Window functions, Wl, for Planck mission[edit]

Beam Window functions, Wl, for LFI channels

Beam Window functions, Wl, for HFI channels

File Names[edit]

The effective beams are stored as unformatted files in directories with the frequency channel's name, e.g., 100GHz, each subdirectory contains N unformatted files with names beams_###.unf, a beam_index.fits and a beams_run.log. For 100GHz and 143GHz: N=160, for 30, 44, 70 217 and 353GHz: N=128; for 545GHz: N=40; and 857GHz: N=32.

beam_index.fits
beams_run.log

Retrieval of effective beam information from the PLA interface[edit]

In order to retrieve the effective beam information, the user should first launch the Java interface from this page: http://www.sciops.esa.int/index.php?project=planck&page=Planck_Legacy_Archive

One should click on "Sky maps" and then open the "Effective beams" area. There is the possibility to either retrieve one beam nearest to the input source (name or coordinates), or to retrieve a set of beams in a grid defined by the Nside and the size of the region around a source (name or coordinates). The resolution of this grid is defined by the Nside parameter. The size of the region is defined by the "Radius of ROI" parameter.

Once the user proceeds with querying the beams, the PLA software retrieves the appropriate set of effective beams from the database and delivers it in a FITS file which can be directly downloaded.

Meta data[edit]

The data format of the effective beams is unformatted.

References[edit]

References

</biblio>

Revision as of 15:05, 17 March 2013 (view source) Agregori (talk \| contribs) (→‎Sky variation of effective beams solid angle and ellipticity of the best-fit Gaussian) ← Older edit		Revision as of 16:45, 17 March 2013 (view source) Agregori (talk \| contribs) (→‎Inputs) Newer edit →
Line 247:		Line 247:


−	[[File:PlanckFocalPlane.png \| 500px\| thumb \| center\| "'Planck Focal Plane''']]	+	[[File:PlanckFocalPlane.png \| 500px\| thumb \| center\|'''Planck Focal Plane''']]