BasicVector.cc

Go to the documentation of this file.

#include <utl/config.h>
#include <utl/BasicVector.h>
#include <utl/TransformationMatrix.h>
#include <CLHEP/Geometry/Point3D.h>
#include <CLHEP/Geometry/Vector3D.h>
#include <cmath>

using std::sin;
using std::cos;


namespace utl {

  template<class HepVector>
  HepVector
  BasicVector<HepVector>::TransformedVector(const CoordinateSystemPtr& newCoordinateSystem)
    const
  {
    CoordinateSystemValid(newCoordinateSystem);
    return HepVector(fCoordinateSystem->GetTransformationTo(newCoordinateSystem).fTransform * fVector);
  }


  template<class HepVector>
  double
  BasicVector<HepVector>::CCylindrical::X(const double rho,
                                          const double phi,
                                          const double /*z*/)
    const
  {
    return rho * cos(phi);
  }


  template<class HepVector>
  double
  BasicVector<HepVector>::CCylindrical::Y(const double rho,
                                          const double phi,
                                          const double /*z*/)
    const
  {
    return rho * sin(phi);
  }


  template<class HepVector>
  double
  BasicVector<HepVector>::CCylindrical::Z(const double /*rho*/,
                                          const double /*phi*/,
                                          const double z)
    const
  {
    return z;
  }


  template<class HepVector>
  double
  BasicVector<HepVector>::CSpherical::X(const double r,
                                        const double theta,
                                        const double phi)
    const
  {
    return r * cos(phi) * sin(theta);
  }


  template<class HepVector>
  double
  BasicVector<HepVector>::CSpherical::Y(const double r,
                                        const double theta,
                                        const double phi)
    const
  {
    return r * sin(phi) * sin(theta);
  }


  template<class HepVector>
  double
  BasicVector<HepVector>::CSpherical::Z(const double r,
                                        const double theta,
                                        const double /*phi*/)
    const
  {
    return r * cos(theta);
  }


  template<class HepVector>
  const typename BasicVector<HepVector>::CCartesian BasicVector<HepVector>::kCartesian;


  template<class HepVector>
  const typename BasicVector<HepVector>::CCylindrical BasicVector<HepVector>::kCylindrical;


  template<class HepVector>
  const typename BasicVector<HepVector>::CSpherical BasicVector<HepVector>::kSpherical;


  template<class HepVector>
  const double BasicVector<HepVector>::fgEpsilon = 1e-9;


  template<class HepVector>
  const char* const BasicVector<HepVector>::CoordinateType::fgKindTags[] = { "cartesian", "cylindrical", "spherical" };


  // instantiations
  template class BasicVector<HepGeom::Point3D<double>>;
  template class BasicVector<HepGeom::Vector3D<double>>;
  template class BasicVector<HepGeom::Normal3D<double>>;

}