MTDSectorForwardLayer Class Reference

#include <MTDSectorForwardLayer.h>

Inheritance diagram for MTDSectorForwardLayer:

Public Member Functions
const std::vector< const GeomDet * > &	basicComponents () 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
	MTDSectorForwardLayer (const std::vector< const MTDDetSector * > &sectors)
	Constructor, takes ownership of pointers. More...
virtual const std::vector< const MTDDetSector * > &	sectors () const
	Return the vector of sectors. More...
SubDetector	subDetector () const override
	The type of detector (PixelBarrel, PixelEndcap, TIB, TOB, TID, TEC, CSC, DT, RPCBarrel, RPCEndcap) More...
	~MTDSectorForwardLayer () override
Public Member Functions inherited from ForwardDetLayer
std::pair< bool, TrajectoryStateOnSurface >	compatible (const TrajectoryStateOnSurface &, const Propagator &, const MeasurementEstimator &) const override
bool	contains (const Local3DPoint &p) const
	ForwardDetLayer (bool doHaveGroups)
Location	location () const final
	Which part of the detector (barrel, endcap) More...
virtual const BoundDisk &	specificSurface () const final
const BoundSurface &	surface () const final
	The surface of the GeometricSearchDet. More...
	~ForwardDetLayer () override
Public Member Functions inherited from DetLayer
	DetLayer (bool doHaveGroup, bool ibar)
bool	isBarrel () const
bool	isForward () const
int	seqNum () const
void	setSeqNum (int sq)
	~DetLayer () 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 Attributes
std::vector< const GeomDet * >	theBasicComps
std::vector< const GeometricSearchDet * >	theComponents
std::vector< const MTDDetSector * >	theSectors

Additional Inherited Members
Public Types inherited from DetLayer
typedef GeomDetEnumerators::Location	Location
typedef GeomDetEnumerators::SubDetector	SubDetector
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 ForwardDetLayer
SimpleDiskBounds const &	bounds () const
virtual BoundDisk *	computeSurface ()
virtual void	initialize ()
float	rmax () const
float	rmin () const
void	setSurface (BoundDisk *cp)
float	zmax () const
float	zmin () const
Protected Attributes inherited from GeometricSearchDet
bool	haveGroups
GeomDetCompatibilityChecker	theCompatibilityChecker

Detailed Description

Definition at line 9 of file MTDSectorForwardLayer.h.

Constructor & Destructor Documentation

MTDSectorForwardLayer()

MTDSectorForwardLayer::MTDSectorForwardLayer

(

const std::vector< const MTDDetSector * > &

sectors

)

Constructor, takes ownership of pointers.

Definition at line 18 of file MTDSectorForwardLayer.cc.

    : ForwardDetLayer(false), theSectors(sectors), theComponents(theSectors.begin(), theSectors.end()) {
  // Initial values for R, Z and Phi bounds
  float theRmin = sectors.front()->basicComponents().front()->position().perp();
  float theRmax = theRmin;
  float theZmin = sectors.front()->position().z();
  float theZmax = theZmin;
 
  // Cache chamber pointers (the basic components_)
  // and find extension in R and Z
  for (const auto& isect : sectors) {
    vector<const GeomDet*> tmp2 = isect->basicComponents();
    theBasicComps.insert(theBasicComps.end(), tmp2.begin(), tmp2.end());
 
    theRmin = min(theRmin, isect->specificSurface().innerRadius());
    theRmax = max(theRmax, isect->specificSurface().outerRadius());
    float halfThick = isect->surface().bounds().thickness() / 2.;
    float zCenter = isect->surface().position().z();
    theZmin = min(theZmin, zCenter - halfThick);
    theZmax = max(theZmax, zCenter + halfThick);
  }
 
  // Build surface
 
  float zPos = (theZmax + theZmin) / 2.;
  PositionType pos(0., 0., zPos);
  RotationType rot;
 
  setSurface(new BoundDisk(pos, rot, new SimpleDiskBounds(theRmin, theRmax, theZmin - zPos, theZmax - zPos)));
 
  LogTrace("MTDDetLayers") << "Constructing MTDSectorForwardLayer: " << std::fixed << std::setw(14)
                           << basicComponents().size() << " Dets, " << std::setw(14) << theSectors.size()
                           << " Sectors, "
                           << " Z: " << std::setw(14) << specificSurface().position().z() << " R1: " << std::setw(14)
                           << specificSurface().innerRadius() << " R2: " << std::setw(14)
                           << specificSurface().outerRadius();
}

References basicComponents(), alignBH_cfg::fixed, LogTrace, SiStripPI::max, min(), makeMuonMisalignmentScenario::rot, sectors(), ForwardDetLayer::setSurface(), ForwardDetLayer::specificSurface(), theBasicComps, and theSectors.

~MTDSectorForwardLayer()

MTDSectorForwardLayer::~MTDSectorForwardLayer ( )

override

Definition at line 56 of file MTDSectorForwardLayer.cc.

                                              {
  for (auto& i : theSectors) {
    delete i;
  }
}

References mps_fire::i, and theSectors.

Member Function Documentation

basicComponents()

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

inlineoverridevirtual

Implements GeometricSearchDet.

Definition at line 18 of file MTDSectorForwardLayer.h.

{ return theBasicComps; }

References theBasicComps.

Referenced by MTDSectorForwardLayer(), and MTDSectorForwardDoubleLayer::selfTest().

compatibleDets()

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

                                                                                                                  {
  vector<DetWithState> result;
 
  LogTrace("MTDDetLayers") << "MTDSectorForwardLayer::compatibleDets,"
                           << " R1 " << std::fixed << std::setw(14) << specificSurface().innerRadius()
                           << " R2: " << std::setw(14) << specificSurface().outerRadius()
                           << " FTS at R: " << std::setw(14) << startingState.globalPosition().perp();
 
  pair<bool, TrajectoryStateOnSurface> compat = compatible(startingState, prop, est);
 
  if (!compat.first) {
    LogTrace("MTDDetLayers") << "     MTDSectorForwardLayer::compatibleDets: not compatible"
                             << " (should not have been selected!)";
    return result;
  }
 
  TrajectoryStateOnSurface& tsos = compat.second;
 
  // as there are either two or four sectors only, avoid complex logic and just loop on all of them
 
  // Use state on layer surface. Note that local coordinates and errors
  // are the same on the layer and on all sectors surfaces, since
  // all BoundDisks are centered in 0,0 and have the same rotation.
  // CAVEAT: if the sectors are not at the same Z, the local position and error
  // will be "Z-projected" to the sectors. This is a fairly good approximation.
  // However in this case additional propagation will be done when calling
  // compatibleDets.
  GlobalPoint startPos = tsos.globalPosition();
 
  for (unsigned int isect = 0; isect < theSectors.size(); isect++) {
    LocalPoint nextPos(theSectors[isect]->specificSurface().toLocal(startPos));
    LogDebug("MTDDetLayers") << "Global point = " << std::fixed << startPos << " local point = " << nextPos
                             << " global sector ref pos = " << theSectors[isect]->specificSurface().position();
    bool inside = false;
    if (tsos.hasError()) {
      inside = theSectors[isect]->specificSurface().bounds().inside(nextPos, tsos.localError().positionError(), 1.);
    } else {
      inside = theSectors[isect]->specificSurface().bounds().inside(nextPos);
    }
    if (inside) {
#ifdef EDM_ML_DEBUG
      LogTrace("MTDDetLayers") << "     MTDSectorForwardLayer::fastCompatibleDets:NextSector " << std::fixed
                               << std::setw(14) << isect << "\n"
                               << (*theSectors[isect]) << "\n FTS at Z,R,phi: " << std::setw(14)
                               << tsos.globalPosition().z() << " , " << std::setw(14) << tsos.globalPosition().perp()
                               << "," << std::setw(14) << tsos.globalPosition().phi();
      if (tsos.hasError()) {
        LogTrace("MTDDetLayers") << " sR: " << sqrt(tsos.localError().positionError().yy())
                                 << " sX: " << sqrt(tsos.localError().positionError().xx());
      }
#endif
      vector<DetWithState> nextRodDets(theSectors[isect]->compatibleDets(tsos, prop, est));
      if (!nextRodDets.empty()) {
        result.insert(result.end(), nextRodDets.begin(), nextRodDets.end());
      } else {
        break;
      }
    }
  }
 
  LogTrace("MTDDetLayers") << "     MTDSectorForwardLayer::fastCompatibleDets: found: " << result.size();
 
  return result;
}

References ForwardDetLayer::compatible(), alignBH_cfg::fixed, TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::hasError(), TrajectoryStateOnSurface::localError(), LogDebug, LogTrace, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), LocalTrajectoryError::positionError(), mps_fire::result, ForwardDetLayer::specificSurface(), mathSSE::sqrt(), theSectors, toLocal(), LocalError::xx(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by MTDSectorForwardDoubleLayer::groupedCompatibleDets().

components()

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

overridevirtual

Returns basic components, if any.

Returns direct components, if any

Implements GeometricSearchDet.

Definition at line 140 of file MTDSectorForwardLayer.cc.

{ return theComponents; }

References theComponents.

groupedCompatibleDets()

vector< DetGroup > MTDSectorForwardLayer::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 128 of file MTDSectorForwardLayer.cc.

                                                                                                     {
  // FIXME should return only 1 group
  edm::LogInfo("MTDDetLayers") << "dummy implementation of MTDSectorForwardLayer::groupedCompatibleDets()";
  return vector<DetGroup>();
}

sectors()

virtual const std::vector<const MTDDetSector*>& MTDSectorForwardLayer::sectors ( ) const

inlinevirtual

Return the vector of sectors.

Definition at line 36 of file MTDSectorForwardLayer.h.

{ return theSectors; }

References theSectors.

Referenced by MTDSectorForwardLayer().

subDetector()

GeomDetEnumerators::SubDetector MTDSectorForwardLayer::subDetector ( ) const

overridevirtual

The type of detector (PixelBarrel, PixelEndcap, TIB, TOB, TID, TEC, CSC, DT, RPCBarrel, RPCEndcap)

Implements DetLayer.

Definition at line 136 of file MTDSectorForwardLayer.cc.

                                                                       {
  return theBasicComps.front()->subDetector();
}

References theBasicComps.

Referenced by MTDSectorForwardDoubleLayer::subDetector().

Member Data Documentation

theBasicComps

std::vector<const GeomDet*> MTDSectorForwardLayer::theBasicComps

private

Definition at line 41 of file MTDSectorForwardLayer.h.

Referenced by basicComponents(), MTDSectorForwardLayer(), and subDetector().

theComponents

std::vector<const GeometricSearchDet*> MTDSectorForwardLayer::theComponents

private

Definition at line 40 of file MTDSectorForwardLayer.h.

Referenced by components().

theSectors

std::vector<const MTDDetSector*> MTDSectorForwardLayer::theSectors

private

Definition at line 39 of file MTDSectorForwardLayer.h.

Referenced by compatibleDets(), MTDSectorForwardLayer(), sectors(), and ~MTDSectorForwardLayer().