StackedDensityProfile.h

Go to the documentation of this file.

/********************************************************************
 *  Copyright (C) 2020 by Giorgio Lesci and Federico Marulli        *
 *  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 __STACKPROFILE__
#define __STACKPROFILE__
 
#include "Catalogue.h"
#include "Cosmology.h"
#include "Measure.h"
 
 
// ===================================================================================================
 
 
namespace cbl {
 
  namespace measure {
 
    namespace stackprofile {
 
      class StackedDensityProfile : public Measure {
 
      protected :     
 
 
    
    void m_set_multiplicative_calib();
    
    bool m_colourSelection_gri (const int i_gal, const int clu_zbin);
    
    bool m_photzSelection (const int i_clu, const int i_gal, const int clu_zbin);
    
    void m_check_catalogue_variables();
    
    void m_resize(const double rad_min, const double rad_max, const int nRad, const bool log_rad);
    
    void m_linked_list();
    
    void m_add_galaxy(const int i_gal, const int clu_index, const double coscludec, const double clu_dist);
    
    void m_profile(const int clu_index);
    
    void m_stacker();
 
    std::shared_ptr<data::Data> m_make_bootstrap(const std::vector<int> z_proxy_bin);
    
    bool m_check_file(const std::string checked_file, const std::vector<int> z_proxy_bin);
    
    void m_write(const std::string output_dir, const std::string output_file);
    
    bool m_write_background;
    
    std::shared_ptr<cosmology::Cosmology> m_cosm;
    
    std::shared_ptr<catalogue::Catalogue> m_galData;
    
    std::shared_ptr<catalogue::Catalogue> m_cluData;
    
    std::vector<bool (StackedDensityProfile::*)(const int, const int)> m_colourSel; // std::vector<std::function<bool(const int)>> m_colourSel;
    
    std::vector<bool (StackedDensityProfile::*)(const int, const int, const int)> m_photzSel;
    
    std::vector<std::function<bool(const std::vector<bool>)>> m_logicSel;
    
    std::vector<bool> m_isSet_colourSel;
    
    std::vector<bool> m_isSet_photzSel;
    
    std::vector<bool> m_isSet_logicSel;
    
    std::vector<std::vector<int>> m_background_idx_z;
    
    std::vector<std::vector<int>> m_background_idx_colour;
    
    bool m_measure_is_read;
    
    double m_delta_redshift;
    
    double m_SN_min;
    
    double m_pix_size;
    
    std::vector<std::vector<double>> m_m_calib;
    
    std::vector<std::vector<double>> m_m_calib_zEdges;
    
    double m_n_resampling;
    
    std::vector<std::string> m_inputs_to_str;
    
    double m_rad_alpha;
    
    double m_obs_gamma;
    
    std::vector<std::vector<double>> m_colour_sel_pars;
    
    std::vector<std::vector<double>> m_zphot_sel_pars;
    
    std::vector<std::vector<std::string>> m_logic_sel_par;
    
    std::vector<double> m_rad_arr;
    
    std::vector<double> m_z_binEdges;
    
    std::vector<std::vector<double>> m_proxy_binEdges;
    
    std::vector<std::vector<int>> m_first;
    
    std::vector<std::vector<int>> m_last;
    
    std::vector<int> m_next;
    
    std::vector<std::vector<double>> m_nClu_inBin;
    
    std::vector<int> m_ngal_arr;
    
    std::vector<std::vector<std::vector<std::vector<int>>>> m_single_ngal_arr;
    
    std::vector<std::vector<std::vector<int>>> m_stacked_ngal_arr;
    
    std::vector<double> m_deltasigma_t;
    
    std::vector<std::vector<std::vector<std::vector<double>>>> m_single_deltasigma_t;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_deltasigma_t;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_deltasigma_t_err;
    
    std::vector<double> m_deltasigma_x;
    
    std::vector<std::vector<std::vector<std::vector<double>>>> m_single_deltasigma_x;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_deltasigma_x;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_deltasigma_x_err;
    
    std::vector<double> m_deltasigma_err;
    
    std::vector<std::vector<std::vector<std::vector<double>>>> m_single_deltasigma_err;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_deltasigma_err;
    
    std::vector<double> m_wetasquareSum;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_wetasquareSum;
    
    std::vector<double> m_wSum;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_wSum;
    
    std::vector<double> m_wetaSum;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_wetaSum;
    
    std::vector<double> m_deltasigmaSum;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_deltasigmaSum;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_deltasigma_wei;
    
    std::vector<double> m_rad_eff_arr; 
    
    std::vector<std::vector<std::vector<double>>> m_stacked_rad_eff_arr;
    
    std::vector<double> m_rad_sigma_arr;
    
    std::vector<std::vector<std::vector<double>>> m_stacked_rad_sigma_arr;
    
    double m_sigma_fac;
    
    double m_ra_start;
    
    double m_dec_start;
    
    std::vector<int> m_nPix;
    
    std::vector<std::vector<double>> m_proxy_eff;
    
    std::vector<std::vector<double>> m_proxy_sigma;
 
    std::vector<std::vector<double>> m_z_eff;
    
    std::vector<std::vector<double>> m_z_sigma;
    
    std::vector<std::vector<std::vector<std::vector<double>>>> m_deltasigma_cov_matr;
    
    
 
 
      public:
 
 
    StackedDensityProfile () = default;
 
    virtual ~StackedDensityProfile () = default;
 
 
    StackedDensityProfile (cosmology::Cosmology cosm, const catalogue::Catalogue gal_cat, const catalogue::Catalogue clu_cat, const double delta_redshift, std::vector<double> z_binEdges, std::vector<std::vector<double>> proxy_binEdges, const double rad_min, const double rad_max, const int nRad, const bool log_rad, const double SN_min, const double pix_size, const std::vector<std::vector<double>> multiplicative_calibration_stats={}, const std::vector<std::vector<double>> multiplicative_calibration_zEdges={}, const double rad_alpha=1., const double obs_gamma=1.);
    
    
    void set_colour_selection(const std::string colour_space, const int z_bin, const double C1, const double C2, const double C3, const double C4, const double z_min);
    
    void set_colour_selection(const std::string colour_space, const double C1, const double C2, const double C3, const double C4, const double z_min);
    
    void set_zphot_selection(const int z_bin, const double deltaz, const double zgal_min, const double zgal_max, const double ODDS_min);
    
    void set_zphot_selection(const double deltaz, const double zgal_min, const double zgal_max, const double ODDS_min);
    
    void set_logic_selection(const int z_bin, const std::string sel);
    
    void set_logic_selection(const std::string sel);
    
    void measure(const std::vector<int> z_proxy_bin, const std::string output_dir, const std::string output_file_root, const ErrorType errorType=ErrorType::_Bootstrap_, const int n_resampling=10000, const bool write_background=false);
    
    
    
    void write(const std::string dir, const std::string file);
    
    
      };
    }
  }
}
 
#endif