Effective Beams
Contents
 1 Product description
 2 Production process
 3 Inputs
 4 Related products
 5 File Names
 6 File format
 7 Previous Releases: (2015) and (2013) effective beams
 8 References
Product description
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 images of the sky. The effective beams are, by definition, the objects whose convolution with the true CMBCosmic Microwave background sky produce the observed sky map.
Details of the beam processing are given in the respective sections for HFI and LFI.
The full algebra involving the effective beams for temperature and polarization was presented in ^{[1]}, and a discussion of its application to Planck data is given in the appropriate LFI(Planck) Low Frequency Instrument (Planck2013IV^{[2]}, Planck2015A05^{[3]}) and HFI(Planck) High Frequency Instrument (Planck2013VII^{[4]}) 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
We show here a comparison of the FEBeCoPderived effective beams, and associated point spread functions, PSF (the transpose of the beam matrix), to the actual images for a few example compact sources observed by Planck, for all LFI(Planck) Low Frequency Instrument and HFI(Planck) High Frequency Instrument frequency channels. We show below a few panels of source images organized as follows:
 Row #1  DX9 images of four ERCSCEarly Release Compact Source Catalog objects with their Galactic (l,b) coordinates shown under the colour bar;
 Row #2  linear scale FEBeCoP PSFs computed using input scanning beams, Grasp Beams (GBs) for LFI(Planck) Low Frequency Instrument and BSpline beams (BSs) of Mars12 apodized for the HFI(Planck) High Frequency Instrument CMBCosmic Microwave background channels and Mars12 for the submm channels (see section "Inputs" below);
 Row #3  log scale of #2, PSF isocontours being shown as solid lines, elliptical Gaussian fit isocontours shown as broken lines.
Histograms of the effective beam parameters
Here we present histograms of the three fit parameters, i.e., beam FWHMFullWidthatHalfMaximum, ellipticity, and orientation with respect to the local meridian, and of the beam solid angle. The sky is uniformly sampled (although fairly sparsely) in 3072 directions that were chosen as HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). N_{side}=16 pixel centres for HFI(Planck) High Frequency Instrument and in 768 directions that were chosen as HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). N_{side}=8 pixel centres for LFI(Planck) Low Frequency Instrument.
Here the beam solid angle is estimated according to the definition: "4π Σ(effbeam)/max(effbeam)" i.e., 4π Σ(B_{ij}) / max(B_{ij}).
Sky variation of effective beams solid angle and ellipticity of the bestfit Gaussian
 The discontinuities at the HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). domain edges in the maps are a visual artefact due to the interplay of the discretized effective beam and the HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). pixel grid.
Statistics of the effective beams computed using FEBeCoP
Below we tabulate the simple statistics of FWHMFullWidthatHalfMaximum, ellipticity (e), orientation (ψ) and beam solid angle, (Ω), for a sample of 3072 and 768 directions on the sky for HFI(Planck) High Frequency Instrument and LFI(Planck) Low Frequency Instrument data, respectively. Statistics shown in the Table are derived from the histograms shown above.
 The derived beam parameters are representative of the DPCData Processing Center N_{side} 1024 and 2048 HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). maps (including the pixel window function).
 The reported FWHMFullWidthatHalfMaximum_{eff} values are derived from the beam solid angles, under a Gaussian approximation. These are appropriate for flux determination, while the the Gaussian fits to the effective beam maps are more suited for source identification.
Frequency [GHz]  Mean(FWHMFullWidthatHalfMaximum) [arcmin]  σ(FWHMFullWidthatHalfMaximum) [arcmin]  Mean(e)  σ(e)  Mean(ψ) [deg]  σ(ψ) [deg]  Mean(Ω) [arcmin^{2}]  σ(Ω) [arcmin^{2}]  FWHMFullWidthatHalfMaximum_{eff} [arcmin] 

30  32.239  0.013  1.320  0.031  0.304  55.349  1189.513  0.842  32.34 
44  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
For the Planck Catalogue of Compact Sources we also need the following definitions.
 Ω_{eff}  is the mean beam solid angle of the effective beam, where the solid angle is estimated according to the definition "4 π Σ(effbeam)/max(effbeam)", i.e., as an integral over the full extent of the effective beam, 4π Σ (B_{ij})/max(B_{ij}).
 From Ω_{eff} we estimate FWHMFullWidthatHalfMaximum_{eff}, under a Gaussian approximation  these are tabulated above.
 Ω^{(1)}_{eff} is the beam solid angle estimated up to a radius equal to 1 FWHMFullWidthatHalfMaximum_{eff} and Ω^{(2)}_{eff} up to a radius equal to twice the FWHMFullWidthatHalfMaximum_{eff}.
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 [GHz]  Ω_{eff}[arcmin^{2}]  Spatial variation [arcmin^{2}]  Ω^{(1)}_{eff} [arcmin^{2}]  Spatial variation1 [arcmin^{2}]  Ω^{(2)}_{eff} [arcmin^{2}]  Spatial variation2 [arcmin^{2}] 
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
FEBeCoP, the Fast Effective Beam Convolution in Pixel space^{[1]}, is an approach for 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 the observations of planets, but in general, an arbitrary mathematical representation of the beam can be used on input);
 exactly follow 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 FWHMFullWidthatHalfMaximum in 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} times 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 pixelbased effective beams using moderate computational resources. The algorithm used different mathematical and computational techniques to bring down the computation cost to a practical level, so that several estimations of the effective beams were possible for all Planck detectors for different scan and beam models, as well as different lengths of data set.
Pixelordered detector angles (PODA)
The main challenge in computing the effective beams is to go through the roughly trillion samples, which becomes severely limited by I/O. In the first stage, for a given data set, ordered lists of pointing angles for each pixel (the pixelordered detector angles, PODA) are made. This is an onetime process for each data set. We used computers with large memory and tedious memory management and bookkeeping to make this step efficient.
Effective beam
The "effBeam" part makes use of the precomputed PODA and unsynchronized reading from the disk to compute the beam. Here we have tried to made sure that no repetition occurs in evaluating any trigonometric quantity.
One important reason for separating the two steps is that they use different parallel computing schemes. The PODA part requires parallelization over timeordered data samples, while the effBeam part requires distribution of pixels among different computers.
Computational Cost
The computation of the effective beams has been performed at the NERSC Supercomputing Center. The table below shows the computational cost for FEBeCoP processing of the nominal mission.
Channel [GHz]  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 hrs]  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 disk and so it depends on the overall usage of the cluster by other users.
Inputs
In order to describe the convention of presentation of the scanning and effective beams, we show the classic view of the Planck focal plane as seen by incoming CMBCosmic Microwave background photons. The scan direction is marked, and the centre of the focal plane is at an angle of 85° with respect to the spin axis, pointing upward in the picture.
The focal plane database (FPDB)
The FPDB contains information on each detector, e.g., the orientation of the polarization axis, and different weight factors. Details are provided in the instrument RIMOs):
 for HFI(Planck) High Frequency Instrument  the HFI RIMO;
 for LFI(Planck) Low Frequency Instrument  the LFI RIMO.
The scanning strategy
The scanning strategy description consists of the three pointing angles for each detector for each sample. Detector pointings for the nominal mission cover about 15 months of observation from Operational Day (ODOperation Day definition is geometric visibility driven as it runs from the start of a DTCP (satellite Acquisition Of Signal) to the start of the next DTCP. Given the different ground stations and spacecraft will takes which station for how long, the OD duration varies but it is basically once a day.) 91 to ODOperation Day definition is geometric visibility driven as it runs from the start of a DTCP (satellite Acquisition Of Signal) to the start of the next DTCP. Given the different ground stations and spacecraft will takes which station for how long, the OD duration varies but it is basically once a day. 563, covering three and a half surveys.
The scanning beam
The "scanbeam" is modelled for each detector through the observation of planets. This was assumed to be constant over the whole mission, although FEBeCoP could be used for a few sets of scanbeams too.
 For LFI(Planck) Low Frequency Instrument the GRASP scanning beams are used; these are based on the radiofrequency tuned model (RFTM), smeared to simulate the inflight optical response.
 For HFI(Planck) High Frequency Instrument the BSpline, BS beams are used based on two observations of Mars.
Beam cutoff radii
The cutoffs are set to N times the geometric mean of the FWHMFullWidthatHalfMaximum of all detectors in a channel, where N is given below.
Channel [GHz]  Cutoff radii in units of FWHMFullWidthatHalfMaximum  FWHMFullWidthatHalfMaximum of full beam extent [arcmin] 
30, 44, 70  2.5  
100  2.25  23.703699 
143  3  21.057402 
217, 353  4  18.782754 
545  4  18.327635 
857  4  17.093706 
Map resolution for the derived beam data object
The resolution is
 N_{side} = 1024 for LFI(Planck) Low Frequency Instrument frequency channels, and
 N_{side} = 2048 for HFI(Planck) High Frequency Instrument frequency channels.
Related products
Monte Carlo simulations
FEBeCoP software enables fast, fullsky convolutions of the sky signals with the effective beams in the pixel domain. Hence, a large number of Monte Carlo simulations of the sky signal maps can be easily generated, which are convolved with realistically rendered, spatially varying, asymmetric Planck beams. We performed the following steps:
 generate the effective beams with FEBeCoP for all frequencies for "dDX9" data and the Nominal Mission;
 generate 100 realizations of maps from a fiducial CMBCosmic Microwave background 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, and 1 σ errors.
Since FEBeCoP enables fast convolutions of the input signal sky with the effective beam, thousands of simulations can be generated. These Monte Carlo simulations of the sky signal (which might be CMBCosmic Microwave background or a foreground, e.g., dust) along with LevelS+Madam noise simulations were used widely for the analysis of Planck data. A suite of simulations were rendered during the mission and tagged as full focalplane simulations, FFP#
(for example the FFP6).
Beam Window Functions
The "Transfer Function" or the "Beam Window Function", W_{ℓ} relates the true angular power spectra C_{ℓ} with the observed angular power spectra Ĉ_{ℓ}:
Note that, the window function can contain a pixel window function (depending on the definition) and it is not the angular power spectra of the scanbeams, although in principle one may be able to connect them via some complicated algebra.
The window functions are estimated by performing Monte Carlo simulations. We generate several random realizations of the CMBCosmic Microwave background sky, starting from a given fiducial C_{ℓ}, convolve the maps with the precomputed effective beams, compute the convolved power spectra C_{ℓ}^{conv}, divide by the power spectra of the unconvolved map C_{ℓ}^{in} and average over their ratio. Thus, the estimated window function is
For subtle reasons, we perform a more rigorous estimation of the window function by comparing C_{ℓ}^{conv} with convolved power spectra of the input maps convolved with a symmetric Gaussian beam of comparable (but not necessarily exactly the same) size and then scaling the estimated window function accordingly.
Beam window functions are provided in the RIMO.
Beam window functions, W_{ℓ}, for the Planck mission.
File Names
The effective beams are provided in the PLAPlanck Legacy Archive as FITSFlexible Image Transfer Specification files containg HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). maps of the beams. For the file names the following convention is used.
 Single beam query, "beams_FFF_PixelNumber.fits" :
 FFF is the channel frequency (one of "30", "44", "70", "100", "143", "217", "353", "545", "857");
 "PixelNumber" is the number of the pixel to which the beam corresponds, 0–12 × N_{side}^{2}  1, with N_{side}=1024 for the LFI(Planck) Low Frequency Instrument and N_{side}=2048 for the HFI(Planck) High Frequency Instrument.
 Multiple beam query, "beams_FFF_FirstPixelNumberLastPixelNumber.zip" .
 The compressed files contains a set of files with the beams for the pixels covereing the selected region, with "FFF" as for sinle beam query and the naming convention for the beam files contained in the .zip file being the same as for single beam queries;
 FirstPixelNumber is the lowest pixel number for the area covered by the request;
 LastPixelNumber is the highest pixel number for the area covered by the request.
File format
The FITSFlexible Image Transfer Specification files provided in the PLAPlanck Legacy Archive contain HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). maps of the beams.
Previous Releases: (2015) and (2013) effective beams
2015 Release of effective beams
Product description 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. The effective beams are, by definition, the objects whose convolution with the true CMBCosmic Microwave background 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 ^{[1]}, and a discussion of its application to Planck data is given in the appropriate LFI(Planck) Low Frequency Instrument Planck2013IV^{[2]}, Planck2015A05^{[3]} and HFI(Planck) High Frequency Instrument Planck2013VII^{[4]} 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
We show here a comparison of the FEBeCoPderived 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(Planck) Low Frequency Instrument and HFI(Planck) High Frequency Instrument frequency channels, as an example. We show below a few panels of source images organized as follows:
 Row #1 DX9 images of four ERCSCEarly Release Compact Source Catalog 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(Planck) Low Frequency Instrument and BSpline beams,BS, Mars12 apodized for the CMBCosmic Microwave background channels and the BS Mars12 for the submm channels, for HFI(Planck) High Frequency Instrument (see section Inputs below).
 Row #3 log scale of #2; PSF isocontours shown in solid line, elliptical Gaussian fit isocontours shown in broken line
Histograms of the effective beam parameters Here we present histograms of the three fit parameters  beam FWHMFullWidthatHalfMaximum, ellipticity, and orientation with respect to the local meridian and of the beam solid angle. The sky is sampled (pretty sparsely) at 3072 directions which were chosen as HEALpix nside=16 pixel centers for HFI(Planck) High Frequency Instrument and at 768 directions which were chosen as HEALpix nside=8 pixel centers for LFI(Planck) Low Frequency Instrument. These uniformly sample the sky.
Where beam solid angle is estimated according to the definition:
(effbeam)/max(effbeam) i.e.,
Sky variation of effective beams solid angle and ellipticity of the bestfit Gaussian
 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.
We tabulate the simple statistics of FWHMFullWidthatHalfMaximum, ellipticity (e), orientation ( ) and beam solid angle, ( ), for a sample of 3072 and 768 directions on the sky for HFI(Planck) High Frequency Instrument and LFI(Planck) Low Frequency Instrument data respectively. Statistics shown in the Table are derived from the histograms shown above.
 The derived beam parameters are representative of the DPCData Processing Center NSIDE 1024 and 2048 healpix maps (they include the pixel window function).
 The reported FWHMFullWidthatHalfMaximum_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.
frequency  mean(fwhm) [arcmin]  sd(fwhm) [arcmin]  mean(e)  sd(e)  mean(  ) [degree]sd(  ) [degree]mean(  ) [arcmin ]sd(  ) [arcmin ]FWHMFullWidthatHalfMaximum_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 
  is the mean beam solid angle of the effective beam, where beam solid angle is estimated according to the definition: (effbeam)/max(effbeam), i.e. as an integral over the full extent of the effective beam, i.e. .
 from

 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).
is the beam solid angle estimated up to a radius equal to one and up to a radius equal to twice the .
we estimate the , under a Gaussian approximation  these are tabulated above

Band  [arcmin ]  spatial variation [arcmin  ][arcmin ]  spatial variation1 [arcmin  ][arcmin ]  spatial variation2 [arcmin  ]
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 
FEBeCoP, or Fast Effective Beam Convolution in Pixel space^{[1]}, 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 FWHMFullWidthatHalfMaximum 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 scan and beam models, as well as different lengths of datasets.
Pixel Ordered Detector Angles (PODA)
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 pixel  the Pixel Ordered Detector Angles (PODA) are made. This is an onetime process for each dataset. We used computers with large memory and used tedious memory management bookkeeping to make this step efficient.
effBeam 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 timeorderdata samples, while the effBeam part requires distribution of pixels among different computers.
Computational Cost
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.
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 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 CMBCosmic Microwave background 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.
The Focal Plane DataBase (FPDB)
The FPDB contains information on each detector, e.g., the orientation of the polarisation axis, different weight factors, (see the instrument RIMOs):
 HFI(Planck) High Frequency Instrument  The HFI RIMO
 LFI(Planck) Low Frequency Instrument  The LFI RIMO
The scanning strategy 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 (ODOperation Day definition is geometric visibility driven as it runs from the start of a DTCP (satellite Acquisition Of Signal) to the start of the next DTCP. Given the different ground stations and spacecraft will takes which station for how long, the OD duration varies but it is basically once a day.) 91 to ODOperation Day definition is geometric visibility driven as it runs from the start of a DTCP (satellite Acquisition Of Signal) to the start of the next DTCP. Given the different ground stations and spacecraft will takes which station for how long, the OD duration varies but it is basically once a day. 563 covering 3 surveys and half.
The scanbeam 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(Planck) Low Frequency Instrument: GRASP scanning beam  the scanning beams used are based on Radio Frequency Tuned Model (RFTM) smeared to simulate the inflight optical response.
 HFI(Planck) High Frequency Instrument: BSpline, BS based on 2 observations of Mars.
(see the instrument RIMOs).
Beam cutoff radii
N times geometric mean of FWHMFullWidthatHalfMaximum of all detectors in a channel, where N
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 
217353  4  18.782754 
submm  4  18.327635(545GHz) ; 17.093706(857GHz) 
 LFI(Planck) Low Frequency Instrument frequency channels for
 HFI(Planck) High Frequency Instrument frequency channels for
Related products Monte Carlo simulations
FEBeCoP software enables fast, fullsky 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 CMBCosmic Microwave background 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, and 1 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 CMBCosmic Microwave background 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 The Transfer Function or the Beam Window Function
relates the true angular power spectra with the observed angular power spectra :
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 MonteCarlo simulations. We generate several random realisations of the CMBCosmic Microwave background sky starting from a given fiducial , convolve the maps with the precomputed effective beams, compute the convolved power spectra , divide by the power spectra of the unconvolved map 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
File Names The effective beams are provided by the PLAPlanck Legacy Archive as FITSFlexible Image Transfer Specification files containg HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). maps of the beams. For the file names the following convention is used:
 Single beam query: beams_FFF_PixelNumber.fits
 FFF is the channel frequency (one of 30, 44, 70, 100, 143, 217, 353, 545, 857);
 PixelNumber is the number of the pixel to which the beam corresponds (LFI(Planck) Low Frequency Instrument , for the HFI(Planck) High Frequency Instrument ;  ). For the
 Multiple beam query: beams_FFF_FirstPixelNumberLastPixelNumber.zip
The compressed files contains a set of files with the beams for the pixels covereing the selected region. FFF as for sinle beam query. The naming convention for the beam files contained in the .zip file is the same as for single beam queries. FirstPixelNumber is the lowest pixel number for the area covered by the request;
 LastPixelNumber is the highest pixel number for the area covered by the request;
File format The FITSFlexible Image Transfer Specification files provided by the PLAPlanck Legacy Archive contain HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). maps of the beams.
2013 Release of effective beams
Product description 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. The effective beams are, by definition, the objects whose convolution with the true CMBCosmic Microwave background 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 ^{[1]}, and a discussion of its application to Planck data is given in the appropriate LFI(Planck) Low Frequency Instrument Planck2013IV^{[2]}, Planck2015A05^{[3]} and HFI(Planck) High Frequency Instrument Planck2013VII^{[4]} 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
We show here a comparison of the FEBeCoPderived 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(Planck) Low Frequency Instrument and HFI(Planck) High Frequency Instrument frequency channels, as an example. We show below a few panels of source images organized as follows:
 Row #1 DX9 images of four ERCSCEarly Release Compact Source Catalog 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(Planck) Low Frequency Instrument and BSpline beams,BS, Mars12 apodized for the CMBCosmic Microwave background channels and the BS Mars12 for the submm channels, for HFI(Planck) High Frequency Instrument (see section Inputs below).
 Row #3 log scale of #2; PSF isocontours shown in solid line, elliptical Gaussian fit isocontours shown in broken line
Histograms of the effective beam parameters Here we present histograms of the three fit parameters  beam FWHMFullWidthatHalfMaximum, ellipticity, and orientation with respect to the local meridian and of the beam solid angle. The sky is sampled (pretty sparsely) at 3072 directions which were chosen as HEALpix nside=16 pixel centers for HFI(Planck) High Frequency Instrument and at 768 directions which were chosen as HEALpix nside=8 pixel centers for LFI(Planck) Low Frequency Instrument. These uniformly sample the sky.
Where beam solid angle is estimated according to the definition:
(effbeam)/max(effbeam) i.e.,
Sky variation of effective beams solid angle and ellipticity of the bestfit Gaussian
 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.
We tabulate the simple statistics of FWHMFullWidthatHalfMaximum, ellipticity (e), orientation ( ) and beam solid angle, ( ), for a sample of 3072 and 768 directions on the sky for HFI(Planck) High Frequency Instrument and LFI(Planck) Low Frequency Instrument data respectively. Statistics shown in the Table are derived from the histograms shown above.
 The derived beam parameters are representative of the DPCData Processing Center NSIDE 1024 and 2048 healpix maps (they include the pixel window function).
 The reported FWHMFullWidthatHalfMaximum_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.
frequency  mean(fwhm) [arcmin]  sd(fwhm) [arcmin]  mean(e)  sd(e)  mean(  ) [degree]sd(  ) [degree]mean(  ) [arcmin ]sd(  ) [arcmin ]FWHMFullWidthatHalfMaximum_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 
  is the mean beam solid angle of the effective beam, where beam solid angle is estimated according to the definition: (effbeam)/max(effbeam), i.e. as an integral over the full extent of the effective beam, i.e. .
 from

 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).
is the beam solid angle estimated up to a radius equal to one and up to a radius equal to twice the .
we estimate the , under a Gaussian approximation  these are tabulated above

Band  [arcmin ]  spatial variation [arcmin  ][arcmin ]  spatial variation1 [arcmin  ][arcmin ]  spatial variation2 [arcmin  ]
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 
FEBeCoP, or Fast Effective Beam Convolution in Pixel space^{[1]}, 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 FWHMFullWidthatHalfMaximum 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 scan and beam models, as well as different lengths of datasets.
Pixel Ordered Detector Angles (PODA)
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 pixel  the Pixel Ordered Detector Angles (PODA) are made. This is an onetime process for each dataset. We used computers with large memory and used tedious memory management bookkeeping to make this step efficient.
effBeam 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 timeorderdata samples, while the effBeam part requires distribution of pixels among different computers.
Computational Cost
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.
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 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 CMBCosmic Microwave background 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.
The Focal Plane DataBase (FPDB)
The FPDB contains information on each detector, e.g., the orientation of the polarisation axis, different weight factors, (see the instrument RIMOs):
 HFI(Planck) High Frequency Instrument  The HFI RIMO
 LFI(Planck) Low Frequency Instrument  The LFI RIMO
The scanning strategy 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 (ODOperation Day definition is geometric visibility driven as it runs from the start of a DTCP (satellite Acquisition Of Signal) to the start of the next DTCP. Given the different ground stations and spacecraft will takes which station for how long, the OD duration varies but it is basically once a day.) 91 to ODOperation Day definition is geometric visibility driven as it runs from the start of a DTCP (satellite Acquisition Of Signal) to the start of the next DTCP. Given the different ground stations and spacecraft will takes which station for how long, the OD duration varies but it is basically once a day. 563 covering 3 surveys and half.
The scanbeam 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(Planck) Low Frequency Instrument: GRASP scanning beam  the scanning beams used are based on Radio Frequency Tuned Model (RFTM) smeared to simulate the inflight optical response.
 HFI(Planck) High Frequency Instrument: BSpline, BS based on 2 observations of Mars.
(see the instrument RIMOs).
Beam cutoff radii
N times geometric mean of FWHMFullWidthatHalfMaximum of all detectors in a channel, where N
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 
217353  4  18.782754 
submm  4  18.327635(545GHz) ; 17.093706(857GHz) 
 LFI(Planck) Low Frequency Instrument frequency channels for
 HFI(Planck) High Frequency Instrument frequency channels for
Related products Monte Carlo simulations
FEBeCoP software enables fast, fullsky 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 CMBCosmic Microwave background 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, and 1 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 CMBCosmic Microwave background 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 The Transfer Function or the Beam Window Function
relates the true angular power spectra with the observed angular power spectra :
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 MonteCarlo simulations. We generate several random realisations of the CMBCosmic Microwave background sky starting from a given fiducial , convolve the maps with the precomputed effective beams, compute the convolved power spectra , divide by the power spectra of the unconvolved map 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
File Names The effective beams are provided by the PLAPlanck Legacy Archive as FITSFlexible Image Transfer Specification files containg HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). maps of the beams. For the file names the following convention is used:
 Single beam query: beams_FFF_PixelNumber.fits
 FFF is the channel frequency (one of 30, 44, 70, 100, 143, 217, 353, 545, 857);
 PixelNumber is the number of the pixel to which the beam corresponds (LFI(Planck) Low Frequency Instrument , for the HFI(Planck) High Frequency Instrument ;  ). For the
 Multiple beam query: beams_FFF_FirstPixelNumberLastPixelNumber.zip
The compressed files contains a set of files with the beams for the pixels covereing the selected region. FFF as for sinle beam query. The naming convention for the beam files contained in the .zip file is the same as for single beam queries. FirstPixelNumber is the lowest pixel number for the area covered by the request;
 LastPixelNumber is the highest pixel number for the area covered by the request;
File format The FITSFlexible Image Transfer Specification files provided by the PLAPlanck Legacy Archive contain HEALPix([http://healpix.sourceforge.net Hierarchical Equal Area isoLatitude Pixelation of a sphere], <ref name="[[:Template:Gorski2005]]">[http://iopscience.iop.org/0004637X/622/2/759/pdf/0004637X_622_2_759.pdf '''HEALPix: A Framework for HighResolution 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''', 759771, (2005). maps of the beams.
References
 ↑ ^{1.0}^{1.1}^{1.2}^{1.3}^{1.4}^{1.5} Fast Pixel Space Convolution for Cosmic Microwave Background Surveys with Asymmetric Beams and Complex Scan Strategies: FEBeCoP, S. Mitra, G. Rocha, K. M. Górski, K. M. Huffenberger, H. K. Eriksen, M. A. J. Ashdown, C. R. Lawrence, ApJS, 193, 5+, (2011).
 ↑ ^{2.0}^{2.1}^{2.2} Planck 2013 results. IV. Low Frequency Instrument beams and window functions, Planck Collaboration, 2014, A&A, 571, A4.
 ↑ ^{3.0}^{3.1}^{3.2} Planck 2015 results. IV. LFI beams and window functions, Planck Collaboration, 2016, A&A, 594, A4.
 ↑ ^{4.0}^{4.1}^{4.2} Planck 2013 results. VII. HFI time response and beams, Planck Collaboration, 2014, A&A, 571, A7.