Modelling_NumberCounts.h

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.       *
 ********************************************************************/
 
#ifndef __MODELLINGNC__
#define __MODELLINGNC__
 
 
#include "NumberCounts.h"
#include "Modelling.h"
#include "ModelFunction_NumberCounts.h"
 
 
// ===================================================================================================
 
 
namespace cbl {
 
  namespace modelling {
 
    namespace numbercounts {
    
      class Modelling_NumberCounts : public Modelling {
      
      protected:
 
    glob::HistogramType m_HistogramType;
 
    double m_fact;
 
    modelling::numbercounts::STR_NC_data_model m_data_model;
 
    modelling::numbercounts::STR_NCSF_data_model m_data_model_SF;
 
      public:
 
 
    Modelling_NumberCounts () = default;
    
    Modelling_NumberCounts (const std::shared_ptr<cbl::measure::numbercounts::NumberCounts> nc)
    { m_HistogramType = nc->HistogramType(); m_fact = nc->fact(); }
    
    Modelling_NumberCounts (glob::HistogramType hist_type, double fact)
    { m_HistogramType = hist_type; m_fact = fact;}
    
    virtual ~Modelling_NumberCounts () = default;
 
 
 
    modelling::numbercounts::STR_NC_data_model data_model () { return m_data_model; }
 
 
 
    modelling::numbercounts::STR_NCSF_data_model data_model_SF () { return m_data_model_SF; }
    
 
    
 
    void set_data_model (const cbl::cosmology::Cosmology cosmology={}, const double redshift=0., const std::string method_Pk="CAMB", const double k_min=1.e-4, const double k_max=100., const int step=500, const bool store_output=true, const int norm=-1, const double Delta=200., const bool isDelta_critical=true, const std::string model_MF="Tinker", const std::string selection_function_file=par::defaultString, const std::vector<int> selection_function_column={}, const double z_min=par::defaultDouble, const double z_max=par::defaultDouble, const int z_step=50, const double Mass_min=par::defaultDouble, const double Mass_max=par::defaultDouble, const int Mass_step=100, const double area_degrees=par::defaultDouble, const double prec=1.e-4);
    
    void set_data_model (const cbl::cosmology::Cosmology cosmology, const std::vector<double> SF_weights, const double z_pivot, const double proxy_pivot, const double mass_pivot, const double log_base, const std::string method_Pk, const bool store_output=true, const int norm=-1, const double Delta=200., const bool isDelta_critical=true, const std::string model_MF="Tinker", const std::string model_bias="Tinker", const double z_min=par::defaultDouble, const double z_max=par::defaultDouble, const double area_degrees=par::defaultDouble, const double prec=1.e-4);
 
        
    void set_data_model_SF (const cosmology::Cosmology cosmology, const std::vector<double> radii, const double redshift, const std::string model, const double b_eff, double slope=0.854, double offset=0.420, const double deltav_NL=-0.795, const double del_c=1.69, const std::string method_Pk="EisensteinHu", const double k_Pk_ratio=-1., const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true); 
 
 
    virtual void set_data_model_snapshot (const cbl::cosmology::Cosmology cosmology={}, const double redshift=0., const std::string method_Pk="CAMB", const double k_min=1.e-4, const double k_max=100., const int step=500, const bool store_output=true, const int norm=-1, const double Delta=200., const bool isDelta_critical=true, const std::string model_MF="Tinker", const double Volume=par::defaultDouble, const double Mass_min=par::defaultDouble, const double Mass_max=par::defaultDouble, const int Mass_step=100, const double prec=1.e-4) 
    { 
      (void)cosmology; (void)redshift; (void)method_Pk; (void)k_min; (void)k_max; (void)step;
      (void)store_output; (void)norm; (void)Delta; (void)isDelta_critical; (void)model_MF;
      (void)Volume; (void)prec; (void)Mass_min; (void)Mass_max; (void)Mass_step;
      cbl::ErrorCBL("", "set_data_model_snapshot", "Modelling_NumberCounts.h");
    }
 
 
    
      };
    }
  }
}
 
#endif