---

AlignableSiStripDet Class Reference

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

#include <Alignment/TrackerAlignment/interface/AlignableSiStripDet.h>

Inheritance diagram for AlignableSiStripDet:

List of all members.

Public Member Functions
	AlignableSiStripDet (const GluedGeomDet *geomDet)
	Constructor.
virtual AlignmentErrors *	alignmentErrors () const
	first consistify with component detunits, then call method from AlignableDet
virtual Alignments *	alignments () const
	first consistify with component detunits, then call method from AlignableDet
virtual	~AlignableSiStripDet ()
	reduntantly make destructor virtual
Private Member Functions
void	consistifyAlignmentErrors ()
	make alignment errors consistent with daughters
void	consistifyAlignments ()
	make alignments consistent with daughters
const AlignTransformError &	errorFromId (const std::vector< AlignTransformError > &trafoErrs, align::ID id) const
	AlignTransformError with 'id'.
Private Attributes
const Bounds *	theMonoBounds
	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.
StripGeomDetType &	theMonoType
const Bounds *	theStereoBounds
StripGeomDetType &	theStereoType

---

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 26 of file AlignableSiStripDet.cc.

References 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 53 of file AlignableSiStripDet.cc.

References theMonoBounds, and theStereoBounds.

{
  delete theMonoBounds;
  delete theStereoBounds;
}

---

Member Function Documentation

AlignmentErrors * AlignableSiStripDet::alignmentErrors

(

void

)

const [virtual]

first consistify with component detunits, then call method from AlignableDet

Reimplemented from AlignableDet.

Definition at line 68 of file AlignableSiStripDet.cc.

References AlignableDet::alignmentErrors(), and consistifyAlignmentErrors().

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

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

Alignments * AlignableSiStripDet::alignments

(

void

)

const [virtual]

first consistify with component detunits, then call method from AlignableDet

Reimplemented from AlignableDet.

Definition at line 60 of file AlignableSiStripDet.cc.

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

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

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

void AlignableSiStripDet::consistifyAlignmentErrors

(

)

[private]

make alignment errors consistent with daughters

Definition at line 131 of file AlignableSiStripDet.cc.

References AlignableDet::alignmentErrors(), errorFromId(), AlignmentErrors::m_alignError, AlignableDet::setAlignmentPositionError(), and funct::sqrt().

Referenced by alignmentErrors().

{
  // make alignment errors consistent with daughters

  AlignmentErrors *oldErrs = this->AlignableDet::alignmentErrors();

  const Alignables units(this->components()); // order mono==0, stereo==1 does not matter here

  const AlignTransformError &gluedErr  = this->errorFromId(oldErrs->m_alignError,
                                                           this->geomDetId());
  const AlignTransformError &monoErr   = this->errorFromId(oldErrs->m_alignError,
                                                           units[0]->geomDetId());
  const AlignTransformError &stereoErr = this->errorFromId(oldErrs->m_alignError,
                                                           units[1]->geomDetId());
  const GlobalError errGlued (gluedErr.matrix());
  const GlobalError errMono  (monoErr.matrix());
  const GlobalError errStereo(stereoErr.matrix());

  //  const GlobalError newErrGlued((errMono + errStereo - errGlued).matrix_new() /= 4.);
  // The above line would be error propagation assuming:
  //   - Glued position is just the mean of its components.
  //   - Components APE is square sum of what has been set to glued and to components itself.
  // But this can be too small, e.g. 
  // - Simply by factor sqrt(4)=2 smaller than the old 'errGlued' in case APE (from position)
  //   was only set to glued (and 1-to-1 propagated to components).
  // - Factor sqrt(2) smaller than components APE in case APE (from position) was only
  //   directly applied to both components.
  //
  // So I choose the max of all three, ignoring correlations (which we do not have?), sigh!
  // And in this way it is safe against repetetive calls to this method!
  double maxX2 = (errMono.cxx() > errStereo.cxx() ? errMono.cxx() : errStereo.cxx());
  maxX2 = (maxX2 > errGlued.cxx() ? maxX2 : errGlued.cxx());
  double maxY2 = (errMono.cyy() > errStereo.cyy() ? errMono.cyy() : errStereo.cyy());
  maxY2 = (maxY2 > errGlued.cyy() ? maxY2 : errGlued.cyy());
  double maxZ2 = (errMono.czz() > errStereo.czz() ? errMono.czz() : errStereo.czz());
  maxZ2 = (maxZ2 > errGlued.czz() ? maxZ2 : errGlued.czz());
  const AlignmentPositionError newApeGlued(sqrt(maxX2), sqrt(maxY2), sqrt(maxZ2));

  // Now set new errors - and reset those of the components, since they get overwritten...
  this->setAlignmentPositionError(newApeGlued);
  units[0]->setAlignmentPositionError(AlignmentPositionError(errMono));
  units[1]->setAlignmentPositionError(AlignmentPositionError(errStereo));

//   edm::LogWarning("Alignment") << "@SUB=consistifyAlignmentErrors" 
//                             << "End Id " << this->geomDetId();
//   AlignmentErrors *newErrs = this->AlignableDet::alignmentErrors();
//   this->dumpCompareAPE(oldErrs->m_alignError, newErrs->m_alignError);
//   delete newErrs;

  delete oldErrs;
}

void AlignableSiStripDet::consistifyAlignments

(

)

[private]

make alignments consistent with daughters

Definition at line 76 of file AlignableSiStripDet.cc.

References Alignable::addDisplacement(), Alignable::addRotation(), PV3DBase< T, PVType, FrameType >::basicVector(), BoundPlane::build(), Alignable::globalPosition(), Alignable::globalRotation(), GloballyPositioned< T >::position(), Alignable::rotation(), GloballyPositioned< T >::rotation(), theMonoBounds, theMonoType, theStereoBounds, theStereoType, and Alignable::theSurface.

Referenced by alignments().

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

  // The aim of all this gymnastics is to have the alignments calculated by
  // 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 GeomDetUnits
  // that are passed to that routine.

  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);
  // Fortunately we do not seem to need a GeometricDet pointer and can use 0:
  const StripGeomDetUnit monoDet(&(*monoPlane), &theMonoType, 0);

  BoundPlane::BoundPlanePointer stereoPlane
    = BoundPlane::build(aliUnits[1]->globalPosition(), aliUnits[1]->globalRotation(),
                        *theStereoBounds);
  // Fortunately we do not seem to need a GeometricDet pointer and can use 0:
  const StripGeomDetUnit stereoDet(&(*stereoPlane), &theStereoType, 0);

  std::vector<const GeomDetUnit*> detComps;
  detComps.push_back(&monoDet);   // order mono first, stereo second should be as in...
  detComps.push_back(&stereoDet); // ...TrackerGeomBuilderFromGeometricDet::buildGeomDet

  // Now we have all to calculate new position and rotation via PlaneBuilderForGluedDet.
  const PositionType oldPos(theSurface.position()); // From old surface for keeping...
  const RotationType oldRot(theSurface.rotation()); // ...track of changes.

  PlaneBuilderForGluedDet planeBuilder;
  theSurface = AlignableSurface(*planeBuilder.plane(detComps));

  // But do not forget to keep track of movements/rotations:
  const GlobalVector movement(theSurface.position().basicVector() - oldPos.basicVector());
  // Seems to be correct down to delta angles 4.*1e-8:
  const RotationType rotation(oldRot.multiplyInverse(theSurface.rotation()));
  this->addDisplacement(movement);
  this->addRotation(rotation);

//   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();
//   }

}

const AlignTransformError & AlignableSiStripDet::errorFromId	(	const std::vector< AlignTransformError > &	trafoErrs,
		align::ID	id
	)			const [private]

AlignTransformError with 'id'.

Definition at line 185 of file AlignableSiStripDet.cc.

References Exception, and i.

Referenced by consistifyAlignmentErrors().

{
  for (unsigned int i = 0; i < trafoErrs.size(); ++i) {
    if (trafoErrs[i].rawId() ==  id) return trafoErrs[i];
  }

  throw cms::Exception("Mismatch") << "[AlignableSiStripDet::indexFromId] "
                                   << id << " not found.";

  return trafoErrs.front(); // never reached due to exception (but pleasing the compiler)
}

---

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 54 of file AlignableSiStripDet.h.

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

StripGeomDetType& AlignableSiStripDet::theMonoType [private]

Definition at line 56 of file AlignableSiStripDet.h.

Referenced by consistifyAlignments().

const Bounds* AlignableSiStripDet::theStereoBounds [private]

Definition at line 55 of file AlignableSiStripDet.h.

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

StripGeomDetType& AlignableSiStripDet::theStereoType [private]

Definition at line 57 of file AlignableSiStripDet.h.

Referenced by consistifyAlignments().

---

The documentation for this class was generated from the following files:

• Alignment/TrackerAlignment/interface/AlignableSiStripDet.h
• Alignment/TrackerAlignment/src/AlignableSiStripDet.cc

---