# kappa.py - plots weak lensing spectra
# by Bjoern Malte Schaefer

import numpy as np
import pylab as plt

plt.close()

data = np.loadtxt('spectrum.data')
l = data[:,0]
c = data[:,1]
aux = l * (l+1.0)/2.0/np.pi

plt.loglog(l,aux*c,'g-')

plt.xlabel('multipole $\ell$')
plt.ylabel('spectrum $\ell^2/(2\pi)C_\kappa(\ell)$')

plt.show()