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.

def main():
    """ 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.
    """
    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')