/cvmfs/atlas.cern.ch/repo/sw/ASG/AnalysisBase/2.4.31/xAODTrigMinBias/xAODTrigMinBias/versions/TrigHisto2D_v1.h

#ifndef xAODTrigMinBias_TrigHisto2D_V1_H
#define xAODTrigMinBias_TrigHisto2D_V1_H

#include "AthContainers/AuxElement.h"

namespace xAOD {
  
  namespace TrigHistoCutType {
    enum CutTypeEnum {
      BELOW_X,
      ABOVE_X,
      BELOW_X_BELOW_Y,
      ABOVE_X_BELOW_Y,
      BELOW_X_ABOVE_Y,
      ABOVE_X_ABOVE_Y};
  }
//   class MsgStream;
  
  class TrigHisto2D_v1 : public SG::AuxElement {
  public:
    TrigHisto2D_v1();
    
    ~TrigHisto2D_v1();
    
    TrigHisto2D_v1(unsigned int nbins_x, float min_x, float max_x, unsigned int nbins_y, float min_y, float max_y);
    const std::vector<float>& contents() const;
    
    unsigned int nbinsX() const;
    
    float minX()const ; 
    
    float maxX() const; 
    
    unsigned int nbinsY() const;
    
    float minY() const;
    
    float maxY() const;
    
    void initialize(unsigned int nbins_x, float min_x, float max_x, unsigned int nbins_y, float min_y, float max_y);
    
    void fill(float value_x, float value_y, float weight);
    
    //double sumEntries(float value_x, float value_y, int cutType);
    
    // it's not very efficient ;/ 
    std::vector<float> profileX() const;
    
    std::vector<float> profileY() const;
    
    double sumEntries(float value_x, float value_y, int cutType) const;
    
    void clear();
    
    void dump();
   
//     void setMsgStream(MsgStream&);
    
    void setContents(const std::vector<float> &cont);
    
    void setMinX(float val);
    
    void setMaxX(float val);
    
    void setNbinsX(unsigned int nx);
    
    void setMinY(float val);
    
    void setMaxY(float val);
    
    void setNbinsY(unsigned int ny);
    
      private:
    unsigned int findBinX(float val) const;
    
    unsigned int findBinY(float val) const;
    
    static const int EXTRA_BINS=2;
    
    /*width of x bins*/
    float m_binWidthX;
   
    /*width of y bins*/
    float m_binWidthY;
  };
  
}

#endif

---