18 #ifndef SOURCEXTRACTORPLUSPLUS_KAPPASIGMABINNING_H
19 #define SOURCEXTRACTORPLUSPLUS_KAPPASIGMABINNING_H
24 namespace SourceXtractor {
39 template <
typename VarType>
53 KappaSigmaBinning(
float kappa1 = 2.,
float kappa2 = 5.,
size_t min_pixels = 4,
size_t max_size = 4096)
68 template<
typename Iterator>
74 for (
auto i = begin; i != end; ++i) {
80 sigma = std::abs(*begin - mean);
87 auto lcut = mean - sigma;
88 auto hcut = mean + sigma;
92 for (
auto i = begin; i != end; ++i) {
93 if (*i >= lcut && *i <= hcut)
131 return std::min(nbins, static_cast<size_t>(4096));
167 #endif //SOURCEXTRACTORPLUSPLUS_KAPPASIGMABINNING_H