code: def main():
      
      dataset = 'RD'
  
      observables = ['log(<\sigma v>)', '\Omega_{\chi}h^2', 'log(\sigma_p^{SI})']
      var  = ['log(m_{\chi})']
      var += ['log(C_1)', 'log(C_2)', 'log(C_3)', 'log(C_4)', 'log(C_5)', 'log(C_6)']
      var += observables
  
      plot_vs_mass(dataset, observables, 'mass_vs_observables.png')
      plot_oneD(dataset, var, 'oneD.png')
      pairplot(dataset, var, 'pairplot.png')
docs: We here demostrate the basic functionality of barrett. We use a global scan
      of scalar dark matter as an example. The details aren't really important.