#include <Phase2EndcapRing.h>

Inheritance diagram for Phase2EndcapRing:

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))

	Phase2EndcapRing (std::vector< const GeomDet * > &innerDets, std::vector< const GeomDet * > &outerDets, const std::vector< const GeomDet * > &innerDetBrothers=std::vector< const GeomDet * >(), const std::vector< const GeomDet * > &outerDetBrothers=std::vector< const GeomDet * >())

virtual const BoundDisk &	specificSurface () const

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

	~Phase2EndcapRing () 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, std::vector< DetGroup > &brotherresult) 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, std::vector< DetGroup > &brotherresult, bool checkClosest) const __attribute__((hot))

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

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

Private Attributes
BinFinderType	theBackBinFinder

std::vector< const GeomDet * >	theBackDetBrothers

std::vector< const GeomDet * >	theBackDets

ReferenceCountingPointer < BoundDisk >	theBackDisk

std::vector< const GeomDet * >	theDets

ReferenceCountingPointer < BoundDisk >	theDisk

BinFinderType	theFrontBinFinder

std::vector< const GeomDet * >	theFrontDetBrothers

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 Phase2EndcapRing.h.

Member Typedef Documentation

typedef PeriodicBinFinderInPhi<float> Phase2EndcapRing::BinFinderType

private

Definition at line 79 of file Phase2EndcapRing.h.

Constructor & Destructor Documentation

Phase2EndcapRing::Phase2EndcapRing	(	std::vector< const GeomDet * > &	innerDets,
		std::vector< const GeomDet * > &	outerDets,
		const std::vector< const GeomDet * > &	innerDetBrothers = `std::vector<const GeomDet*>()`,
		const std::vector< const GeomDet * > &	outerDetBrothers = `std::vector<const GeomDet*>()`
	)

Definition at line 29 of file Phase2EndcapRing.cc.

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

     : GeometricSearchDet(true),
       theFrontDets(innerDets.begin(), innerDets.end()),
       theBackDets(outerDets.begin(), outerDets.end()),
       theFrontDetBrothers(innerDetBrothers.begin(), innerDetBrothers.end()),
       theBackDetBrothers(outerDetBrothers.begin(), outerDetBrothers.end()) {
   theDets.assign(theFrontDets.begin(), theFrontDets.end());
   theDets.insert(theDets.end(), theBackDets.begin(), theBackDets.end());
   theDets.insert(theDets.end(), theFrontDetBrothers.begin(), theFrontDetBrothers.end());
   theDets.insert(theDets.end(), theBackDetBrothers.begin(), theBackDetBrothers.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());
 
 #ifdef EDM_ML_DEBUG
   LogDebug("TkDetLayers") << "DEBUG INFO for Phase2EndcapRing";
   for (vector<const GeomDet*>::const_iterator it = theFrontDets.begin(); it != theFrontDets.end(); it++) {
     LogDebug("TkDetLayers") << "frontDet detId,phi,z,r: " << (*it)->geographicalId().rawId() << " , "
                             << (*it)->surface().position().phi() << " , " << (*it)->surface().position().z() << " , "
                             << (*it)->surface().position().perp();
   }
 
   if (!theFrontDetBrothers.empty()) {
     for (vector<const GeomDet*>::const_iterator it = theFrontDetBrothers.begin(); it != theFrontDetBrothers.end();
          it++) {
       LogDebug("TkDetLayers") << "frontDet brothers detId,phi,z,r: " << (*it)->geographicalId().rawId() << " , "
                               << (*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 detId,phi,z,r: " << (*it)->geographicalId().rawId() << " , "
                             << (*it)->surface().position().phi() << " , " << (*it)->surface().position().z() << " , "
                             << (*it)->surface().position().perp();
   }
 
   if (!theBackDetBrothers.empty()) {
     for (vector<const GeomDet*>::const_iterator it = theBackDetBrothers.begin(); it != theBackDetBrothers.end(); it++) {
       LogDebug("TkDetLayers") << "backDet brothers detId,phi,z,r: " << (*it)->geographicalId().rawId() << " , "
                               << (*it)->surface().position().phi() << " , " << (*it)->surface().position().z() << " , "
                               << (*it)->surface().position().perp();
     }
   }
 #endif
 }

Phase2EndcapRing::~Phase2EndcapRing ( )

override

Definition at line 88 of file Phase2EndcapRing.cc.

88 {}

Member Function Documentation

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

private

Definition at line 191 of file Phase2EndcapRing.cc.

References CompatibleDetToGroupAdder::add(), SubLayerCrossing::closestDetIndex(), subLayer(), subLayerBrothers(), SubLayerCrossing::subLayerIndex(), theBackDetBrothers, and theFrontDetBrothers.

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);
   if (theFrontDetBrothers.empty() && theBackDetBrothers.empty())
     return firstgroup;
   // it assumes that the closestDetIndex is ok also for the brother detectors: the crossing is NOT recomputed
   const vector<const GeomDet*>& subBrothers(subLayerBrothers(crossing.subLayerIndex()));
   const GeomDet* detBrother(subBrothers[crossing.closestDetIndex()]);
   bool brothergroup = CompatibleDetToGroupAdder::add(*detBrother, tsos, prop, est, brotherresult);
   return firstgroup || brothergroup;
 }

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

inlineoverridevirtual

Implements GeometricSearchDet.

Definition at line 24 of file Phase2EndcapRing.h.

References theDets.

24 { return theDets; }

Phase2EndcapRing::theDets

std::vector< const GeomDet * > theDets

Definition: Phase2EndcapRing.h:69

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

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

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

overridevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 90 of file Phase2EndcapRing.cc.

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

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

private

Definition at line 155 of file Phase2EndcapRing.cc.

References funct::abs(), PeriodicBinFinderInPhi< T >::binIndex(), srCondWrite_cfg::deltaPhi, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), rho, theBackBinFinder, theBackDets, theBackDisk, theFrontBinFinder, theFrontDets, theFrontDisk, 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(*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);
   }
 }

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

overridevirtual

Reimplemented from GeometricSearchDet.

Definition at line 101 of file Phase2EndcapRing.cc.

References addClosest(), SubLayerCrossings::closest(), SubLayerCrossings::closestIndex(), computeCrossings(), tkDetUtil::computeWindowSize(), LayerCrossingSide::endcapSide(), SubLayerCrossings::isValid(), LogTrace, eostools::move(), DetGroupMerger::orderAndMergeTwoLevels(), SubLayerCrossings::other(), Propagator::propagationDirection(), searchNeighbors(), theBackDetBrothers, and theFrontDetBrothers.

                                                                                  {
   SubLayerCrossings crossings;
   crossings = computeCrossings(tsos, prop.propagationDirection());
   if (!crossings.isValid())
     return;
 
   std::vector<DetGroup> closestResult;
   std::vector<DetGroup> closestBrotherResult;
   addClosest(tsos, prop, est, crossings.closest(), closestResult, closestBrotherResult);
   if (closestResult.empty())
     return;
 
   DetGroupElement closestGel(closestResult.front().front());
   int crossingSide = LayerCrossingSide().endcapSide(closestGel.trajectoryState(), prop);
   float phiWindow = tkDetUtil::computeWindowSize(closestGel.det(), closestGel.trajectoryState(), est);
   searchNeighbors(tsos, prop, est, crossings.closest(), phiWindow, closestResult, closestBrotherResult, false);
 
   vector<DetGroup> closestCompleteResult;
   DetGroupMerger::orderAndMergeTwoLevels(
       std::move(closestResult), std::move(closestBrotherResult), closestCompleteResult, 0, crossingSide);
 
   vector<DetGroup> nextResult;
   vector<DetGroup> nextBrotherResult;
   searchNeighbors(tsos, prop, est, crossings.other(), phiWindow, nextResult, nextBrotherResult, true);
 
   vector<DetGroup> nextCompleteResult;
   DetGroupMerger::orderAndMergeTwoLevels(
       std::move(nextResult), std::move(nextBrotherResult), nextCompleteResult, 0, crossingSide);
 
   DetGroupMerger::orderAndMergeTwoLevels(
       std::move(closestCompleteResult), std::move(nextCompleteResult), result, crossings.closestIndex(), crossingSide);
 
   //due to propagator problems, when we add single pt sub modules, we should order them in z (endcap)
   if (!theFrontDetBrothers.empty() && !theBackDetBrothers.empty())
     sort(result.begin(), result.end(), DetGroupElementZLess());
 
 #ifdef EDM_ML_DEBUG
   LogTrace("TkDetLayers") << "Number of groups : " << result.size() << std::endl;
   for (auto& grp : result) {
     if (grp.empty())
       continue;
     LogTrace("TkDetLayers") << "New group in Phase2EndcapRing made by : " << std::endl;
     for (auto const& det : grp) {
       LogTrace("TkDetLayers") << " geom det at r: " << det.det()->position().perp()
                               << " id:" << det.det()->geographicalId().rawId()
                               << " tsos at:" << det.trajectoryState().globalPosition() << std::endl;
     }
   }
 #endif
 }

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

private

Definition at line 209 of file Phase2EndcapRing.cc.

References PVValHelper::add(), PV3DBase< T, PVType, FrameType >::barePhi(), PeriodicBinFinderInPhi< T >::binIndex(), SubLayerCrossing::closestDetIndex(), tkDetUtil::overlapInPhi(), Geom::phiLess(), SubLayerCrossing::position(), subLayer(), subLayerBrothers(), SubLayerCrossing::subLayerIndex(), theBackBinFinder, theBackDetBrothers, theFrontBinFinder, and theFrontDetBrothers.

Referenced by groupedCompatibleDetsV().

                                                                 {
   const GlobalPoint& gCrossingPos = crossing.position();
 
   const vector<const GeomDet*>& sLayer(subLayer(crossing.subLayerIndex()));
   // It assumes that what is ok for the front modules in the pt modules is ok also for the back module
   const vector<const GeomDet*>& sBrotherLayer(subLayerBrothers(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;
     if (theFrontDetBrothers.empty() && theBackDetBrothers.empty())
       break;
     // If the two above checks are passed also the brother module will be added with no further checks
     const GeomDet& neighborBrotherDet = *sBrotherLayer[binFinder.binIndex(idet)];
     Adder::add(neighborBrotherDet, tsos, prop, est, brotherresult);
     // 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;
     if (theFrontDetBrothers.empty() && theBackDetBrothers.empty())
       break;
     // If the two above checks are passed also the brother module will be added with no further checks
     const GeomDet& neighborBrotherDet = *sBrotherLayer[binFinder.binIndex(idet)];
     Adder::add(neighborBrotherDet, tsos, prop, est, brotherresult);
     // maybe also add shallow crossing angle test here???
   }
 }

virtual const BoundDisk& Phase2EndcapRing::specificSurface ( ) const

inlinevirtual

Definition at line 38 of file Phase2EndcapRing.h.

References theDisk.

38 { return *theDisk; }

Phase2EndcapRing::theDisk

ReferenceCountingPointer< BoundDisk > theDisk

Definition: Phase2EndcapRing.h:75

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

inlineprivate

Definition at line 62 of file Phase2EndcapRing.h.

References theBackDets, and theFrontDets.

Referenced by addClosest(), and searchNeighbors().

62 { return (ind == 0 ? theFrontDets : theBackDets); }

Phase2EndcapRing::theFrontDets

std::vector< const GeomDet * > theFrontDets

Definition: Phase2EndcapRing.h:70

Phase2EndcapRing::theBackDets

std::vector< const GeomDet * > theBackDets

Definition: Phase2EndcapRing.h:71

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

inlineprivate

Definition at line 64 of file Phase2EndcapRing.h.

References theBackDetBrothers, and theFrontDetBrothers.

Referenced by addClosest(), and searchNeighbors().

                                                                    {
     return (ind == 0 ? theFrontDetBrothers : theBackDetBrothers);
   }

const BoundSurface& Phase2EndcapRing::surface ( ) const

inlineoverridevirtual

The surface of the GeometricSearchDet.

Implements GeometricSearchDet.

Definition at line 22 of file Phase2EndcapRing.h.

References theDisk.

22 { return *theDisk; }

Phase2EndcapRing::theDisk

ReferenceCountingPointer< BoundDisk > theDisk

Definition: Phase2EndcapRing.h:75

Member Data Documentation

BinFinderType Phase2EndcapRing::theBackBinFinder

private

Definition at line 82 of file Phase2EndcapRing.h.

Referenced by computeCrossings(), Phase2EndcapRing(), and searchNeighbors().

std::vector<const GeomDet*> Phase2EndcapRing::theBackDetBrothers

private

Definition at line 73 of file Phase2EndcapRing.h.

Referenced by addClosest(), groupedCompatibleDetsV(), Phase2EndcapRing(), searchNeighbors(), and subLayerBrothers().

std::vector<const GeomDet*> Phase2EndcapRing::theBackDets

private

Definition at line 71 of file Phase2EndcapRing.h.

Referenced by computeCrossings(), Phase2EndcapRing(), and subLayer().

ReferenceCountingPointer<BoundDisk> Phase2EndcapRing::theBackDisk

private

Definition at line 77 of file Phase2EndcapRing.h.

Referenced by computeCrossings(), and Phase2EndcapRing().

std::vector<const GeomDet*> Phase2EndcapRing::theDets

private

Definition at line 69 of file Phase2EndcapRing.h.

Referenced by basicComponents(), and Phase2EndcapRing().

ReferenceCountingPointer<BoundDisk> Phase2EndcapRing::theDisk

private

Definition at line 75 of file Phase2EndcapRing.h.

Referenced by Phase2EndcapRing(), specificSurface(), and surface().

BinFinderType Phase2EndcapRing::theFrontBinFinder

private

Definition at line 81 of file Phase2EndcapRing.h.

Referenced by computeCrossings(), Phase2EndcapRing(), and searchNeighbors().

std::vector<const GeomDet*> Phase2EndcapRing::theFrontDetBrothers

private

Definition at line 72 of file Phase2EndcapRing.h.

Referenced by addClosest(), groupedCompatibleDetsV(), Phase2EndcapRing(), searchNeighbors(), and subLayerBrothers().

std::vector<const GeomDet*> Phase2EndcapRing::theFrontDets

private

Definition at line 70 of file Phase2EndcapRing.h.

Referenced by computeCrossings(), Phase2EndcapRing(), and subLayer().

ReferenceCountingPointer<BoundDisk> Phase2EndcapRing::theFrontDisk

private

Definition at line 76 of file Phase2EndcapRing.h.

Referenced by computeCrossings(), and Phase2EndcapRing().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation