LSS Mocks - User contributions [en-ca]

CIB lensing

2018-08-24T19:35:26Z

Remim:

Final picture, the CIB is divided in slices,each slice is lensed by the appropriate kappa shell. Then all the slices are added up

[[File:Cib totalfield powerspectrum err norm.png|200px]]

and smoothed :

[[File:Cib totalfield powerspectrum sm.png|200px]]

This was done with field only, Cal_sky was not working with halos cutout.

Here is the fractional difference :

[[File:Cib totalfield powerspectrum reldiff.png|200px]]

With smoothing and Planck data :

[[File:Cib totalfield powerspectrum reldiff err.png|200px]]

File:Cib totalfield powerspectrum sm.png

2018-08-24T19:35:04Z

Remim: Lensed-unlensed CIB slice by slice, then added up

Lensed-unlensed CIB slice by slice, then added up

CIB lensing

2018-08-24T19:32:11Z

Remim:

Final picture, the CIB is divided in slices,each slice is lensed by the appropriate kappa shell. Then all the slices are added up

[[File:Cib totalfield powerspectrum err norm.png|200px]]

This was done with field only, Cal_sky was not working with halos cutout.

Here is the fractional difference :

[[File:Cib totalfield powerspectrum reldiff.png|200px]]

With smoothing and Planck data :

[[File:Cib totalfield powerspectrum reldiff err.png|200px]]

File:Cib totalfield powerspectrum reldiff err.png

2018-08-24T19:31:33Z

Remim: Relative difference of lensed-unlensed CIB + Planck data

Relative difference of lensed-unlensed CIB + Planck data

CIB lensing

2018-08-24T19:05:33Z

Remim:

Final picture, the CIB is divided in slices,each slice is lensed by the appropriate kappa shell. Then all the slices are added up

[[File:Cib totalfield powerspectrum err norm.png|200px]]

This was done with field only, Cal_sky was not working with halos cutout.

Here is the fractional difference :

[[File:Cib totalfield powerspectrum reldiff.png|200px]]

File:Cib totalfield powerspectrum reldiff.png

2018-08-24T19:05:15Z

Remim: Fractional difference between lensed CIB and unlensed

Fractional difference between lensed CIB and unlensed

CIB lensing

2018-08-24T19:00:14Z

Remim:

Final picture, the CIB is divided in slices,each slice is lensed by the appropriate kappa shell. Then all the slices are added up

[[File:Cib totalfield powerspectrum err norm.png|200px]]

File:Cib powerspectrum.png

2018-08-24T18:59:47Z

Remim:

Final CIB, divided in slices, each slice lensed by a kappa shell, then the whole is added up

[[File:Cib totalfield powerspectrum err norm.png|200px]]

File:Cib totalfield powerspectrum err norm.png

2018-08-24T18:59:37Z

Remim: CIB lensed slice by slice by an appropriate lensing kernel.

CIB lensed slice by slice by an appropriate lensing kernel.

CIB lensing

2018-08-24T18:49:04Z

Remim:

Final picture, the CIB is divided in slices,each slice is lensed by the appropriate kappa shell. Then all the slices are added up

[[File:Cib powerspectrum.png|200px]]

CIB lensing

2018-08-24T18:48:09Z

Remim: Created page with "File:Cib powerspectrum.png"

[[File:Cib powerspectrum.png]]

File:Cib powerspectrum.png

2018-08-24T18:47:57Z

Remim: Final CIB, divided in slices, each slice lensed by a kappa shell, then the whole is added up

Final CIB, divided in slices, each slice lensed by a kappa shell, then the whole is added up

Logbook

2018-08-24T18:46:08Z

Remim:

Here is the research log for Remi's CIB lensing project. Links to research notes are below.

[[Final Work - Remi-links]]

[[CIB lensing]]

[[August 17 - Cal-Sky]]

[[August 9 - Cal-Sky - CIB 2]]

[[July 20 - Cal-Sky - CIB]]

[[June 11 - CIB - Lensing]]

[[June 06 - CMB, Websky and LensPix]]

[[June 05 - LensPix parameters]]

[[May 14 - Issue with lensing]]

[[Apr 26 - First look at lensing sims]]

Final Work - Remi-links

2018-08-24T18:45:41Z

Remim:

Lenspix, the code to lens healpix maps given an input kappa map, can be found here:

https://github.com/louis-pham/modified_lenspix

cal-sky, the code to create "field" maps of kappa, cib, etc, can be found here, where scripts to break up maps into redshift slices can be found in George Stein's fork:

https://github.com/marcelo-alvarez/cal-sky

lenspower, the code to lens power spectrum theoretically, compute the change in correlation function. Can be found in Remi Mourgues's fork

https://github.com/CITA/lenspower

Final Work - Remi-links

2018-08-24T18:45:23Z

Remim:

August 17 - Cal-Sky

2018-08-22T20:53:35Z

Remim:

Implemented Halos + field for CIB and Kappa, slices of delta z = 0.2

Lensed with LensPix vs Theory

[[File:Cib total reldiff log nf.png|200px]]

Kappa generated vs CAMB

[[File:Kappa split halofits+lin log nf.png|200px]]

A closer look :

[[File:Kappa one halofits+lin log nf.png|200px]]

If we take a look at the maps : 2 deg side
Zmin = 1.4
[[File:1deg un-len zmin1.4 nosmth.png|200px]]

Zmin = 3.4
[[File:1deg un-len zmin3.4 nosmth.png|200px]]

The points on the 1.4 slice are halos.

August 17 - Cal-Sky

2018-08-22T20:53:11Z

Remim:

Implemented Halos + field for CIB and Kappa, slices of delta z = 0.2

Lensed with LensPix vs Theory

[[File:Cib total reldiff log nf.png|200px]]

Kappa generated vs CAMB

[[File:Kappa split halofits+lin log nf.png|200px]]

A closer look :

[[File:Kappa one halofits+lin log nf]]

If we take a look at the maps : 2 deg side
Zmin = 1.4
[[File:1deg un-len zmin1.4 nosmth.png|200px]]

Zmin = 3.4
[[File:1deg un-len zmin3.4 nosmth.png|200px]]

The points on the 1.4 slice are halos.

File:Kappa one halofits+lin log nf.png

2018-08-22T20:52:39Z

Remim: Kappa power spectrum for one redshift Linear theory + non linear corrections

Kappa power spectrum for one redshift
Linear theory + non linear corrections

August 17 - Cal-Sky

2018-08-22T20:50:56Z

Remim:

Implemented Halos + field for CIB and Kappa, slices of delta z = 0.2

Lensed with LensPix vs Theory

[[File:Cib total reldiff log nf.png|200px]]

Kappa generated vs CAMB

[[File:Kappa split halofits+lin log nf.png|200px]]

If we take a look at the maps : 2 deg side
Zmin = 1.4
[[File:1deg un-len zmin1.4 nosmth.png|200px]]

Zmin = 3.4
[[File:1deg un-len zmin3.4 nosmth.png|200px]]

The points on the 1.4 slice are halos.

File:Kappa split halofits+lin log nf.png

2018-08-22T20:50:29Z

Remim: Kappa power spectrum for different slices Linear theory + non-linear corrections

Kappa power spectrum for different slices
Linear theory + non-linear corrections

August 17 - Cal-Sky

2018-08-22T16:50:49Z

Remim:

Implemented Halos + field for CIB and Kappa, slices of delta z = 0.2

Lensed with LensPix vs Theory

[[File:Cib total reldiff log nf.png|200px]]

Kappa generated vs CAMB

[[File:Kappa split log nf.png|200px]]

If we take a look at the maps : 2 deg side
Zmin = 1.4
[[File:1deg un-len zmin1.4 nosmth.png|200px]]

Zmin = 3.4
[[File:1deg un-len zmin3.4 nosmth.png|200px]]

The points on the 1.4 slice are halos.

File:Kappa split log nf.png

2018-08-22T16:50:14Z

Remim: Kappa shells: simulation vs theory

Kappa shells: simulation vs theory

August 17 - Cal-Sky

2018-08-17T21:18:03Z

Remim:

Implemented Halos + field for CIB and Kappa, slices of delta z = 0.2

Lensed with LensPix vs Theory

[[File:Cib total reldiff log nf.png|200px]]

Kappa generated vs CAMB

[[File:Kappa total log nf.png|200px]]

If we take a look at the maps : 2 deg side
Zmin = 1.4
[[File:1deg un-len zmin1.4 nosmth.png|200px]]

Zmin = 3.4
[[File:1deg un-len zmin3.4 nosmth.png|200px]]

The points on the 1.4 slice are halos.

TO DO next : smoothe the maps

File:1deg un-len zmin3.4 nosmth.png

2018-08-17T21:15:08Z

Remim: Unlensed CIB - Lensed CIB Zmin 3.4 2 deg side

Unlensed CIB - Lensed CIB
Zmin 3.4
2 deg side

File:1deg un-len zmin1.4 nosmth.png

2018-08-17T21:13:44Z

Remim: Unlensed CIB minus Lensed CIB Slice 1.4-1.6 2 deg side

Unlensed CIB minus Lensed CIB
Slice 1.4-1.6
2 deg side

August 17 - Cal-Sky

2018-08-17T21:03:26Z

Remim: Created page with "Implemented Halos + field for CIB and Kappa, slices of delta z = 0.2 Lensed with LensPix vs Theory 200px Kappa generated vs CAMB F..."

File:Kappa total log nf.png

2018-08-17T21:02:28Z

Remim: Kappa shells : cal-sky vs camb

Kappa shells : cal-sky vs camb

File:Cib total reldiff log nf.png

2018-08-17T21:01:18Z

Remim: Lensed CIB (halos+field) by slices

Lensed CIB (halos+field) by slices

Logbook

2018-08-17T21:00:12Z

Remim:

Here is a place for a research log. Please put links to research notes below.

[[August 17 - Cal-Sky]]

[[August 9 - Cal-Sky - CIB 2]]

[[July 20 - Cal-Sky - CIB]]

[[June 11 - CIB - Lensing]]

[[June 06 - CMB, Websky and LensPix]]

[[June 05 - LensPix parameters]]

[[May 14 - Issue with lensing]]

[[Apr 26 - First look at lensing sims]]

August 9 - Cal-Sky - CIB 2

2018-08-13T19:35:28Z

Remim: Created page with "Did the same as previous for a larger set of slices : zmin = 0.2 -> 4.8 with delta z = 0.2 200px Computed also the kappa shells zith zsource = zm..."

Did the same as previous for a larger set of slices : zmin = 0.2 -> 4.8 with delta z = 0.2

[[File:Reldiff log nf.png|200px]]

Computed also the kappa shells zith zsource = zmax-0.1, and compared them to simulation by CAMB

[[File:Kappa theory log nf.png|200px]]

File:Kappa theory log nf.png

2018-08-13T19:30:00Z

Remim: Kappa power spectrum for different slices of Cal-Sky simulation, compared to CAMB

Kappa power spectrum for different slices of Cal-Sky simulation, compared to CAMB

File:Reldiff log nf.png

2018-08-13T19:28:05Z

Remim: Fractional difference of the lensed CIB power spectrum for simulation and theory (lenspower)

Fractional difference of the lensed CIB power spectrum for simulation and theory (lenspower)

Logbook

2018-08-13T19:19:37Z

Remim:

Here is a place for a research log. Please put links to research notes below.

[[August 9 - Cal-Sky - CIB 2]]

[[July 20 - Cal-Sky - CIB]]

[[June 11 - CIB - Lensing]]

[[June 06 - CMB, Websky and LensPix]]

[[June 05 - LensPix parameters]]

[[May 14 - Issue with lensing]]

[[Apr 26 - First look at lensing sims]]

July 20 - Cal-Sky - CIB

2018-07-30T20:44:46Z

Remim:

Ran Cal-Sky with fields, with a box of redshift depth equal to 0.5.
Computed slices of CIB, with a delta z of 0.1.
Lensed those slices by kappa shells from redshift 0 to zi, where CIB is from zi to zi+0.1.

[[File:Calsky lenspix theory.png|400px]]

There is a mismatch between theoretical computation (lenspower) and simulation (lenspix).

Ran the same thing but with both now GRF, simulated from the Cal-Sky power spectrum of CIB/kappa with Healpy

[[File:Calsky lenspix theory bGRF.png|400px]]

There is still a mismatch.

Here is the statistics of kappa fields

[[File:Stats kappa.png|400px]]

and the same with GRF

[[File:Stats kappa GRF.png|400px]]

Finally, theoretical computation of kappa power spectrum (pycamb) vs cal-sky simulation

[[File:Kap cambvscalsky.png|400px]]

UPDATE :

Ran lensing with zero kappa, GRF CIB

[[File:GRF kapzero.png|400px]]

NEXT : run for more slices and deeper redshifts, halos included

File:GRF kapzero.png

2018-07-30T20:43:59Z

Remim: Lensing of GRF CIB by zero kappa

Lensing of GRF CIB by zero kappa

July 20 - Cal-Sky - CIB

2018-07-23T18:48:44Z

Remim: Created page with "Ran Cal-Sky with fields, with a box of redshift depth equal to 0.5. Computed slices of CIB, with a delta z of 0.1. Lensed those slices by kappa shells from redshift 0 to zi, w..."

File:Kap cambvscalsky.png

2018-07-23T18:37:05Z

Remim: Theoretical computation of kappa power spectrum by pycamb vs simulation

Theoretical computation of kappa power spectrum by pycamb vs simulation

File:Stats kappa GRF.png

2018-07-23T18:35:47Z

Remim: Statistics of kappa field from Cal-Sky : histogram - qq plot - values from d'Agostino and Pearson's test

Statistics of kappa field from Cal-Sky : histogram - qq plot - values from d'Agostino and Pearson's test

File:Stats kappa.png

2018-07-23T18:34:25Z

Remim: Statistics of kappa field : histogram - qqplot and some statistical values from D'Agostino and Person's test

Statistics of kappa field : histogram - qqplot and some statistical values from D'Agostino and Person's test

File:Calsky lenspix theory bGRF.png

2018-07-23T18:32:04Z

Remim: Lensing of CIB slices by kappa shells, both are GRF

Lensing of CIB slices by kappa shells, both are GRF

File:Calsky lenspix theory.png

2018-07-23T18:29:51Z

Remim: Lensing of CIB slices by kappa shells.

Lensing of CIB slices by kappa shells.

Logbook

2018-07-23T18:27:00Z

Remim:

Here is a place for a research log. Please put links to research notes below.

[[July 20 - Cal-Sky - CIB]]

[[June 11 - CIB - Lensing]]

[[June 06 - CMB, Websky and LensPix]]

[[June 05 - LensPix parameters]]

[[May 14 - Issue with lensing]]

[[Apr 26 - First look at lensing sims]]

June 11 - CIB - Lensing

2018-06-12T15:22:35Z

Remim: Created page with "First results of CIB lensing with Lenspix and Lenspower File:Cib ps reldiff kappa-cmbhalf 2.png"

First results of CIB lensing with Lenspix and Lenspower

[[File:Cib ps reldiff kappa-cmbhalf 2.png]]

File:Cib ps reldiff kappa-cmbhalf 2.png

2018-06-12T15:22:22Z

Remim: Relative difference of power spectrum for different lensing methods

Relative difference of power spectrum for different lensing methods

Logbook

2018-06-12T15:21:23Z

Remim:

Here is a place for a research log. Please put links to research notes below.

[[June 11 - CIB - Lensing]]

[[June 06 - CMB, Websky and LensPix]]

[[June 05 - LensPix parameters]]

[[May 14 - Issue with lensing]]

[[Apr 26 - First look at lensing sims]]

June 06 - CMB, Websky and LensPix

2018-06-07T18:38:16Z

Remim:

I took the map from the Websky simulation (those on this website) and did :

- lensing of unlensed CMB with kappa from Websky, with two different parameters for LensPix: the old and one of the best obtained according to the results of yesterday;

- lensing of unlensed CMB with kappa from Websky, using Joel's & Alex's code.

Then I plotted the relative difference for each lensed map (with filter), the result is below.

There is a difference between the Websky map and the other, if the theory calculus is the right one, with no surprise, the best parameter for LensPix (interpolation factor = 6) is near the theory's computation.

UPDATE : I included CAMB simulation, see below.

[[File:Reldiff web theo lenspix 2.png]]

File:Reldiff web theo lenspix 2.png

2018-06-07T18:37:36Z

Remim:

June 06 - CMB, Websky and LensPix

2018-06-06T21:17:45Z

Remim:

I took the map from the Websky simulation (those on this website) and did :

- lensing of unlensed CMB with kappa from Websky, with two different parameters for LensPix: the old and one of the best obtained according to the results of yesterday;

- lensing of unlensed CMB with kappa from Websky, using Joel's & Alex's code.

Then I plotted the relative difference for each lensed map (with filter), the result is below.

There is a difference between the Websky map and the other, if the theory calculus is the right one, with no surprise, the best parameter for LensPix (interpolation factor = 6) is near the theory's computation.

I will compare those results with CAMB later, I just need to have the exact parameters used in the Peak Patch simulation, according to CAMB nomenclature.

[[File:Reldiff web theo lenspix.png]]

June 05 - LensPix parameters

2018-06-06T21:16:49Z

Remim:

I changed the parameters for LensPix in order to see what changed in the power spectrum.
Those plots have been obtained with a zero field for kappa.
The 'theory - code' on the first plot is a code from Joel/Alex, similar to CAMB.

Loss of energy here means the sum of Cl over all l.

I will check the mocks on this site, and see if we can increase the agreement with theoretical plots.

[[File:Interp factor comparison n2048.png|200px|thumb|left|Comparison of time execution / loss of energy when lensing with LensPix for different parameters : interp_factors and MPI division method.]]

[[File:Reldiff interp comp.png|200px|thumb|left|Relative difference of power spectrum for different values of LensPix's parameters.]]

[[File:Reldiff interp comp lin.png|200px|thumb|left|Relative difference of power spectrum for different values of LensPix's parameters : same scale and high l.]]

File:Reldiff interp comp lin.png

2018-06-06T21:16:13Z

Remim: Relative difference of power spectrum for different values of LensPix's parameters, same scale and high l.

Relative difference of power spectrum for different values of LensPix's parameters, same scale and high l.