#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 ( const std::vector< const MTDDetSector * > & sectors )

Constructor, takes ownership of pointers.

Definition at line 18 of file MTDSectorForwardLayer.cc.

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

     : 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();
 }

MTDSectorForwardLayer::~MTDSectorForwardLayer ( )

override

Definition at line 56 of file MTDSectorForwardLayer.cc.

References mps_fire::i, and theSectors.

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

Member Function Documentation

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

inlineoverridevirtual

Implements GeometricSearchDet.

Definition at line 18 of file MTDSectorForwardLayer.h.

References theBasicComps.

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

18 { return theBasicComps; }

MTDSectorForwardLayer::theBasicComps

std::vector< const GeomDet * > theBasicComps

Definition: MTDSectorForwardLayer.h:41

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.

References ForwardDetLayer::compatible(), 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().

                                                                                                                   {
   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;
 }

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.

References theComponents.

140 { return theComponents; }

MTDSectorForwardLayer::theComponents

std::vector< const GeometricSearchDet * > theComponents

Definition: MTDSectorForwardLayer.h:40

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>();
 }

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

inlinevirtual

Return the vector of sectors.

Definition at line 36 of file MTDSectorForwardLayer.h.

References theSectors.

36 { return theSectors; }

MTDSectorForwardLayer::theSectors

std::vector< const MTDDetSector * > theSectors

Definition: MTDSectorForwardLayer.h:39

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.

References theBasicComps.

Referenced by MTDSectorForwardDoubleLayer::subDetector().

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

Member Data Documentation

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

private

Definition at line 41 of file MTDSectorForwardLayer.h.

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

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

private

Definition at line 40 of file MTDSectorForwardLayer.h.

Referenced by components().

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

private

Definition at line 39 of file MTDSectorForwardLayer.h.

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

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation