====== Beam sidelobes simulations with Conviqt ======


===== Software versions =====

See [[https://github.com/zonca/pico-simulations/tree/master/conviqt_sim/runs/201810_boresight_1pix4det_conviqt_realbeams_mask5deg_1y|Simulation page on Github]]

===== Configuration =====

  * 1 pixel in band 12 at 154.80 GHz
  * 4 detectors, 0,90,45,135 degrees
  * realistic beams with center masked (''center_150GHz_mask400px5d_l720_alm.fits'')
  * boresight pointing
  * PySM input sky
  * 1 year

===== Discussion and plots =====

  * [[https://github.com/zonca/pico-simulations/issues/24|Github issue with discussion]]
  * [[https://github.com/zonca/pico-simulations/pull/26/files|Github pull request with implementation]]
  * [[https://gist.github.com/zonca/9f195e82e3e887f770b4ede52ad0e0fc|Notebook with map plots]]
  * Also scattered among daily updates on the systematic effects [[systematicswg|telecon page]]


===== Summary plots =====

  * 4 pi beams, main beam is at the plot center: {{::4pi_beam_moll_unclipped.png?direct&100| full sky beam}}
    * These beams had the central 10 deg peak removed and filled with the value of the neighboring pixels. This clipped beam was then used in conviqt. Original and clipped x-pol beam: {{::fsl_slice_xpol_w-clipped_beam.png?direct&100|}} y-pol is similar.
  * Input Galaxy maps, smoothed to 10 deg:{{::conviqt_in_201810_smooth10deg_iqu.png?direct&100|input to Conviqt, smoothed 10deg gaussian}}
  * Ouput maps from Conviqt nside=512:{{::conviqt_out_4detect_20181025_smooth0.1d_iqu.png?direct&100|output of conviqt, 4 detectors}}
  * Power spectrum of far sidelobe map masked with 43% Planck mask:{{::fsl_cls_41prct.png?direct&100|}}
    * Mask:{{::43prct_planck_mask.png?direct&100|}}




===== Partial cleaning =====

To give some idea of what level of removal is necessary I multiplied the sidelobe by a scaling factor, alpha, and added noise then recalculated the Cl's.  I based the scaling factor on the map noise level for these 155 GHz detectors.

Possible noise levels to compare sidelobes to:
  * Depth for 1 year for 1 wafer is 0.1 uK deg.  
  * Depth for 5 year for 1 wafer is 0.05 uK deg.  
  * Depth for 5 year for all 6 wafers is 0.02 uK deg.  
  * Depth for the full mission is 0.01 uK deg.

Intensity of sidelobes, all units are uK, Q and U are the median and max of absolute value of Q, U:
|                     |^  I    ^  Q     ^  U     ^
|  Masked map  |  Mean  | 0.10  | 0.018  | 0.015  |
| :::          |  Max  | 0.21  | 0.12   | 0.11   |
|  Full map    |  Mean  | 0.28  | 0.022  | 0.014  |
| :::          |  Max  | 2.6   | 0.58   | 0.14   |

When looking at the masked map, which is what makes sense for B-modes as this is the cleanest 40% of the sky, the average polarized sidelobe level is less than the map noise.  This takes away the simplest option of reducing the sidelobe power until its average is equal to the noise level.  Instead I reduced the sidelobe power until the max in the maxed region was approximately at the noise level.  This gives the following power spectra.

  * Power spectra for masked FSL, no noise added. {{::fsl_cls_41prct_sky_no_noise_20181206.png?direct&100|FSL, no noise}}
  * Power spectra for masked FSL*0.6 + 0.05 uK deg noise: {{::fsl_cls_41prct_sky_a0.6_0.05ukdeg_noise_20181206.png?direct&100|}}  (1 wafer 5 years)
  * Power spectra for masked FSL*0.25 + 0.02 uK deg noise: {{::fsl_cls_41prct_sky_a0.25_0.02ukdeg_noise_20181206.png?direct&100|}}  (6 wafers 5 years)
  * Power spectra for masked FSL*0.10 + 0.01 uK deg noise: {{::fsl_cls_41prct_sky_a0.10_0.01ukdeg_noise_20181206.png?direct&100|}} (Full CBE mission)

For these power spectra I multiplied by the FSL map by alpha, added white noise where each map pixel value is drawn from a gaussian with sigma equal to the quoted noise level (scaled to appropriate pixel size, nside=512), and applied the galaxy mask. The power spectra are calculated using the healpy anafast function.

**Bottom line:**  We need to remove about 60% of the sidelobes to get to r of 10^-3 and about 90% to reach r = 10^-4
===== Results location =====

NERSC

''/global/cscratch1/sd/zonca/pico/conviqt_sim/201810_boresight_1pix4det_conviqt_realbeams_mask5deg_1y''

tape

''~zonca/pico/conviqt_sim/201810_boresight_1pix4det_conviqt_realbeams_mask5deg_1y''