SuperSampleCovariance.h

Go to the documentation of this file.

/********************************************************************
 *  Copyright (C) 2021 by Federico Marulli and Giorgio Lesci        *
 *  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 __SSC__
#define __SSC__
 
 
#include "Cosmology.h"
 
 
namespace cbl {
 
  namespace cosmology {
 
     
    class SuperSampleCovariance
    {
      
    protected:
    
      std::vector<std::shared_ptr<statistics::Model>> m_response_func;
    
      std::shared_ptr<cosmology::Cosmology> m_cosmo;
      
      std::vector<cbl::cosmology::CosmologicalParameter> m_cosmo_param;
      
      std::string m_method_Pk;
      
      bool m_NL;
      
      bool m_store_output;
      
      int m_nbins;
      
      int m_nsteps;
      
      std::vector<double> m_redshifts;
      
      double m_area;
      
      double m_precision;
      
      std::vector<std::vector<double>> m_windows;
 
      void m_compute_topHat_window (const double delta_z, const std::vector<double> redshift_edges);
      
      void m_compute_gaussian_window (const double delta_z, const std::vector<double> W_mean, const std::vector<double> W_std);
      
      std::vector<std::vector<double>> m_compute_Sij (cbl::cosmology::Cosmology cosmo) const;
 
 
    public:
 
 
      SuperSampleCovariance (cbl::cosmology::Cosmology cosm, const std::vector<cbl::cosmology::CosmologicalParameter> cosmo_param, const std::vector<double> redshift_edges, const double area, const std::string method_Pk="EisensteinHu", const double delta_z=0.001, const double precision=10, const bool NL=false, const bool store_output=false);
      
      SuperSampleCovariance (cbl::cosmology::Cosmology cosm, const std::vector<cbl::cosmology::CosmologicalParameter> cosmo_param, const double area, const std::vector<double> W_mean, const std::vector<double> W_std, const std::string method_Pk="EisensteinHu", const double delta_z=0.001, const double precision=10, const bool NL=false, const bool store_output=false);
 
      virtual ~SuperSampleCovariance () = default;
 
      
      
      std::vector<std::vector<double>> operator () (std::vector<double> &parameter) const;
       
       
       
      
      
       std::vector<std::vector<double>> get_window_function ();
       
       int Sij_dimension ();
       
       
       
      
      
       void write_window_function (const std::string dir, const std::string file);
       
       void write_Sij (const std::string dir, const std::string file);
       
       
      
    };
 
 
  }
}
 
#endif