TransverseBoundPlaneFactory.cc

Go to the documentation of this file.

//#include "Utilities/Configuration/interface/Architecture.h"
 
#include "DataFormats/GeometryCommonDetAlgo/interface/TransverseBoundPlaneFactory.h"
 
BoundPlane* TransverseBoundPlaneFactory::operator()(const Surface::GlobalPoint& origin,
                                                    const Surface::GlobalVector& dir) const {
  // z axis perpendicular to transverse momentum
  GlobalVector zAxis = GlobalVector(dir.x(), dir.y(), 0.).unit();
 
  // y axis coincides with global z
  GlobalVector yAxis(0., 0., 1.);
 
  // x axis obtained by cross product
  GlobalVector xAxis = (yAxis.cross(zAxis)).unit();
 
  Surface::RotationType rot(
      xAxis.x(), xAxis.y(), xAxis.z(), yAxis.x(), yAxis.y(), yAxis.z(), zAxis.x(), zAxis.y(), zAxis.z());
 
  //  Surface::RotationType rot(yAxis, zAxis);
 
  return new BoundPlane(origin, rot);
}