---

GeometryAligner Class Reference

Class to update a given geometry with a set of alignments. More...

#include <Geometry/TrackingGeometryAligner/interface/GeometryAligner.h>

Inheritance diagram for GeometryAligner:

List of all members.

Public Member Functions
template<class C>
void	applyAlignments (C *geometry, const Alignments *alignments, const AlignmentErrors *alignmentErrors, const AlignTransform &globalCoordinates)

---

Detailed Description

Class to update a given geometry with a set of alignments.

Definition at line 28 of file GeometryAligner.h.

---

Member Function Documentation

template<class C>

void GeometryAligner::applyAlignments	(	C *	geometry,
		const Alignments *	alignments,
		const AlignmentErrors *	alignmentErrors,
		const AlignTransform &	globalCoordinates
	)			[inline]

Definition at line 41 of file GeometryAligner.h.

References edmNew::copy(), error, Alignments::m_align, AlignmentErrors::m_alignError, AlignTransform::rotation(), DetPositioner::setAlignmentPositionError(), DetPositioner::setGeomDetPosition(), and AlignTransform::translation().

Referenced by TrackerGeometryIntoNtuples::beginJob(), MuonGeometryIntoNtuples::beginJob(), AlignmentMonitorAsAnalyzer::beginJob(), AlignmentProducer::beginOfJob(), TrackerGeometryCompare::createROOTGeometry(), MuonAlignmentInputDB::newAlignableMuon(), TrackerDigiGeometryESModule::produce(), CSCGeometryESModule::produce(), DTGeometryESModule::produce(), MisalignedMuonESProducer::produce(), MisalignedTrackerESProducer::produce(), and AlignmentProducer::startingNewLoop().

{

  edm::LogInfo("Starting") << "Starting to apply alignments";

  // Preliminary checks (or we can't loop!)
  if ( alignments->m_align.size() != geometry->theMap.size() )
        throw cms::Exception("GeometryMismatch") 
          << "Size mismatch between geometry (size=" << geometry->theMap.size() 
          << ") and alignments (size=" << alignments->m_align.size() << ")";
  if ( alignments->m_align.size() != alignmentErrors->m_alignError.size() )
        throw cms::Exception("GeometryMismatch") 
          << "Size mismatch between geometry (size=" << geometry->theMap.size() 
          << ") and alignment errors (size=" << alignmentErrors->m_alignError.size() << ")";

  const AlignTransform::Translation &globalShift = globalCoordinates.translation();
  const AlignTransform::Rotation globalRotation = globalCoordinates.rotation(); // by value!

  // Parallel loop on alignments, alignment errors and geomdets
  std::vector<AlignTransform>::const_iterator iAlign = alignments->m_align.begin();
  std::vector<AlignTransformError>::const_iterator 
        iAlignError = alignmentErrors->m_alignError.begin();
  //copy  geometry->theMap to a real map to order it....
  std::map<unsigned int, GeomDet*> theMap;
  std::copy(geometry->theMap.begin(), geometry->theMap.end(), std::inserter(theMap,theMap.begin()));
  for ( std::map<unsigned int, GeomDet*>::const_iterator iPair = theMap.begin(); 
                iPair != theMap.end(); ++iPair, ++iAlign, ++iAlignError )
        {
          // Check DetIds
          if ( (*iPair).first != (*iAlign).rawId() )
            throw cms::Exception("GeometryMismatch") 
              << "DetId mismatch between geometry (rawId=" << (*iPair).first
              << ") and alignments (rawId=" << (*iAlign).rawId();
          
          if ( (*iPair).first != (*iAlignError).rawId() )
            throw cms::Exception("GeometryMismatch") 
              << "DetId mismatch between geometry (rawId=" << (*iPair).first
              << ") and alignment errors (rawId=" << (*iAlignError).rawId();

          // Apply global correction
          CLHEP::Hep3Vector positionHep = globalRotation * CLHEP::Hep3Vector( (*iAlign).translation() ) + globalShift;
          CLHEP::HepRotation rotationHep = globalRotation * CLHEP::HepRotation( (*iAlign).rotation() );

          // Define new quantities
          Surface::PositionType position( positionHep.x(), positionHep.y(), positionHep.z() );
          Surface::RotationType rotation( rotationHep.xx(), rotationHep.xy(), rotationHep.xz(), 
                                          rotationHep.yx(), rotationHep.yy(), rotationHep.yz(), 
                                          rotationHep.zx(), rotationHep.zy(), rotationHep.zz() );

          // Alignment Position Error
          GlobalError error( (*iAlignError).matrix() );
          AlignmentPositionError ape( error );

          // Apply new quantities
          GeomDet* iGeomDet = (*iPair).second;
          DetPositioner::setGeomDetPosition( *iGeomDet, position, rotation );
          DetPositioner::setAlignmentPositionError( *iGeomDet, ape );
        }

  edm::LogInfo("Done") << "Finished to apply alignments";
}

---

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

• Geometry/TrackingGeometryAligner/interface/GeometryAligner.h

---