Since an arbitrary beam is related to its SH Transform via
a circular beam has a radial profile (as returned by bl2beam)
where Pl is Legendre Polynomial and
is the beam window (or transfer) function, expected as input to bl2beam.
bl = gaussbeam(15.d0, 4000, 1) |
theta = dindgen(3000)/100. |
beam = bl2beam(bl, theta, /arcmin) |
plot, theta, beam |
the example above generates a beam window function (defined for all l in ) for a 15arcmin-FWHM gaussian beam, computes the beam profile for angles in [0,30] arcmin and then plots it.
Version 3.83, 2024-11-13