ModelFunction_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 __MODFUNCNC__
#define __MODFUNCNC__
 
#include "Cosmology.h"
#include "NumberCounts.h"
 
 
// ============================================================================
 
 
namespace cbl {
 
  namespace modelling {
 
    namespace numbercounts {
      
      struct STR_NC_data_model {
 
    bool isSnapshot;
    
    std::vector<double> edges_x;
    
    std::vector<double> edges_y;
 
    std::shared_ptr<cosmology::Cosmology> cosmology;
 
    std::vector<cosmology::CosmologicalParameter> Cpar;
    
    std::function<double(const double, const double, const std::shared_ptr<void>)> fz;
    
    double (*response_fact)(const double, const double, const double, const double, const std::string, const double, const std::string, std::shared_ptr<void>);
    
    std::function<double(const double, const double)> z_error;
    
    std::function<double(const double, const double)> proxy_error;
 
    double redshift;
 
    std::string method_Pk;
      
    bool NL;
 
    double k_min;
 
    double k_max;
 
    int step;
 
    std::vector<double> kk;
 
    bool store_output;
      
    std::string output_root;
 
    int norm;
 
    std::string file_par;
 
    double prec;
 
    double Delta;
 
    bool isDelta_critical;
 
    std::string model_MF;
    
    std::string model_bias;
 
    double z_min; 
 
    double z_max;
 
    int z_step;
 
    std::vector<double> z_vector;
 
    double Mass_min; 
 
    double Mass_max; 
 
    int Mass_step;
 
    std::vector<double> Mass_vector;
 
    double area_rad;
 
    double Volume;
 
        bool use_SF;
 
    bool is_sigma8_free;
    
    std::shared_ptr<glob::FuncGrid2D> interp_SelectionFunction;
    
    std::vector<double> SF_weights;
    
    double z_pivot;
    
    double proxy_pivot;
    
    double mass_pivot;
    
    double log_base;
 
    STR_NC_data_model () = default;
      };
 
 
   
      struct STR_NCSF_data_model {
 
    std::vector<double> radii;
 
    std::shared_ptr<cosmology::Cosmology> cosmology;
 
    std::vector<cosmology::CosmologicalParameter> Cpar;
 
    double redshift;
 
    std::string model_SF;
 
    double b_eff;
 
    double b_slope;
 
    double b_offset;
 
    double deltav_NL;
 
    double delta_c;
    
    std::string method_Pk;
 
    double k_Pk_ratio;
 
    bool store_output;
 
    std::string output_root;
 
    std::string interpType;
 
    double k_max;
 
    std::string input_file;
 
    bool is_parameter_file;
      
    STR_NCSF_data_model () = default;
      };
      
 
      double Filter_sigmaR (const double kk, const double radius); 
 
      double Filter_dsigmaR (const double kk, const double radius); 
 
      void sigmaM_dlnsigmaM (double &sigmaM, double &dlnsigmaM, const double mass, const cbl::glob::FuncGrid interp_Pk, const double kmax, const double rho);
 
      void sigmaM_dlnsigmaM (std::vector<double> &sigmaM, std::vector<double> &dlnsigmaM, const std::vector<double> mass, const std::vector<double> kk, const std::vector<double> Pk, const std::string interpType, const double kmax, const double rho);
 
      std::vector<cbl::glob::FuncGrid> sigmaM_dlnsigmaM (const std::vector<double> mass, cosmology::Cosmology cosmology, const std::vector<double> kk, const std::vector<double> Pk, const std::string interpType, const double kmax);
 
      double mass_function (const double mass, cosmology::Cosmology cosmology, const double redshift, const std::string model_MF, const bool store_output, const double Delta, const bool isDelta_critical, const cbl::glob::FuncGrid interp_Pk, const double kmax);
 
      std::vector<double> mass_function (const std::vector<double> mass, cosmology::Cosmology cosmology, const double redshift, const std::string model_MF, const bool store_output, const double Delta, const bool isDelta_critical, const std::vector<double> kk, const std::vector<double> Pk, const std::string interpType, const double kmax);
 
      std::vector<std::vector<double>> mass_function (const std::vector<double> redshift, const std::vector<double> mass, cosmology::Cosmology cosmology, const std::string model_MF, const bool store_output, const double Delta, const bool isDelta_critical, const std::vector<double> kk, const std::vector<double> Pk, const std::string interpType, const double kmax);
 
      std::vector<double> size_function (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 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);
 
      double number_counts (const double redshift_min, const double redshift_max, const double Mass_min, const double Mass_max, cosmology::Cosmology cosmology, const double Area, const std::string model_MF, const bool store_output, const double Delta, const bool isDelta_critical, const glob::FuncGrid interp_sigmaM, const  glob::FuncGrid interp_DlnsigmaM);
      
      double counts_proxy (const double alpha, const double beta, const double gamma, const double scatter0, const double scatterM, const double scatterM_exp, const double scatterz, const double scatterz_exp, const double z_bias, const double proxy_bias, const double z_err, const double proxy_err, const double Plambda_a, const double Plambda_b, const double Plambda_c, std::function<double(const double, const double, const std::shared_ptr<void>)> fz, std::function<double(const double, const double)> z_error, std::function<double(const double, const double)> proxy_error, double (*response_fact)(const double, const double, const double, const double, const std::string, const double, const std::string, std::shared_ptr<void>), const double redshift_min, const double redshift_max, const double proxy_min, const double proxy_max, cbl::cosmology::Cosmology cosmology, const double Area, const std::string model_MF, const std::string model_bias, const bool store_output, const double Delta, const bool isDelta_critical, const cbl::glob::FuncGrid interp_sigmaM, const cbl::glob::FuncGrid interp_DlnsigmaM, const cbl::glob::FuncGrid interp_DN, const double proxy_pivot, const double z_pivot, const double mass_pivot, const double log_base, const double weight);
      
    }
  }
}
 
#endif