Difference between revisions of "Sep 26, 2018 - Excercise1"
Jump to navigation
Jump to search
| (2 intermediate revisions by the same user not shown) | |||
| Line 7: | Line 7: | ||
(ex. wget https://mocks.cita.utoronto.ca/data/lensing/cmb_temperature_lensed.fits) | (ex. wget https://mocks.cita.utoronto.ca/data/lensing/cmb_temperature_lensed.fits) | ||
| − | 2. | + | 2. Lens CMB to get a CMB lensed map |
Used the command listed in | Used the command listed in | ||
| Line 19: | Line 19: | ||
[[File:Lensed_vs._unlensed.png|300x225px|]] [[File:Lensed vs. unlensed log.png|300x225px|]] | [[File:Lensed_vs._unlensed.png|300x225px|]] [[File:Lensed vs. unlensed log.png|300x225px|]] | ||
| − | [[File:Delta T.png|300x225px|]] [[File:Delta T log.png|300x225px|]] | + | [[File:Delta T.png|300x225px|]] [[File:Delta T log.png|300x225px|(log)]] |
| − | + | (log) | |
| − | + | '''Kappa map''' | |
[[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]] | ||
view of a patch from 5º~10ºN, | view of a patch from 5º~10ºN, | ||
| Line 30: | Line 30: | ||
4. Take the power spectrum of CMB lensed and unlensed; plot this, and plot the fractional difference of the spectrums. | 4. Take the power spectrum of CMB lensed and unlensed; plot this, and plot the fractional difference of the spectrums. | ||
| − | [[File:Fd.png|300x225px|]] [[File:Fd log.png|300x225px|]] | + | [[File:Fd.png|300x225px|]] [[File:Fd log.png|300x225px|(log)]] |
| − | + | (log) | |
Latest revision as of 04:08, 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) 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)
