ECLLocalRunCalibAcc.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2019 - Belle II Collaboration                             *
 *                                                                        *
 * ECLLocalRunCalibAcc                                                    *
 *                                                                        *
 * Feature (mean value, standard deviation and number of accepted         *
 * events) accumulator.                                                   *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Sergei Gribanov (S.S.Gribanov@inp.nsk.su),               *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
// STL
#include <vector>
// BOOST
#include <boost/accumulators/accumulators.hpp>
#include <boost/accumulators/statistics.hpp>
namespace bacc = boost::accumulators;
namespace Belle2 {
  class ECLLocalRunCalibAcc {
  public:
    ECLLocalRunCalibAcc(
      const float& min_value,
      const float& max_value,
      const int* const ndevs);
    ~ECLLocalRunCalibAcc();
    int getNOfEvents() const;
    int getCount() const;
    float getMean() const;
    float getStdDev() const;
    void add(const float& value);
    void calc();
  private:
    bool isValueInRange(const float& value) const;
    float calcStdDev(const float& variance,
                     const int& count) const;
    void updateLimits(const float& mean,
                      const float& stddev);
    float m_leftLimit;
    float m_rightLimit;
    const int* const c_ndevs;
    int m_nevents;
    int m_count;
    float m_mean;
    float m_stddev;
    std::vector<int> m_data;
  };
}