#include <MTDDetTray.h>

Inheritance diagram for MTDDetTray:

Public Member Functions
std::pair< bool, TrajectoryStateOnSurface >	compatible (const TrajectoryStateOnSurface &ts, const Propagator &prop, const MeasurementEstimator &est) const override

std::vector< DetWithState >	compatibleDets (const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const override

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

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

	MTDDetTray (std::vector< const GeomDet >::const_iterator first, std::vector< const GeomDet >::const_iterator last)
	Construct from iterators on GeomDet*. More...

	MTDDetTray (const std::vector< const GeomDet *> &dets)
	Construct from a std::vector of GeomDet*. More...

	~MTDDetTray () override
	Destructor. More...

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

	DetRodOneR (std::vector< const GeomDet >::const_iterator first, std::vector< const GeomDet >::const_iterator last)
	Construct from iterators on GeomDet*. More...

	DetRodOneR (const std::vector< const GeomDet *> &dets)
	Construct from a std::vector of GeomDet*. More...

	~DetRodOneR () override

Public Member Functions inherited from DetRod
	GeometricSearchDet (bool doHaveGroups)

virtual const Plane &	specificSurface () const final
	Return the rod surface as a Plane. More...

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

	~DetRod () override

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

	GeometricSearchDet (bool doHaveGroups)

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

bool	hasGroups () const

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

virtual	~GeometricSearchDet ()

Private Types
typedef GenericBinFinderInZ< float, GeomDet >	BinFinderType

Private Member Functions
void	init ()

Private Attributes
BinFinderType	theBinFinder

Additional Inherited Members
Public Types inherited from DetRodOneR
typedef std::vector< GeometricSearchDet * >	DetContainer

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

typedef BoundSurface::PositionType	PositionType

typedef BoundSurface::RotationType	RotationType

typedef TrajectoryStateOnSurface	TrajectoryState

Protected Member Functions inherited from DetRodOneR
bool	add (int idet, std::vector< DetWithState > &result, const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
	Query detector idet for compatible and add the output to result. More...

void	initialize ()

Protected Member Functions inherited from DetRod
void	setPlane (Plane *plane)
	Set the rod's plane. More...

Protected Attributes inherited from DetRodOneR
std::vector< const GeomDet * >	theDets

Protected Attributes inherited from GeometricSearchDet
bool	haveGroups

GeomDetCompatibilityChecker	theCompatibilityChecker

Detailed Description

A tray of aligned equal-sized non-overlapping detectors.
Designed for barrel timing layer.

Author: L. Gray - FNAL

Definition at line 17 of file MTDDetTray.h.

Member Typedef Documentation

◆ BinFinderType

typedef GenericBinFinderInZ<float, GeomDet> MTDDetTray::BinFinderType

private

Definition at line 45 of file MTDDetTray.h.

Constructor & Destructor Documentation

◆ MTDDetTray() [1/2]

MTDDetTray::MTDDetTray	(	std::vector< const GeomDet *>::const_iterator	first,
		std::vector< const GeomDet *>::const_iterator	last
	)

Construct from iterators on GeomDet*.

Definition at line 18 of file MTDDetTray.cc.

References init().

     : DetRodOneR(first, last) {
   init();
 }

◆ MTDDetTray() [2/2]

MTDDetTray::MTDDetTray ( const std::vector< const GeomDet *> & dets )

Construct from a std::vector of GeomDet*.

Definition at line 23 of file MTDDetTray.cc.

References init().

23 : DetRodOneR(vdets) { init(); }

MTDDetTray::init

void init()

Definition: MTDDetTray.cc:25

DetRodOneR::DetRodOneR

DetRodOneR(std::vector< const GeomDet *>::const_iterator first, std::vector< const GeomDet *>::const_iterator last)

Construct from iterators on GeomDet*.

Definition: DetRodOneR.cc:15

◆ ~MTDDetTray()

MTDDetTray::~MTDDetTray ( )

override

Destructor.

Definition at line 27 of file MTDDetTray.cc.

27 {}

Member Function Documentation

◆ compatible()

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

References MeasurementEstimator::estimate(), TrajectoryStateOnSurface::isValid(), Propagator::propagate(), and DetRod::specificSurface().

Referenced by compatibleDets().

                                                                                                    {
   TrajectoryStateOnSurface ms = prop.propagate(ts, specificSurface());
   return make_pair(ms.isValid() and est.estimate(ms, specificSurface()) != 0, ms);
 }

◆ compatibleDets()

vector< GeometricSearchDet::DetWithState > MTDDetTray::compatibleDets	(	const TrajectoryStateOnSurface &	startingState,
		const Propagator &	prop,
		const MeasurementEstimator &	est
	)		const

overridevirtual

Returns all Dets compatible with a trajectory state according to the estimator est. The startingState should be propagated to the surface of each compatible Det using the Propagator passed as an argument. The default implementation should be overridden in dets with specific surface types to avoid propagation to a generic Surface

Reimplemented from GeometricSearchDet.

Definition at line 43 of file MTDDetTray.cc.

References DetRodOneR::add(), DetRodOneR::basicComponents(), GenericBinFinderInZ< T, G >::binIndex(), pv::closest(), compatible(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::localError(), LogTrace, particleFlowClusterHGC_cfi::maxDistance, MeasurementEstimator::maximalLocalDisplacement(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GloballyPositioned< T >::position(), LocalTrajectoryError::positionError(), mps_fire::result, mathSSE::sqrt(), DetRod::surface(), theBinFinder, toLocal(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by groupedCompatibleDets().

                                                                                                            {
   LogTrace("MTDDetLayers") << "MTDDetTray::compatibleDets, Surface at R,phi: " << surface().position().perp() << ","
                            << surface().position().phi() << "     DetRod pos."
                            << " TS at R,phi: " << startingState.globalPosition().perp() << ","
                            << startingState.globalPosition().phi() << " TSOS dx/dy "
                            << std::sqrt(startingState.localError().positionError().xx()) << " "
                            << std::sqrt(startingState.localError().positionError().yy());
 
   vector<DetWithState> result;
 
   // Propagate and check that the result is within bounds
   pair<bool, TrajectoryStateOnSurface> compat = compatible(startingState, prop, est);
 
   if (!compat.first) {
     LogTrace("MTDDetLayers") << "    MTDDetTray::compatibleDets: not compatible"
                              << "    (should not have been selected!)";
     return result;
   }
 
   // Find the most probable destination component
   TrajectoryStateOnSurface& tsos = compat.second;
   GlobalPoint startPos = tsos.globalPosition();
   int closest = theBinFinder.binIndex(startPos.z());
   const vector<const GeomDet*> dets = basicComponents();
   LogTrace("MTDDetLayers") << "     MTDDetTray::compatibleDets, closest det: " << closest
                            << " pos: " << dets[closest]->surface().position() << " impact " << startPos;
 
   // Add this detector, if it is compatible
   // NOTE: add performs a null propagation
   add(closest, result, tsos, prop, est);
 
 #ifdef EDM_ML_DEBUG
   int nclosest = result.size();
   int nnextdet = 0;  // just DEBUG counters
 #endif
 
   // Try the neighbors on each side until no more compatible.
   // If closest is not compatible the next cannot be either
   if (!result.empty()) {
     const BoundPlane& closestPlane(dets[closest]->surface());
     MeasurementEstimator::Local2DVector maxDistance = est.maximalLocalDisplacement(result.front().second, closestPlane);
 
     // detHalfLen is assumed to be the same for all detectors.
     float detHalfLen = closestPlane.bounds().length() / 2.;
 
     for (unsigned int idet = closest + 1; idet < dets.size(); idet++) {
       LocalPoint nextPos(dets[idet]->toLocal(startPos));
       if (fabs(nextPos.y()) < detHalfLen + maxDistance.y()) {
 #ifdef EDM_ML_DEBUG
         LogTrace("MTDDetLayers") << "     negativeZ: det:" << idet << " pos " << nextPos.y() << " maxDistance "
                                  << maxDistance.y();
         nnextdet++;
 #endif
         if (!add(idet, result, tsos, prop, est))
           break;
       } else {
         break;
       }
     }
 
     for (int idet = closest - 1; idet >= 0; idet--) {
       LocalPoint nextPos(dets[idet]->toLocal(startPos));
       if (fabs(nextPos.y()) < detHalfLen + maxDistance.y()) {
 #ifdef EDM_ML_DEBUG
         LogTrace("MTDDetLayers") << "     positiveZ: det:" << idet << " pos " << nextPos.y() << " maxDistance "
                                  << maxDistance.y();
         nnextdet++;
 #endif
         if (!add(idet, result, tsos, prop, est))
           break;
       } else {
         break;
       }
     }
   }
 
 #ifdef EDM_ML_DEBUG
   LogTrace("MTDDetLayers") << "     MTDDetTray::compatibleDets, size: " << result.size() << " on closest: " << nclosest
                            << " # checked dets: " << nnextdet + 1;
   if (result.empty()) {
     LogTrace("MTDDetLayers") << "   ***Rod not compatible---should have been discarded before!!!";
   }
 #endif
   return result;
 }

◆ components()

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

overridevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 29 of file MTDDetTray.cc.

References mps_fire::result.

                                                                       {
   // FIXME dummy impl.
   edm::LogError("MTDDetLayers") << "temporary dummy implementation of MTDDetTray::components()!!";
   static const vector<const GeometricSearchDet*> result;
   return result;
 }

◆ groupedCompatibleDets()

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

overridevirtual

Similar to compatibleDets(), but the compatible Dets are grouped in one or more groups. Dets are put in the same group if they are mutually exclusive for track crossing, i.e. a reconstructible track cannot cross more than one Det from a group. Pathological tracks (spirals etc.) can of course violate this rule.
The DetGroups are sorted in the sequence of crossing by a track. In order to define the direction of crossing the Propagator used in this method should have a defined direction() : either "alongMomentum" or "oppositeToMomentum" but not "anyDirection".
The three signatures of this method differ by the input trajectory state arguments: the starting state can be a TrajectoryStateOnSurface or a FreeTrajectoryState, and the state on this CompositeDet may be already known or not. The last two arguments are as for the method compatibleDets().
First signature: The first argument is a TrajectoryStateOnSurface, usually not on the surface of this CompositeDet.

Reimplemented from GeometricSearchDet.

Definition at line 131 of file MTDDetTray.cc.

References compatibleDets(), LogTrace, and mps_fire::result.

                                                                                           {
   vector<GeometricSearchDet::DetWithState> detWithStates;
 
   detWithStates = compatibleDets(startingState, prop, est);
 
   vector<DetGroup> result;
   if (!detWithStates.empty()) {
     result.push_back(DetGroup(detWithStates));
   }
   LogTrace("MTDDetLayers") << "MTDdetTray Compatible modules: " << result.size();
   return result;
 }