VectorXYZT.h

Go to the documentation of this file.

// AUTOMATICALLY GENERATED FILE - DO NOT EDIT
 
#ifndef EICD_VectorXYZT_H
#define EICD_VectorXYZT_H
 
#include <cmath>
#include <tuple>
#include <ostream>
 
namespace eic {
 
 
class VectorXYZT {
public:
  double x{}; ///< [mm] or [GeV]
  double y{};
  double z{};
  double t{}; ///< [ns] or [GeV]
 
 VectorXYZT() : x{0}, y{0}, z{0}, t{0} {}
 VectorXYZT(double xx, double yy, double zz, double tt) : x{xx}, y{yy}, z{zz}, t{tt} {}
 double& operator[](unsigned i) {return *(&x + i);}
 const double& operator[](unsigned i) const {return *(&x + i);}
 double mag() const {return std::hypot(x, y, z);}
 double r() const {return mag();}
 double theta() const {return acos(z/mag());}
 double phi() const {return atan2(y,x);}
 double eta() const {return -log(tan(0.5*theta()));} double energy() const {return t;}
 double mass() const {return sqrt(t*t - x*x - y*y - z*z);}
 operator std::tuple<double, double, double, double>() {return {x, y, z, t};}
 double dot(const VectorXYZT& rhs) const {return t*rhs.t - x*rhs.x - y*rhs.y - z*rhs.z;}
 VectorXYZT add(const VectorXYZT& rhs) const {return {x+rhs.x, y+rhs.y, z+rhs.z, t+rhs.t};}
 VectorXYZT subtract(const VectorXYZT& rhs) const {return {x-rhs.x, y-rhs.y, z-rhs.z, t-rhs.t};}
 VectorXYZT scale(double f) const {return {f*x, f*y, f*z, f*t};}
 
 
};
 
inline std::ostream& operator<<(std::ostream& o, const eic::VectorXYZT& value) {
  o << value.x << " ";
  o << value.y << " ";
  o << value.z << " ";
  o << value.t << " ";
 
  return o;
}
 
} // namespace eic
 
 
 
#endif