ValidationUtil.h

Go to the documentation of this file.

#ifndef _ValidationUtil_h_
#define _ValidationUtil_h_

#include <cmath>
#include <TVector3.h>


inline
double
Angle(const double theta1, const double phi1,
      const double theta2, const double phi2)
{
  TVector3 a;
  a.SetMagThetaPhi(1, theta1, phi1);
  TVector3 b;
  b.SetMagThetaPhi(1, theta2, phi2);
  return 2*std::asin(0.5 * (a - b).Mag());
}


inline
double
RelativeError(const double a, const double b)
{
  return std::fabs(2 * (a - b) / (a + b));
}


#endif