BinnedHistogram.hh

// -*- C++ -*-
#ifndef RIVET_BINNEDHISTOGRAM_HH
#define RIVET_BINNEDHISTOGRAM_HH

#include "Rivet/Config/RivetCommon.hh"
#include "Rivet/Tools/RivetYODA.hh"

namespace Rivet {

  class Analysis;


  template<typename T>
  class BinnedHistogram {
  public:

    BinnedHistogram() {    }

    BinnedHistogram(const vector<T>& edges, const vector<Histo1DPtr>& histos) {
      assert(edges.size() == histos.size()+1);
      for (size_t i = 0; i < histos.size(); ++i)
        addHistogram(edges[i], edges[i+1], histos[i]);
    }



    const BinnedHistogram<T>& add(const T& binMin, const T& binMax, Histo1DPtr histo);
    const BinnedHistogram<T>& addHistogram(const T& binMin, const T& binMax, Histo1DPtr histo) {
      return add(binMin, binMax, histo);
    }


    void fill(const T& binval, double val, double weight);


    const Histo1DPtr histo(const T& binval) const;
    Histo1DPtr histo(const T& binval);

    const vector<Histo1DPtr>& histos() const { return _histos; }
    vector<Histo1DPtr>& histos() { return _histos; }
    const vector<Histo1DPtr>& getHistograms() const { return _histos; }
    vector<Histo1DPtr>& getHistograms() { return _histos; }


    void scale(const T& scale, Analysis* ana);


  private:

    map<T, Histo1DPtr> _histosByUpperBound, _histosByLowerBound;
    vector<Histo1DPtr> _histos;
    map<Histo1DPtr, T> _binWidths;

  };


}

#endif