#include <BarrelDetLayer.h>

Inheritance diagram for BarrelDetLayer:

Public Member Functions
	BarrelDetLayer (bool doHaveGroup)

std::pair< bool, TrajectoryStateOnSurface >	compatible (const TrajectoryStateOnSurface &ts, const Propagator &, const MeasurementEstimator &) const final

bool	contains (const Local3DPoint &p) const

Location	location () const final
	DetLayer interface. More...

virtual const BoundCylinder &	specificSurface () const final
	Extension of the interface. More...

const BoundSurface &	surface () const final
	GeometricSearchDet interface. More...

	~BarrelDetLayer () 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)

virtual SubDetector	subDetector () const =0
	The type of detector (PixelBarrel, PixelEndcap, TIB, TOB, TID, TEC, CSC, DT, RPCBarrel, RPCEndcap) More...

	~DetLayer () override

Public Member Functions inherited from GeometricSearchDet
virtual const std::vector< const GeomDet * > &	basicComponents () const =0

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

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

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

	GeometricSearchDet (bool doHaveGroups)

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

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
SimpleCylinderBounds const &	bounds () const

virtual BoundCylinder *	computeSurface ()

virtual void	initialize ()

void	setSurface (BoundCylinder *cp)

Private Attributes
ReferenceCountingPointer< BoundCylinder >	theCylinder

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 Attributes inherited from GeometricSearchDet
bool	haveGroups

GeomDetCompatibilityChecker	theCompatibilityChecker

Detailed Description

A specialization of the DetLayer interface for barrel layers. Barrel layers are cylinders with their axes parallel to the global Z axis. The methods that have a common implementation for all BarrelDetLayers are implemented in this class, but some methods are left abstract.

Definition at line 22 of file BarrelDetLayer.h.

Constructor & Destructor Documentation

◆ BarrelDetLayer()

BarrelDetLayer::BarrelDetLayer ( bool doHaveGroup )

inline

Definition at line 24 of file BarrelDetLayer.h.

24 : DetLayer(doHaveGroup, true), theCylinder(nullptr) {}

◆ ~BarrelDetLayer()

BarrelDetLayer::~BarrelDetLayer ( )

override

Definition at line 11 of file BarrelDetLayer.cc.

11 {}

Member Function Documentation

◆ bounds()

SimpleCylinderBounds const& BarrelDetLayer::bounds ( ) const

inlineprotected

Definition at line 49 of file BarrelDetLayer.h.

49 { return static_cast<SimpleCylinderBounds const&>(theCylinder->bounds()); }

References theCylinder.

◆ compatible()

pair< bool, TrajectoryStateOnSurface > BarrelDetLayer::compatible	(	const TrajectoryStateOnSurface &	ts,
		const Propagator &	,
		const MeasurementEstimator &
	)		const

finalvirtual

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 58 of file BarrelDetLayer.cc.

                                                                                                    {
   if
     UNLIKELY(theCylinder == nullptr)
   edm::LogError("DetLayers") << "ERROR: BarrelDetLayer::compatible() is used before the layer surface is initialized";
   // throw an exception? which one?
  
   TrajectoryStateOnSurface myState = prop.propagate(ts, specificSurface());
   if
     UNLIKELY(!myState.isValid()) return make_pair(false, myState);
  
   // check z assuming symmetric bounds around position().z()
   auto z0 = std::abs(myState.globalPosition().z() - specificSurface().position().z());
   auto deltaZ = 0.5f * bounds().length();
   if (z0 < deltaZ)
     return make_pair(true, myState);
  
   // take into account the thickness of the layer
   deltaZ += 0.5f * bounds().thickness() * std::abs(myState.globalDirection().z()) / myState.globalDirection().perp();
  
   // take also into account the error on the predicted state
   const float nSigma = 3.;
   if (myState.hasError())
     deltaZ += nSigma * sqrt(myState.cartesianError().position().czz());
   //
   //  check z again
   return make_pair(z0 < deltaZ, myState);
 }

References funct::abs(), TrajectoryStateOnSurface::cartesianError(), GlobalErrorBase< T, ErrorWeightType >::czz(), TrajectoryStateOnSurface::globalDirection(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::hasError(), TrajectoryStateOnSurface::isValid(), HLTSiStripMonitoring_cff::nSigma, PV3DBase< T, PVType, FrameType >::perp(), CartesianTrajectoryError::position(), Propagator::propagate(), mathSSE::sqrt(), UNLIKELY, PV3DBase< T, PVType, FrameType >::z(), and HLTMuonOfflineAnalyzer_cfi::z0.

◆ computeSurface()

BoundCylinder * BarrelDetLayer::computeSurface ( )

protectedvirtual

Definition at line 21 of file BarrelDetLayer.cc.

                                               {
   vector<const GeomDet*> comps = basicComponents();
  
   // Find extension in Z
   float theRmin = comps.front()->position().perp();
   float theRmax = theRmin;
   float theZmin = comps.front()->position().z();
   float theZmax = theZmin;
   for (vector<const GeomDet*>::const_iterator deti = comps.begin(); deti != comps.end(); deti++) {
     vector<GlobalPoint> corners = BoundingBox().corners(dynamic_cast<const Plane&>((*deti)->surface()));
     for (vector<GlobalPoint>::const_iterator ic = corners.begin(); ic != corners.end(); ic++) {
       float r = ic->perp();
       float z = ic->z();
       theRmin = min(theRmin, r);
       theRmax = max(theRmax, r);
       theZmin = min(theZmin, z);
       theZmax = max(theZmax, z);
     }
     // in addition to the corners we have to check the middle of the
     // det +/- thickness/2
     // , since the min  radius for some barrel dets is reached there
     float rdet = (**deti).position().perp();
     float thick = (**deti).surface().bounds().thickness();
     theRmin = min(theRmin, rdet - thick / 2.F);
     theRmax = max(theRmax, rdet + thick / 2.F);
   }
  
   // By default the barrel layers are positioned at the center of the
   // global frame, and the axes of their local frame coincide with
   // those of the global grame (z along the cylinder axis)
   PositionType pos(0., 0., 0.);
   RotationType rot;
  
   auto scp = new SimpleCylinderBounds(theRmin, theRmax, theZmin, theZmax);
   return new Cylinder(Cylinder::computeRadius(*scp), pos, rot, scp);
 }