VariableBinwidthHistogram.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
//   VariableBinwidthHistogram.h
//   Header file for class VariableBinwidthHistogram
// (c) ATLAS Detector software
//
 
#ifndef VARIABLEBINWIDTHHISTOGRAM_H
#define VARIABLEBINWIDTHHISTOGRAM_H
 
// c - c++
#include <algorithm>
#include <iostream>
#include <vector>
 
// ROOT
//#include "TH1.h"
//#include "TGraph.h"
 
class TH1F;
class TGraph;
 
// this
#include "VariableBinwidthHistogramBin.h"
 
namespace MuonCalib {

    class VariableBinwidthHistogram {
    public:
        //----------------------------constructors--------------------------------------
        inline VariableBinwidthHistogram() : m_binc(0), m_max_bin_width(0), m_error(false), m_sorted(false) {}

        inline ~VariableBinwidthHistogram() {
            for (unsigned int i = 0; i < m_bins.size(); i++) { delete m_bins[i]; }
        }
 
        //-----------------------------public member functinos--------------------------

        bool Initialize(TH1F *hist, double binc_r, double max_bin_width, double min_x = -9e9, double max_x = 9e9);

        bool Smooth(double width);

        inline const VariableBinwidthHistogramBin &GetBin(unsigned int bin) const { return *(m_bins[bin]); }

        inline const VariableBinwidthHistogramBin &GetSortedBin(unsigned int bin) {
            if (!m_sorted) {
                sort(m_sort_bins.begin(), m_sort_bins.end());
                m_sorted = true;
            }
            return m_sort_bins[bin].Bin();
        }

        inline unsigned int GetNumberOfBins() const { return m_bins.size(); }

        TGraph *DenistyGraph() const;

        TGraph *BinWidthGraph() const;

        TGraph *BinContentGraph() const;

        TGraph *DiffDensityGraph() const;

        TGraph *DiffBinwidthGraph() const;
 
    private:
        //----------------------------private data members------------------------------
        std::vector<VariableBinwidthHistogramBin *> m_bins;

        std::vector<VBHBinPtrSrt> m_sort_bins;

        unsigned int m_binc;

        double m_max_bin_width;

        bool m_error;
 
        bool m_sorted;
 
        //--------------------------privite member functions----------------------------
        inline double sign(double val) const {
            if (val > 0.0) return 1.0;
            if (val < 0.0) return -1.0;
            if (val == 0.0) return 0;
            return 0;
        }
    };
 
}  // namespace MuonCalib
 
#endif