# BK-style processing of DC1 maps to spectra

— J. Willmert, C. Bischoff

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

## Final pager

 Figure 1 Final Pager

### Legend for first half of the "Figure types":

• 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

### Legend for second half of the "Figure types":

#### spectra

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

#### bp devs

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

#### chi2

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

#### sum bp devs

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

#### uncertainties

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