Modelling_NumberCounts1D_Redshift.cpp

Go to the documentation of this file.

/********************************************************************
 *  Copyright (C) 2016 by Federico Marulli and Alfonso Veropalumbo  *
 *  federico.marulli3@unibo.it                                      *
 *                                                                  *
 *  This program is free software; you can redistribute it and/or   * 
 *  modify it under the terms of the GNU General Public License as  *
 *  published by the Free Software Foundation; either version 2 of  *
 *  the License, or (at your option) any later version.             *
 *                                                                  *
 *  This program is distributed in the hope that it will be useful, *
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of  *
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the   *
 *  GNU General Public License for more details.                    *
 *                                                                  *
 *  You should have received a copy of the GNU General Public       *
 *  License along with this program; if not, write to the Free      *
 *  Software Foundation, Inc.,                                      *
 *  59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.       *
 ********************************************************************/
 
#include "ModelFunction_NumberCounts1D_Redshift.h"
#include "Modelling_NumberCounts1D_Redshift.h"
 
using namespace std;
 
using namespace cbl;
 
using namespace modelling::numbercounts;
 
 
// ===========================================================================================
 
 
void cbl::modelling::numbercounts::Modelling_NumberCounts1D_Redshift::set_model_NumberCounts_cosmology (const std::vector<cbl::cosmology::CosmologicalParameter> cosmo_param, const std::vector<statistics::PriorDistribution> cosmo_param_prior)
{
  m_data_model.Cpar = cosmo_param;
 
  const size_t nParams = cosmo_param.size();
  checkDim(cosmo_param, nParams, "cosmoPar_prior");
 
  vector<statistics::ParameterType> cosmoPar_type(nParams, statistics::ParameterType::_Base_);
  vector<string> cosmoPar_string(nParams);
 
  for (size_t i=0; i<nParams; i++)
    cosmoPar_string[i] = CosmologicalParameter_name(cosmo_param[i]);
 
  // input data used to construct the model
  auto inputs = make_shared<STR_NC_data_model>(m_data_model);
 
  // set prior
  m_set_prior(cosmo_param_prior);
 
  // construct the model
  switch (m_HistogramType) {
  case (glob::HistogramType::_N_V_):
    m_model = make_shared<statistics::Model1D>(statistics::Model1D(&number_counts_redshift, nParams, cosmoPar_type, cosmoPar_string, inputs));
    break;
  case (glob::HistogramType::_n_V_):
    m_model = make_shared<statistics::Model1D>(statistics::Model1D(&number_density_redshift, nParams, cosmoPar_type, cosmoPar_string, inputs));
    break;
  default:
    ErrorCBL("no such a variable in the list!", "set_model_NumberCounts_cosmology", "Modelling_NumberCounts1D_Redshift.cpp");
  }
}