# plot_cmb.py - plots the output from CAMB
# by Bjoern Malte Schaefer, bjoern.malte.schaefer@uni-heidelberg.de

import numpy as np
import pylab as plt

lmin = 0
lmax = 1800

plt.close()

data = np.loadtxt('camb_cmb.data')

l = data[lmin:lmax,0]
cu = data[lmin:lmax,1]

su = cu / l / (l+1.0) * 2.0 * np.pi

plt.plot(l,cu,'b-',label='unlensed')

data = np.loadtxt('camb_lensed.data')

l = data[lmin:lmax,0]
cl = data[lmin:lmax,1]

sl = cl / l / (l+1.0) * 2.0 * np.pi

plt.plot(l,cl,'r-',label='lensed')
plt.legend(loc='upper right')

plt.xlabel('multipole $\ell$')
plt.ylabel('temperature spectrum $\ell(\ell+1) C(\ell)/(2\pi)$')

plt.show()