#include <MTDDetSector.h>

Inheritance diagram for MTDDetSector:

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

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

void	compatibleDetsLine (const size_t idetMin, std::vector< DetWithState > &result, const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est) const

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

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

	GeometricSearchDet (bool doHaveGroups)

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

size_t	hshift (const uint32_t detid, const int horizontalShift) const

	MTDDetSector (const std::vector< const GeomDet * > &dets, const MTDTopology &topo)
	Construct from a vector of GeomDet*. More...

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

const BoundDiskSector &	specificSurface () const

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

size_t	vshift (const uint32_t detid, const int verticalShift, size_t &closest) const

	~MTDDetSector () override

Public Member Functions inherited from GeometricSearchDet
	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 ()

Protected Member Functions
bool	add (size_t idet, std::vector< DetWithState > &result, const TrajectoryStateOnSurface &tsos, const Propagator &prop, const MeasurementEstimator &est) const

void	setDisk (BoundDiskSector *diskS)

Private Member Functions
void	init ()

Private Attributes
std::vector< const GeomDet * >	theDets

ReferenceCountingPointer< BoundDiskSector >	theDiskS

const MTDTopology *	topo_

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

Definition at line 12 of file MTDDetSector.h.

Constructor & Destructor Documentation

◆ MTDDetSector() [1/2]

MTDDetSector::MTDDetSector	(	std::vector< const GeomDet * >::const_iterator	first,
		std::vector< const GeomDet * >::const_iterator	last,
		const MTDTopology &	topo
	)

Construct from iterators on GeomDet*.

Definition at line 18 of file MTDDetSector.cc.

     : GeometricSearchDet(false), theDets(first, last), topo_(&topo) {
   init();
 }

References init().

◆ MTDDetSector() [2/2]

MTDDetSector::MTDDetSector	(	const std::vector< const GeomDet * > &	dets,
		const MTDTopology &	topo
	)

Construct from a vector of GeomDet*.

Definition at line 25 of file MTDDetSector.cc.

     : GeometricSearchDet(false), theDets(vdets), topo_(&topo) {
   init();
 }

References init().

◆ ~MTDDetSector()

MTDDetSector::~MTDDetSector ( )

inlineoverride

Definition at line 24 of file MTDDetSector.h.

24 {};

Member Function Documentation

◆ add()

bool MTDDetSector::add	(	size_t	idet,
		std::vector< DetWithState > &	result,
		const TrajectoryStateOnSurface &	tsos,
		const Propagator &	prop,
		const MeasurementEstimator &	est
	)		const

protected

Definition at line 162 of file MTDDetSector.cc.

                                                               {
   pair<bool, TrajectoryStateOnSurface> compat = theCompatibilityChecker.isCompatible(theDets[idet], tsos, prop, est);
  
   if (compat.first) {
     result.push_back(DetWithState(theDets[idet], compat.second));
     LogTrace("MTDDetLayers") << "MTDDetSector::compatibleDets found compatible det idetMin " << idet
                              << " detId = " << theDets[idet]->geographicalId().rawId() << " at "
                              << theDets[idet]->position()
                              << " dist = " << std::sqrt((tsos.globalPosition() - theDets[idet]->position()).mag2());
   }
  
   return compat.first;
 }

References TrajectoryStateOnSurface::globalPosition(), GeomDetCompatibilityChecker::isCompatible(), LogTrace, mag2(), mps_fire::result, mathSSE::sqrt(), GeometricSearchDet::theCompatibilityChecker, and theDets.

Referenced by compatibleDets(), compatibleDetsLine(), and counter.Counter::register().

◆ basicComponents()

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

inlineoverridevirtual

Implements GeometricSearchDet.

Definition at line 28 of file MTDDetSector.h.

28 { return theDets; }

References theDets.

Referenced by compatibleDets().

◆ compatible()

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

                                                                                                      {
   TrajectoryStateOnSurface ms = prop.propagate(ts, specificSurface());
  
 #ifdef EDM_ML_DEBUG
   LogTrace("MTDDetLayers") << "MTDDetSector::compatible, sector: \n"
                            << (*this) << "\n  TS at Z,R,phi: " << std::fixed << std::setw(14) << ts.globalPosition().z()
                            << " , " << std::setw(14) << ts.globalPosition().perp() << " , " << std::setw(14)
                            << ts.globalPosition().phi();
   if (ms.isValid()) {
     LogTrace("MTDDetLayers") << " DEST at Z,R,phi: " << std::fixed << std::setw(14) << ms.globalPosition().z() << " , "
                              << std::setw(14) << ms.globalPosition().perp() << " , " << std::setw(14)
                              << ms.globalPosition().phi() << " local Z: " << std::setw(14) << ms.localPosition().z();
   } else {
     LogTrace("MTDDetLayers") << " DEST: not valid";
   }
 #endif
  
   return make_pair(ms.isValid() and est.estimate(ms, specificSurface()) != 0, ms);
 }

References MeasurementEstimator::estimate(), alignBH_cfg::fixed, TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localPosition(), LogTrace, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), Propagator::propagate(), specificSurface(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by compatibleDets().

◆ compatibleDets()

vector< GeometricSearchDet::DetWithState > MTDDetSector::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 66 of file MTDDetSector.cc.

                                                                                                              {
   LogTrace("MTDDetLayers") << "MTDDetSector::compatibleDets, sector: \n"
                            << (*this) << "\n  TS at Z,R,phi: " << std::fixed << std::setw(14)
                            << startingState.globalPosition().z() << " , " << std::setw(14)
                            << startingState.globalPosition().perp() << " , " << std::setw(14)
                            << startingState.globalPosition().phi();
  
   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") << "    MTDDetSector::compatibleDets: not compatible"
                              << "    (should not have been selected!)";
     return result;
   }
  
   TrajectoryStateOnSurface& tsos = compat.second;
   GlobalPoint startPos = tsos.globalPosition();
  
   LogTrace("MTDDetLayers") << "Starting position: " << startPos << " starting p/pT: " << tsos.globalMomentum().mag()
                            << " / " << tsos.globalMomentum().perp();
  
   // determine distance of det center from extrapolation on the surface, sort dets accordingly
  
   size_t idetMin = basicComponents().size();
   double dist2Min = std::numeric_limits<double>::max();
   std::vector<std::pair<double, size_t> > tmpDets;
   tmpDets.reserve(basicComponents().size());
  
   for (size_t idet = 0; idet < basicComponents().size(); idet++) {
     double dist2 = (startPos - theDets[idet]->position()).mag2();
     tmpDets.emplace_back(dist2, idet);
     if (dist2 < dist2Min) {
       dist2Min = dist2;
       idetMin = idet;
     }
   }
  
   //look for the compatibledets considering each line of the sector
  
   if (add(idetMin, result, tsos, prop, est)) {
     compatibleDetsLine(idetMin, result, tsos, prop, est);
  
     for (int iside = -1; iside <= 1; iside += 2) {
       bool isCompatible(true);
       size_t idetNew(idetMin);
       size_t closest = theDets.size();
  
       while (isCompatible) {
         idetNew = vshift(theDets[idetNew]->geographicalId().rawId(), iside, closest);
         if (idetNew >= theDets.size()) {
           if (closest < theDets.size()) {
             idetNew = closest;
           } else {
             break;
           }
         }
         isCompatible = add(idetNew, result, tsos, prop, est);
         if (isCompatible) {
           compatibleDetsLine(idetNew, result, tsos, prop, est);
         }
       }
     }
   }
  
 #ifdef EDM_ML_DEBUG
   if (result.empty()) {
     LogTrace("MTDDetLayers") << "MTDDetSector::compatibleDets, closest not compatible!";
   } else {
     LogTrace("MTDDetLayers") << "MTDDetSector::compatibleDets, found " << result.size() << " compatible dets";
   }
 #endif
  
   return result;
 }

References add(), basicComponents(), compatible(), compatibleDetsLine(), alignBH_cfg::fixed, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), LogTrace, PV3DBase< T, PVType, FrameType >::mag(), mag2(), SiStripPI::max, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), mps_fire::result, findQualityFiles::size, theDets, vshift(), and PV3DBase< T, PVType, FrameType >::z().

◆ compatibleDetsLine()

void MTDDetSector::compatibleDetsLine	(	const size_t	idetMin,
		std::vector< DetWithState > &	result,
		const TrajectoryStateOnSurface &	tsos,
		const Propagator &	prop,
		const MeasurementEstimator &	est
	)		const

Definition at line 192 of file MTDDetSector.cc.

                                                                              {
   for (int iside = -1; iside <= 1; iside += 2) {
     bool isCompatible(true);
     size_t idetNew(idetMin);
  
     while (isCompatible) {
       idetNew = hshift(theDets[idetNew]->geographicalId().rawId(), iside);
       if (idetNew >= theDets.size()) {
         break;
       }
       isCompatible = add(idetNew, result, tsos, prop, est);
     }
   }
  
   return;
 }

References add(), hshift(), mps_fire::result, and theDets.

Referenced by compatibleDets().

◆ compatibleDetsV()

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

overridevirtual

Reimplemented from GeometricSearchDet.

Definition at line 145 of file MTDDetSector.cc.

                                                                    {
   edm::LogError("MTDDetLayers") << "At the moment not a real implementation";
 }

◆ components()

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

overridevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 37 of file MTDDetSector.cc.

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

References mps_fire::result.

◆ GeometricSearchDet()

GeometricSearchDet::GeometricSearchDet

inline

Definition at line 24 of file GeometricSearchDet.h.

24 : haveGroups(doHaveGroups) {}

◆ groupedCompatibleDets()

vector< DetGroup > MTDDetSector::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 152 of file MTDDetSector.cc.

                                                                                             {
   // FIXME should be implemented to allow returning  overlapping chambers
   // as separate groups!
   edm::LogInfo("MTDDetLayers") << "dummy implementation of MTDDetSector::groupedCompatibleDets()";
   vector<DetGroup> result;
   return result;
 }

References mps_fire::result.

◆ hshift()

size_t MTDDetSector::hshift	(	const uint32_t	detid,
		const int	horizontalShift
	)		const

Definition at line 213 of file MTDDetSector.cc.

                                                                                  {
   return topo_->hshiftETL(detid, horizontalShift);
 }

References MTDTopology::hshiftETL(), and topo_.

Referenced by compatibleDetsLine().

◆ init()

void MTDDetSector::init ( void )

private

Definition at line 30 of file MTDDetSector.cc.

                         {
   // Add here the sector build based on a collection of GeomDets, mimic what done in ForwardDetRingOneZ
   // using the code from tracker BladeShapeBuilderFromDet
   // simple initial version, no sorting for the time being
   setDisk(MTDDiskSectorBuilderFromDet()(theDets));
 }

References setDisk(), and theDets.

Referenced by MTDDetSector().

◆ setDisk()

void MTDDetSector::setDisk ( BoundDiskSector * diskS )

inlineprotected

Definition at line 65 of file MTDDetSector.h.

65 { theDiskS = diskS; }

References theDiskS.

Referenced by init().

◆ specificSurface()

const BoundDiskSector& MTDDetSector::specificSurface ( ) const

inline

Definition at line 53 of file MTDDetSector.h.

53 { return *theDiskS; }

References theDiskS.

Referenced by compatible().

◆ surface()

const BoundSurface& MTDDetSector::surface ( ) const

inlinefinalvirtual

The surface of the GeometricSearchDet.

Implements GeometricSearchDet.

Definition at line 30 of file MTDDetSector.h.

30 { return *theDiskS; }

References theDiskS.

◆ vshift()

size_t MTDDetSector::vshift	(	const uint32_t	detid,
		const int	verticalShift,
		size_t &	closest
	)		const

Definition at line 217 of file MTDDetSector.cc.

                                                                                                 {
   return topo_->vshiftETL(detid, verticalShift, closest);
 }

References topo_, and MTDTopology::vshiftETL().

Referenced by compatibleDets().

Member Data Documentation

◆ theDets

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

private

Definition at line 75 of file MTDDetSector.h.

Referenced by add(), basicComponents(), compatibleDets(), compatibleDetsLine(), and init().

◆ theDiskS

ReferenceCountingPointer<BoundDiskSector> MTDDetSector::theDiskS

private

Definition at line 74 of file MTDDetSector.h.

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

◆ topo_

const MTDTopology* MTDDetSector::topo_

private

Definition at line 77 of file MTDDetSector.h.

Referenced by hshift(), and vshift().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MTDDetSector() [1/2]

◆ MTDDetSector() [2/2]

◆ ~MTDDetSector()

Member Function Documentation

◆ add()

◆ basicComponents()

◆ compatible()

◆ compatibleDets()

◆ compatibleDetsLine()

◆ compatibleDetsV()

◆ components()

◆ GeometricSearchDet()

◆ groupedCompatibleDets()

◆ hshift()

◆ init()

◆ setDisk()

◆ specificSurface()

◆ surface()

◆ vshift()

Member Data Documentation

◆ theDets

◆ theDiskS

◆ topo_