PlaneBuilderForGluedDet Class Reference

#include <PlaneBuilderForGluedDet.h>

Public Types
typedef ReferenceCountingPointer < BoundPlane >	ResultType

Public Member Functions
std::pair < RectangularPlaneBounds, GlobalVector >	computeRectBounds (const std::vector< const GeomDetUnit * > &dets, const BoundPlane &plane) const
Surface::RotationType	computeRotation (const std::vector< GeomDetUnit * > &dets, const Surface::PositionType &meanPos) const
ResultType	plane (const std::vector< const GeomDetUnit * > &dets) const

Detailed Description

Builds the minimal rectangular box that contains all input GeomDetUnits fully.

Definition at line 16 of file PlaneBuilderForGluedDet.h.

Member Typedef Documentation

typedef ReferenceCountingPointer<BoundPlane> PlaneBuilderForGluedDet::ResultType

Warning, remember to assign this pointer to a ReferenceCountingPointer! Should be changed to return a ReferenceCountingPointer<BoundPlane>

Definition at line 21 of file PlaneBuilderForGluedDet.h.

Member Function Documentation

std::pair< RectangularPlaneBounds, GlobalVector > PlaneBuilderForGluedDet::computeRectBounds	(	const std::vector< const GeomDetUnit * > &	dets,
		const BoundPlane &	plane
	)		const

Definition at line 34 of file PlaneBuilderForGluedDet.cc.

References BoundingBox::corners(), i, evf::evtn::offset(), L1TEmulatorMonitor_cff::p, query::result, Surface::toGlobal(), GloballyPositioned< T >::toLocal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by plane().

                                                                                                                                                                  {
  // go over all corners and compute maximum deviations from mean pos.
  std::vector<GlobalPoint> corners;
  for (std::vector<const GeomDetUnit*>::const_iterator idet=dets.begin();
       idet != dets.end(); idet++) {
    const BoundPlane& plane = dynamic_cast<const BoundPlane&>((*idet)->surface());
    std::vector<GlobalPoint> dc = BoundingBox().corners(plane);
    corners.insert( corners.end(), dc.begin(), dc.end());
  }
  
  float xmin(0), xmax(0), ymin(0), ymax(0), zmin(0), zmax(0);
  for (std::vector<GlobalPoint>::const_iterator i=corners.begin();
       i!=corners.end(); i++) {
    LocalPoint p = plane.toLocal(*i);
    if (p.x() < xmin) xmin = p.x();
    if (p.x() > xmax) xmax = p.x();
    if (p.y() < ymin) ymin = p.y();
    if (p.y() > ymax) ymax = p.y();
    if (p.z() < zmin) zmin = p.z();
    if (p.z() > zmax) zmax = p.z();
  }

  LocalVector localOffset( (xmin+xmax)/2., (ymin+ymax)/2., (zmin+zmax)/2.);
  GlobalVector offset( plane.toGlobal(localOffset));

  std::pair<RectangularPlaneBounds, GlobalVector> result(RectangularPlaneBounds((xmax-xmin)/2, (ymax-ymin)/2, (zmax-zmin)/2), offset);

  return result;
}

Surface::RotationType PlaneBuilderForGluedDet::computeRotation	(	const std::vector< GeomDetUnit * > &	dets,
		const Surface::PositionType &	meanPos
	)		const

Definition at line 64 of file PlaneBuilderForGluedDet.cc.

References Vector3DBase< T, FrameTag >::cross(), n, plane(), Surface::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                                                                          {

  // choose first mono out-pointing rotation
  // the rotations of GluedDets coincide with the mono part
  // Simply take the x,y of the first Det if z points out,
  // or -x, y if it doesn't
  const BoundPlane& plane = dynamic_cast<const BoundPlane&>(dets.front()->surface());
  //GlobalVector n = plane.normalVector();

  GlobalVector xAxis;
  GlobalVector yAxis;
  GlobalVector planeYAxis = plane.toGlobal( LocalVector( 0, 1, 0));
  if (planeYAxis.z() < 0) yAxis = -planeYAxis;
  else                    yAxis =  planeYAxis;

  GlobalVector planeXAxis = plane.toGlobal( LocalVector( 1, 0, 0));
  GlobalVector n = planeXAxis.cross( planeYAxis);

  if (n.x() * meanPos.x() + n.y() * meanPos.y() > 0) {
    xAxis = planeXAxis;
  }
  else {
    xAxis = -planeXAxis;
  }

  return Surface::RotationType( xAxis, yAxis);
}

PlaneBuilderForGluedDet::ResultType PlaneBuilderForGluedDet::plane

(

const std::vector< const GeomDetUnit * > &

dets

)

const

Definition at line 11 of file PlaneBuilderForGluedDet.cc.

References BoundPlane::build(), computeRectBounds(), i, MediumProperties::radLen(), and MediumProperties::xi().

Referenced by TrackerGeomBuilderFromGeometricDet::buildGeomDet(), and computeRotation().

                                                                                                                 {
  // find mean position
  typedef Surface::PositionType::BasicVectorType Vector;
  Vector posSum(0,0,0);
  for (std::vector<const GeomDetUnit*>::const_iterator i=dets.begin(); i!=dets.end(); i++) {
    posSum += (**i).surface().position().basicVector();
  }
  Surface::PositionType meanPos( posSum/float(dets.size()));
  
  Surface::RotationType rotation =  dets.front()->surface().rotation();
  //  Surface::RotationType rotation = computeRotation( dets, meanPos);
  BoundPlane::BoundPlanePointer tmpPlane = BoundPlane::build( meanPos, rotation, OpenBounds());

  // Take the medium properties from the first DetUnit 
  const MediumProperties* mp = dets.front()->surface().mediumProperties();
  MediumProperties newmp(0,0);
  if (mp != 0) newmp = MediumProperties( mp->radLen()*2.0,mp->xi()*2.0);

  std::pair<RectangularPlaneBounds,GlobalVector> bo = computeRectBounds( dets, *tmpPlane);
  return new BoundPlane( meanPos+bo.second, rotation, bo.first, &newmp);
}