35 #ifndef __TWOPOINT2DPOL__
36 #define __TWOPOINT2DPOL__
145 void measureJackknife (
const std::string dir_output_pairs,
const std::vector<std::string> dir_input_pairs={},
const std::string dir_output_resample =
"NULL",
const bool count_dd=
true,
const bool count_rr=
true,
const bool count_dr=
true,
const bool tcount=
true,
const Estimator estimator=
Estimator::_LandySzalay_,
const double fact=0.1)
override;
192 void measureBootstrap (
const int nMocks,
const std::string dir_output_pairs,
const std::vector<std::string> dir_input_pairs={},
const std::string dir_output_resample=
"NULL",
const bool count_dd=
true,
const bool count_rr=
true,
const bool count_dr=
true,
const bool tcount=
true,
const Estimator estimator=
Estimator::_LandySzalay_,
const double fact=0.1,
const int seed=3213)
override;
240 TwoPointCorrelation2D_polar (
catalogue::Catalogue data,
catalogue::Catalogue random,
const BinType binType_rad,
const double rMin,
const double rMax,
const int nbins_rad,
const double shift_rad,
const BinType binType_mu,
const double muMin,
const double muMax,
const int nbins_mu,
const double shift_mu,
const CoordinateUnits angularUnits=
CoordinateUnits::_radians_, std::function<
double(
double)> angularWeight=
nullptr,
const bool compute_extra_info=
false,
const double random_dilution_fraction=1.)
241 :
TwoPointCorrelation2D(
data,
random,
compute_extra_info,
random_dilution_fraction) {
m_twoPType =
TwoPType::_2D_polar_;
set_parameters(binType_rad, rMin, rMax, nbins_rad, shift_rad, binType_mu, muMin, muMax, nbins_mu, shift_mu, angularUnits, angularWeight,
compute_extra_info); }
275 TwoPointCorrelation2D_polar (
catalogue::Catalogue data,
catalogue::Catalogue random,
const BinType binType_rad,
const double rMin,
const double rMax,
const double binSize_rad,
const double shift_rad,
const BinType binType_mu,
const double muMin,
const double muMax,
const double binSize_mu,
const double shift_mu,
const CoordinateUnits angularUnits=
CoordinateUnits::_radians_, std::function<
double(
double)> angularWeight=
nullptr,
const bool compute_extra_info=
false,
const double random_dilution_fraction=1.)
276 :
TwoPointCorrelation2D(
data,
random,
compute_extra_info,
random_dilution_fraction) {
m_twoPType =
TwoPType::_2D_polar_;
set_parameters(binType_rad, rMin, rMax, binSize_rad, shift_rad, binType_mu, muMin, muMax, binSize_mu, shift_mu, angularUnits, angularWeight,
compute_extra_info); }
317 void set_parameters (
const BinType binType_rad,
const double rMin,
const double rMax,
const int nbins_rad,
const double shift_rad,
const BinType binType_mu,
const double muMin,
const double muMax,
const int nbins_mu,
const double shift_mu,
const CoordinateUnits angularUnits=
CoordinateUnits::_radians_, std::function<
double(
double)> angularWeight=
nullptr,
const bool compute_extra_info=
false);
344 void set_parameters (
const BinType binType_rad,
const double rMin,
const double rMax,
const double binSize_rad,
const double shift_rad,
const BinType binType_mu,
const double muMin,
const double muMax,
const double binSize_mu,
const double shift_mu,
const CoordinateUnits angularUnits=
CoordinateUnits::_radians_, std::function<
double(
double)> angularWeight=
nullptr,
const bool compute_extra_info=
false);
397 void measure (
const ErrorType errorType=
ErrorType::_Poisson_,
const std::string dir_output_pairs=
par::defaultString,
const std::vector<std::string> dir_input_pairs={},
const std::string dir_output_resample=
par::defaultString,
const int nMocks=0,
const bool count_dd=
true,
const bool count_rr=
true,
const bool count_dr=
true,
const bool tcount=
true,
const Estimator estimator=
Estimator::_LandySzalay_,
const double fact=0.1,
const int seed=3213)
override;
413 void read (
const std::string dir,
const std::string file)
override;
425 void write (
const std::string dir,
const std::string file,
const bool full,
const int rank=0)
const override;
435 {
write(dir, file,
true, rank); }
The class TwoPointCorrelation2D.
The class TwoPointCorrelation2D_polar.
void measure(const ErrorType errorType=ErrorType::_Poisson_, const std::string dir_output_pairs=par::defaultString, const std::vector< std::string > dir_input_pairs={}, const std::string dir_output_resample=par::defaultString, const int nMocks=0, const bool count_dd=true, const bool count_rr=true, const bool count_dr=true, const bool tcount=true, const Estimator estimator=Estimator::_LandySzalay_, const double fact=0.1, const int seed=3213) override
measure the 2D two-point correlation function in polar coordinates
void measurePoisson(const std::string dir_output_pairs=par::defaultString, const std::vector< std::string > dir_input_pairs={}, const bool count_dd=true, const bool count_rr=true, const bool count_dr=true, const bool tcount=true, const Estimator estimator=Estimator::_LandySzalay_, const double fact=0.1) override
measure the 2D two-point correlation function in polar coordinates, with Poisson errors
void write(const std::string dir, const std::string file, const bool full, const int rank=0) const override
write the 2D two-point correlation function
void compute_covariance(const std::vector< std::string > file, const bool JK) override
compute the covariance matrix
void write(const std::string dir=par::defaultString, const std::string file=par::defaultString, const int rank=0) const override
write the 2D two-point correlation function
void measureJackknife(const std::string dir_output_pairs, const std::vector< std::string > dir_input_pairs={}, const std::string dir_output_resample="NULL", const bool count_dd=true, const bool count_rr=true, const bool count_dr=true, const bool tcount=true, const Estimator estimator=Estimator::_LandySzalay_, const double fact=0.1) override
measure the 2D two-point correlation function in polar coordinates, estimating the covariance with Ja...
~TwoPointCorrelation2D_polar()=default
default destructor
void compute_covariance(const std::vector< std::shared_ptr< data::Data >> xi, const bool JK) override
compute the covariance matrix
TwoPointCorrelation2D_polar(catalogue::Catalogue data, catalogue::Catalogue random, const BinType binType_rad, const double rMin, const double rMax, const double binSize_rad, const double shift_rad, const BinType binType_mu, const double muMin, const double muMax, const double binSize_mu, const double shift_mu, const CoordinateUnits angularUnits=CoordinateUnits::_radians_, std::function< double(double)> angularWeight=nullptr, const bool compute_extra_info=false, const double random_dilution_fraction=1.)
constructor
void measureBootstrap(const int nMocks, const std::string dir_output_pairs, const std::vector< std::string > dir_input_pairs={}, const std::string dir_output_resample="NULL", const bool count_dd=true, const bool count_rr=true, const bool count_dr=true, const bool tcount=true, const Estimator estimator=Estimator::_LandySzalay_, const double fact=0.1, const int seed=3213) override
measure the 2D two-point correlation function in polar coordinates, estimating the covariance with Bo...
void write_covariance(const std::string dir, const std::string file) const override
write the measured two-point correlation
TwoPointCorrelation2D_polar()
default constructor 2D_polar
TwoPointCorrelation2D_polar(catalogue::Catalogue data, catalogue::Catalogue random, const BinType binType_rad, const double rMin, const double rMax, const int nbins_rad, const double shift_rad, const BinType binType_mu, const double muMin, const double muMax, const int nbins_mu, const double shift_mu, const CoordinateUnits angularUnits=CoordinateUnits::_radians_, std::function< double(double)> angularWeight=nullptr, const bool compute_extra_info=false, const double random_dilution_fraction=1.)
constructor
void set_parameters(const BinType binType_rad, const double rMin, const double rMax, const int nbins_rad, const double shift_rad, const BinType binType_mu, const double muMin, const double muMax, const int nbins_mu, const double shift_mu, const CoordinateUnits angularUnits=CoordinateUnits::_radians_, std::function< double(double)> angularWeight=nullptr, const bool compute_extra_info=false)
set the binning parameters
void read_covariance(const std::string dir, const std::string file) override
read the measured covariance matrix
void read(const std::string dir, const std::string file) override
read the 2D two-point correlation function
The class TwoPointCorrelation2D.
std::shared_ptr< catalogue::Catalogue > random() const
get the protected member m_random
bool compute_extra_info() const
get the protected member m_compute_extra_info
bool random_dilution_fraction() const
get the protected member m_random_dilution_fraction
std::shared_ptr< catalogue::Catalogue > data() const
get the protected member m_data
TwoPType m_twoPType
two-point correlation function type
static const std::string defaultString
default std::string value
@ _workInProgress_
error due to work in progress
Estimator
the two-point correlation estimator
@ _LandySzalay_
Landy&Szalay estimator.
@ _2D_polar_
2D two-point correlation function in polar coordinates, ξ(r,μ)
ErrorType
the two-point correlation function error type
@ _Poisson_
Poissonian error.
The global namespace of the CosmoBolognaLib
int ErrorCBL(const std::string msg, const std::string functionCBL, const std::string fileCBL, const cbl::glob::ExitCode exitCode=cbl::glob::ExitCode::_error_)
throw an exception: it is used for handling exceptions inside the CosmoBolognaLib
CoordinateUnits
the coordinate units
@ _radians_
angle in radians