Sep 26, 2018 - Excercise1

This is an exercise George and Pavel gave me to get an idea of working with lenspix and healpy.

1. Get CMB unlensed + kappa map from mocks.cita.utoronto.ca

  Used the 'wget' command 
  (ex. wget https://mocks.cita.utoronto.ca/data/lensing/cmb_temperature_lensed.fits)

2. Lens CMB to get a CMB lensed map

  Used the command listed in 
  (python lens.py ../cmb_temperature_unlensed.fits ../cmb_temperature_unlensed_alm.fits ../p_180921.fits 
  ../cmb_temperature_lensed.fits -np 8)

3. Plot CMB unlensed + cmb lensed with unit labels, also ΔT = lens-unlensed, also zoom in on 5º, 10º patch on Kappa map.

  Used 'healpy', 'anafast', 'matplotlib.pyplot' etc.

       
                                                                   (log)
   Kappa map
      
      view of a patch from 5º~10ºN,
      5º~10ºE using 'cartview'
  

4. Take the power spectrum of CMB lensed and unlensed; plot this, and plot the fractional difference of the spectrums.

      
                                                                   (log)