geo_vectors.h

Go to the documentation of this file.

#ifndef LARCOREOBJ_SIMPLETYPESANDCONSTANTS_GEO_VECTORS_H
#define LARCOREOBJ_SIMPLETYPESANDCONSTANTS_GEO_VECTORS_H

// ROOT libraries
#include "Math/GenVector/CoordinateSystemTags.h"
#include "Math/GenVector/Cartesian3D.h"
#include "Math/GenVector/PositionVector3D.h"
#include "Math/GenVector/DisplacementVector3D.h"
#include "Math/GenVector/Rotation3D.h"


// BEGIN -- GENVECTOR_CONSTEXPR issue ------------------------------------------
#include "RVersion.h"
#if ROOT_VERSION_CODE <= ROOT_VERSION(6,12,6)
#define GENVECTOR_CONSTEXPR
#else
# error ("Please check if new ROOT supports constexpr GenVector objects.")
# define GENVECTOR_CONSTEXPR constexpr // ok, this is useless!
#endif
// END -- GENVECTOR_CONSTEXPR issue --------------------------------------------


// BEGIN Geometry group --------------------------------------------------------

namespace geo {
  
  
  template <typename T, typename C = ROOT::Math::GlobalCoordinateSystemTag>
  using GenVector3DBase_t
    = ROOT::Math::DisplacementVector3D<ROOT::Math::Cartesian3D<T>, C>;
  
  template <typename T, typename C = ROOT::Math::GlobalCoordinateSystemTag>
  using GenPoint3DBase_t
    = ROOT::Math::PositionVector3D<ROOT::Math::Cartesian3D<T>, C>;
  
  template <typename C>
  using Vector3DBase_t = GenVector3DBase_t<double, C>;
  
  template <typename C>
  using Point3DBase_t = GenPoint3DBase_t<double, C>;
  
  

  using Length_t = double;
  
  using GlobalCoords = ROOT::Math::GlobalCoordinateSystemTag;
  
  // the actual definition commented out is not understood by GenReflex
//  using Vector_t = Vector3DBase_t<GlobalCoords>;
  using Vector_t = ROOT::Math::DisplacementVector3D
    <ROOT::Math::Cartesian3D<double>, ROOT::Math::GlobalCoordinateSystemTag>;
  
  // the actual definition commented out is not understood by GenReflex
//  using Point_t = Point3DBase_t<GlobalCoords>;
   using Point_t = ROOT::Math::PositionVector3D
     <ROOT::Math::Cartesian3D<double>, ROOT::Math::GlobalCoordinateSystemTag>;
  
  
  template <typename CoordSystemTag>
  using VectorIn_t = Vector3DBase_t<CoordSystemTag>;
  
  template <typename CoordSystemTag>
  using PointIn_t = Point3DBase_t<CoordSystemTag>;
  
  
  using Rotation_t = ROOT::Math::Rotation3D;
  
  
  //--------------------------------------------------------------------------
  template <typename Vector = Vector_t>
  constexpr Vector Xaxis() { return { 1.0, 0.0, 0.0 }; }
  
  template <typename Vector = Vector_t>
  constexpr Vector Yaxis() { return { 0.0, 1.0, 0.0 }; }
  
  template <typename Vector = Vector_t>
  constexpr Vector Zaxis() { return { 0.0, 0.0, 1.0 }; }
  
  template <typename Point = Point_t>
  constexpr Point origin() { return { 0.0, 0.0, 0.0 }; }
  
  
  
  //----------------------------------------------------------------------------

} // namespace geo
// END Geometry group ----------------------------------------------------------


#endif // LARCOREOBJ_SIMPLETYPESANDCONSTANTS_GEO_VECTORS_H