closeWithinTolerance.h

Go to the documentation of this file.

#ifndef NPSTAT_CLOSEWITHINTOLERANCE_HH_
#define NPSTAT_CLOSEWITHINTOLERANCE_HH_

#include <cmath>
#include <algorithm>
#include "JetMETCorrections/InterpolationTables/interface/NpstatException.h"

namespace npstat {
  inline bool closeWithinTolerance(const double& a, const double& b, const double& tol) {
    if (tol < 0.0)
      throw npstat::NpstatInvalidArgument(
          "In npstat::closeWithinTolerance: "
          "negative tolerance is not allowed");
    if (a == b)
      return true;
    else
      return fabs(a - b) / std::max(fabs(a), fabs(b)) <= tol;
  }
}  // namespace npstat

#endif  // NPSTAT_CLOSEWITHINTOLERANCE_HH_