CoordinateSystem_ROOT.cc

Go to the documentation of this file.

#include <io/CoordinateSystem_ROOT.h>

#include <det/Detector.h>

#include <utl/Point.h>
#include <utl/Vector.h>
#include <utl/CoordinateSystemPtr.h>
#include <utl/CoordinateSystem.h>
#include <utl/TransformationMatrix.h>

#include <fwk/LocalCoordinateSystem.h>

#include <iostream>
using namespace std;

using namespace io;
using namespace utl;


CoordinateSystem_ROOT::CoordinateSystem_ROOT(const utl::CoordinateSystemPtr& cs)
{
  const utl::Point p(0, 0, 0, cs);
  const utl::CoordinateSystemPtr augerCS = fwk::LocalCoordinateSystem::Create(p);
  const utl::CoordinateSystemPtr siteCS = det::Detector::GetInstance().GetSiteCoordinateSystem();
  const utl::Point::Triple t = p.GetCoordinates(siteCS);
  fOrigin.SetXYZ(t.get<0>(), t.get<1>(), t.get<2>());
  const utl::TransformationMatrix M = cs->GetTransformationTo(augerCS);    
  fxx = M.fTransform.xx();
  fxy = M.fTransform.xy();
  fxz = M.fTransform.xz();
  fyx = M.fTransform.yx();
  fyy = M.fTransform.yy();
  fyz = M.fTransform.yz();
  fzx = M.fTransform.zx();
  fzy = M.fTransform.zy();
  fzz = M.fTransform.zz();
  fdx = M.fTransform.dx();
  fdy = M.fTransform.dy();
  fdz = M.fTransform.dz();
}


void
CoordinateSystem_ROOT::operator>>(utl::CoordinateSystemPtr& cs)
  const
{
  utl::Point p(fOrigin.x(), fOrigin.y(), fOrigin.z(),
               det::Detector::GetInstance().GetSiteCoordinateSystem());
  
  const utl::CoordinateSystemPtr augerCS = fwk::LocalCoordinateSystem::Create(p);
  
  class TC : public HepGeom::Transform3D {
  public:
    TC(const double xx, const double xy, const double xz, const double dx,
       const double yx, const double yy, const double yz, const double dy,
       const double zx, const double zy, const double zz, const double dz)
      : HepGeom::Transform3D(xx, xy, xz, dx,
                             yx, yy, yz, dy,
                             zx, zy, zz, dz) { }
  };

  const TransformationMatrix TM(TC(fxx, fxy, fxz, fdx,
                                   fyx, fyy, fyz, fdy,
                                   fzx, fzy, fzz, fdz));
  
  cs = utl::CoordinateSystem::FromTransformation(TM, augerCS);  
}