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...") |
|||
| (3 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 29: | 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)
