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 */