AlignableSiStripDet Class Reference

#include <AlignableSiStripDet.h>

Inheritance diagram for AlignableSiStripDet:

Public Member Functions
	AlignableSiStripDet (const GluedGeomDet *geomDet)
	Constructor. More...
virtual Alignments *	alignments () const
	first consistify with component detunits, then call method from AlignableDet More...
virtual	~AlignableSiStripDet ()
	reduntantly make destructor virtual More...
Public Member Functions inherited from AlignableDet
virtual void	addAlignmentPositionError (const AlignmentPositionError &ape, bool propagateDown)
virtual void	addAlignmentPositionErrorFromRotation (const RotationType &rot, bool propagateDown)
	AlignableDet (const GeomDet *geomDet, bool addComponents=true)
virtual AlignmentErrors *	alignmentErrors () const
	Return vector of alignment errors. More...
const AlignmentPositionError *	alignmentPositionError () const
	alignment position error - for checking only, otherwise use alignmentErrors() above! More...
virtual void	setAlignmentPositionError (const AlignmentPositionError &ape, bool propagateDown)
	Set the AlignmentPositionError and, if (propagateDown), to all components. More...
virtual	~AlignableDet ()
	Destructor. More...
Public Member Functions inherited from AlignableComposite
virtual void	addAlignmentPositionErrorFromLocalRotation (const RotationType &rotation, bool propagateDown)
virtual void	addComponent (Alignable *component)
virtual void	addSurfaceDeformation (const SurfaceDeformation *deformation, bool propagateDown)
	AlignableComposite (align::ID id, StructureType aType, const RotationType &rot=RotationType())
virtual StructureType	alignableObjectId () const
	Return the alignable type identifier. More...
virtual Alignables	components () const
	Return vector of direct components. More...
virtual void	dump () const
	Recursive printout of alignable structure. More...
virtual void	move (const GlobalVector &displacement)
	Move with respect to the global reference frame. More...
virtual void	moveComponentLocal (const int i, const LocalVector &localDisplacement)
	Move a single component with respect to the local reference frame. More...
virtual void	moveComponentsLocal (const LocalVector &localDisplacement)
	Move with respect to the local reference frame. More...
virtual void	recursiveComponents (Alignables &result) const
	Provide all components, subcomponents etc. (cf. description in base class) More...
virtual void	rotateInGlobalFrame (const RotationType &rotation)
	Rotation interpreted in global reference frame. More...
virtual void	setSurfaceDeformation (const SurfaceDeformation *deformation, bool propagateDown)
	Set the surface deformation parameters - if (!propagateDown) do not affect daughters. More...
virtual int	surfaceDeformationIdPairs (std::vector< std::pair< int, SurfaceDeformation * > > &) const
	Return surface deformations. More...
virtual	~AlignableComposite ()
	deleting its components More...
Public Member Functions inherited from Alignable
	Alignable (align::ID, const AlignableSurface &)
	Alignable (align::ID, const RotationType &)
AlignmentParameters *	alignmentParameters () const
	Get the AlignmentParameters. More...
const Alignables &	deepComponents () const
const GlobalVector &	displacement () const
	Return change of the global position since the creation of the object. More...
bool	firstCompsWithParams (Alignables &paramComps) const
const DetId &	geomDetId () const
const PositionType &	globalPosition () const
	Return the global position of the object. More...
const RotationType &	globalRotation () const
	Return the global orientation of the object. More...
align::ID	id () const
	Return the ID of Alignable, i.e. DetId of 'first' component GeomDet(Unit). More...
Alignable *	mother () const
	Return pointer to container alignable (if any) More...
virtual void	rotateAroundGlobalAxis (const GlobalVector &axis, Scalar radians)
	Rotation around arbitratry global axis. More...
virtual void	rotateAroundGlobalX (Scalar radians)
	Rotation around global x-axis. More...
virtual void	rotateAroundGlobalY (Scalar radians)
	Rotation around global y-axis. More...
virtual void	rotateAroundGlobalZ (Scalar radians)
	Rotation around global z-axis. More...
virtual void	rotateAroundLocalAxis (const LocalVector &axis, Scalar radians)
	Rotation around arbitratry local axis. More...
virtual void	rotateAroundLocalX (Scalar radians)
	Rotation around local x-axis. More...
virtual void	rotateAroundLocalY (Scalar radians)
	Rotation around local y-axis. More...
virtual void	rotateAroundLocalZ (Scalar radians)
	Rotation around local z-axis. More...
virtual void	rotateInLocalFrame (const RotationType &rotation)
	Rotation intepreted in the local reference frame. More...
const RotationType &	rotation () const
	Return change of orientation since the creation of the object. More...
void	setAlignmentParameters (AlignmentParameters *dap)
	Set the AlignmentParameters. More...
void	setMother (Alignable *mother)
	Assign mother to alignable. More...
void	setSurvey (const SurveyDet *)
	Set survey info. More...
const int	size () const
	Return number of direct components. More...
const AlignableSurface &	surface () const
	Return the Surface (global position and orientation) of the object. More...
AlignmentSurfaceDeformations *	surfaceDeformations () const
	Return surface deformations, sorted by DetId. More...
const SurveyDet *	survey () const
	Return survey info. More...
virtual	~Alignable ()
	Destructor. More...

Private Member Functions
void	consistifyAlignments ()
	make alignments consistent with daughters More...

Private Attributes
const Bounds *	theMonoBounds
StripGeomDetType &	theMonoType
const Bounds *	theStereoBounds
StripGeomDetType &	theStereoType

Additional Inherited Members
Public Types inherited from Alignable
typedef align::Alignables	Alignables
typedef align::GlobalVector	GlobalVector
typedef align::LocalVector	LocalVector
typedef align::PositionType	PositionType
typedef align::RotationType	RotationType
typedef align::Scalar	Scalar
typedef align::StructureType	StructureType
Protected Member Functions inherited from AlignableComposite
	AlignableComposite (const GeomDet *geomDet)
	Constructor from GeomDet, only for use in AlignableDet. More...
void	setSurface (const AlignableSurface &s)
Protected Member Functions inherited from Alignable
void	addDisplacement (const GlobalVector &displacement)
void	addRotation (const RotationType &rotation)
Protected Attributes inherited from AlignableComposite
StructureType	theStructureType
Protected Attributes inherited from Alignable
Alignables	theDeepComponents
DetId	theDetId
GlobalVector	theDisplacement
align::ID	theId
RotationType	theRotation
AlignableSurface	theSurface

Detailed Description

An alignable for GluedDets in Strip tracker, taking care of consistency with AlignableDet components.

First implementation April/May 2008

Author

Gero Flucke, Hamburg University

Date:

2008/07/13 12:26:41

Revision:

1.4

Definition at line 23 of file AlignableSiStripDet.h.

Constructor & Destructor Documentation

AlignableSiStripDet::AlignableSiStripDet

(

const GluedGeomDet *

geomDet

)

Constructor.

Definition at line 29 of file AlignableSiStripDet.cc.

References AlignableComposite::components(), edm::hlt::Exception, GeomDet::geographicalId(), GluedGeomDet::monoDet(), and GluedGeomDet::stereoDet().

  : AlignableDet(gluedDet, true), // true: adding DetUnits
    theMonoBounds  (gluedDet->monoDet()  ->surface().bounds().clone()),
    theStereoBounds(gluedDet->stereoDet()->surface().bounds().clone()),
    theMonoType  (static_cast<const StripGeomDetUnit*>(gluedDet->monoDet())  ->specificType()),
    theStereoType(static_cast<const StripGeomDetUnit*>(gluedDet->stereoDet())->specificType())
{
  // It is not allowed to store a pointer to GeomDet within objects with a life time
  // longer than an Event: 
  // GeomDet comes from TrackerGeometry that is created from GeometricDet that depends on
  // IdealGeometryRecord from EventSetup, so it could change in next event!
  //  ==> Need to store directly what I need from it. 
  // Unfortunately the current way with references for the type does not solve that,
  // either. But currently no way out, see header file.

  // check order mono/stereo
  const Alignables units(this->components());
  if (units.size() != 2
      || gluedDet->monoDet()->geographicalId() != units[0]->geomDetId()
      || gluedDet->stereoDet()->geographicalId() != units[1]->geomDetId()) {
    throw cms::Exception("LogicError")
      << "[AlignableSiStripDet] " << "Either != 2 components or "
      << "mono/stereo in wrong order for consistifyAlignments.";
  }
}

AlignableSiStripDet::~AlignableSiStripDet ( )

virtual

reduntantly make destructor virtual

Definition at line 56 of file AlignableSiStripDet.cc.

References theMonoBounds, and theStereoBounds.

{
  delete theMonoBounds;
  delete theStereoBounds;
}

Member Function Documentation

Alignments * AlignableSiStripDet::alignments ( void  ) const

virtual

first consistify with component detunits, then call method from AlignableDet

Reimplemented from AlignableDet.

Definition at line 63 of file AlignableSiStripDet.cc.

References AlignableDet::alignments(), and consistifyAlignments().

{
  const_cast<AlignableSiStripDet*>(this)->consistifyAlignments();

  return this->AlignableDet::alignments();
}

void AlignableSiStripDet::consistifyAlignments ( )

private

make alignments consistent with daughters

Definition at line 71 of file AlignableSiStripDet.cc.

References Alignable::addDisplacement(), Alignable::addRotation(), BoundPlane::build(), AlignableComposite::components(), BoundingBox::corners(), Alignable::globalPosition(), Alignable::globalRotation(), i, evf::evtn::offset(), L1TEmulatorMonitor_cff::p, GloballyPositioned< T >::position(), MediumProperties::radLen(), GloballyPositioned< T >::rotation(), Alignable::rotation(), theMonoBounds, Alignable::theRotation, theStereoBounds, Alignable::theSurface, PV3DBase< T, PVType, FrameType >::x(), MediumProperties::xi(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by alignments().

{
  // make alignments consistent with daughters, calling method from geometry

  // The aim of all this gymnastics is to have the alignments calculated the same way as
  // in PlaneBuilderForGluedDet::plane(const std::vector<const GeomDetUnit*> &detComps);
  // 
  // So we take the (new) position and orientation from the AligableDetUnits,
  // but bounds and GeomDetType from original GeomDetUnits to create new a new glued
  // surface.

  const PositionType oldPos(theSurface.position()); // From old surface for keeping...
  const RotationType oldRot(theSurface.rotation()); // ...track of changes.

  const Alignables aliUnits(this->components()); // order mono==0, stereo==1 checked in ctr.
  
  BoundPlane::BoundPlanePointer monoPlane
    = BoundPlane::build(aliUnits[0]->globalPosition(), aliUnits[0]->globalRotation(),
            *theMonoBounds);

  BoundPlane::BoundPlanePointer stereoPlane
    = BoundPlane::build(aliUnits[1]->globalPosition(), aliUnits[1]->globalRotation(),
            *theStereoBounds);
  
  Surface::PositionType::BasicVectorType posSum = monoPlane->position().basicVector();
  posSum += stereoPlane->position().basicVector();
  Surface::PositionType meanPos(posSum/(float)aliUnits.size());
  Surface::RotationType rotation =  monoPlane->rotation();
  
  BoundPlane::BoundPlanePointer tmpPlane = BoundPlane::build(meanPos, rotation, OpenBounds());

  const MediumProperties* mp = monoPlane->mediumProperties();
  MediumProperties newmp(0,0);
  if (mp != 0) newmp = MediumProperties(mp->radLen()*2.0,mp->xi()*2.0);

  std::vector<GlobalPoint> corners;
  std::vector<GlobalPoint> monoDC = BoundingBox().corners(*monoPlane);
  corners.insert(corners.end(), monoDC.begin(), monoDC.end());
  std::vector<GlobalPoint> stereoDC = BoundingBox().corners(*stereoPlane);
  corners.insert(corners.end(), stereoDC.begin(), stereoDC.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 = tmpPlane->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(tmpPlane->toGlobal(localOffset));

  theSurface = AlignableSurface(BoundPlane(meanPos+offset, rotation, 
                       RectangularPlaneBounds((xmax-xmin)/2, (ymax-ymin)/2, (zmax-zmin)/2),
                       &newmp));
  
  // But do not forget to keep track of movements/rotations:
  const GlobalVector theMovement(theSurface.position().basicVector() - oldPos.basicVector());
  // Seems to be correct down to delta angles 4.*1e-8:
  const RotationType theRotation(oldRot.multiplyInverse(theSurface.rotation()));
  this->addDisplacement(theMovement);
  this->addRotation(theRotation);
  
//   this->dumpCompareEuler(oldRot, theSurface.rotation());

//   if (movement.mag2()) { // > 1.e-10) { 
//     edm::LogWarning("Alignment") << "@SUB=consistifyAlignments" 
//                   << "Delta: " << movement.x() << " " << movement.y() << " " << movement.z()
//               << "\nPos: " << oldPos.perp() << " " << oldPos.phi() << " " << oldPos.z();
//   }
}

Member Data Documentation

const Bounds* AlignableSiStripDet::theMonoBounds

private

The following four members are needed to recalculate the surface in consistifyAlignments, to get rid of a GluedDet* which is disregarded since it could become an invalid pointer in the next event (theoretically...). But this solution is not better, the references for the types would become invalid together with the GeomDets they are taken from. StripGeomDetType has neither clone() and nor a decent copy constructor, so I cannot go the the same way as for the bounds. Sigh!

Definition at line 45 of file AlignableSiStripDet.h.

Referenced by consistifyAlignments(), and ~AlignableSiStripDet().

StripGeomDetType& AlignableSiStripDet::theMonoType

private

Definition at line 47 of file AlignableSiStripDet.h.

const Bounds* AlignableSiStripDet::theStereoBounds

private

Definition at line 46 of file AlignableSiStripDet.h.

Referenced by consistifyAlignments(), and ~AlignableSiStripDet().

StripGeomDetType& AlignableSiStripDet::theStereoType

private

Definition at line 48 of file AlignableSiStripDet.h.