#include <TIDRing.h>

Inheritance diagram for TIDRing:

Public Member Functions
const std::vector< const GeomDet * > &	basicComponents () const override

std::pair< bool, TrajectoryStateOnSurface >	compatible (const TrajectoryStateOnSurface &, const Propagator &, const MeasurementEstimator &) const override

const std::vector< const GeometricSearchDet * > &	components () const override __attribute__((cold))
	Returns basic components, if any. More...

void	groupedCompatibleDetsV (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetGroup > &result) const override __attribute__((hot))

virtual const BoundDisk &	specificSurface () const

const BoundSurface &	surface () const override
	The surface of the GeometricSearchDet. More...

	TIDRing (std::vector< const GeomDet * > &innerDets, std::vector< const GeomDet * > &outerDets)

	~TIDRing () override

Public Member Functions inherited from GeometricSearchDet
virtual std::vector< DetWithState >	compatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const

virtual void	compatibleDetsV (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est, std::vector< DetWithState > &result) const

	GeometricSearchDet (bool doHaveGroups)

virtual std::vector< DetGroup >	groupedCompatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const

bool	hasGroups () const

virtual const Surface::PositionType &	position () const
	Returns position of the surface. More...

virtual	~GeometricSearchDet ()

Private Types
typedef PeriodicBinFinderInPhi< float >	BinFinderType

Private Member Functions
bool	addClosest (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, const SubLayerCrossing &crossing, std::vector< DetGroup > &result) const __attribute__((hot))

SubLayerCrossings	computeCrossings (const TrajectoryStateOnSurface &tsos, PropagationDirection propDir) const __attribute__((hot))

void	searchNeighbors (const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est, const SubLayerCrossing &crossing, float window, std::vector< DetGroup > &result, bool checkClosest) const __attribute__((hot))

const std::vector< const GeomDet * > &	subLayer (int ind) const

Private Attributes
BinFinderType	theBackBinFinder

std::vector< const GeomDet * >	theBackDets

ReferenceCountingPointer< BoundDisk >	theBackDisk

std::vector< const GeomDet * >	theDets

ReferenceCountingPointer< BoundDisk >	theDisk

BinFinderType	theFrontBinFinder

std::vector< const GeomDet * >	theFrontDets

ReferenceCountingPointer< BoundDisk >	theFrontDisk

Additional Inherited Members
Public Types inherited from GeometricSearchDet
typedef std::pair< const GeomDet *, TrajectoryStateOnSurface >	DetWithState

typedef BoundSurface::PositionType	PositionType

typedef BoundSurface::RotationType	RotationType

typedef TrajectoryStateOnSurface	TrajectoryState

Protected Attributes inherited from GeometricSearchDet
bool	haveGroups

GeomDetCompatibilityChecker	theCompatibilityChecker

Detailed Description

A concrete implementation for TID rings

Definition at line 13 of file TIDRing.h.

Member Typedef Documentation

◆ BinFinderType

typedef PeriodicBinFinderInPhi<float> TIDRing::BinFinderType

private

Definition at line 68 of file TIDRing.h.

Constructor & Destructor Documentation

◆ TIDRing()

TIDRing::TIDRing	(	std::vector< const GeomDet * > &	innerDets,
		std::vector< const GeomDet * > &	outerDets
	)

Definition at line 24 of file TIDRing.cc.

     : GeometricSearchDet(true),
       theFrontDets(innerDets.begin(), innerDets.end()),
       theBackDets(outerDets.begin(), outerDets.end()) {
   theDets.assign(theFrontDets.begin(), theFrontDets.end());
   theDets.insert(theDets.end(), theBackDets.begin(), theBackDets.end());
  
   // the dets should be already phi-ordered. TO BE CHECKED
   //sort( theFrontDets.begin(), theFrontDets.end(), DetLessPhi() );
   //sort( theBackDets.begin(), theBackDets.end(), DetLessPhi() );
  
   theDisk = ForwardRingDiskBuilderFromDet()(theDets);
  
   theFrontDisk = ForwardRingDiskBuilderFromDet()(theFrontDets);
   theBackDisk = ForwardRingDiskBuilderFromDet()(theBackDets);
  
   theFrontBinFinder = BinFinderType(theFrontDets.front()->surface().position().phi(), theFrontDets.size());
   theBackBinFinder = BinFinderType(theBackDets.front()->surface().position().phi(), theBackDets.size());
  
   LogDebug("TkDetLayers") << "DEBUG INFO for TIDRing";
   for (vector<const GeomDet*>::const_iterator it = theFrontDets.begin(); it != theFrontDets.end(); it++) {
     LogDebug("TkDetLayers") << "frontDet phi,z,r: " << (*it)->surface().position().phi() << " , "
                             << (*it)->surface().position().z() << " , " << (*it)->surface().position().perp();
   }
  
   for (vector<const GeomDet*>::const_iterator it = theBackDets.begin(); it != theBackDets.end(); it++) {
     LogDebug("TkDetLayers") << "backDet phi,z,r: " << (*it)->surface().position().phi() << " , "
                             << (*it)->surface().position().z() << " , " << (*it)->surface().position().perp();
   }
 }

References LogDebug, theBackBinFinder, theBackDets, theBackDisk, theDets, theDisk, theFrontBinFinder, theFrontDets, and theFrontDisk.

◆ ~TIDRing()

TIDRing::~TIDRing ( )

override

Definition at line 55 of file TIDRing.cc.

55 {}

Member Function Documentation

◆ addClosest()

bool TIDRing::addClosest	(	const TrajectoryStateOnSurface &	tsos,
		const Propagator &	prop,
		const MeasurementEstimator &	est,
		const SubLayerCrossing &	crossing,
		std::vector< DetGroup > &	result
	)		const

private

Definition at line 132 of file TIDRing.cc.

                                                          {
   const vector<const GeomDet*>& sub(subLayer(crossing.subLayerIndex()));
   const GeomDet* det(sub[crossing.closestDetIndex()]);
   return CompatibleDetToGroupAdder::add(*det, tsos, prop, est, result);
 }

References CompatibleDetToGroupAdder::add(), SubLayerCrossing::closestDetIndex(), mps_fire::result, subLayer(), and SubLayerCrossing::subLayerIndex().

Referenced by groupedCompatibleDetsV().

◆ basicComponents()

const std::vector<const GeomDet*>& TIDRing::basicComponents ( ) const

inlineoverridevirtual

Implements GeometricSearchDet.

Definition at line 21 of file TIDRing.h.

21 { return theDets; }

References theDets.

◆ compatible()

pair< bool, TrajectoryStateOnSurface > TIDRing::compatible	(	const TrajectoryStateOnSurface &	ts,
		const Propagator &	,
		const MeasurementEstimator &
	)		const

overridevirtual

tests the geometrical compatibility of the Det with the predicted state. The FreeTrajectoryState argument is propagated to the Det surface using the Propagator argument. The resulting TrajectoryStateOnSurface is tested for compatibility with the surface bounds. If compatible, a std::pair< true, propagatedState> is returned. If the propagation fails, or if the state is not compatible, a std::pair< false, propagatedState> is returned.

Implements GeometricSearchDet.

Definition at line 61 of file TIDRing.cc.

                                                                                             {
   edm::LogError("TkDetLayers") << "temporary dummy implementation of TIDRing::compatible()!!";
   return pair<bool, TrajectoryStateOnSurface>();
 }

◆ components()

const vector< const GeometricSearchDet * > & TIDRing::components ( ) const

overridevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 57 of file TIDRing.cc.

                                                                    {
   throw DetLayerException("TIDRing doesn't have GeometricSearchDet components");
 }

◆ computeCrossings()

SubLayerCrossings TIDRing::computeCrossings	(	const TrajectoryStateOnSurface &	tsos,
		PropagationDirection	propDir
	)		const

private

Definition at line 95 of file TIDRing.cc.

                                                                                 {
   HelixPlaneCrossing::PositionType startPos(startingState.globalPosition());
   HelixPlaneCrossing::DirectionType startDir(startingState.globalMomentum());
  
   auto rho = startingState.transverseCurvature();
  
   HelixForwardPlaneCrossing crossing(startPos, startDir, rho, propDir);
  
   pair<bool, double> frontPath = crossing.pathLength(*theFrontDisk);
   if (!frontPath.first)
     return SubLayerCrossings();
  
   pair<bool, double> backPath = crossing.pathLength(*theBackDisk);
   if (!backPath.first)
     return SubLayerCrossings();
  
   GlobalPoint gFrontPoint(crossing.position(frontPath.second));
   GlobalPoint gBackPoint(crossing.position(backPath.second));
  
   int frontIndex = theFrontBinFinder.binIndex(gFrontPoint.barePhi());
   SubLayerCrossing frontSLC(0, frontIndex, gFrontPoint);
  
   int backIndex = theBackBinFinder.binIndex(gBackPoint.barePhi());
   SubLayerCrossing backSLC(1, backIndex, gBackPoint);
  
   // 0ss: frontDisk has index=0, backDisk has index=1
   float frontDist = std::abs(Geom::deltaPhi(gFrontPoint.barePhi(), theFrontDets[frontIndex]->surface().phi()));
   float backDist = std::abs(Geom::deltaPhi(gBackPoint.barePhi(), theBackDets[backIndex]->surface().phi()));
  
   if (frontDist < backDist) {
     return SubLayerCrossings(frontSLC, backSLC, 0);
   } else {
     return SubLayerCrossings(backSLC, frontSLC, 1);
   }
 }

References funct::abs(), PeriodicBinFinderInPhi< T >::binIndex(), SiPixelRawToDigiRegional_cfi::deltaPhi, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), HelixForwardPlaneCrossing::pathLength(), HelixForwardPlaneCrossing::position(), rho, theBackBinFinder, theBackDets, theBackDisk, theFrontBinFinder, theFrontDets, theFrontDisk, and TrajectoryStateOnSurface::transverseCurvature().

Referenced by groupedCompatibleDetsV().

◆ groupedCompatibleDetsV()

void TIDRing::groupedCompatibleDetsV	(	const TrajectoryStateOnSurface &	tsos,
		const Propagator &	prop,
		const MeasurementEstimator &	est,
		std::vector< DetGroup > &	result
	)		const

overridevirtual

Reimplemented from GeometricSearchDet.

Definition at line 68 of file TIDRing.cc.

                                                                         {
   SubLayerCrossings crossings;
   crossings = computeCrossings(tsos, prop.propagationDirection());
   if (!crossings.isValid())
     return;
  
   std::vector<DetGroup> closestResult;
   addClosest(tsos, prop, est, crossings.closest(), closestResult);
   if (closestResult.empty())
     return;
  
   DetGroupElement closestGel(closestResult.front().front());
   float phiWindow = tkDetUtil::computeWindowSize(closestGel.det(), closestGel.trajectoryState(), est);
   searchNeighbors(tsos, prop, est, crossings.closest(), phiWindow, closestResult, false);
  
   vector<DetGroup> nextResult;
   searchNeighbors(tsos, prop, est, crossings.other(), phiWindow, nextResult, true);
  
   int crossingSide = LayerCrossingSide().endcapSide(closestGel.trajectoryState(), prop);
   DetGroupMerger::orderAndMergeTwoLevels(
       std::move(closestResult), std::move(nextResult), result, crossings.closestIndex(), crossingSide);
 }

References addClosest(), SubLayerCrossings::closest(), SubLayerCrossings::closestIndex(), computeCrossings(), tkDetUtil::computeWindowSize(), LayerCrossingSide::endcapSide(), SubLayerCrossings::isValid(), eostools::move(), DetGroupMerger::orderAndMergeTwoLevels(), SubLayerCrossings::other(), trackingPOGFilters_cfi::phiWindow, Propagator::propagationDirection(), mps_fire::result, and searchNeighbors().

Referenced by TIDLayer::groupedCompatibleDetsV().

◆ searchNeighbors()

void TIDRing::searchNeighbors	(	const TrajectoryStateOnSurface &	tsos,
		const Propagator &	prop,
		const MeasurementEstimator &	est,
		const SubLayerCrossing &	crossing,
		float	window,
		std::vector< DetGroup > &	result,
		bool	checkClosest
	)		const

private

Definition at line 142 of file TIDRing.cc.

                                                        {
   const GlobalPoint& gCrossingPos = crossing.position();
  
   const vector<const GeomDet*>& sLayer(subLayer(crossing.subLayerIndex()));
  
   int closestIndex = crossing.closestDetIndex();
   int negStartIndex = closestIndex - 1;
   int posStartIndex = closestIndex + 1;
  
   if (checkClosest) {  // must decide if the closest is on the neg or pos side
     if (Geom::phiLess(gCrossingPos.barePhi(), sLayer[closestIndex]->surface().phi())) {
       posStartIndex = closestIndex;
     } else {
       negStartIndex = closestIndex;
     }
   }
  
   const BinFinderType& binFinder = (crossing.subLayerIndex() == 0 ? theFrontBinFinder : theBackBinFinder);
  
   typedef CompatibleDetToGroupAdder Adder;
   int half = sLayer.size() / 2;  // to check if dets are called twice....
   for (int idet = negStartIndex; idet >= negStartIndex - half; idet--) {
     const GeomDet& neighborDet = *sLayer[binFinder.binIndex(idet)];
     if (!tkDetUtil::overlapInPhi(gCrossingPos, neighborDet, window))
       break;
     if (!Adder::add(neighborDet, tsos, prop, est, result))
       break;
     // maybe also add shallow crossing angle test here???
   }
   for (int idet = posStartIndex; idet < posStartIndex + half; idet++) {
     const GeomDet& neighborDet = *sLayer[binFinder.binIndex(idet)];
     if (!tkDetUtil::overlapInPhi(gCrossingPos, neighborDet, window))
       break;
     if (!Adder::add(neighborDet, tsos, prop, est, result))
       break;
     // maybe also add shallow crossing angle test here???
   }
 }