NumberCounts.h

Go to the documentation of this file.

/********************************************************************
 *  Copyright (C) 2010 by 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 __NCOUNTS__
#define __NCOUNTS__
 
 
#include "Catalogue.h"
#include "Measure.h"
#include "Histogram.h"
 
 
// ===================================================================================================
 
 
namespace cbl {
 
  namespace measure {
 
    namespace numbercounts {
 
      class NumberCounts : public Measure {
 
      protected :
 
 
    std::shared_ptr<catalogue::Catalogue> m_data;
 
 
      
    glob::HistogramType m_HistogramType;
 
    double m_fact;
 
    std::shared_ptr<glob::Histogram> m_histogram;
 
      
 
 
    virtual std::shared_ptr<data::Data> m_measurePoisson () 
    { ErrorCBL("", "m_measurePoisson", "NumberCounts.h"); return NULL; }
 
    virtual std::shared_ptr<data::Data> m_measureJackknife (const std::string dir_output_resample=par::defaultString)
    { (void)dir_output_resample; ErrorCBL("", "m_measureJackknife", "NumberCounts.h"); return NULL; }
 
    virtual std::shared_ptr<data::Data> m_measureBootstrap (const std::string dir_output_resample=par::defaultString, const int nResamplings=0, const int seed=3213)
    { (void)dir_output_resample; (void)nResamplings; (void)seed; ErrorCBL("", "m_measureBootstrap", "NumberCounts.h");  return NULL; }
 
 
      public:
 
 
    NumberCounts () = default;
 
    virtual ~NumberCounts () = default;
 
 
 
    void set_data (const catalogue::Catalogue data)
    { m_data = std::make_shared<catalogue::Catalogue>(catalogue::Catalogue(std::move(data))); }
 
 
 
    std::shared_ptr<catalogue::Catalogue> catalogue () { return m_data; }
 
 
 
    virtual void measure (const ErrorType errorType=ErrorType::_Poisson_, const std::string dir_output_resample=par::defaultString, const int nResamplings=0, const int seed=3213, const bool conv=false, const double sigma=0.)
    { (void)errorType; (void)dir_output_resample; (void)nResamplings; (void)seed; (void)conv; (void)sigma; ErrorCBL("", "measure", "NumberCounts.h"); }
 
 
 
    virtual void write (const std::string dir=par::defaultString, const std::string file=par::defaultString, const int rank=0) const
    { (void)dir; (void)file; (void)rank; ErrorCBL("", "write", "NumberCounts.h"); }
 
 
 
 
    virtual void write_covariance (const std::string dir, const std::string file) const
    { (void)dir; (void)file; ErrorCBL("", "write_covariance", "NumberCounts.h"); }
 
    virtual void compute_covariance (const std::vector<std::shared_ptr<glob::Histogram>> histo, const bool JK)
    { (void)histo; (void)JK; ErrorCBL("", "compute_covariance", "NumberCounts.h"); }
          
    virtual std::shared_ptr<data::Data> Gaussian_smoothing (const double sigma)
    { (void)sigma; ErrorCBL("", "Gaussian_smoothing", "NumberCounts.h"); return NULL; }
 
 
          
    std::shared_ptr<glob::Histogram> histogram () { return m_histogram; }
 
    glob::HistogramType HistogramType () { return m_HistogramType; }
 
    double fact () { return m_fact; }
 
 
 
      };
    }
  }
}
 
#endif