Difference between revisions of "Sep 26, 2018 - Excercise1"
Jump to navigation
Jump to search
(Created page with "'''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 th...") |
|||
| Line 20: | Line 20: | ||
[[File:Delta T.png|300x225px|]] [[File:Delta T log.png|300x225px|]] | [[File:Delta T.png|300x225px|]] [[File:Delta T log.png|300x225px|]] | ||
| + | (log) | ||
[[File:Kappa.png|250x250px|view of a patch from 5º~10ºN, 5º~10ºE]] | [[File:Kappa.png|250x250px|view of a patch from 5º~10ºN, 5º~10ºE]] | ||
| Line 30: | Line 31: | ||
[[File:Fd.png|300x225px|]] [[File:Fd log.png|300x225px|]] | [[File:Fd.png|300x225px|]] [[File:Fd log.png|300x225px|]] | ||
| + | (log) | ||
Revision as of 03:13, 27 September 2018
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)
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)
