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