Axis.h

// -*- C++ -*-
//
// Axis.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef LWH_Axis_H
#define LWH_Axis_H
//
// This is the declaration of the Axis class representing
//


#include <limits>
#include <cmath>
#include <algorithm>
#include "AIAxis.h"

namespace LWH {

using namespace AIDA;

class Axis: public IAxis {

public:

  Axis(int n, double lo, double up)
    : lower(lo), upper(up), nbins(n) {}

  Axis(const Axis & a)
    : IAxis(a), lower(a.lower), upper(a.upper), nbins(a.nbins) {}

  virtual ~Axis() { }

  bool isFixedBinning() const {return true; }

  double lowerEdge() const { return lower; }

  double upperEdge() const { return upper; }

  int bins() const { return nbins; }

  double binLowerEdge(int index) const {
    return index < 0? -std::numeric_limits<double>::max():
      lower + double(std::min(index, nbins))*binWidth(0);
  }

  double binUpperEdge(int index) const {
    return index >= nbins? std::numeric_limits<double>::max():
      lower + double(std::max(index, -1) + 1)*binWidth(0);
  }

  double binWidth(int) const {
    return (upper - lower)/double(nbins);
  }

  int coordToIndex(double coord) const {
    if ( coord >= upper ) return OVERFLOW_BIN;
    else if ( coord < lower ) return UNDERFLOW_BIN;
    else return int((coord - lower)/binWidth(0));
  }

  double binMidPoint(int index) const {
    return lower + (double(index) + 0.5)*binWidth(0);
  }

private:

  double lower;

  double upper;

  int nbins;

};

}

#endif /* LWH_Axis_H */