LogNormalFull.h

Go to the documentation of this file.

/********************************************************************
 *  Copyright (C) 2010 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 __LOGNORMALF__
#define __LOGNORMALF__
 
#include "Catalogue.h"
#include <fftw3.h>
 
namespace cbl {
    
  namespace lognormal {
 
    class LogNormalFull {
 
    protected:
      
      std::vector<std::shared_ptr<catalogue::Catalogue>> m_data;
 
      std::vector<std::shared_ptr<catalogue::Catalogue>> m_random;
      
      std::shared_ptr<cosmology::Cosmology> m_cosmology;
 
      std::string m_author;
 
      bool m_NL;
 
      std::default_random_engine m_generator;
 
      std::shared_ptr<data::Field3D> m_clustering_signal;
 
      double m_rmin;
 
      int m_nx;
 
      int m_ny;
 
      int m_nz;
 
      int m_nzF;
 
      double m_xMin;
 
      double m_xMax;
 
      double m_yMin;
 
      double m_yMax;
 
      double m_zMin;
 
      double m_zMax;
 
      int m_nCells_eff;
 
      int m_nLN;
 
      double m_rsmooth;
 
      bool m_use_visibility;
 
      std::shared_ptr<data::Field3D> m_density;  
 
      std::shared_ptr<data::Field3D> m_densityG;  
 
      double m_sigma2G;
 
      std::shared_ptr<data::Field3D> m_rsd_displacement;  
 
      std::shared_ptr<data::Field3D> m_displacement; 
 
      std::shared_ptr<data::Field3D> m_potential;  
 
      std::shared_ptr<data::Field3D> m_velocity; 
 
      std::shared_ptr<data::Field3D> m_los_velocity; 
 
      std::shared_ptr<data::Field3D> m_visibility; 
 
      std::vector<std::shared_ptr<data::Field3D>> m_visibility_random; 
 
      std::vector<std::vector<std::shared_ptr<catalogue::Catalogue>>> m_LNCat;  
 
      std::shared_ptr<glob::FuncGrid> m_func_redshift;
 
      std::shared_ptr<glob::FuncGrid> m_func_DC;
 
      std::shared_ptr<glob::FuncGrid> m_func_HH;
 
      std::vector<std::shared_ptr<glob::FuncGrid> > m_func_bias;
 
      std::shared_ptr<glob::FuncGrid> m_func_growth_factor;
 
      std::shared_ptr<glob::FuncGrid> m_func_growth_rate;
 
      std::shared_ptr<glob::FuncGrid> m_func_pk;
 
      void m_set_grid_parameters ();
 
      void m_set_fields (const bool use_random, const bool doRSD);
 
      void m_reset_fields (const bool use_random, const bool doRSD);
 
      void m_set_clustering_signal ();
 
      void m_set_density_field (const double smoothing_radius);
 
      void m_set_potential ();
 
      void m_set_radial_velocity ();
 
      void m_set_visibility ();
 
      void m_set_visibility_from_random ();
 
      void m_set_visibility_from_random_RSD ();
 
      void m_extract_points_lognormal_field (const double nObjects, const bool doRSD, const std::vector<double> redshift, const std::vector<double> bias, const std::shared_ptr<data::Field3D> visibility, const std::string file_out);
 
    
    public:
 
      LogNormalFull () = default;
 
      LogNormalFull (const cosmology::Cosmology cosmology, const double redshift_min=0., const double redshift_max=10., const int n_redshift_bins=500, const std::string author="CAMB");
 
      LogNormalFull (const double rmin, const double xMin, const double xMax, const double yMin, const double yMax, const double zMin, const double zMax, const cosmology::Cosmology cosmology, const double redshift_min=0., const double redshift_max=10., const int n_redshift_bins=500, const std::string author="CAMB");
 
      LogNormalFull (const double rmin, const std::vector<std::shared_ptr<catalogue::Catalogue>> random, const double pad, const cosmology::Cosmology cosmology, const double redshift_min=0., const double redshift_max=10., const int n_redshift_bins=500, const std::string author="CAMB");
 
      ~LogNormalFull () = default;
 
      void set_cosmo_function (const cosmology::Cosmology cosmology, const double redshift_min=0., const double redshift_max=10., const int n_redshift_bins=500, const std::string author="CAMB");
 
      void set_grid_parameters (const double rmin, const double xMin, const double xMax, const double yMin, const double yMax, const double zMin, const double zMax);
 
      void set_grid_parameters (const double rmin, const std::vector<std::shared_ptr<catalogue::Catalogue>> random, const double pad);
 
      void generate_lognormal (const int start, const int stop, const bool doRSD, const double smoothing_radius, const std::vector<double> nObjects, const std::vector<std::vector<double>> redshift, const std::vector<std::vector<double> > bias, const std::string dir, const std::string filename, const int seed, const bool set_fields=true, const bool use_random=true);
 
    };
  }
}
 
#endif