alm2dl

Description

alm2dl reads $a_{{\ell}m}$ and calculates $C_{{\ell}}$ and $D_{{\ell}}$ where

$$C_{{\ell}}=\frac{1}{2{\ell}+1}\sum_{m=-{\ell}}^{{\ell}} \vert a_{{\ell}m}\vert^2.$$ (1)

$$D_{{\ell}}=\frac{ {\ell}({\ell}+1)C_{{\ell}}}{2 \pi T_0^2}.$$ (2)

where $T_0 = 2.726$K, $C_{{\ell}}$ is output in unit of $T^2$ and $D_{{\ell}}$ is unitless. alm2dl can read the input $a_{{\ell}m}$ file either in FITS or in ASCII formats. The FITS binary format should contain 3 fields: index= ${\ell}\times {\ell}+{\ell}+m+1$ (in HEALPix definition), $a_{{\ell}m}$ real part ($\Re$) and imaginary part ($\Im$). The ASCII format should contain 4 entries : L, M, Re, Im, where L is the ${\ell}$-multipole and M is the $m$-multipole. It can also read data from a single 3-$a_{{\ell}m}$ file ($a_{{\ell}m}$, $e_{{\ell}m}$ and $b_{{\ell}m}$) to calculate correlation spectra BB, EE, TE or TB between $a_{{\ell}m}$ and $e,b_{{\ell}m}$-coefficients. If the input file contains only $a_{{\ell}m}$-coefficients, to calculate TE or TB-spectra, a file with $e_{lm}$ or $b_{lm}$-coefficients should be added in the command string.

alm2dl calculates $a_{{\ell}m}$ amplitudes: $\vert a_{{\ell}m}\vert$ and phases defined as $\arctan(\Im/\Re)$ in $[0,2\pi]$.

alm2dl also lists the values of the FITS $a_{{\ell}m}$ coefficients, which can be output either to a file or to the standard output device for display.

Examples

1. Calculating $C_{{\ell}}$ or $D_{{\ell}}$:
   • alm2dl alm.txt -lmax 25
     reads ASCII ' alm.txt' and displays $D_{{\ell}}$ up to 25.

   • alm2dl -cl alm.fits -lmax 1000 -o cl.txt
     calculates $C_{{\ell}}$ up to ${\ell}=1000$ from ' alm.fits' and outputs in 'cl.txt'.

2. Listing $a_{{\ell}m}$:
   • alm2dl -g alm.fits | less
     reads FITS file ' alm.fits' and displays in standard output device the $a_{{\ell}m}$ coefficients in 4 columns: ${\ell}$, $m$, real and imaginary parts.

   • alm2dl -p alm.fits -o nlmalm.txt
     reads ' alm.fits' and outputs $a_{{\ell}m}$ and phases to file ' nlmalm.txt' with 6 entries: nlm index (= ${\ell}\times {\ell}+{\ell}+m+1$), ${\ell}$, $m$, real, imaginary part and the phase.

   • alm2dl 3alm.fits -g -e -lmax 50 -o elm.dat
     reads a file ' 3alm.fits', and extracts data of E-mode writing them in ASCII format to file ' elm.at' up to ${\ell}=50$.

   • alm2dl 3alm.fits -b -lmax 50 -o blm.fts
     reads a file ' 3alm.fits', and extracts data of E-mode writing them in FITS to file ' blm.fts' up to ${\ell}=50$.

3. Calculating amplitudes:
   • alm2dl -amp alm.fits -lmax 10
     reads ' alm.fits' and calculates amplitudes up to ${\ell}=10$.

   • alm2dl -amp -a alm.txt -lmax 50 -o amp.txt
     reads ASCII ' alm.txt' and outputs ASCII amplitudes up to ${\ell}=50$ in ' amp.txt'.

4. Calculating phases:
   • alm2dl -ph alm.fits -lmax 10
     reads ' alm.fits' and displays phases up to ${\ell}=10$.

5. Calculating polarizatiom modes $B_{{\ell}m}$ and $E_{{\ell}m}$ :
   • alm2dl 3alm.fits -tbl -lmax 50 -o tbl.txt
     reads a file ' 3alm.fits', calculates T-B correlation spectrum and writes ASCII coefficients up to ${\ell}=50$ to ' tbl.txt'.

   • alm2dl alm.fits elm.fits -tel -lmax 50 -o tel.txt
     reads files ' alm.fits' and ' elm.fits', calculates T-E correlation spectrum and writes ASCII coefficients up to ${\ell}=50$ to ' tel.txt'.

   • alm2dl 3alm.fits -e -lmax 50 -o elm.fts
     reads a file ' 3alm.fits', and extracts data of E-mode writing them in FITS format to file ' elm.fts' up to ${\ell}=50$.