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;
    }
}

#endif // NPSTAT_CLOSEWITHINTOLERANCE_HH_