#include <MkFitGeometry.h>

Public Member Functions
const std::vector< const DetLayer * > &	detLayers () const

mkfit::LayerNumberConverter const &	layerNumberConverter () const

	MkFitGeometry (const TrackerGeometry &geom, const GeometricSearchTracker &tracker, const TrackerTopology &ttopo, std::unique_ptr< mkfit::TrackerInfo > trackerInfo)

int	mkFitLayerNumber (DetId detId) const

const TrackerTopology *	topology () const

mkfit::TrackerInfo const &	trackerInfo () const

unsigned int	uniqueIdInLayer (int layer, unsigned int detId) const

	~MkFitGeometry ()

Private Attributes
std::vector < std::unordered_map< unsigned int, unsigned int > >	detIdToShortId_

std::vector< const DetLayer * >	dets_

std::unique_ptr < mkfit::LayerNumberConverter >	lnc_

std::unique_ptr < mkfit::TrackerInfo >	trackerInfo_

const TrackerTopology *	ttopo_

Detailed Description

Collection of geometry-related objects for mkFit

Definition at line 23 of file MkFitGeometry.h.

Constructor & Destructor Documentation

MkFitGeometry::MkFitGeometry	(	const TrackerGeometry &	geom,
		const GeometricSearchTracker &	tracker,
		const TrackerTopology &	ttopo,
		std::unique_ptr< mkfit::TrackerInfo >	trackerInfo
	)

explicit

Definition at line 17 of file MkFitGeometry.cc.

References GeometricSearchTracker::allLayers(), cms::cuda::assert(), AlCaHLTBitMon_QueryRunRegistry::comp, TrackerGeometry::detIds(), detIdToShortId_, dets_, Exception, phase1PixelTopology::layer, TrackerTopology::layer(), lnc_, mkFitLayerNumber(), eostools::move(), TrackerGeometry::numberOfLayers(), or, mkfit::phase1, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, trackerInfo_, ttopo_, and UNLIKELY.

     : ttopo_(&ttopo),
       lnc_{std::make_unique<mkfit::LayerNumberConverter>(mkfit::TkLayout::phase1)},
       trackerInfo_(std::move(trackerInfo)) {
   if (geom.numberOfLayers(PixelSubdetector::PixelBarrel) != 4 ||
       geom.numberOfLayers(PixelSubdetector::PixelEndcap) != 3) {
     throw cms::Exception("Assert") << "For now this code works only with phase1 tracker, you have something else";
   }
 
   // Create DetLayer structure
   dets_.resize(lnc_->nLayers(), nullptr);
   auto setDet = [this](const int subdet, const int layer, const int isStereo, const DetId& detId, const DetLayer* lay) {
     const int index = lnc_->convertLayerNumber(subdet, layer, false, isStereo, isPlusSide(*ttopo_, detId));
     if (index < 0 or static_cast<unsigned>(index) >= dets_.size()) {
       throw cms::Exception("LogicError") << "Invalid mkFit layer index " << index << " for DetId " << detId.rawId()
                                          << " subdet " << subdet << " layer " << layer << " isStereo " << isStereo;
     }
     dets_[index] = lay;
   };
   constexpr int monoLayer = 0;
   constexpr int stereoLayer = 1;
   for (const DetLayer* lay : tracker.allLayers()) {
     const auto& comp = lay->basicComponents();
     if (UNLIKELY(comp.empty())) {
       throw cms::Exception("LogicError") << "Got a tracker layer (subdet " << lay->subDetector()
                                          << ") with empty basicComponents.";
     }
     // First component is enough for layer and side information
     const auto& detId = comp.front()->geographicalId();
     const auto subdet = detId.subdetId();
     const auto layer = ttopo.layer(detId);
 
     // TODO: mono/stereo structure is still hardcoded for phase0/1 strip tracker
     setDet(subdet, layer, monoLayer, detId, lay);
     if (((subdet == StripSubdetector::TIB or subdet == StripSubdetector::TOB) and (layer == 1 or layer == 2)) or
         subdet == StripSubdetector::TID or subdet == StripSubdetector::TEC) {
       setDet(subdet, layer, stereoLayer, detId, lay);
     }
   }
 
   // Create "short id" aka "unique id within layer"
   detIdToShortId_.resize(lnc_->nLayers());
   for (const auto& detId : geom.detIds()) {
     const auto ilay = mkFitLayerNumber(detId);
     auto& map = detIdToShortId_[ilay];
     const unsigned int ind = map.size();
     // Make sure the short id fits in the 12 bits...
     assert(ind < (int)1 << 11);
     map[detId.rawId()] = ind;
   }
 }

MkFitGeometry::~MkFitGeometry ( )

Definition at line 74 of file MkFitGeometry.cc.

74 {}

Member Function Documentation

const std::vector<const DetLayer*>& MkFitGeometry::detLayers ( ) const

inline

Definition at line 34 of file MkFitGeometry.h.

References dets_.

34 { return dets_; }

MkFitGeometry::dets_

std::vector< const DetLayer * > dets_

Definition: MkFitGeometry.h:42

mkfit::LayerNumberConverter const& MkFitGeometry::layerNumberConverter ( ) const

inline

Definition at line 32 of file MkFitGeometry.h.

References lnc_.

Referenced by MkFitSeedConverter::convertSeeds().

32 { return *lnc_; }

MkFitGeometry::lnc_

std::unique_ptr< mkfit::LayerNumberConverter > lnc_

Definition: MkFitGeometry.h:40

int MkFitGeometry::mkFitLayerNumber ( DetId detId ) const

Definition at line 76 of file MkFitGeometry.cc.

References TrackerTopology::isStereo(), TrackerTopology::layer(), lnc_, DetId::subdetId(), and ttopo_.

Referenced by mkfit::convertHits(), MkFitEventOfHitsProducer::fill(), and MkFitGeometry().

                                                      {
   return lnc_->convertLayerNumber(
       detId.subdetId(), ttopo_->layer(detId), false, ttopo_->isStereo(detId), isPlusSide(*ttopo_, detId));
 }