Direction.h

Go to the documentation of this file.

// AUTOMATICALLY GENERATED FILE - DO NOT EDIT
 
#ifndef EICD_Direction_H
#define EICD_Direction_H
 
#include <cmath>
#include <tuple>
#include <ostream>
 
namespace eic {
 
 
class Direction {
public:
  float theta{}; ///< [rad, 0->pi]
  float phi{}; ///< [rad, -pi->pi]
 
 Direction() : theta{0}, phi{0} {}
 Direction(double th, double ph) : theta{static_cast<float>(th)}, phi{static_cast<float>(ph)} {}
 Direction(double x, double y, double z)
 : theta{static_cast<float>(acos(z/std::hypot(x,y,z)))}
 , phi{static_cast<float>(atan2(y,x))} {}
 template <class VectorType> Direction(const VectorType& v) : Direction(v.theta(), v.phi()) {}
 operator std::pair<float, float>() const {return {theta, phi};}
 float eta() const {return -log(tan(0.5*theta));}
 Direction add(const Direction& rhs) const {return {theta+rhs.theta, phi+rhs.phi};}
 Direction subtract(const Direction& rhs) const {return {theta-rhs.theta, phi-rhs.phi};}
 
 
};
 
inline std::ostream& operator<<(std::ostream& o, const eic::Direction& value) {
  o << value.theta << " ";
  o << value.phi << " ";
 
  return o;
}
 
} // namespace eic
 
 
 
#endif