CMB from BICEP2/Keck-Array Oct 2018

Synopsis:Likelihood of BICEP2/Keck-Array, October 2018 (2015 data)
Author:BICEP2/Keck team and Antony Lewis

Warning

This is still in beta. Coming soon!

For now, use bicep_keck_2014 instead, which was used in the Planck 2015 analysis.

Note

If you use this likelihood, please cite it as:
Keck Array and BICEP2 Collaborations, BICEP2 / Keck Array x: Constraints on Primordial Gravitational Waves using Planck, WMAP, and New BICEP2/Keck Observations through the 2015 Season (arXiv:1810.05216)

Usage

To use this likelihood, bicep_keck_2015, you simply need to mention it in the likelihood block, or add them using the input generator.

The corresponding nuisance parameters will be added automatically, so you don’t have to care about listing them in the params block.

The nuisance parameters and their default priors are reproduced below, in the defaults yaml file.

You shouldn’t need to modify any of the options of this simple likelihood, but if you really need to, just copy the likelihood block into your input yaml file and modify whatever options you want (you can delete the rest).

# BICEP2/Keck Array September 2018 Data Products
# BICEP2/Keck Array X: Constraints On Primordial Gravitational Waves using Planck, WMAP, and New BICEP2/Keck Observations through the 2015 Season
# http://bicepkeck.org/
#
# File: BK15.ini
# Date: 2018-08-30
#
# BICEP2/Keck Array likelihood, including external data from Planck and WMAP.
#
# Default settings for data selection and foreground nuisance parameters from
# the baseline analysis described in BICEP2/Keck VI.
#
#   * Includes BICEP2+KECK 95, 150, and 220 GHz maps
#   * Additional maps are: Planck 30, 44, 70, 100, 143, 217, and 353 GHz;
#                          WMAP 23 GHz (K-band) and 33 GHz (Ka-band)
#   * By default uses B-modes only.
#   * Uses all nine ell bins.
#   * Gaussian priors on beta_dust and beta_sync.
#   * Dust blackbody temperature fixed to 19.6 K because this parameter
#     is degenerate with A_dust for our frequency range.
#   * Dust and sync decorrelation turned off by default.
#   * EE/BB ratio for dust and sync fixed to 2 (not relevant for analysis
#     with B modes only)
#   * Band center errors turned off by default.

likelihood:
  bicep_keck_2015:
    # Path to the BK15_cosmomc folder under the official data distribution
    path: null
    dataset_file: BK15_cosmomc/data/BK15/BK15_dust.dataset
    dataset_params:
      maps_use: [BK15_95_B, BK15_150_B, BK15_220_B, W023_B, P030_B, W033_B, P044_B, P070_B, P100_B, P143_B, P217_B, P353_B]
      use_min: 1
      use_max: 9
    field_names: [T, E, B, P]
    map_separator: x
    # Actually computed l_max (takes into account the lensing window function)
    l_max: 2700
    # Aliases for automatic covariance matrix
    renames: [BK15]
    # Speed in evaluations/second
    speed: 90

params:
  # dust power at ell=80, nu=353 GHz [uK^2]
  BBdust:
    prior:
      min: 0
      max: 15
    ref:
      dist: norm
      loc: 3
      scale: 0.1
    proposal: 0.1
    latex: A_{B,\mathrm{dust}}
  # sync power at ell=80, nu=23 GHz [uK^2]
  BBsync:
    prior:
      min: 0.0
      max: 50.0
    ref:
      dist: norm
      loc: 1.0
      scale: 1.0
    proposal: 1.0
    latex: A_{B,\mathrm{sync}}
  # dust spatial power spectrum power law index
  BBalphadust:
    prior:
      min: -1.0
      max: 0.0
    ref:
      dist: norm
      loc: -0.42
      scale: 0.01
    proposal: 0.01
    latex: \alpha_{B,\mathrm{dust}}
  # sync spatial power specturm power law index
  BBalphasync:
    prior:
      min: -1.0
      max: 0.0
    ref:
      dist: norm
      loc: -0.6
      scale: 0.01
    proposal: 0.01
    latex: \alpha_{B,\mathrm{sync}}
  # dust SED power law index
  BBbetadust:
    prior:
      dist: norm
      loc: 1.59
      scale: 0.11
    ref:
      dist: norm
      loc: 1.59
      scale: 0.02
    proposal: 0.02
    latex: \beta_{B,\mathrm{dust}}
  # sync SED power law index
  BBbetasync:
    prior:
      dist: norm
      loc: -3.1
      scale: 0.3
    ref:
      dist: norm
      loc: -3.1
      scale: 0.02
    proposal: 0.02
    latex: \beta_{B,\mathrm{sync}}
  # correlation between dust and sync
  # NOTE: For BK15 analysis, we are allowing this parameter to range
  #       from -1 to 1 (allowed range was 0 to 1 for BK14).
  BBdustsynccorr:
    prior:
      min: -1.0
      max: 1.0
    ref:
      dist: norm
      loc: 0.2
      scale: 0.01
    proposal: 0.01
    latex: \epsilon_{\mathrm{dust,sync}}
  # EE/BB ratios -- fixed / only used if E-modes are turned on
  EEtoBB_dust: 2
  #   latex: 2EE_{\mathrm{dust}}/BB_{\mathrm{dust}}
  EEtoBB_sync: 2
  #   latex: 2EE_{\mathrm{sync}}/BB_{\mathrm{sync}}
  # dust blackbody temperature [K] -- fixed / very insensitive to this
  BBTdust: 19.6
  #   latex: T_{\mathrm{dust}}
  # dust correlation ratio between 217 and 353 GHz, ell=80
  delta_dust: 1.0
  # sync correlation ratio between 23 and 33 GHz, ell=80
  delta_sync: 1.0
  # Band center errors, fixed to zero
  gamma_corr: 0.0
  gamma_95: 0.0
  gamma_150: 0.0
  gamma_220: 0.0

Installation

This likelihood can be installed automatically as explained in Installing cosmological codes and data. If are following the instructions there (you should!), you don’t need to read the rest of this section.

Manual installation of the BICEP2/Keck-Array likelihood

Assuming you are installing all your likelihoods under /path/to/likelihoods, simply do

$ cd /path/to/likelihoods
$ mkdir bicep_keck_2015
$ cd bicep_keck_2015
$ wget http://bicepkeck.org/BK15_datarelease/BK15_cosmomc.tgz
$ tar xvf BK15_cosmomc.tgz
$ rm BK15_cosmomc.tgz

After this, mention the path to this likelihood when you include it in an input file as

likelihood:
  bicep_keck_2015:
    path: /path/to/likelihoods/bicep_keck_2015