The power spectra presented here are analyzed from the input HEALPix maps using the S2HAT and PS2HAT libraries (see also arXiv:0903.2350). Custom MATLAB mex bindings interface to the (P)S2HAT library and permit us to finishing processing using the MATLAB pipeline used for standard BICEP/Keck analysis (i.e. our FFT steps are replaced with (P)S2HAT, but debiasing, suppression factor correction, etc follow the same code path as used in BK analysis; see for example arXiv:1403.3985 and arXiv:1603.05976 for more details on the BK pipeline).

Notes and caveats:

- Maps are those presented in [20170112_data_challenge_1].
- Direct BPWFs generated for a no-beam case, and then post-processed to apply a beam (before binning) for each frequency.
- Variable lensing + noise + foregrounds are the
`simd`simset. Some plots may still use the normal`sim`simset which always uses the \(A_L=1\) lensed-ΛCDM + noise sims. (I think this only happens in the first four Figure types made by`reduc_final`.) - “Real” data points are realization 0 of the signal+noise+foregrounds simset.
- Lensing theory line comes from

`aux_data/official_cl/camb_planck2013_r0_rev1_lensing.fits`

which is the updated spectrum, not the one used in our standard BK sims. - For 30GHz and 40GHz, the spectra diverge towards large negative values for \( \ell \gtrsim 280\). The cause is still being investigated. (First guess is that noise debiasing is off a bit, and the large suppression factor correction blows up the offset.)

Figure 1 | Final Pager |

- Thick red : ΛCDM unlensed scalars, zero for jacks
- Thick black : Real data
- Colored : ΛCDM lensed scalars + noise realizations

Note: Noise Debias done with LCDM unlensed scalars realizations

- thick red for all but BB : ΛCDM unlensed scalars, zero for jacks
- thick red for BB : ΛCDM tensors ($r=0.1$), zero for jacks
- thin gray : s+n simulations with ΛCDM lensed scalars as signal (sim type 7), debiased with ΛCDM unlensed scalars (sim type 2) and noise (sim type 6)
- thick gray : mean of sims (thin gray)
- red points with blue error bars : real data, debiased with the ΛCDM unlensed scalars and signflip noise (as s+n sims), uncertainties from the spread of the s+n sims (thin gray)

- thin gray : s+n sims
- thin red : 1 sigma contours
- thin green : 2 sigma contours
- thick black + markers : real data
- red markers : bins that go into chi2 analysis

Note: The contours mark the simulated bandpower in each bin closest to the respective percentile.

- red vertical : real data
- blue : distribution from s+n sims
- green : chi2 probability distribution fitted to blue

- red vertical : real data
- blue : distribution from s+n sims

- red : std of s+n sims
- blue : std of noise-only sims

Note: Remember that bandpower is a variance. Hence bandpower error bar is the std of a variance. This means noise and sample fluctuation components straight add, no quadrature — i.e. the difference between blue and red is the signal fluctatution component.