CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes
MTDDetTray Class Reference

#include <MTDDetTray.h>

Inheritance diagram for MTDDetTray:
DetRodOneR DetRod GeometricSearchDet

Public Member Functions

std::pair< bool, TrajectoryStateOnSurfacecompatible (const TrajectoryStateOnSurface &ts, const Propagator &prop, const MeasurementEstimator &est) const override
 
std::vector< DetWithStatecompatibleDets (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< DetGroupgroupedCompatibleDets (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 PlanespecificSurface () const final
 Return the rod surface as a Plane. More...
 
const BoundSurfacesurface () 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::PositionTypeposition () const
 Returns position of the surface. More...
 
virtual ~GeometricSearchDet ()
 

Private Types

typedef GenericBinFinderInZ< float, GeomDetBinFinderType
 

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 *, TrajectoryStateOnSurfaceDetWithState
 
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

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

19  : DetRodOneR(first, last) {
20  init();
21 }
void init()
Definition: MTDDetTray.cc:25
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() [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(); }
void init()
Definition: MTDDetTray.cc:25
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().

38  {
39  TrajectoryStateOnSurface ms = prop.propagate(ts, specificSurface());
40  return make_pair(ms.isValid() and est.estimate(ms, specificSurface()) != 0, ms);
41 }
virtual const Plane & specificSurface() const final
Return the rod surface as a Plane.
Definition: DetRod.h:24

◆ 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().

45  {
46  LogTrace("MTDDetLayers") << "MTDDetTray::compatibleDets, Surface at R,phi: " << surface().position().perp() << ","
47  << surface().position().phi() << " DetRod pos."
48  << " TS at R,phi: " << startingState.globalPosition().perp() << ","
49  << startingState.globalPosition().phi() << " TSOS dx/dy "
50  << std::sqrt(startingState.localError().positionError().xx()) << " "
51  << std::sqrt(startingState.localError().positionError().yy());
52 
53  vector<DetWithState> result;
54 
55  // Propagate and check that the result is within bounds
56  pair<bool, TrajectoryStateOnSurface> compat = compatible(startingState, prop, est);
57 
58  if (!compat.first) {
59  LogTrace("MTDDetLayers") << " MTDDetTray::compatibleDets: not compatible"
60  << " (should not have been selected!)";
61  return result;
62  }
63 
64  // Find the most probable destination component
65  TrajectoryStateOnSurface& tsos = compat.second;
66  GlobalPoint startPos = tsos.globalPosition();
67  int closest = theBinFinder.binIndex(startPos.z());
68  const vector<const GeomDet*> dets = basicComponents();
69  LogTrace("MTDDetLayers") << " MTDDetTray::compatibleDets, closest det: " << closest
70  << " pos: " << dets[closest]->surface().position() << " impact " << startPos;
71 
72  // Add this detector, if it is compatible
73  // NOTE: add performs a null propagation
74  add(closest, result, tsos, prop, est);
75 
76 #ifdef EDM_ML_DEBUG
77  int nclosest = result.size();
78  int nnextdet = 0; // just DEBUG counters
79 #endif
80 
81  // Try the neighbors on each side until no more compatible.
82  // If closest is not compatible the next cannot be either
83  if (!result.empty()) {
84  const BoundPlane& closestPlane(dets[closest]->surface());
86 
87  // detHalfLen is assumed to be the same for all detectors.
88  float detHalfLen = closestPlane.bounds().length() / 2.;
89 
90  for (unsigned int idet = closest + 1; idet < dets.size(); idet++) {
91  LocalPoint nextPos(dets[idet]->toLocal(startPos));
92  if (fabs(nextPos.y()) < detHalfLen + maxDistance.y()) {
93 #ifdef EDM_ML_DEBUG
94  LogTrace("MTDDetLayers") << " negativeZ: det:" << idet << " pos " << nextPos.y() << " maxDistance "
95  << maxDistance.y();
96  nnextdet++;
97 #endif
98  if (!add(idet, result, tsos, prop, est))
99  break;
100  } else {
101  break;
102  }
103  }
104 
105  for (int idet = closest - 1; idet >= 0; idet--) {
106  LocalPoint nextPos(dets[idet]->toLocal(startPos));
107  if (fabs(nextPos.y()) < detHalfLen + maxDistance.y()) {
108 #ifdef EDM_ML_DEBUG
109  LogTrace("MTDDetLayers") << " positiveZ: det:" << idet << " pos " << nextPos.y() << " maxDistance "
110  << maxDistance.y();
111  nnextdet++;
112 #endif
113  if (!add(idet, result, tsos, prop, est))
114  break;
115  } else {
116  break;
117  }
118  }
119  }
120 
121 #ifdef EDM_ML_DEBUG
122  LogTrace("MTDDetLayers") << " MTDDetTray::compatibleDets, size: " << result.size() << " on closest: " << nclosest
123  << " # checked dets: " << nnextdet + 1;
124  if (result.empty()) {
125  LogTrace("MTDDetLayers") << " ***Rod not compatible---should have been discarded before!!!";
126  }
127 #endif
128  return result;
129 }
const BoundSurface & surface() const final
The surface of the GeometricSearchDet.
Definition: DetRod.h:19
T perp() const
Definition: PV3DBase.h:69
const LocalTrajectoryError & localError() const
T z() const
Definition: PV3DBase.h:61
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
int closest(std::vector< int > const &vec, int value)
virtual Local2DVector maximalLocalDisplacement(const TrajectoryStateOnSurface &ts, const Plane &plane) const =0
LocalError positionError() const
#define LogTrace(id)
float yy() const
Definition: LocalError.h:24
GlobalPoint globalPosition() const
T sqrt(T t)
Definition: SSEVec.h:19
int binIndex(T z) const override
returns an index in the valid range for the bin closest to Z
LocalVector toLocal(const reco::Track::Vector &v, const Surface &s)
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.
Definition: DetRodOneR.cc:32
const PositionType & position() const
const std::vector< const GeomDet * > & basicComponents() const override
Definition: DetRodOneR.h:26
std::pair< bool, TrajectoryStateOnSurface > compatible(const TrajectoryStateOnSurface &ts, const Propagator &prop, const MeasurementEstimator &est) const override
Definition: MTDDetTray.cc:36
BinFinderType theBinFinder
Definition: MTDDetTray.h:46
float xx() const
Definition: LocalError.h:22

◆ 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.

29  {
30  // FIXME dummy impl.
31  edm::LogError("MTDDetLayers") << "temporary dummy implementation of MTDDetTray::components()!!";
32  static const vector<const GeometricSearchDet*> result;
33  return result;
34 }
Log< level::Error, false > LogError

◆ 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.

133  {
134  vector<GeometricSearchDet::DetWithState> detWithStates;
135 
136  detWithStates = compatibleDets(startingState, prop, est);
137 
138  vector<DetGroup> result;
139  if (!detWithStates.empty()) {
140  result.push_back(DetGroup(detWithStates));
141  }
142  LogTrace("MTDDetLayers") << "MTDdetTray Compatible modules: " << result.size();
143  return result;
144 }
std::vector< DetWithState > compatibleDets(const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const override
Definition: MTDDetTray.cc:43
#define LogTrace(id)

◆ init()

void MTDDetTray::init ( void  )
private

Definition at line 25 of file MTDDetTray.cc.

References DetRodOneR::basicComponents(), mps_fire::end, and theBinFinder.

Referenced by MTDDetTray().

GenericBinFinderInZ< float, GeomDet > BinFinderType
Definition: MTDDetTray.h:45
const std::vector< const GeomDet * > & basicComponents() const override
Definition: DetRodOneR.h:26
BinFinderType theBinFinder
Definition: MTDDetTray.h:46

Member Data Documentation

◆ theBinFinder

BinFinderType MTDDetTray::theBinFinder
private

Definition at line 46 of file MTDDetTray.h.

Referenced by compatibleDets(), and init().