Oct 17, 2018 - Lensing the CIB

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.)

  Now there is a difference for sure.

  Again, the two maps look virtually the same and these cannot be differentiated even with blinking.

difference -30º~-60º

  This looks better.

  The same goes for these maps even though the range is smaller.

difference 0º~20º

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

    z = 2.0 to 2.2 

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

  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:
  

       
                                                                   (log)
   Kappa map
      
      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)

      
                                                                   (log)

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