#include <Phase2EndcapSingleRing.h>

Inheritance diagram for Phase2EndcapSingleRing:

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
	Returns basic components, if any. More...

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

	Phase2EndcapSingleRing (std::vector< const GeomDet *> &allDets)

virtual const BoundDisk &	specificSurface () const

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

	~Phase2EndcapSingleRing () 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

SubLayerCrossing	computeCrossing (const TrajectoryStateOnSurface &tsos, PropagationDirection propDir) const

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

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

Private Attributes
BinFinderType	theBinFinder

std::vector< const GeomDet * >	theDets

ReferenceCountingPointer< BoundDisk >	theDisk

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 Phase2 SubDisk rings

Definition at line 13 of file Phase2EndcapSingleRing.h.

Member Typedef Documentation

◆ BinFinderType

typedef PeriodicBinFinderInPhi<float> Phase2EndcapSingleRing::BinFinderType

private

Definition at line 63 of file Phase2EndcapSingleRing.h.

Constructor & Destructor Documentation

◆ Phase2EndcapSingleRing()

Phase2EndcapSingleRing::Phase2EndcapSingleRing ( std::vector< const GeomDet *> & allDets )

Definition at line 24 of file Phase2EndcapSingleRing.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, LogDebug, theBinFinder, theDets, and theDisk.

     : GeometricSearchDet(true), theDets(allDets.begin(), allDets.end()) {
   theDisk = ForwardRingDiskBuilderFromDet()(theDets);
 
   theBinFinder = BinFinderType(theDets.front()->surface().position().phi(), theDets.size());
 
 #ifdef EDM_ML_DEBUG
   LogDebug("TkDetLayers") << "DEBUG INFO for Phase2EndcapSingleRing";
   for (vector<const GeomDet*>::const_iterator it = theDets.begin(); it != theDets.end(); it++) {
     LogDebug("TkDetLayers") << "Det detId,phi,z,r: " << (*it)->geographicalId().rawId() << " , "
                             << (*it)->surface().position().phi() << " , " << (*it)->surface().position().z() << " , "
                             << (*it)->surface().position().perp();
   }
 
 #endif
 }

◆ ~Phase2EndcapSingleRing()

Phase2EndcapSingleRing::~Phase2EndcapSingleRing ( )

overridedefault

Member Function Documentation

◆ addClosest()

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

private

Definition at line 100 of file Phase2EndcapSingleRing.cc.

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

Referenced by groupedCompatibleDetsV().

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

◆ basicComponents()

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

inlineoverridevirtual

Implements GeometricSearchDet.

Definition at line 21 of file Phase2EndcapSingleRing.h.

References theDets.

21 { return theDets; }

Phase2EndcapSingleRing::theDets

std::vector< const GeomDet * > theDets

Definition: Phase2EndcapSingleRing.h:59

◆ compatible()

pair< bool, TrajectoryStateOnSurface > Phase2EndcapSingleRing::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 47 of file Phase2EndcapSingleRing.cc.

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

◆ components()

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

overridevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 43 of file Phase2EndcapSingleRing.cc.

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

◆ computeCrossing()

SubLayerCrossing Phase2EndcapSingleRing::computeCrossing	(	const TrajectoryStateOnSurface &	tsos,
		PropagationDirection	propDir
	)		const

private

Definition at line 80 of file Phase2EndcapSingleRing.cc.

References PeriodicBinFinderInPhi< T >::binIndex(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), rho, theBinFinder, theDisk, and TrajectoryStateOnSurface::transverseCurvature().

Referenced by groupedCompatibleDetsV().

                                                                                              {
   auto rho = startingState.transverseCurvature();
 
   HelixPlaneCrossing::PositionType startPos(startingState.globalPosition());
   HelixPlaneCrossing::DirectionType startDir(startingState.globalMomentum());
   HelixForwardPlaneCrossing crossing(startPos, startDir, rho, propDir);
 
   pair<bool, double> frontPath = crossing.pathLength(*theDisk);
   if (!frontPath.first)
     return SubLayerCrossing();
 
   GlobalPoint gFrontPoint(crossing.position(frontPath.second));  //There is only one path
 
   int frontIndex = theBinFinder.binIndex(gFrontPoint.barePhi());
   SubLayerCrossing frontSLC(0, frontIndex, gFrontPoint);
 
   return frontSLC;
 }

◆ groupedCompatibleDetsV()

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

overridevirtual

Reimplemented from GeometricSearchDet.

Definition at line 54 of file Phase2EndcapSingleRing.cc.

References addClosest(), DetGroupMerger::addSameLevel(), computeCrossing(), tkDetUtil::computeWindowSize(), SubLayerCrossing::isValid(), eostools::move(), trackingPOGFilters_cfi::phiWindow, Propagator::propagationDirection(), mps_fire::result, and searchNeighbors().

                                                                                        {
   SubLayerCrossing crossing;
 
   crossing = computeCrossing(tsos, prop.propagationDirection());
 
   if (!crossing.isValid())
     return;
 
   std::vector<DetGroup> closestResult;
 
   addClosest(tsos, prop, est, crossing, closestResult);
   if (closestResult.empty())
     return;
 
   DetGroupElement closestGel(closestResult.front().front());
 
   float phiWindow = tkDetUtil::computeWindowSize(closestGel.det(), closestGel.trajectoryState(), est);
 
   searchNeighbors(tsos, prop, est, crossing, phiWindow, closestResult, false);
 
   DetGroupMerger::addSameLevel(std::move(closestResult), result);
 }

◆ searchNeighbors()

void Phase2EndcapSingleRing::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 114 of file Phase2EndcapSingleRing.cc.

References PVValHelper::add(), PV3DBase< T, PVType, FrameType >::barePhi(), PeriodicBinFinderInPhi< T >::binIndex(), SubLayerCrossing::closestDetIndex(), tkDetUtil::overlapInPhi(), Geom::phiLess(), SubLayerCrossing::position(), mps_fire::result, subLayer(), SubLayerCrossing::subLayerIndex(), theBinFinder, and svgfig::window().

Referenced by groupedCompatibleDetsV().

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