VariableBinwidthHistogramBin.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef VARIABLEBINWIDTHHISTOGRAMBIN_H
#define VARIABLEBINWIDTHHISTOGRAMBIN_H
 
#include <TString.h>  // for Form
 
#include <cstdlib>
#include <iostream>
#include <vector>
 
namespace MuonCalib {

    class VariableBinwidthHistogramBin {
    public:
        //----------------------------------constructor---------------------------------

        inline VariableBinwidthHistogramBin() : m_center(0.0), m_width(0.0), m_content(0) {}

        inline VariableBinwidthHistogramBin(double center, double width, unsigned int content) :
            m_center(center), m_width(width), m_content(content) {}
 
        //-----------------------------------public member functions--------------------

        inline void Set(double center, double width, unsigned int content) {
            m_center = center;
            m_width = width;
            m_content = content;
        }

        inline void SetContent(unsigned int n) { m_content = n; }

        inline VariableBinwidthHistogramBin operator+=(unsigned int n) {
            m_content += n;
            return *this;
        }

        inline VariableBinwidthHistogramBin operator+(unsigned int n) {
            VariableBinwidthHistogramBin ret;
            ret += n;
            return ret;
        }

        void MoveRight(double new_right) {
            double left = Left();
            if (new_right < left) {
                throw std::runtime_error(
                    Form("File: %s, Line: %d\nVariableBinwidthHistogramBin::MoveRight() - new right is too small!", __FILE__, __LINE__));
            }
            m_width = new_right - left;
            m_center = 0.5 * (left + new_right);
        }

        void MoveLeft(double new_left) {
            double right = Right();
            if (new_left > right) {
                throw std::runtime_error(
                    Form("File: %s, Line: %d\nVariableBinwidthHistogramBin::MoveRight() - new left is too small!", __FILE__, __LINE__));
            }
            m_width = right - new_left;
            m_center = 0.5 * (new_left + right);
        }

        inline double Center() const { return m_center; }

        inline double Width() const { return m_width; }

        inline double Left() const { return m_center - 0.5 * m_width; }

        inline double Right() const { return m_center + 0.5 * m_width; }

        inline unsigned int Entries() const { return m_content; }

        inline double Density() const { return static_cast<double>(m_content) / m_width; }
 
    private:
        //-----------------------------------------private data members-----------------
 
        double m_center, m_width;
        unsigned int m_content;
    };

    class VBHBinPtrSrt {
    public:
        //-----------------------------------constructor--------------------------------
        inline VBHBinPtrSrt() {}

        inline VBHBinPtrSrt(VariableBinwidthHistogramBin *bin) : m_bin(bin) {}
 
        //--------------------------------------public member functions-----------------

        inline void Initialize(VariableBinwidthHistogramBin *bin) { m_bin = bin; }

        inline VariableBinwidthHistogramBin &Bin() { return *m_bin; }

        inline bool operator>(const VBHBinPtrSrt &other) const { return m_bin->Width() > other.m_bin->Width(); }

        inline bool operator<(const VBHBinPtrSrt &other) const { return m_bin->Width() < other.m_bin->Width(); }
 
    private:
        //------------------------------private data members----------------------------

        VariableBinwidthHistogramBin *m_bin;
    };
 
}  // namespace MuonCalib
 
#endif