AlignableNavigator Class Reference

#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=nullptr)
	Constructor from one or two Alignables. More...
	AlignableNavigator (AlignableExtras *extras, Alignable *tracker, Alignable *muon=nullptr)
	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 25 of file AlignableNavigator.h.

Member Typedef Documentation

MapType

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

Definition at line 36 of file AlignableNavigator.h.

PairType

typedef MapType::value_type AlignableNavigator::PairType

Definition at line 37 of file AlignableNavigator.h.

Constructor & Destructor Documentation

AlignableNavigator() [1/3]

AlignableNavigator::AlignableNavigator ( Alignable *  tracker, Alignable *  muon = nullptr  )

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() [2/3]

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

explicit

Constructor from one or two Alignables.

Definition at line 38 of file AlignableNavigator.cc.

References AlignableExtras::components(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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() [3/3]

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

explicit

Constructor from list of Alignbable.

Definition at line 64 of file AlignableNavigator.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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

alignableDetOrUnits()

std::vector< AlignableDetOrUnitPtr > AlignableNavigator::alignableDetOrUnits

(

)

return all AlignableDetOrUnitPtrs

Definition at line 181 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;
}

alignableFromDetId()

AlignableDetOrUnitPtr AlignableNavigator::alignableFromDetId

(

const DetId &

detid

)

Returns AlignableDetOrUnitPtr corresponding to given DetId.

Definition at line 91 of file AlignableNavigator.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid, Exception, position, and theMap.

Referenced by alignableFromGeomDet(), alignablesFromHits(), SurveyInputCSCfromPins::analyze(), MuonGeometrySVGTemplate::analyze(), CSCAlignmentCorrections::applyAlignment(), AlignmentMonitorTemplate::event(), AlignmentMonitorGeneric::event(), ApeSettingAlgorithm::initialize(), MuonAlignment::moveAlignableGlobalCoord(), MuonAlignment::moveAlignableLocalCoord(), 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);
}

alignableFromGeomDet()

AlignableDetOrUnitPtr AlignableNavigator::alignableFromGeomDet

(

const GeomDet *

geomDet

)

Returns AlignableDetOrUnitPtr corresponding to given GeomDet.

Definition at line 86 of file AlignableNavigator.cc.

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

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

alignablesFromHits() [1/2]

std::vector< AlignableDetOrUnitPtr > AlignableNavigator::alignablesFromHits

(

const std::vector< const TransientTrackingRecHit *> &

hitvec

)

Returns vector AlignableDetOrUnitPtr for given vector of Hits.

Definition at line 152 of file AlignableNavigator.cc.

References alignableFromDetId(), and mps_fire::result.

                                                             {
  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;
}

alignablesFromHits() [2/2]

std::vector< AlignableDetOrUnitPtr > AlignableNavigator::alignablesFromHits

(

const TransientTrackingRecHit::ConstRecHitContainer &

hitVec

)

Returns vector of AlignableDetOrUnitPtr for given vector of Hits.

Definition at line 166 of file AlignableNavigator.cc.

References alignableFromDetId(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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;
}

detAndSubdetInMap()

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 193 of file AlignableNavigator.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid, mps_fire::i, theDetAndSubdet, and trackerHitRTTI::vector.

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

                                                                   {
  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;
}

recursiveGetId()

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 104 of file AlignableNavigator.cc.

References align::AlignableDet, Alignable::alignableObjectId(), AlCaHLTBitMon_QueryRunRegistry::comp, Alignable::components(), detAndSubdetInMap(), hcalRecHitTable_cff::detId, Alignable::geomDetId(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, theDetAndSubdet, and theMap.

Referenced by AlignableNavigator().

                                                                    {
  // 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;
}

size()

int AlignableNavigator::size ( void  )

inline

Returns number of elements in map.

Definition at line 55 of file AlignableNavigator.h.

References theMap.

Referenced by ntupleDataFormat._Collection::__iter__(), and ntupleDataFormat._Collection::__len__().

{ return theMap.size(); }

Member Data Documentation

theDetAndSubdet

std::vector<std::pair<int, int> > AlignableNavigator::theDetAndSubdet

private

Definition at line 67 of file AlignableNavigator.h.

Referenced by detAndSubdetInMap(), and recursiveGetId().

theMap

MapType AlignableNavigator::theMap

private

Definition at line 66 of file AlignableNavigator.h.

Referenced by alignableDetOrUnits(), alignableFromDetId(), AlignableNavigator(), recursiveGetId(), and size().