#include <AlignableNavigator.h>

Public Types
typedef std::map< DetId, AlignableDetOrUnitPtr >	MapType

typedef MapType::value_type	PairType

Public Member Functions
std::vector< AlignableDetOrUnitPtr >	alignableDetOrUnits ()
	return all AlignableDetOrUnitPtrs More...

AlignableDetOrUnitPtr	alignableFromDetId (const DetId &detid)
	Returns AlignableDetOrUnitPtr corresponding to given DetId. More...

AlignableDetOrUnitPtr	alignableFromGeomDet (const GeomDet *geomDet)
	Returns AlignableDetOrUnitPtr corresponding to given GeomDet. More...

	AlignableNavigator (Alignable tracker, Alignable muon=0)
	Constructor from one or two Alignables. More...

	AlignableNavigator (AlignableExtras extras, Alignable tracker, Alignable *muon=0)
	Constructor from one or two Alignables. More...

	AlignableNavigator (const align::Alignables &)
	Constructor from list of Alignbable. More...

std::vector< AlignableDetOrUnitPtr >	alignablesFromHits (const std::vector< const TransientTrackingRecHit * > &hitvec)
	Returns vector AlignableDetOrUnitPtr for given vector of Hits. More...

std::vector< AlignableDetOrUnitPtr >	alignablesFromHits (const TransientTrackingRecHit::ConstRecHitContainer &hitVec)
	Returns vector of AlignableDetOrUnitPtr for given vector of Hits. More...

bool	detAndSubdetInMap (const DetId &detid) const
	Given a DetId, returns true if DetIds with this detector and subdetector id are in the map (not necessarily the exact DetId) More...

int	size (void)
	Returns number of elements in map. More...

Private Member Functions
unsigned int	recursiveGetId (Alignable *alignable)

Private Attributes
std::vector< std::pair< int, int > >	theDetAndSubdet

MapType	theMap

Detailed Description

A class to navigate from a DetId to the corresponding AlignableDetOrUnitPtr. A map is created at construction time from all sub-structures of the constructor's argument(s) that are AlignableDet or AlignableDetUnit.

Definition at line 26 of file AlignableNavigator.h.

Member Typedef Documentation

typedef std::map<DetId, AlignableDetOrUnitPtr> AlignableNavigator::MapType

Definition at line 40 of file AlignableNavigator.h.

typedef MapType::value_type AlignableNavigator::PairType

Definition at line 41 of file AlignableNavigator.h.

Constructor & Destructor Documentation

AlignableNavigator::AlignableNavigator	(	Alignable *	tracker,
		Alignable *	muon = `0`
	)

explicit

Constructor from one or two Alignables.

Definition at line 18 of file AlignableNavigator.cc.

References recursiveGetId(), and theMap.

 {
   theMap.clear();
 
   const unsigned int numNonDets = this->recursiveGetId(tracker) + this->recursiveGetId(muon);
 
   if (numNonDets) {
     edm::LogWarning("Alignment") <<"@SUB=AlignableNavigator" << "Created with map of size "
                                  << theMap.size() << ", but found also " << numNonDets 
                                  << " Alignables that have DetId!=0,\nbeing neither "
                  << "AlignableDet nor AlignableDetUnit. This will "
                                  << "lead to an exception in case alignableFromDetId(..) "
                  << "is called for one of these DetIds.\n" 
                                  << "If there is no exception, you can ignore this message.";
   } else {
     edm::LogInfo("Alignment") <<"@SUB=AlignableNavigator" << "Created with map of size "
                               << theMap.size() << ".";
   }
 }

AlignableNavigator::AlignableNavigator	(	AlignableExtras *	extras,
		Alignable *	tracker,
		Alignable *	muon = `0`
	)

explicit

Constructor from one or two Alignables.

Definition at line 39 of file AlignableNavigator.cc.

References AlignableExtras::components(), recursiveGetId(), and theMap.

 {
   theMap.clear();
 
   unsigned int numNonDets = this->recursiveGetId(tracker) + this->recursiveGetId(muon);
 
   if (extras) {
     for (const auto& it: extras->components()) {
       numNonDets += this->recursiveGetId(it);
     }
   }
 
   if (numNonDets) {
     edm::LogWarning("Alignment") <<"@SUB=AlignableNavigator" << "Created with map of size "
                                  << theMap.size() << ", but found also " << numNonDets 
                                  << " Alignables that have DetId!=0,\nbeing neither "
                  << "AlignableDet nor AlignableDetUnit. This will "
                                  << "lead to an exception in case alignableFromDetId(..) "
                  << "is called for one of these DetIds.\n" 
                                  << "If there is no exception, you can ignore this message.";
   } else {
     edm::LogInfo("Alignment") <<"@SUB=AlignableNavigator" << "Created with map of size "
                               << theMap.size() << ".";
   }
 }

AlignableNavigator::AlignableNavigator ( const align::Alignables & alignables )

explicit

Constructor from list of Alignbable.

Definition at line 66 of file AlignableNavigator.cc.

References recursiveGetId(), and theMap.

 {
   theMap.clear();
 
   unsigned int numNonDets = 0;
   for (const auto& it: alignables) {
     numNonDets += this->recursiveGetId(it);
   }
   if (numNonDets) {
     edm::LogWarning("Alignment") <<"@SUB=AlignableNavigator" << "Created with map of size "
                                  << theMap.size() << ", but found also " << numNonDets 
                                  << " Alignables that have DetId!=0,\nbeing neither "
                  << "AlignableDet nor AlignableDetUnit. This will "
                                  << "lead to an exception in case alignableFromDetId(..) "
                  << "is called for one of these DetIds.\n" 
                                  << "If there is no exception, you can ignore this message.";
   } else {
     edm::LogInfo("Alignment") <<"@SUB=AlignableNavigator" << "created with map of size "
                               << theMap.size() << ".";
   }
 }

Member Function Documentation

std::vector< AlignableDetOrUnitPtr > AlignableNavigator::alignableDetOrUnits ( )

return all AlignableDetOrUnitPtrs

Definition at line 194 of file AlignableNavigator.cc.

References mps_fire::result, and theMap.

 {
   std::vector<AlignableDetOrUnitPtr> result;
   result.reserve(theMap.size());
 
   for (MapType::const_iterator iIdAli = theMap.begin(); iIdAli != theMap.end(); ++iIdAli) {
     result.push_back(iIdAli->second);
   }
 
   return result;
 }

AlignableDetOrUnitPtr AlignableNavigator::alignableFromDetId ( const DetId & detid )

Returns AlignableDetOrUnitPtr corresponding to given DetId.

Definition at line 95 of file AlignableNavigator.cc.

References Exception, position, DetId::rawId(), and theMap.

Referenced by alignableFromGeomDet(), alignablesFromHits(), SurveyInputCSCfromPins::analyze(), MuonGeometrySVGTemplate::analyze(), CSCAlignmentCorrections::applyAlignment(), AlignmentMonitorTemplate::event(), AlignmentMonitorGeneric::event(), ApeSettingAlgorithm::initialize(), MuonAlignment::moveAlignableGlobalCoord(), MuonAlignment::moveAlignableLocalCoord(), MuonResidualsFromTrack::MuonResidualsFromTrack(), CSCChamberFitter::radiusCorrection(), and ApeSettingAlgorithm::terminate().

 {
   MapType::iterator position = theMap.find( detid );
   if ( position != theMap.end() ) {
     return position->second;
   }
    
   throw cms::Exception("BadLogic") 
     << "[AlignableNavigator::alignableDetFromDetId] DetId " << detid.rawId() << " not found";
 
   return static_cast<AlignableDet*>(nullptr);
 }

AlignableDetOrUnitPtr AlignableNavigator::alignableFromGeomDet ( const GeomDet * geomDet )

Returns AlignableDetOrUnitPtr corresponding to given GeomDet.

Definition at line 89 of file AlignableNavigator.cc.

References alignableFromDetId(), and GeomDet::geographicalId().

Referenced by MuonMillepedeAlgorithm::run().

 {
   return alignableFromDetId( geomDet->geographicalId() );
 }

std::vector< AlignableDetOrUnitPtr > AlignableNavigator::alignablesFromHits ( const std::vector< const TransientTrackingRecHit * > & hitvec )

Returns vector AlignableDetOrUnitPtr for given vector of Hits.

Definition at line 163 of file AlignableNavigator.cc.

References alignableFromDetId(), and mps_fire::result.

Referenced by MuonMillepedeAlgorithm::run().

 {
   std::vector<AlignableDetOrUnitPtr> result;
   result.reserve(hitvec.size());
 
   for(std::vector<const TransientTrackingRecHit*>::const_iterator ih
         = hitvec.begin(), iEnd = hitvec.end(); ih != iEnd; ++ih) {
     result.push_back(this->alignableFromDetId((*ih)->geographicalId()));
   }
 
   return result;
 }

std::vector< AlignableDetOrUnitPtr > AlignableNavigator::alignablesFromHits ( const TransientTrackingRecHit::ConstRecHitContainer & hitVec )

Returns vector of AlignableDetOrUnitPtr for given vector of Hits.

Definition at line 179 of file AlignableNavigator.cc.

References alignableFromDetId(), and mps_fire::result.

 {
 
   std::vector<AlignableDetOrUnitPtr> result;
   result.reserve(hitVec.size());
   
   for (TransientTrackingRecHit::ConstRecHitContainer::const_iterator it
          = hitVec.begin(), iEnd = hitVec.end(); it != iEnd; ++it) {
     result.push_back(this->alignableFromDetId((*it)->geographicalId()));
   }
 
   return result;
 }

bool AlignableNavigator::detAndSubdetInMap ( const DetId & detid ) const

Given a DetId, returns true if DetIds with this detector and subdetector id are in the map (not necessarily the exact DetId)

Definition at line 207 of file AlignableNavigator.cc.

References DetId::det(), mps_fire::i, DetId::subdetId(), and theDetAndSubdet.

Referenced by AlignmentMonitorTemplate::event(), AlignmentMonitorMuonResiduals::event(), recursiveGetId(), MuonMillepedeAlgorithm::run(), and size().

 {
    int det = detid.det();
    int subdet = detid.subdetId();
    for (std::vector<std::pair<int, int> >::const_iterator i = theDetAndSubdet.begin();  i != theDetAndSubdet.end();  ++i) {
       if (det == i->first  &&  subdet == i->second) return true;
    }
    return false;
 }

unsigned int AlignableNavigator::recursiveGetId ( Alignable * alignable )

private

Add recursively DetId-AlignableDetOrUnitPtr pairs to map. Returns number of Alignables with DetId!=0 which are (illegaly) neither AlignableDet nor AlignableDetUnit and are thus not added to the map.

Definition at line 109 of file AlignableNavigator.cc.

References align::AlignableDet, Alignable::alignableObjectId(), AlCaHLTBitMon_QueryRunRegistry::comp, Alignable::components(), detAndSubdetInMap(), Alignable::geomDetId(), theDetAndSubdet, and theMap.

Referenced by AlignableNavigator(), and size().

 {
   // Recursive method to get the detIds of an alignable and its childs
   // and add them to the map.
   // Returns number of Alignables with DetId which are neither AlignableDet
   // nor AlignableDetUnit and are thus not added to the map.
 
   if (!alignable) return 0;
 
   unsigned int nProblem = 0;
   const DetId detId(alignable->geomDetId());
   if (detId.rawId()) {
     
     AlignableDet *aliDet;
     AlignableDetUnit *aliDetUnit;
     AlignableBeamSpot *aliBeamSpot;
 
     if ((aliDet = dynamic_cast<AlignableDet*>(alignable))) {
       
       theMap.insert( PairType( detId, aliDet ) );
     
     } else if ((aliDetUnit = dynamic_cast<AlignableDetUnit*>(alignable))) {
 
       theMap.insert( PairType( detId, aliDetUnit ) );
 
     } else if ((aliBeamSpot = dynamic_cast<AlignableBeamSpot*>(alignable))) {
 
       theMap.insert( PairType( detId, aliBeamSpot ) );
 
     } else {
       
       nProblem = 1; // equivalent to '++nProblem;' which could confuse to be ina loop...
       // Cannot be an exception since it happens (illegaly) in Muon DT hierarchy:
       //         throw cms::Exception("BadLogic") 
       //           << "[AlignableNavigator::recursiveGetId] Alignable with DetId " << detId.rawId() 
       //           << " neither AlignableDet nor AlignableDetUnit";
     }
     
     if (!nProblem && !this->detAndSubdetInMap(detId)) {
       theDetAndSubdet.push_back(std::pair<int, int>( detId.det(), detId.subdetId() ));
     }
   }
   const auto& comp = alignable->components();
   if ( alignable->alignableObjectId() != align::AlignableDet
        || comp.size() > 1 ) { // Non-glued AlignableDets contain themselves
     for (const auto& it: comp) {
       nProblem += this->recursiveGetId(it);
     }
   }
   return nProblem;
 }