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/Cartesian3D.h"
#include "Math/GenVector/CoordinateSystemTags.h"
#include "Math/GenVector/DisplacementVector3D.h"
#include "Math/GenVector/PositionVector3D.h"
#include "Math/GenVector/Rotation3D.h"

// BEGIN -- GENVECTOR_CONSTEXPR issue ------------------------------------------
#include "RVersion.h"
#define GENVECTOR_CONSTEXPR
// use the following definition if the issue is ever resolved
// # define GENVECTOR_CONSTEXPR constexpr
// 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