Oct 17, 2018 - Lensing the CIB

From LSS Mocks
Revision as of 19:13, 17 October 2018 by Jlee (talk | contribs)
Jump to navigation Jump to search

CIB Lensing of Individual Shells and Analysis

Here is what I was able to do.

1. Lens the 23 CIB maps with the 23 kappa maps, creating a "cib_i^lensed" and a "cib_i^unlensed" set of 23 maps.

  I only used the field maps for both the kappa and the CIB maps. There were only 16 maps for the field (from z = 0 to z = 3.2). 
  Used the bash script Remi left. The following is the part of my script that did the lensing.
  It turned out this was where I made a mistake before and got such high values for the fractional difference between the lensed 
  and unlensed maps. I had been lensing a CIB map with the kappa map of the same range of z which does not make sense since a CIB shell 
  will only be lensed by the kappa before it. For example, CIB_z_1.0_1.2 has to be lensed by Kappa_z_0.8_1.0_zsource_0.9 but I was 
  lensing it with Kappa_z_1.0_1.2_zsource_1.1. 
     for i in seq 0.2 0.2 3.2;
  do
      echo "-----------$i----------"
      zmid=`awk "BEGIN {print $i-0.1}"`
      zfield=`awk "BEGIN {print $i+0.2}"`
      echo "$zmid"
      echo "$zfield"
      python lens.py ../jason_cib_lensing/maps/kappa/field/kappa_field_nside2048_zmin0.0_zmax${i}_zsource${zmid}_kap.fits 
      ../jason_cib_lensing/maps/cib/field/2048/unlensed_maps/cib_field_nside2048_zmin${i}_zmax${zfield}_cib.fits  
      ../jason_cib_lensing/maps/lensed/phi/run2/kap_zmax${i}_phi_field.fits 
      ../jason_cib_lensing/maps/lensed/field/run2/cib_field_nside2048_zmin${i}_zmax${zfield}_cib_lensed.fits -np 40

2. Plot cartview "zoom ins" of a few cib maps, lensed and unlensed, to make sure the lensing makes sense.

  I tried this out for z = 1.0 to 1.2 and z = 2.0 to 2.2, the whole map and some cartviews. The scales at the bottom were synchronized between
  lensed and unlensed maps.


   z = 1.0 to 1.2 
unlensed total
lensed total
  There isn't a clear difference between the two maps so I subtracted the lensed map from the unlensed (linearly) and made a map out of it.
  I will refer to these as 'difference maps.' (It may be different from the difference maps mentioned in the last meeting.)
difference total
                                  difference total
  Now there is a difference for sure.

unlensed -30º~-60ºlensed -30º~-60º

                      unlensed -30º~ -60º                                               lensed -30º~ -60º
  Again, the two maps look virtually the same and these cannot be differentiated even with blinking.

difference -30º~-60º

                                  difference -30º~ -60º
  This looks better.

unlensed 0º~20ºlensed 0º~20º

  The same goes for these maps even though the range is smaller.
                      unlensed 0º~ 20º                                               lensed 0º~ 20º

difference 0º~20º

                                  difference 0º~ 20º

unlensed 35º~40ºlensed 35º~40º

  When the maps are zoomed into 5º patches, the two maps certainly display a difference even though details are blurred due to the resolution.
                      unlensed 35º~ 40º                                               lensed 35º~ 40º


    z = 2.0 to 2.2 

unlensed totallensed total

                           unlensed total                                                          lensed total

difference total

                                  difference total

unlensed 0º~20ºlensed 0º~20º

                      unlensed 0º~ 20º                                               lensed 0º~ 20º

difference 0º~20º

                                  difference 0º~ 20º

unlensed 35º~40ºlensed 35º~40º

                      unlensed 35º~ 40º                                               lensed 35º~ 40º
  These two maps seem much different from each other compared to the 35º~40º for z = 1.0 to 1.2.

difference 35º~40º

                                  difference 35º~ 40º
  What luck!

3. add these 23 maps together to get a cib_total^lensed = \sum_i cib_i^lensed, and a cib_total^unlensed = \sum_i cib_i^unlensed

  Used another bash script to add the maps up. It goes:
  
  Lensed vs. unlensed.png    Lensed vs. unlensed log.png 
  Delta T.png    (log) 
                                                                   (log)
   Kappa map
     view of a patch from 5º~10ºN, 5º~10ºE 
      view of a patch from 5º~10ºN,
      5º~10ºE using 'cartview'
  


4. plot these, and also take their power spectrums and plot those just like you did for the CMB lensing (getting a delta Cl/Cl too)

  Fd.png    (log)
                                                                   (log)



5. power spectra of each of the 23 slices lensed and unlensed seperately to see whats going on