VectorXYZ.h

Go to the documentation of this file.

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