MatrixHelpers.h

Go to the documentation of this file.

//==========================================================================
//  AIDA Detector description implementation 
//--------------------------------------------------------------------------
// Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
// All rights reserved.
//
// For the licensing terms see $DD4hepINSTALL/LICENSE.
// For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
//
// Author     : M.Frank
//
//==========================================================================
#ifndef DD4HEP_MATRIXHELPERS_H
#define DD4HEP_MATRIXHELPERS_H
 
// Framework include files
#include <DD4hep/Objects.h>
 
// Forward declarations
class TGeoTranslation;
class TGeoRotation;
class TGeoHMatrix;
class TGeoMatrix;
 
namespace dd4hep {
 
  namespace detail   {
 
    namespace matrix {
 
      TGeoIdentity*    _identity();
      ROOT::Math::XYZVector _scale(const TGeoMatrix* matrix);
      ROOT::Math::XYZVector _scale(const TGeoMatrix& matrix);
 
      TGeoTranslation* _translation(const Position& pos);
      Position         _translation(const TGeoMatrix* matrix);
      Position         _translation(const TGeoMatrix& matrix);
 
      TGeoRotation*    _rotationZYX(const RotationZYX& rot);
      TGeoRotation*    _rotation3D(const Rotation3D& rot);
      Rotation3D       _rotation3D(const TGeoMatrix* matrix);
      Rotation3D       _rotation3D(const TGeoMatrix& matrix);
      
      TGeoHMatrix*     _transform(const Transform3D& trans);
      TGeoHMatrix*     _transform(const Position& pos);
      TGeoHMatrix*     _transform(const RotationZYX& rot);
      TGeoHMatrix*     _transform(const Rotation3D& rot3D);
      TGeoHMatrix*     _transform(const Position& pos, const RotationZYX& rot);
 
      TGeoHMatrix&     _transform(TGeoHMatrix& mat, const Transform3D& trans);
      TGeoHMatrix&     _transform(TGeoHMatrix& mat, const Position& pos);
      TGeoHMatrix&     _transform(TGeoHMatrix& mat, const RotationZYX& rot);
      TGeoHMatrix&     _transform(TGeoHMatrix& mat, const Rotation3D& rot3D);
      TGeoHMatrix&     _transform(TGeoHMatrix& mat, const Position& pos, const RotationZYX& rot);
 
      Transform3D      _transform(const TGeoMatrix* matrix);
      Transform3D      _transform(const TGeoMatrix& matrix);
      void             _transform(const TGeoMatrix* matrix, Transform3D& trafo);
      void             _transform(const TGeoMatrix& matrix, Transform3D& trafo);
 
      void _decompose(const TGeoMatrix& matrix, Position& pos, Rotation3D& rot);
      void _decompose(const TGeoMatrix& matrix, Position& pos, RotationZYX& rot);
      void _decompose(const TGeoMatrix& matrix, Position& pos, XYZAngles& rot);
      void _decompose(const Transform3D& matrix, Position& pos, Rotation3D& rot);
      void _decompose(const Transform3D& matrix, Position& pos, RotationZYX& rot);
      void _decompose(const Transform3D& matrix, Position& pos, XYZAngles& rot);
      void _decompose(const Transform3D& matrix, Translation3D& pos, RotationZYX& rot);
      void _decompose(const Transform3D& matrix, Translation3D& pos, XYZAngles& rot);
      void _decompose(const Rotation3D&  rot, Position& x, Position& y, Position& z);
 
      XYZAngles _xyzAngles(const double* matrix);
      XYZAngles _xyzAngles(const TGeoMatrix* matrix);
 
      double determinant(const TGeoMatrix* matrix);
      double determinant(const TGeoMatrix& matrix);
      double determinant(const Transform3D& trafo);
      double determinant(const Rotation3D& rot);
 
      enum MatrixEqualityCode  {
        MATRICES_EQUAL = 1,
        MATRICES_DIFFER_TRANSLATION = 2,
        MATRICES_DIFFER_ROTATION = 4
      };
 
      int _matrixEqual(const TGeoMatrix& left, const TGeoMatrix& right);
    }
  }
}       /* End namespace dd4hep            */
 
#endif // DD4HEP_MATRIXHELPERS_H