#include <CylinderLayer.h>

Inheritance diagram for Trk::CylinderLayer:

Collaboration diagram for Trk::CylinderLayer:

Public Types
using	ChargedTrackParametersUniquePtr = std::unique_ptr< ParametersBase< 5, Trk::Charged > >
	Unique ptr types. More...

using	NeutralTrackParametersUniquePtr = std::unique_ptr< ParametersBase< 5, Trk::Neutral > >

Public Member Functions
	CylinderLayer ()
	Default Constructor. More...

	CylinderLayer (const Amg::Transform3D &transform, CylinderBounds *cbounds, const LayerMaterialProperties &laymatprop, double thickness=0., std::unique_ptr< OverlapDescriptor > od=nullptr, int laytyp=int(Trk::active))
	Constructor with CylinderSurface components and MaterialProperties. More...

	CylinderLayer (CylinderSurface *cyl, const LayerMaterialProperties &laymatprop, double thickness=0., std::unique_ptr< OverlapDescriptor > od=nullptr, int laytyp=int(Trk::active))
	Constructor with CylinderSurface and MaterialProperties. More...

	CylinderLayer (const Amg::Transform3D &transform, CylinderBounds cbounds, std::unique_ptr< SurfaceArray > surfaceArray, double thickness=0., std::unique_ptr< OverlapDescriptor > od=nullptr, IApproachDescriptor ad=nullptr, int laytyp=int(Trk::active))
	Constructor with CylinderSurface components and pointer to SurfaceArray (passing ownership) More...

	CylinderLayer (const Amg::Transform3D &transform, CylinderBounds cbounds, std::unique_ptr< SurfaceArray > surfaceArray, const LayerMaterialProperties &laymatprop, double thickness=0., std::unique_ptr< OverlapDescriptor > od=nullptr, IApproachDescriptor ad=nullptr, int laytyp=int(Trk::active))
	Constructor with CylinderSurface components, MaterialProperties and pointer SurfaceArray (passing ownership) More...

	CylinderLayer (CylinderBounds *cbounds, const LayerMaterialProperties &laymatprop, double thickness=0., std::unique_ptr< OverlapDescriptor > od=nullptr, int laytyp=int(Trk::active))
	Concentric Layer: Constructor with CylinderSurface components and MaterialProperties. More...

	CylinderLayer (CylinderBounds cbounds, std::unique_ptr< SurfaceArray > surfaceArray, double thickness=0., std::unique_ptr< OverlapDescriptor > od=nullptr, IApproachDescriptor ad=nullptr, int laytyp=int(Trk::active))
	Concentric Layer: Constructor with CylinderSurface components and pointer to SurfaceArray (passing ownership) More...

	CylinderLayer (CylinderBounds cbounds, std::unique_ptr< SurfaceArray > surfaceArray, const LayerMaterialProperties &laymatprop, double thickness=0., std::unique_ptr< OverlapDescriptor > od=nullptr, IApproachDescriptor ad=nullptr, int laytyp=int(Trk::active))
	Concentric Layer: Constructor with CylinderSurface components, MaterialProperties and pointer SurfaceArray (passing ownership) More...

	CylinderLayer (const CylinderLayer &cla)
	Copy constructor of CylinderLayer. More...

	CylinderLayer (const CylinderLayer &cla, const Amg::Transform3D &tr)
	Copy constructor with shift. More...

CylinderLayer &	operator= (const CylinderLayer &)
	Assignment operator for CylinderLayers. More...

virtual	~CylinderLayer ()=default
	Destructor. More...

virtual const CylinderSurface &	surfaceRepresentation () const override final
	Transforms the layer into a Surface representation for extrapolation. More...

virtual CylinderSurface &	surfaceRepresentation () override final

virtual double	preUpdateMaterialFactor (const Trk::TrackParameters &par, Trk::PropDirection dir) const override final
	getting the MaterialProperties back - for pre-update More...

virtual double	postUpdateMaterialFactor (const Trk::TrackParameters &par, Trk::PropDirection dir) const override final
	getting the MaterialProperties back - for post-update More...

virtual const Surface &	surfaceOnApproach (const Amg::Vector3D &pos, const Amg::Vector3D &dir, PropDirection pdir, const BoundaryCheck &bcheck, bool resolveSubSurfaces=0, const ICompatibilityEstimator *ice=nullptr) const override final
	Surface seen on approach - if not defined differently, it is the surfaceRepresentation() More...

virtual void	moveLayer (Amg::Transform3D &shift) override final
	move the Layer More...

virtual void	resizeLayer (const VolumeBounds &vBounds, double envelope) override final
	Resize the layer to the tracking volume - only works for CylinderVolumeBouns. More...

virtual void	resizeAndRepositionLayer (const VolumeBounds &vBounds, const Amg::Vector3D &cCenter, double envelope) override final
	Resize the layer to the tracking volume. More...

virtual bool	operator== (const Surface &sf) const override
	Equality operator. More...

bool	operator== (const CylinderSurface &cf) const

virtual CylinderSurface *	clone () const override
	Implicit Constructor. More...

virtual Surface::ChargedTrackParametersUniquePtr	createUniqueTrackParameters (double l1, double l2, double phi, double theta, double qop, std::optional< AmgSymMatrix(5)> cov=std::nullopt) const override final
	Use the Surface as a ParametersBase constructor, from local parameters - charged. More...

virtual Surface::ChargedTrackParametersUniquePtr	createUniqueTrackParameters (const Amg::Vector3D &position, const Amg::Vector3D &momentum, double charge, std::optional< AmgSymMatrix(5)> cov=std::nullopt) const override final
	Use the Surface as a ParametersBase constructor, from global parameters - charged. More...

virtual NeutralTrackParametersUniquePtr	createUniqueNeutralParameters (double l1, double l2, double phi, double theta, double qop, std::optional< AmgSymMatrix(5)> cov=std::nullopt) const override final
	Use the Surface as a ParametersBase constructor, from local parameters - neutral. More...

virtual NeutralTrackParametersUniquePtr	createUniqueNeutralParameters (const Amg::Vector3D &position, const Amg::Vector3D &momentum, double charge, std::optional< AmgSymMatrix(5)> cov=std::nullopt) const override final
	Use the Surface as a ParametersBase constructor, from global parameters - neutral. More...

template<int DIM, class T >
std::unique_ptr< ParametersT< DIM, T, CylinderSurface > >	createUniqueParameters (double l1, double l2, double phi, double theta, double qop, std::optional< AmgSymMatrix(DIM)> cov=std::nullopt) const
	Use the Surface as a ParametersBase constructor, from local parameters. More...

template<int DIM, class T >
std::unique_ptr< ParametersT< DIM, T, CylinderSurface > >	createUniqueParameters (const Amg::Vector3D &position, const Amg::Vector3D &momentum, double charge, std::optional< AmgSymMatrix(DIM)> cov=std::nullopt) const
	Use the Surface as a ParametersBase constructor, from global parameters. More...

virtual Amg::RotationMatrix3D	measurementFrame (const Amg::Vector3D &glopos, const Amg::Vector3D &glomom) const override final
	Return the measurement frame - this is needed for alignment, in particular for StraightLine and Perigee Surface. More...

constexpr virtual SurfaceType	type () const override final
	Return the surface type. More...

virtual const Amg::Vector3D &	globalReferencePoint () const override final
	Returns a global reference point: For the Cylinder this is \( (Rcos(\phi), Rsin(\phi),0)*transform() \) Where \( \phi \) denotes the averagePhi() of the cylinderBounds. More...

virtual Amg::Vector3D	normal (const Amg::Vector2D &locpo) const override final
	Return method for surface normal information at a given local point, overwrites the normal() from base class. More...

virtual const Amg::Vector3D &	normal () const
	Returns the normal vector of the Surface (i.e. More...

virtual const Amg::Vector3D &	rotSymmetryAxis () const
	Return method for the rotational symmetry axis - the z-Axis of the HepTransform. More...

virtual const CylinderBounds &	bounds () const override final
	This method returns the CylinderBounds by reference (NoBounds is not possible for cylinder) More...

bool	hasBounds () const

virtual bool	insideBounds (const Amg::Vector2D &locpos, double tol1=0., double tol2=0.) const override
	This method calls the inside method of CylinderBounds. More...

virtual bool	insideBoundsCheck (const Amg::Vector2D &locpos, const BoundaryCheck &bchk) const override final

virtual Amg::Vector2D	localParametersToPosition (const LocalParameters &locpars) const override final
	Specialized for CylinderSurface : LocalParameters to Vector2D. More...

virtual void	localToGlobal (const Amg::Vector2D &locp, const Amg::Vector3D &mom, Amg::Vector3D &glob) const override
	Specialized for CylinderSurface : LocalToGlobal method without dynamic memory allocation. More...

Amg::Vector3D	localToGlobal (const Amg::Vector2D &locpos) const
	This method returns the GlobalPosition from a LocalPosition uses the per surface localToGlobal. More...

Amg::Vector3D	localToGlobal (const Amg::Vector2D &locpos, const Amg::Vector3D &glomom) const
	This method returns the GlobalPosition from a LocalPosition The LocalPosition can be outside Surface bounds - for generality with momentum. More...

Amg::Vector3D	localToGlobal (const LocalParameters &locpars) const
	This method returns the GlobalPosition from LocalParameters The LocalParameters can be outside Surface bounds. More...

Amg::Vector3D	localToGlobal (const LocalParameters &locpars, const Amg::Vector3D &glomom) const
	This method returns the GlobalPosition from LocalParameters The LocalParameters can be outside Surface bounds - for generality with momentum. More...

virtual bool	globalToLocal (const Amg::Vector3D &glob, const Amg::Vector3D &mom, Amg::Vector2D &loc) const override
	Specialized for CylinderSurface : GlobalToLocal method without dynamic memory allocation - boolean checks if on surface. More...

std::optional< Amg::Vector2D >	globalToLocal (const Amg::Vector3D &glopos, double tol=0.) const
	This method returns the LocalPosition from a provided GlobalPosition. More...

std::optional< Amg::Vector2D >	globalToLocal (const Amg::Vector3D &glopos, const Amg::Vector3D &glomom) const
	This method returns the LocalPosition from a provided GlobalPosition. More...

virtual bool	isOnSurface (const Amg::Vector3D &glopo, const BoundaryCheck &bchk=true, double tol1=0., double tol2=0.) const override
	This method returns true if the GlobalPosition is on the Surface for both, within or without check of whether the local position is inside boundaries or not. More...

virtual Intersection	straightLineIntersection (const Amg::Vector3D &pos, const Amg::Vector3D &dir, bool forceDir=false, Trk::BoundaryCheck bchk=false) const override final
	fast straight line intersection schema - provides closest intersection and (signed) path length More...

template<class T >
Intersection	straightLineIntersection (const T &pars, bool forceDir=false, const Trk::BoundaryCheck &bchk=false) const
	fst straight line intersection schema - templated for charged and neutral parameters More...

virtual DistanceSolution	straightLineDistanceEstimate (const Amg::Vector3D &pos, const Amg::Vector3D &dir) const override
	fast distance to Surface More...

virtual DistanceSolution	straightLineDistanceEstimate (const Amg::Vector3D &pos, const Amg::Vector3D &dir, bool bound) const override
	fast distance to Surface - with bounds directive More...

virtual double	pathCorrection (const Amg::Vector3D &pos, const Amg::Vector3D &mom) const override
	the pathCorrection for derived classes with thickness More...

virtual std::string	name () const override
	Return properly formatted class name for screen output. More...

bool	operator!= (const Surface &sf) const
	Non-equality operator. More...

std::unique_ptr< Surface >	uniqueClone () const
	NVI method returning unique_ptr clone. More...

const Amg::Transform3D *	cachedTransform () const
	Return the cached transformation directly. More...

const Amg::Transform3D &	transform () const
	Returns HepGeom::Transform3D by reference. More...

const Amg::Vector3D &	center () const
	Returns the center position of the Surface. More...

const TrkDetElementBase *	associatedDetectorElement () const
	return associated Detector Element More...

Identifier	associatedDetectorElementIdentifier () const
	return Identifier of the associated Detector Element More...

const Trk::Layer *	associatedLayer () const
	return the associated Layer More...

const Trk::Layer *	materialLayer () const
	return the material Layer More...

Trk::Layer *	materialLayer ()

virtual const Trk::Surface *	baseSurface () const
	return the base surface (simplified for persistification) More...

std::optional< Amg::Vector2D >	positionOnSurface (const Amg::Vector3D &glopo, const BoundaryCheck &bchk=true, double tol1=0., double tol2=0.) const
	positionOnSurface() returns the LocalPosition on the Surface, If BoundaryCheck==false it just returns the value of globalToLocal (including nullptr possibility), if BoundaryCheck==true it checks whether the point is inside bounds or not (returns std::nullopt in this case). More...

template<class T >
bool	onSurface (const T &parameters, const BoundaryCheck &bchk=BoundaryCheck(true)) const
	The templated Parameters OnSurface method - checks on surface pointer first. More...

bool	isFree () const
	Returns 'true' if this surface is 'free', i.e. More...

bool	isActive () const
	Return 'true' if this surface is owned by the detector element. More...

void	setTransform (const Amg::Transform3D &trans)
	Set the transform updates center and normal. More...

void	setOwner (SurfaceOwner x)
	set Ownership More...

SurfaceOwner	owner () const
	return ownership More...

void	setMaterialLayer (std::shared_ptr< Layer > mlay)
	set material layer More...

virtual MsgStream &	dump (MsgStream &sl) const
	Output Method for MsgStream, to be overloaded by child classes. More...

virtual std::ostream &	dump (std::ostream &sl) const
	Output Method for std::ostream, to be overloaded by child classes. More...

void	associateLayer (const Layer &lay)
	method to associate the associated Trk::Layer which is alreay owned More...

const SurfaceArray *	surfaceArray () const
	Return the entire SurfaceArray, returns nullptr if no SurfaceArray. More...

SurfaceArray *	surfaceArray ()
	Return the entire SurfaceArray, returns nullptr if no SurfaceArray. More...

const Surface *	subSurface (const Amg::Vector3D &gp) const
	If no subSurface array is defined or no subSurface can be found to the given Amg::Vector3D, it would return 0. More...

const Surface *	subSurface (const Amg::Vector2D &lp) const
	If no subSurface array is defined or no subSurface can be found to the given Amg::Vector2D, it would return 0. More...

const Surface *	subSurfaceReference (unsigned int idx=0) const
	Return a reference sub surface of the layer, usually the first one in the array. More...

double	thickness () const
	Return the Thickness of the Layer. More...

template<class T >
bool	onLayer (const T &parameters, const BoundaryCheck &bcheck=BoundaryCheck(true)) const
	templated on Layer method More...

virtual bool	isOnLayer (const Amg::Vector3D &gp, const BoundaryCheck &bcheck=BoundaryCheck(true)) const
	isOnLayer() method, using isOnSurface() with Layer specific tolerance More...

const LayerMaterialProperties *	layerMaterialProperties () const
	getting the LayerMaterialProperties including full/pre/post update More...

const MaterialProperties *	fullUpdateMaterialProperties (const TrackParameters &par) const
	getting the MaterialProperties back - for full update More...

const OverlapDescriptor *	overlapDescriptor () const
	gettint hte overlap descriptor More...

const Layer *	previousLayer (bool skipNavLayer=false) const
	getting what's stored to be the previous Layer, boolean to skip navigation layers More...

void	setPreviousLayer (const Layer *)
	set the previous Layer More...

const Layer *	nextLayer (const Amg::Vector3D &gp, const Amg::Vector3D &udir) const
	getting the next/previous Layer if registered - unit for direction vector required More...

const Layer *	nextLayer (bool skipNavLayer=false) const
	getting what's stored to be the next Layer, boolean to skip navigation layers More...

void	setNextLayer (const Layer *)
	set the next Layer More...

const BinUtility *	binUtility () const
	access the BinUtility More...

void	setBinUtility (const BinUtility *)
	set the BinUtility More...

size_t	compatibleSurfaces (std::vector< SurfaceIntersection > &cSurfaces, const TrackParameters &pars, PropDirection pdir, const BoundaryCheck &bcheck, bool materialSurfacesOnly=true, const Surface startSurface=nullptr, const Surface endSurface=nullptr, const ICompatibilityEstimator *ice=nullptr) const
	get compatible surfaces starting from charged parameters More...

size_t	compatibleSurfaces (std::vector< SurfaceIntersection > &cSurfaces, const NeutralParameters &pars, PropDirection pdir, const BoundaryCheck &bcheck, bool materialSurfacesOnly=true, const Surface startSurface=nullptr, const Surface endSurface=nullptr, const ICompatibilityEstimator *ice=nullptr) const
	get compatible surfaces starting from neutral parameters More...

virtual bool	hasSubStructure (bool resolveSensitive=false) const
	Has sub-structure method: More...

const TrackingVolume *	enclosingTrackingVolume () const
	get the confining TrackingVolume More...

const DetachedTrackingVolume *	enclosingDetachedTrackingVolume () const
	get the confining DetachedTrackingVolume More...

const LayerIndex &	layerIndex () const
	get the layerIndex More...

int	layerType () const
	get the Layer coding More...

void	setLayerType (int identifier)
	set the Layer coding More...

void	assignMaterialProperties (const LayerMaterialProperties &, double scale=1.0)
	assignMaterialPropeties More...

void	registerRepresentingVolume (const Volume *theVol)
	register Volume associated to the layer More...

const Volume *	representingVolume () const
	get the Volume associated to the layer More...

void	setRef (double)
	set the reference measure More...

double	getRef () const
	get the reference measure More...

void	encloseTrackingVolume (const TrackingVolume &tvol)
	private method to set the enclosed detached TV More...

void	encloseDetachedTrackingVolume (const DetachedTrackingVolume &tvol)

template<class T >
size_t	getCompatibleSurfaces (std::vector< SurfaceIntersection > &cSurfaces, const T &pars, PropDirection pdir, const BoundaryCheck &bcheck, bool materialSurfacesOnly=true, const Surface startSurface=nullptr, const Surface endSurface=nullptr, const ICompatibilityEstimator *ice=nullptr) const
	get compatible surfaces starting from charged parameters More...

void	compactify (size_t &cSurfaces, size_t &tSurfaces)
	register layer index for material map registration More...

void	registerLayerIndex (const LayerIndex &lIdx)

Static Public Member Functions
static std::size_t	numberOfInstantiations ()

Static Public Attributes
static constexpr SurfaceType	staticType = SurfaceType::Cylinder
	The surface type static constexpr. More...

static std::atomic_size_t	s_numberOfInstantiations

Protected Member Functions
Amg::Transform3D	inverseTransformHelper () const
	Helper method to factorize in one place common operations calculate inverse transofrm and multiply with position. More...

Amg::Vector3D	inverseTransformMultHelper (const Amg::Vector3D &glopos) const

Protected Attributes
std::unique_ptr< IApproachDescriptor >	m_approachDescriptor
	surfaces on approach to the layer More...

SharedObject< const CylinderBounds >	m_bounds
	The global reference point (== a point on the surface) More...

CxxUtils::CachedUniquePtr< Amg::Vector3D >	m_referencePoint
	The rotational symmetry axis. More...

CxxUtils::CachedUniquePtr< Amg::Vector3D >	m_rotSymmetryAxis

std::unique_ptr< Transforms >	m_transforms {}

const TrkDetElementBase *	m_associatedDetElement {}
	Not owning Pointer to the TrkDetElementBase. More...

Identifier	m_associatedDetElementId {}
	Identifier for the TrkDetElementBase. More...

const Layer *	m_associatedLayer {}
	The associated layer Trk::Layer. More...

std::shared_ptr< Layer >	m_materialLayer {}
	Possibility to attach a material descrption. More...

SurfaceOwner	m_owner {SurfaceOwner::noOwn}
	enum for surface owner : 0 free surface More...

std::unique_ptr< SurfaceArray >	m_surfaceArray
	MaterialPoperties of this layer Surface. More...

std::unique_ptr< LayerMaterialProperties >	m_layerMaterialProperties
	thickness of the Layer More...

double	m_layerThickness
	descriptor for overlap/next surface (owning ptr) More...

std::unique_ptr< OverlapDescriptor >	m_overlapDescriptor

const Layer *	m_previousLayer
	< the previous Layer according to BinGenUtils More...

const Layer *	m_nextLayer
	BinUtility for next/previous decision. More...

const BinUtility *	m_binUtility
	Enclosing TrackingVolume. More...

const TrackingVolume *	m_enclosingTrackingVolume
	Enclosing DetachedTrackingVolume. More...

const DetachedTrackingVolume *	m_enclosingDetachedTrackingVolume

LayerIndex	m_index
	LayerIndex. More...

int	m_layerType
	active passive layer More...

double	m_ref
	reference measure for local coordinate convertors More...

Static Protected Attributes
static constexpr double	s_onSurfaceTolerance = 10e-5
	Tolerance for being on Surface. More...

Private Member Functions
const Surface &	approachSurface (const Amg::Vector3D &pos, const Amg::Vector3D &dir, const BoundaryCheck &bcheck) const
	Surface seen on approach - if not defined differently, it is the surfaceRepresentation() More...

void	buildApproachDescriptor ()
	build approach surfaces More...

Detailed Description

Class to describe a cylindrical detector layer for tracking, it inhertis from both, Layer base class and CylinderSurface class

Author: Andre.nosp@m.as.S.nosp@m.alzbu.nosp@m.rger.nosp@m.@cern.nosp@m..ch

Definition at line 40 of file CylinderLayer.h.

Member Typedef Documentation

◆ ChargedTrackParametersUniquePtr

using Trk::Surface::ChargedTrackParametersUniquePtr = std::unique_ptr<ParametersBase<5, Trk::Charged> >

inherited

Unique ptr types.

Definition at line 124 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

◆ NeutralTrackParametersUniquePtr

using Trk::Surface::NeutralTrackParametersUniquePtr = std::unique_ptr<ParametersBase<5, Trk::Neutral> >

inherited

Definition at line 126 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

Constructor & Destructor Documentation

◆ CylinderLayer() [1/10]

Trk::CylinderLayer::CylinderLayer ( )

inline

Default Constructor.

Definition at line 46 of file CylinderLayer.h.

46 {}

◆ CylinderLayer() [2/10]

Trk::CylinderLayer::CylinderLayer	(	const Amg::Transform3D &	transform,
		Trk::CylinderBounds *	cbounds,
		const LayerMaterialProperties &	laymatprop,
		double	thickness = `0.`,
		std::unique_ptr< OverlapDescriptor >	od = `nullptr`,
		int	laytyp = `int(Trk::active)`
	)

Constructor with CylinderSurface components and MaterialProperties.

Definition at line 22 of file CylinderLayer.cxx.

   : CylinderSurface(transform, cbounds)
   , Layer(laymatprop, thickness, std::move(olap), laytyp)
   , m_approachDescriptor(nullptr)
 {
   CylinderSurface::associateLayer(*this);
 }

◆ CylinderLayer() [3/10]

Trk::CylinderLayer::CylinderLayer	(	Trk::CylinderSurface *	cyl,
		const LayerMaterialProperties &	laymatprop,
		double	thickness = `0.`,
		std::unique_ptr< OverlapDescriptor >	od = `nullptr`,
		int	laytyp = `int(Trk::active)`
	)

Constructor with CylinderSurface and MaterialProperties.

Definition at line 36 of file CylinderLayer.cxx.

   : CylinderSurface(*cyl)
   , Layer(laymatprop, thickness, std::move(olap), laytyp)
   , m_approachDescriptor(nullptr)
 {
   CylinderSurface::associateLayer(*this);
 }

◆ CylinderLayer() [4/10]

Trk::CylinderLayer::CylinderLayer	(	const Amg::Transform3D &	transform,
		Trk::CylinderBounds *	cbounds,
		std::unique_ptr< SurfaceArray >	surfaceArray,
		double	thickness = `0.`,
		std::unique_ptr< OverlapDescriptor >	od = `nullptr`,
		Trk::IApproachDescriptor *	ades = `nullptr`,
		int	laytyp = `int(Trk::active)`
	)

Constructor with CylinderSurface components and pointer to SurfaceArray (passing ownership)

has optionally an OverlapDescriptor and ApproachDescriptor for the surfaceArray

Definition at line 49 of file CylinderLayer.cxx.

   : CylinderSurface(transform, cbounds)
   , Layer(std::move(surfaceArray), thickness, std::move(olap), laytyp)
   , m_approachDescriptor(ades)
 {
   CylinderSurface::associateLayer(*this);
   if (!ades && m_surfaceArray)
     buildApproachDescriptor();
   // register the layer
   if (ades)
     m_approachDescriptor->registerLayer(*this);
 }

◆ CylinderLayer() [5/10]

Trk::CylinderLayer::CylinderLayer	(	const Amg::Transform3D &	transform,
		Trk::CylinderBounds *	cbounds,
		std::unique_ptr< SurfaceArray >	surfaceArray,
		const LayerMaterialProperties &	laymatprop,
		double	thickness = `0.`,
		std::unique_ptr< OverlapDescriptor >	od = `nullptr`,
		Trk::IApproachDescriptor *	ades = `nullptr`,
		int	laytyp = `int(Trk::active)`
	)

Constructor with CylinderSurface components, MaterialProperties and pointer SurfaceArray (passing ownership)

Definition at line 68 of file CylinderLayer.cxx.

   : CylinderSurface(transform, cbounds)
   , Layer(std::move(surfaceArray), laymatprop, thickness, std::move(olap), laytyp)
   , m_approachDescriptor(ades)
 {
   CylinderSurface::associateLayer(*this);
   if (!ades && m_surfaceArray)
     buildApproachDescriptor();
   // register the layer
   if (ades)
     m_approachDescriptor->registerLayer(*this);
 }

◆ CylinderLayer() [6/10]

Trk::CylinderLayer::CylinderLayer	(	Trk::CylinderBounds *	cbounds,
		const LayerMaterialProperties &	laymatprop,
		double	thickness = `0.`,
		std::unique_ptr< OverlapDescriptor >	od = `nullptr`,
		int	laytyp = `int(Trk::active)`
	)

Concentric Layer: Constructor with CylinderSurface components and MaterialProperties.

Definition at line 89 of file CylinderLayer.cxx.

     : CylinderSurface(cbounds),
       Layer(laymatprop, thickness, std::move(olap), laytyp),
       m_approachDescriptor(nullptr) {
   CylinderSurface::associateLayer(*this);
 }

◆ CylinderLayer() [7/10]

Trk::CylinderLayer::CylinderLayer	(	Trk::CylinderBounds *	cbounds,
		std::unique_ptr< SurfaceArray >	surfaceArray,
		double	thickness = `0.`,
		std::unique_ptr< OverlapDescriptor >	od = `nullptr`,
		Trk::IApproachDescriptor *	ades = `nullptr`,
		int	laytyp = `int(Trk::active)`
	)

Concentric Layer: Constructor with CylinderSurface components and pointer to SurfaceArray (passing ownership)

Definition at line 99 of file CylinderLayer.cxx.

     : CylinderSurface(cbounds),
       Layer(std::move(surfaceArray), thickness, std::move(olap), laytyp),
       m_approachDescriptor(ades) {
   CylinderSurface::associateLayer(*this);
   if (!ades && m_surfaceArray) buildApproachDescriptor();
   // register the layer
   if (ades) m_approachDescriptor->registerLayer(*this);
 }

◆ CylinderLayer() [8/10]

Trk::CylinderLayer::CylinderLayer	(	Trk::CylinderBounds *	cbounds,
		std::unique_ptr< SurfaceArray >	surfaceArray,
		const LayerMaterialProperties &	laymatprop,
		double	thickness = `0.`,
		std::unique_ptr< OverlapDescriptor >	od = `nullptr`,
		Trk::IApproachDescriptor *	ades = `nullptr`,
		int	laytyp = `int(Trk::active)`
	)

Concentric Layer: Constructor with CylinderSurface components, MaterialProperties and pointer SurfaceArray (passing ownership)

Definition at line 113 of file CylinderLayer.cxx.

     : CylinderSurface(cbounds),
       Layer(std::move(surfaceArray), laymatprop, thickness, std::move(olap), laytyp),
       m_approachDescriptor(ades) {
   CylinderSurface::associateLayer(*this);
   if (!ades && m_surfaceArray) buildApproachDescriptor();
   // register the layer
   if (ades) m_approachDescriptor->registerLayer(*this);
 }

◆ CylinderLayer() [9/10]

Trk::CylinderLayer::CylinderLayer ( const CylinderLayer & cla )

Copy constructor of CylinderLayer.

< TODO use clone when exists

Definition at line 126 of file CylinderLayer.cxx.

     : CylinderSurface(clay), Layer(clay), m_approachDescriptor(nullptr) {
   CylinderSurface::associateLayer(*this);
   m_approachDescriptor.reset();
   if (m_surfaceArray)
     buildApproachDescriptor();  
 }

◆ CylinderLayer() [10/10]

Trk::CylinderLayer::CylinderLayer	(	const CylinderLayer &	cla,
		const Amg::Transform3D &	tr
	)

Copy constructor with shift.

Definition at line 134 of file CylinderLayer.cxx.

     : CylinderSurface(clay, transf),
       Layer(clay),
       m_approachDescriptor(nullptr) {
   if (m_surfaceArray) buildApproachDescriptor();
 }

◆ ~CylinderLayer()

virtual Trk::CylinderLayer::~CylinderLayer ( )

virtualdefault

Destructor.

Member Function Documentation

◆ approachSurface()

const Trk::Surface & Trk::CylinderLayer::approachSurface	(	const Amg::Vector3D &	pos,
		const Amg::Vector3D &	dir,
		const BoundaryCheck &	bcheck
	)		const

private

Surface seen on approach - if not defined differently, it is the surfaceRepresentation()

Definition at line 237 of file CylinderLayer.cxx.

                                           {
   if (m_approachDescriptor) {
     // get the test surfaces from the approach Descriptor
     const Trk::ApproachSurfaces* surfacesOnApproach =
         m_approachDescriptor->approachSurfaces(pos, dir);
     if (surfacesOnApproach) {
       // test the intersections and go
       std::vector<Trk::Intersection> sfIntersections;
       const Trk::Surface* aSurface = nullptr;
       double aPathLength = 10e10;
       // get the surfaces
       for (const auto& sfIter : (*surfacesOnApproach)) {
         // get the intersection with the surface
         Trk::Intersection sIntersection =
             sfIter->straightLineIntersection(pos, dir, true, bcheck);
         // validation
         if (sIntersection.valid && sIntersection.pathLength < aPathLength) {
           aPathLength = sIntersection.pathLength;
           aSurface = sfIter;
         }
       }
       if (aSurface) return (*aSurface);
     }
   }
   return surfaceRepresentation();
 }

◆ assignMaterialProperties()

void Trk::Layer::assignMaterialProperties	(	const LayerMaterialProperties &	prop,
		double	scale = `1.0`
	)

inherited

assignMaterialPropeties

Definition at line 191 of file Layer.cxx.

                                                         {
   m_layerMaterialProperties.reset(prop.clone());
   if (scale != 1.0) {
     (*(m_layerMaterialProperties)) *= scale;
   }
 }

◆ associatedDetectorElement()

const TrkDetElementBase* Trk::Surface::associatedDetectorElement ( ) const

inherited

return associated Detector Element

◆ associatedDetectorElementIdentifier()

Identifier Trk::Surface::associatedDetectorElementIdentifier ( ) const

inherited

return Identifier of the associated Detector Element

◆ associatedLayer()

const Trk::Layer* Trk::Surface::associatedLayer ( ) const

inherited

return the associated Layer

◆ associateLayer()

void Trk::Surface::associateLayer ( const Layer & lay )

inherited

method to associate the associated Trk::Layer which is alreay owned

only allowed by LayerBuilder
only done if no Layer is set already

◆ baseSurface()

virtual const Trk::Surface* Trk::Surface::baseSurface ( ) const

virtualinherited

return the base surface (simplified for persistification)

◆ binUtility()

const BinUtility* Trk::Layer::binUtility ( ) const

inherited

access the BinUtility

◆ bounds()

virtual const CylinderBounds& Trk::CylinderSurface::bounds ( ) const

finaloverridevirtualinherited

This method returns the CylinderBounds by reference (NoBounds is not possible for cylinder)

Implements Trk::Surface.

◆ buildApproachDescriptor()

void Trk::CylinderLayer::buildApproachDescriptor ( )

private

build approach surfaces

Definition at line 281 of file CylinderLayer.cxx.

                                                {
   // delete it
   // delete the surfaces
   auto aSurfaces = std::make_unique<Trk::ApproachSurfaces>();
   // create the new surfaces
   if (m_transforms) {
     Amg::Transform3D asTransform = m_transforms->transform;
     aSurfaces->push_back(new Trk::CylinderSurface(
       asTransform, m_bounds->r() - 0.5 * thickness(), m_bounds->halflengthZ()));
     aSurfaces->push_back(new Trk::CylinderSurface(
       asTransform, m_bounds->r() + 0.5 * thickness(), m_bounds->halflengthZ()));
   } else {
     aSurfaces->push_back(new Trk::CylinderSurface(
       m_bounds->r() - 0.5 * thickness(), m_bounds->halflengthZ()));
     aSurfaces->push_back(new Trk::CylinderSurface(
       m_bounds->r() + 0.5 * thickness(), m_bounds->halflengthZ()));
   }
   // set the layer and make TGOwn
   for (auto& sIter : (*aSurfaces)) {
     sIter->associateLayer(*this);
     sIter->setOwner(Trk::TGOwn);
   }
   m_approachDescriptor =
       std::make_unique<Trk::ApproachDescriptor>(std::move(aSurfaces));
 }

◆ cachedTransform()

const Amg::Transform3D* Trk::Surface::cachedTransform ( ) const

inherited

Return the cached transformation directly.

Don't try to make a new transform if it's not here.

◆ center()

const Amg::Vector3D& Trk::Surface::center ( ) const

inherited

Returns the center position of the Surface.

◆ clone()

virtual CylinderSurface* Trk::CylinderSurface::clone ( ) const

overridevirtualinherited

Implicit Constructor.

Implements Trk::Surface.

◆ compactify()

void Trk::Layer::compactify	(	size_t &	cSurfaces,
		size_t &	tSurfaces
	)

inherited

register layer index for material map registration

Definition at line 200 of file Layer.cxx.

 {
   // set the layer surface representation
   ++tSurfaces;
   if (surfaceRepresentation().owner() == Trk::noOwn) {
     surfaceRepresentation().setOwner(Trk::TGOwn);
     ++cSurfaces;
   }
   // set the subsurface representation, usually already owned by DetElement
   if (m_surfaceArray) {
     BinnedArraySpan<Trk::Surface * const> surfaces = m_surfaceArray->arrayObjects();
     for (Trk::Surface * const sIter : surfaces) {
       if (sIter && (*sIter).owner() == Trk::noOwn) {
         (*sIter).setOwner(Trk::TGOwn);
         ++cSurfaces;
       }
       ++tSurfaces;
     }
   }
 }

◆ compatibleSurfaces() [1/2]

size_t Trk::Layer::compatibleSurfaces	(	std::vector< SurfaceIntersection > &	cSurfaces,
		const NeutralParameters &	pars,
		PropDirection	pdir,
		const BoundaryCheck &	bcheck,
		bool	materialSurfacesOnly = `true`,
		const Surface *	startSurface = `nullptr`,
		const Surface *	endSurface = `nullptr`,
		const ICompatibilityEstimator *	ice = `nullptr`
	)		const

inherited

get compatible surfaces starting from neutral parameters

◆ compatibleSurfaces() [2/2]

size_t Trk::Layer::compatibleSurfaces	(	std::vector< SurfaceIntersection > &	cSurfaces,
		const TrackParameters &	pars,
		PropDirection	pdir,
		const BoundaryCheck &	bcheck,
		bool	materialSurfacesOnly = `true`,
		const Surface *	startSurface = `nullptr`,
		const Surface *	endSurface = `nullptr`,
		const ICompatibilityEstimator *	ice = `nullptr`
	)		const

inherited

get compatible surfaces starting from charged parameters

◆ createUniqueNeutralParameters() [1/2]

Trk::Surface::NeutralTrackParametersUniquePtr Trk::CylinderSurface::createUniqueNeutralParameters	(	const Amg::Vector3D &	position,
		const Amg::Vector3D &	momentum,
		double	charge,
		std::optional< AmgSymMatrix(5)>	cov = `std::nullopt`
	)		const

finaloverridevirtualinherited

Use the Surface as a ParametersBase constructor, from global parameters - neutral.

Implements Trk::Surface.

Definition at line 159 of file CylinderSurface.cxx.

                                             {
   return std::make_unique<ParametersT<5, Neutral, CylinderSurface>>(
       position, momentum, charge, *this, std::move(cov));
 }

◆ createUniqueNeutralParameters() [2/2]

Trk::Surface::NeutralTrackParametersUniquePtr Trk::CylinderSurface::createUniqueNeutralParameters	(	double	l1,
		double	l2,
		double	phi,
		double	theta,
		double	qop,
		std::optional< AmgSymMatrix(5)>	cov = `std::nullopt`
	)		const

finaloverridevirtualinherited

Use the Surface as a ParametersBase constructor, from local parameters - neutral.

Implements Trk::Surface.

Definition at line 149 of file CylinderSurface.cxx.

                                             {
   return std::make_unique<ParametersT<5, Neutral, CylinderSurface>>(
       l1, l2, phi, theta, qop, *this, std::move(cov));
 }

◆ createUniqueParameters() [1/2]

template<int DIM, class T >

std::unique_ptr<ParametersT<DIM, T, CylinderSurface> > Trk::CylinderSurface::createUniqueParameters	(	const Amg::Vector3D &	position,
		const Amg::Vector3D &	momentum,
		double	charge,
		std::optional< AmgSymMatrix(DIM)>	cov = `std::nullopt`
	)		const

inherited

Use the Surface as a ParametersBase constructor, from global parameters.

◆ createUniqueParameters() [2/2]

template<int DIM, class T >

std::unique_ptr<ParametersT<DIM, T, CylinderSurface> > Trk::CylinderSurface::createUniqueParameters	(	double	l1,
		double	l2,
		double	phi,
		double	theta,
		double	qop,
		std::optional< AmgSymMatrix(DIM)>	cov = `std::nullopt`
	)		const

inherited

Use the Surface as a ParametersBase constructor, from local parameters.

◆ createUniqueTrackParameters() [1/2]

Trk::Surface::ChargedTrackParametersUniquePtr Trk::CylinderSurface::createUniqueTrackParameters	(	const Amg::Vector3D &	position,
		const Amg::Vector3D &	momentum,
		double	charge,
		std::optional< AmgSymMatrix(5)>	cov = `std::nullopt`
	)		const

finaloverridevirtualinherited

Use the Surface as a ParametersBase constructor, from global parameters - charged.

Implements Trk::Surface.

Definition at line 139 of file CylinderSurface.cxx.

                                             {
   return std::make_unique<ParametersT<5, Charged, CylinderSurface>>(
       position, momentum, charge, *this, std::move(cov));
 }

◆ createUniqueTrackParameters() [2/2]

Trk::Surface::ChargedTrackParametersUniquePtr Trk::CylinderSurface::createUniqueTrackParameters	(	double	l1,
		double	l2,
		double	phi,
		double	theta,
		double	qop,
		std::optional< AmgSymMatrix(5)>	cov = `std::nullopt`
	)		const

finaloverridevirtualinherited

Use the Surface as a ParametersBase constructor, from local parameters - charged.

Implements Trk::Surface.

Definition at line 129 of file CylinderSurface.cxx.

                                             {
   return std::make_unique<ParametersT<5, Charged, CylinderSurface>>(
       l1, l2, phi, theta, qop, *this, std::move(cov));
 }

◆ dump() [1/2]

MsgStream & Surface::dump ( MsgStream & sl ) const

virtualinherited

Output Method for MsgStream, to be overloaded by child classes.

Reimplemented in Trk::PerigeeSurface, and Trk::DetElementSurface.

Definition at line 157 of file Surface.cxx.

 {
   sl << std::setiosflags(std::ios::fixed);
   sl << std::setprecision(4);
   sl << name() << std::endl;
   if (associatedDetectorElement()!=nullptr){
     sl<<"     Detector Type              = "<<associatedDetectorElement()->detectorTypeString()<<std::endl;
   }
   sl << "     Center position  (x, y, z) = (" << center().x() << ", " << center().y() << ", " << center().z() << ")"
      << std::endl;
   Amg::RotationMatrix3D rot(transform().linear());
   Amg::Vector3D rotX(rot.col(0));
   Amg::Vector3D rotY(rot.col(1));
   Amg::Vector3D rotZ(rot.col(2));
   sl << std::setprecision(6);
   sl << "     Rotation:             colX = (" << rotX(0) << ", " << rotX(1) << ", " << rotX(2) << ")" << std::endl;
   sl << "                           colY = (" << rotY(0) << ", " << rotY(1) << ", " << rotY(2) << ")" << std::endl;
   sl << "                           colZ = (" << rotZ(0) << ", " << rotZ(1) << ", " << rotZ(2) << ")" << std::endl;
   sl << "     Bounds  : " << bounds();
   if (!checkTransform(*this)) {
      sl << std::endl << "     NOT a strict rotation matrix." << std::endl;
   }
   sl << std::setprecision(-1);
   return sl;
 }

◆ dump() [2/2]

std::ostream & Surface::dump ( std::ostream & sl ) const

virtualinherited

Output Method for std::ostream, to be overloaded by child classes.

Reimplemented in Trk::PerigeeSurface, and Trk::DetElementSurface.

Definition at line 185 of file Surface.cxx.

 {
   sl << std::setiosflags(std::ios::fixed);
   sl << std::setprecision(4);
   sl << name() << std::endl;
   if (associatedDetectorElement()!=nullptr){
     sl<<"     Detector Type              = "<<associatedDetectorElement()->detectorTypeString()<<std::endl;
   }
   sl << "     Center position  (x, y, z) = (" << center().x() << ", " << center().y() << ", " << center().z() << ")"
      << std::endl;
   Amg::RotationMatrix3D rot(transform().linear());
   Amg::Vector3D rotX(rot.col(0));
   Amg::Vector3D rotY(rot.col(1));
   Amg::Vector3D rotZ(rot.col(2));
   sl << std::setprecision(6);
   sl << "     Rotation:             colX = (" << rotX(0) << ", " << rotX(1) << ", " << rotX(2) << ")" << std::endl;
   sl << "                           colY = (" << rotY(0) << ", " << rotY(1) << ", " << rotY(2) << ")" << std::endl;
   sl << "                           colZ = (" << rotZ(0) << ", " << rotZ(1) << ", " << rotZ(2) << ")" << std::endl;
   sl << "     Bounds  : " << bounds();
   if (!checkTransform(*this)) {
      sl << std::endl << "     NOT a strict rotation matrix." << std::endl;
   }
   sl << std::setprecision(-1);
   return sl;
 }

◆ encloseDetachedTrackingVolume()

void Trk::Layer::encloseDetachedTrackingVolume ( const DetachedTrackingVolume & tvol )

inherited

◆ encloseTrackingVolume()

void Trk::Layer::encloseTrackingVolume ( const TrackingVolume & tvol )

inherited

private method to set the enclosed detached TV

◆ enclosingDetachedTrackingVolume()

const DetachedTrackingVolume* Trk::Layer::enclosingDetachedTrackingVolume ( ) const

inherited

get the confining DetachedTrackingVolume

◆ enclosingTrackingVolume()

const TrackingVolume* Trk::Layer::enclosingTrackingVolume ( ) const

inherited

get the confining TrackingVolume

◆ fullUpdateMaterialProperties()

const Trk::MaterialProperties * Trk::Layer::fullUpdateMaterialProperties ( const TrackParameters & par ) const

inherited

getting the MaterialProperties back - for full update

Definition at line 183 of file Layer.cxx.

                                           {
   if (m_layerMaterialProperties){
     return m_layerMaterialProperties->fullMaterial(parm.position());
   }
   return nullptr;
 }

◆ getCompatibleSurfaces()

template<class T >

size_t Trk::Layer::getCompatibleSurfaces	(	std::vector< SurfaceIntersection > &	cSurfaces,
		const T &	pars,
		PropDirection	pdir,
		const BoundaryCheck &	bcheck,
		bool	materialSurfacesOnly = `true`,
		const Surface *	startSurface = `nullptr`,
		const Surface *	endSurface = `nullptr`,
		const ICompatibilityEstimator *	ice = `nullptr`
	)		const

inherited

get compatible surfaces starting from charged parameters

propagate TrackingGeometry owner downstream

◆ getRef()

double Trk::Layer::getRef ( ) const

inherited

get the reference measure

◆ globalReferencePoint()

const Amg::Vector3D & Trk::CylinderSurface::globalReferencePoint ( ) const

finaloverridevirtualinherited

Returns a global reference point: For the Cylinder this is \( (R*cos(\phi), R*sin(\phi),0)*transform() \) Where \( \phi \) denotes the averagePhi() of the cylinderBounds.

Reimplemented from Trk::Surface.

Definition at line 167 of file CylinderSurface.cxx.

 {
   if (!m_referencePoint) {
     double rMedium = bounds().r();
     double phi = bounds().averagePhi();
     Amg::Vector3D gp(rMedium * cos(phi), rMedium * sin(phi), 0.);
     m_referencePoint.set(std::make_unique<Amg::Vector3D>(transform() * gp));
   }
   return (*m_referencePoint);
 }

◆ globalToLocal() [1/3]

bool Trk::CylinderSurface::globalToLocal	(	const Amg::Vector3D &	glob,
		const Amg::Vector3D &	mom,
		Amg::Vector2D &	loc
	)		const

overridevirtualinherited

Specialized for CylinderSurface : GlobalToLocal method without dynamic memory allocation - boolean checks if on surface.

Implements Trk::Surface.

Reimplemented in Trk::SlidingCylinderSurface.

Definition at line 238 of file CylinderSurface.cxx.

 {
   // get the transform & transform global position into cylinder frame
   // transform it to the globalframe: CylinderSurfaces are allowed to have 0
   // pointer transform
   double radius = 0.;
   double inttol = bounds().r() * 0.0001;
   if (inttol < 0.01)
     inttol = 0.01;
   // do the transformation or not
   if (Trk::Surface::m_transforms) {
     Amg::Vector3D loc3Dframe(inverseTransformMultHelper(glopos));
     locpos = Amg::Vector2D(bounds().r() * loc3Dframe.phi(), loc3Dframe.z());
     radius = loc3Dframe.perp();
   } else {
     locpos = Amg::Vector2D(bounds().r() * glopos.phi(), glopos.z());
     radius = glopos.perp();
   }
   // return true or false
   return (fabs(radius - bounds().r()) <= inttol);
 }

◆ globalToLocal() [2/3]

std::optional<Amg::Vector2D> Trk::Surface::globalToLocal	(	const Amg::Vector3D &	glopos,
		const Amg::Vector3D &	glomom
	)		const

inherited

This method returns the LocalPosition from a provided GlobalPosition.

If the GlobalPosition is not on the Surface, it returns a nullopt This method does not check if the calculated LocalPosition is inside surface bounds. If this check is needed, use positionOnSurface - for generality with momentum

◆ globalToLocal() [3/3]

std::optional<Amg::Vector2D> Trk::Surface::globalToLocal	(	const Amg::Vector3D &	glopos,
		double	tol = `0.`
	)		const

inherited

This method returns the LocalPosition from a provided GlobalPosition.

If the GlobalPosition is not on the Surface, it returns nullopt This method does not check if the calculated LocalPosition is inside surface bounds. If this check is needed, use positionOnSurface - only for planar, cylinder surface fully defined

◆ hasBounds()

bool Trk::CylinderSurface::hasBounds ( ) const

inherited

◆ hasSubStructure()

bool Trk::Layer::hasSubStructure ( bool resolveSensitive = false ) const

inlinevirtualinherited

Has sub-structure method:

sub-structure depending on : (a) only when required to resolve sub surfaces for sensitive hits (b) also material is ordered with sub structure

Definition at line 252 of file Layer.cxx.

                                                                  {
   return resolveSensitive && m_surfaceArray;
 }

◆ insideBounds()

virtual bool Trk::CylinderSurface::insideBounds	(	const Amg::Vector2D &	locpos,
		double	tol1 = `0.`,
		double	tol2 = `0.`
	)		const

overridevirtualinherited

This method calls the inside method of CylinderBounds.

Implements Trk::Surface.

Reimplemented in Trk::SubtractedCylinderSurface.

◆ insideBoundsCheck()

virtual bool Trk::CylinderSurface::insideBoundsCheck	(	const Amg::Vector2D &	locpos,
		const BoundaryCheck &	bchk
	)		const

finaloverridevirtualinherited

Implements Trk::Surface.

◆ inverseTransformHelper()

Amg::Transform3D Trk::Surface::inverseTransformHelper ( ) const

protectedinherited

Helper method to factorize in one place common operations calculate inverse transofrm and multiply with position.

◆ inverseTransformMultHelper()

Amg::Vector3D Trk::Surface::inverseTransformMultHelper ( const Amg::Vector3D & glopos ) const

protectedinherited

◆ isActive()

bool Trk::Surface::isActive ( ) const

inherited

Return 'true' if this surface is owned by the detector element.

◆ isFree()

bool Trk::Surface::isFree ( ) const

inherited

Returns 'true' if this surface is 'free', i.e.

it does not belong to a detector element (and returns false otherwise

◆ isOnLayer()

bool Trk::Layer::isOnLayer	(	const Amg::Vector3D &	gp,
		const BoundaryCheck &	bcheck = `BoundaryCheck(true)`
	)		const

virtualinherited

isOnLayer() method, using isOnSurface() with Layer specific tolerance

Reimplemented in Trk::NavigationLayer, and Trk::MaterialLayer.

Definition at line 149 of file Layer.cxx.

                                                             {
   return (surfaceRepresentation())
       .isOnSurface(gp, bchk, 0.5 * m_layerThickness);
 }

◆ isOnSurface()

bool Trk::CylinderSurface::isOnSurface	(	const Amg::Vector3D &	glopo,
		const BoundaryCheck &	bchk = `true`,
		double	tol1 = `0.`,
		double	tol2 = `0.`
	)		const

overridevirtualinherited

This method returns true if the GlobalPosition is on the Surface for both, within or without check of whether the local position is inside boundaries or not.

Reimplemented from Trk::Surface.

Reimplemented in Trk::SlidingCylinderSurface.

Definition at line 263 of file CylinderSurface.cxx.

 {
   Amg::Vector3D loc3Dframe =
     Trk::Surface::m_transforms ? inverseTransformMultHelper(glopo) : glopo;
   return (bchk ? bounds().inside3D(loc3Dframe,
                                    tol1 + s_onSurfaceTolerance,
                                    tol2 + s_onSurfaceTolerance)
                : true);
 }

◆ layerIndex()

const LayerIndex& Trk::Layer::layerIndex ( ) const

inherited

get the layerIndex

◆ layerMaterialProperties()

const LayerMaterialProperties* Trk::Layer::layerMaterialProperties ( ) const

inherited

getting the LayerMaterialProperties including full/pre/post update

◆ layerType()

int Trk::Layer::layerType ( ) const

inherited

get the Layer coding

◆ localParametersToPosition()

virtual Amg::Vector2D Trk::CylinderSurface::localParametersToPosition ( const LocalParameters & locpars ) const

finaloverridevirtualinherited

Specialized for CylinderSurface : LocalParameters to Vector2D.

Reimplemented from Trk::Surface.

◆ localToGlobal() [1/5]

ATH_FLATTEN void Trk::CylinderSurface::localToGlobal	(	const Amg::Vector2D &	locp,
		const Amg::Vector3D &	mom,
		Amg::Vector3D &	glob
	)		const

overridevirtualinherited

Specialized for CylinderSurface : LocalToGlobal method without dynamic memory allocation.

Implements Trk::Surface.

Reimplemented in Trk::SlidingCylinderSurface.

Definition at line 223 of file CylinderSurface.cxx.

 {
   // create the position in the local 3d frame
   double r = bounds().r();
   double phi = locpos[Trk::locRPhi] / r;
   glopos = Amg::Vector3D(r * cos(phi), r * sin(phi), locpos[Trk::locZ]);
   // transform it to the globalframe: CylinderSurfaces are allowed to have 0
   // pointer transform
   if (Trk::Surface::m_transforms)
     glopos = transform() * glopos;
 }

◆ localToGlobal() [2/5]

Amg::Vector3D Trk::Surface::localToGlobal ( const Amg::Vector2D & locpos ) const

inherited

This method returns the GlobalPosition from a LocalPosition uses the per surface localToGlobal.

◆ localToGlobal() [3/5]

Amg::Vector3D Trk::Surface::localToGlobal	(	const Amg::Vector2D &	locpos,
		const Amg::Vector3D &	glomom
	)		const

inherited

This method returns the GlobalPosition from a LocalPosition The LocalPosition can be outside Surface bounds - for generality with momentum.

◆ localToGlobal() [4/5]

Amg::Vector3D Trk::Surface::localToGlobal ( const LocalParameters & locpars ) const

inherited

This method returns the GlobalPosition from LocalParameters The LocalParameters can be outside Surface bounds.

◆ localToGlobal() [5/5]

Amg::Vector3D Trk::Surface::localToGlobal	(	const LocalParameters &	locpars,
		const Amg::Vector3D &	glomom
	)		const

inherited

This method returns the GlobalPosition from LocalParameters The LocalParameters can be outside Surface bounds - for generality with momentum.

◆ materialLayer() [1/2]

Trk::Layer* Trk::Surface::materialLayer ( )

inherited

◆ materialLayer() [2/2]

const Trk::Layer* Trk::Surface::materialLayer ( ) const

inherited

return the material Layer

◆ measurementFrame()

Amg::RotationMatrix3D Trk::CylinderSurface::measurementFrame	(	const Amg::Vector3D &	glopos,
		const Amg::Vector3D &	glomom
	)		const

finaloverridevirtualinherited

Return the measurement frame - this is needed for alignment, in particular for StraightLine and Perigee Surface.

the default implementation is the the RotationMatrix3D of the transform

Reimplemented from Trk::Surface.

Definition at line 190 of file CylinderSurface.cxx.

 {
   Amg::RotationMatrix3D mFrame;
   // construct the measurement frame
   Amg::Vector3D measY(
     transform().rotation().col(2)); // measured Y is the z axis
   Amg::Vector3D measDepth =
     Amg::Vector3D(pos.x(), pos.y(), 0.)
       .unit(); // measured z is the position transverse normalized
   Amg::Vector3D measX(
     measY.cross(measDepth).unit()); // measured X is what comoes out of it
   // the columnes
   mFrame.col(0) = measX;
   mFrame.col(1) = measY;
   mFrame.col(2) = measDepth;
   // return the rotation matrix
   return mFrame;
 }

◆ moveLayer()

void Trk::CylinderLayer::moveLayer ( Amg::Transform3D & shift )

finaloverridevirtual

move the Layer

Implements Trk::Layer.

Definition at line 186 of file CylinderLayer.cxx.

                                                       {
   Trk::CylinderSurface::m_transforms =
       std::make_unique<Transforms>(shift * (m_transforms->transform));
  
   if (m_approachDescriptor && m_approachDescriptor->rebuild()) {
     // build the new approach descriptor - deletes the current one
     buildApproachDescriptor();
   }
 }

◆ name()

virtual std::string Trk::CylinderSurface::name ( ) const

overridevirtualinherited

Return properly formatted class name for screen output.

Implements Trk::Surface.

Reimplemented in Trk::SubtractedCylinderSurface, and Trk::SlidingCylinderSurface.

◆ nextLayer() [1/2]

const Trk::Layer * Trk::Layer::nextLayer ( bool skipNavLayer = false ) const

inherited

getting what's stored to be the next Layer, boolean to skip navigation layers

Definition at line 165 of file Layer.cxx.

                                                            {
   if (!skipNavLayer) return m_nextLayer;
   const Trk::Layer* nextMatLayer = m_nextLayer;
   // get the next Material layer
   while (nextMatLayer && !nextMatLayer->layerMaterialProperties() &&
          !nextMatLayer->surfaceArray())
     nextMatLayer = nextMatLayer->nextLayer();
   return nextMatLayer;
 }

◆ nextLayer() [2/2]

const Trk::Layer * Trk::Layer::nextLayer	(	const Amg::Vector3D &	gp,
		const Amg::Vector3D &	udir
	)		const

inherited

getting the next/previous Layer if registered - unit for direction vector required

Definition at line 175 of file Layer.cxx.

                                                                       {
   // no binutility -> no chance to find out the direction
   if (!m_binUtility) return nullptr;
   return (m_binUtility->orderDirection(gp, mom) == Trk::next) ? m_nextLayer
                                                               : m_previousLayer;
 }

◆ normal() [1/2]

virtual const Amg::Vector3D& Trk::Surface::normal ( ) const

virtualinherited

Returns the normal vector of the Surface (i.e.

in generall z-axis of rotation)

Reimplemented in Trk::PerigeeSurface.

◆ normal() [2/2]

virtual Amg::Vector3D Trk::CylinderSurface::normal ( const Amg::Vector2D & locpo ) const

finaloverridevirtualinherited

Return method for surface normal information at a given local point, overwrites the normal() from base class.

Reimplemented from Trk::Surface.

◆ numberOfInstantiations()

static std::size_t Trk::ObjectCounter< Trk::Surface >::numberOfInstantiations ( )

inlinestaticinherited

Definition at line 25 of file TrkObjectCounter.h.

   {
 #ifndef NDEBUG
     return s_numberOfInstantiations.load();
 #endif
     return 0;
   }

◆ onLayer()

template<class T >

bool Trk::Layer::onLayer	(	const T &	parameters,
		const BoundaryCheck &	bcheck = `BoundaryCheck(true)`
	)		const

inherited

templated on Layer method

◆ onSurface()

template<class T >

bool Trk::Surface::onSurface	(	const T &	parameters,
		const BoundaryCheck &	bchk = `BoundaryCheck(true)`
	)		const

inherited

The templated Parameters OnSurface method - checks on surface pointer first.

◆ operator!=()

bool Trk::Surface::operator!= ( const Surface & sf ) const

inherited

Non-equality operator.

◆ operator=()

Trk::CylinderLayer & Trk::CylinderLayer::operator= ( const CylinderLayer & clay )

Assignment operator for CylinderLayers.

Definition at line 142 of file CylinderLayer.cxx.

                                                                        {
   if (this != &clay) {
     // call the assignments of the base classes
     Trk::CylinderSurface::operator=(clay);
     Trk::Layer::operator=(clay);
     if (m_surfaceArray) buildApproachDescriptor();
     CylinderSurface::associateLayer(*this);
   }
   return (*this);
 }

◆ operator==() [1/2]

bool Trk::CylinderSurface::operator== ( const CylinderSurface & cf ) const

inherited

◆ operator==() [2/2]

bool Trk::CylinderSurface::operator== ( const Surface & sf ) const

overridevirtualinherited

Equality operator.

Implements Trk::Surface.

Reimplemented in Trk::SubtractedCylinderSurface, and Trk::SlidingCylinderSurface.

Definition at line 179 of file CylinderSurface.cxx.

 {
   // first check the type not to compare apples with oranges
   if (sf.type()!=Trk::SurfaceType::Cylinder){
       return false;
   }
   return (*this) == static_cast<const Trk::CylinderSurface&>(sf);
 }

◆ overlapDescriptor()

const OverlapDescriptor* Trk::Layer::overlapDescriptor ( ) const

inherited

gettint hte overlap descriptor

◆ owner()

SurfaceOwner Trk::Surface::owner ( ) const

inherited

return ownership

◆ pathCorrection()

virtual double Trk::CylinderSurface::pathCorrection	(	const Amg::Vector3D &	pos,
		const Amg::Vector3D &	mom
	)		const

overridevirtualinherited

the pathCorrection for derived classes with thickness

Reimplemented from Trk::Surface.

◆ positionOnSurface()

std::optional< Amg::Vector2D > Surface::positionOnSurface	(	const Amg::Vector3D &	glopo,
		const BoundaryCheck &	bchk = `true`,
		double	tol1 = `0.`,
		double	tol2 = `0.`
	)		const

inherited

positionOnSurface() returns the LocalPosition on the Surface,
If BoundaryCheck==false it just returns the value of globalToLocal (including nullptr possibility), if BoundaryCheck==true it checks whether the point is inside bounds or not (returns std::nullopt in this case).

Definition at line 106 of file Surface.cxx.

 {
   std::optional<Amg::Vector2D> posOnSurface = globalToLocal(glopo, tol1);
   if (!bchk){
     return posOnSurface;
   }
   if (posOnSurface && insideBounds(*posOnSurface, tol1, tol2)){
     return posOnSurface;
   }
   return std::nullopt;
 }

◆ postUpdateMaterialFactor()

double Trk::CylinderLayer::postUpdateMaterialFactor	(	const Trk::TrackParameters &	par,
		Trk::PropDirection	dir
	)		const

finaloverridevirtual

getting the MaterialProperties back - for post-update

Reimplemented from Trk::Layer.

Definition at line 176 of file CylinderLayer.cxx.

                                                                 {
   if (!Trk::Layer::m_layerMaterialProperties) return 0;
   const Amg::Vector3D& parmPos = parm.position();
   Amg::Vector3D pastStep(parmPos + dir * parm.momentum().normalized());
   if (pastStep.perp() > parm.position().perp())
     return Trk::Layer::m_layerMaterialProperties->alongPostFactor();
   return Trk::Layer::m_layerMaterialProperties->oppositePostFactor();
 }

◆ preUpdateMaterialFactor()

double Trk::CylinderLayer::preUpdateMaterialFactor	(	const Trk::TrackParameters &	par,
		Trk::PropDirection	dir
	)		const

finaloverridevirtual

getting the MaterialProperties back - for pre-update

Reimplemented from Trk::Layer.

Definition at line 165 of file CylinderLayer.cxx.

                                                                 {
   if (!Trk::Layer::m_layerMaterialProperties) return 0.;
   // calculate the direction to the normal
   const Amg::Vector3D& parmPos = parm.position();
   Amg::Vector3D pastStep(parmPos + dir * parm.momentum().normalized());
   if (pastStep.perp() > parm.position().perp())
     return Trk::Layer::m_layerMaterialProperties->alongPreFactor();
   return Trk::Layer::m_layerMaterialProperties->oppositePreFactor();
 }

◆ previousLayer()

const Trk::Layer * Trk::Layer::previousLayer ( bool skipNavLayer = false ) const

inherited

getting what's stored to be the previous Layer, boolean to skip navigation layers

Definition at line 155 of file Layer.cxx.

                                                                {
   if (!skipNavLayer) return m_previousLayer;
   const Trk::Layer* prevMatLayer = m_previousLayer;
   // get the previoys Material layer
   while (prevMatLayer && !prevMatLayer->layerMaterialProperties() &&
          !prevMatLayer->surfaceArray())
     prevMatLayer = prevMatLayer->previousLayer();
   return prevMatLayer;
 }

◆ registerLayerIndex()

void Trk::Layer::registerLayerIndex ( const LayerIndex & lIdx )

inherited

◆ registerRepresentingVolume()

void Trk::Layer::registerRepresentingVolume ( const Volume * theVol )

inherited

register Volume associated to the layer

◆ representingVolume()

const Volume* Trk::Layer::representingVolume ( ) const

inherited

get the Volume associated to the layer

◆ resizeAndRepositionLayer()

void Trk::CylinderLayer::resizeAndRepositionLayer	(	const VolumeBounds &	vBounds,
		const Amg::Vector3D &	cCenter,
		double	envelope
	)

finaloverridevirtual

Resize the layer to the tracking volume.

Implements Trk::Layer.

Definition at line 307 of file CylinderLayer.cxx.

                                                                    {
   // resize first of all
   resizeLayer(vBounds, envelope);
   // now reposition to the potentially center if necessary, do not change layers
   // with no transform
   if (Trk::CylinderSurface::m_transforms || center().isApprox(vCenter)) {
     return;
   }
   Amg::Transform3D transf(vCenter);
   Trk::CylinderSurface::m_transforms =
       std::make_unique<Transforms>(Amg::Transform3D(vCenter), vCenter);
   // rebuild approaching layers if needed
   if (m_approachDescriptor && m_approachDescriptor->rebuild())
     buildApproachDescriptor();
 }

◆ resizeLayer()

void Trk::CylinderLayer::resizeLayer	(	const VolumeBounds &	vBounds,
		double	envelope
	)

finaloverridevirtual

Resize the layer to the tracking volume - only works for CylinderVolumeBouns.

Implements Trk::Layer.

Definition at line 196 of file CylinderLayer.cxx.

                                                       {
   // only do this if the volume bounds a CylinderVolumeBounds
   const Trk::CylinderVolumeBounds* cvb =
       dynamic_cast<const Trk::CylinderVolumeBounds*>(&bounds);
   if (cvb) {
     // get the dimensions
     double hLengthZ = cvb->halflengthZ();
     double r = surfaceRepresentation().bounds().r();
     // (0) first, resize the layer itself
     Trk::CylinderBounds* rCylinderBounds =
         new Trk::CylinderBounds(r, hLengthZ - envelope);
     Trk::CylinderSurface::m_bounds =
         Trk::SharedObject<const Trk::CylinderBounds>(rCylinderBounds);
     // (1) resize the material properties by updating the BinUtility, assuming
     // rphi/z binning
     if (Trk::Layer::m_layerMaterialProperties) {
       const BinUtility* layerMaterialBU =
           Trk::Layer::m_layerMaterialProperties->binUtility();
       if (layerMaterialBU && layerMaterialBU->dimensions() > 1) {
         size_t binsRPhi = layerMaterialBU->max(0) + 1;
         size_t binsZ = layerMaterialBU->max(1) + 1;
         // create a new binning with the new dimensions
         Trk::BinUtility* rBinUtility = new Trk::BinUtility(
             binsRPhi, -r * M_PI, r * M_PI, Trk::closed, Trk::binRPhi);
         (*rBinUtility) +=
             Trk::BinUtility(binsZ, -hLengthZ + envelope, hLengthZ - envelope,
                             Trk::open, Trk::binZ);
         Trk::Layer::m_layerMaterialProperties->updateBinning(rBinUtility);
       }
     }
   }
  
   if (m_approachDescriptor && m_approachDescriptor->rebuild()) {
     // build the approach descriptor - delete the current approach descriptor
     buildApproachDescriptor();
   }
 }

◆ rotSymmetryAxis()

const Amg::Vector3D & Trk::CylinderSurface::rotSymmetryAxis ( ) const

virtualinherited

Return method for the rotational symmetry axis - the z-Axis of the HepTransform.

Definition at line 211 of file CylinderSurface.cxx.

 {
   if (!m_rotSymmetryAxis) {
     Amg::Vector3D zAxis(transform().rotation().col(2));
     m_rotSymmetryAxis.set(std::make_unique<Amg::Vector3D>(zAxis));
   }
   return (*m_rotSymmetryAxis);
 }

◆ setBinUtility()

void Trk::Layer::setBinUtility ( const BinUtility * )

inherited

set the BinUtility

◆ setLayerType()

void Trk::Layer::setLayerType ( int identifier )

inherited

set the Layer coding

◆ setMaterialLayer()

void Trk::Surface::setMaterialLayer ( std::shared_ptr< Layer > mlay )

inherited

set material layer

◆ setNextLayer()

void Trk::Layer::setNextLayer ( const Layer * )

inherited

set the next Layer

◆ setOwner()

void Trk::Surface::setOwner ( SurfaceOwner x )

inherited

set Ownership

◆ setPreviousLayer()

void Trk::Layer::setPreviousLayer ( const Layer * )

inherited

set the previous Layer

◆ setRef()

void Trk::Layer::setRef ( double )

inherited

set the reference measure

◆ setTransform()

void Trk::Surface::setTransform ( const Amg::Transform3D & trans )

inherited

Set the transform updates center and normal.

◆ straightLineDistanceEstimate() [1/2]

ATH_FLATTEN Trk::DistanceSolution Trk::CylinderSurface::straightLineDistanceEstimate	(	const Amg::Vector3D &	pos,
		const Amg::Vector3D &	dir
	)		const

overridevirtualinherited

fast distance to Surface

distance to surface

Implements Trk::Surface.

Reimplemented in Trk::SlidingCylinderSurface.

Definition at line 377 of file CylinderSurface.cxx.

 {
   double tol = 0.001;
  
   const Amg::Vector3D& X = center(); // point
   const Amg::Vector3D& S = normal(); // vector
  
   double radius = bounds().r();
   double sp = pos.dot(S);
   double sc = X.dot(S);
   double dp = dir.dot(S);
   Amg::Vector3D dx = X - pos - (sc - sp) * S; // vector
   Amg::Vector3D ax = dir - dp * S;            // vector
  
   double A = ax.dot(ax); // size of projected direction (squared)
   double B = ax.dot(dx); // dot product (->cos angle)
   double C = dx.dot(dx); // distance to axis (squared)
   double currDist = radius - sqrt(C);
  
   if (A == 0.) { // direction parallel to cylinder axis
     if (fabs(currDist) < tol) {
       return {1, 0., true, 0.}; // solution at surface
     }
     return {
       0, currDist, true, 0.}; // point of closest approach without intersection
   }
  
   // minimal distance to cylinder axis
   // The [[maybe_unused]] declaration is to suppress redundant division checking
   // here. Even a tiny change in rmin (~1e-13) can cause huge changes in the
   // reconstructed output, so don't change how it's evaluated.
   [[maybe_unused]] const double rmin_tmp = B * B / A;
   const double rmin2 = C - rmin_tmp;
   const double rmin = rmin2 < 0 ? 0 : sqrt(rmin2);
  
   if (rmin > radius) { // no intersection
     double first = B / A;
     return {
       0,
       currDist,
       true,
       first}; // point of closest approach without intersection
   }
   if (fabs(rmin - radius) <
       tol) { // tangential 'intersection' - return double solution
     double first = B / A;
     return {2, currDist, true, first, first};
   }
   // The [[maybe_unused]] declaration here suppresses redundant division
   // checking. We don't want to rewrite how this is evaluated due to
   // instabilities.
   [[maybe_unused]] const double b_a = B / A;
   const double x = sqrt((radius - rmin) * (radius + rmin) / A);
   double first = b_a - x;
   double second = b_a + x;
   if (first >= 0.) {
     return {2, currDist, true, first, second};
   }
   if (second <= 0.) {
     return {2, currDist, true, second, first};
   } // inside cylinder
   return {2, currDist, true, second, first};
 }

◆ straightLineDistanceEstimate() [2/2]

Trk::DistanceSolution Trk::CylinderSurface::straightLineDistanceEstimate	(	const Amg::Vector3D &	pos,
		const Amg::Vector3D &	dir,
		bool	bound
	)		const

overridevirtualinherited

fast distance to Surface - with bounds directive

Implements Trk::Surface.

Reimplemented in Trk::SlidingCylinderSurface.

Definition at line 444 of file CylinderSurface.cxx.

 {
   const double tolb = .01;
  
   const Amg::Transform3D& T = transform();
   //  double                Ax[3] = {T.xx(),T.yx(),T.zx()};
   //  double                Ay[3] = {T.xy(),T.yy(),T.zy()};
   //  double                Az[3] = {T.xz(),T.yz(),T.zz()};
  
   // Transformation to cylinder system coordinates
   //
  
   // BEGIN here is what i guess this might mean: BEGIN
   Amg::Vector3D Ax = T.rotation().col(0);
   Amg::Vector3D Ay = T.rotation().col(1);
   Amg::Vector3D Az = T.rotation().col(2);
  
   Amg::Vector3D dxyz = pos - T.translation();
   double x = dxyz.dot(Ax);
   double y = dxyz.dot(Ay);
   double z = dxyz.dot(Az);
   double ax = dir.dot(Ax);
   double ay = dir.dot(Ay);
   double at = ax * ax + ay * ay;
   double r = sqrt(x * x + y * y);
   double R = bounds().r();
  
   // END here is what i guessed this means END
  
   //  double dx  = pos[0]-T.dx()                         ;
   //  double dy  = pos[1]-T.dy()                         ;
   //  double dz  = pos[2]-T.dz()                         ;
   //  double x   = dx*Ax[0]+dy*Ax[1]+dz*Ax[2]            ;
   //  double y   = dx*Ay[0]+dy*Ay[1]+dz*Ay[2]            ;
   //  double z   = dx*Az[0]+dy*Az[1]+dz*Az[2]            ;
   //  double ax  = dir[0]*Ax[0]+dir[1]*Ax[1]+dir[2]*Ax[2];
   //  double ay  = dir[0]*Ay[0]+dir[1]*Ay[1]+dir[2]*Ay[2];
   //  double at  = ax*ax+ay*ay                           ;
   //  double r   = sqrt(x*x+y*y)                         ;
   //  double R   = bounds().r()                          ;
  
   // Step to surface
   //
   int ns = 0;
   double s1 = 0.;
   double s2 = 0.;
  
   if (at != 0.) {
  
     const double inv_at = 1. / at;
     double A = -(ax * x + ay * y) * inv_at;
     double B = A * A + (R - r) * (R - r) * inv_at;
  
     if (B >= 0.) {
  
       B = sqrt(B);
       if (B > tolb) {
         if (A > 0.) {
           s1 = A - B;
           s2 = A + B;
         } else {
           s1 = A + B;
           s2 = A - B;
         }
         ns = 2;
       } else {
         s1 = A;
         ns = 1;
       }
     }
   }
   double sr = r - R;
   if (!bound)
     return {ns, fabs(sr), true, s1, s2};
  
   // Min distance to surface
   //
   Amg::Vector2D lp(atan2(y, x) * R, 0.);
  
   double d = bounds().minDistance(lp);
   double sz = fabs(z) - bounds().halflengthZ();
   if (sz <= 0.)
     sz = 0.;
   double dist = sr * sr + sz * sz;
   if (d > 0.)
     dist += ((d * d) * (sr / R + 1.));
  
   return {ns, sqrt(dist), true, s1, s2};
 }

◆ straightLineIntersection() [1/2]

Trk::Intersection Trk::CylinderSurface::straightLineIntersection	(	const Amg::Vector3D &	pos,
		const Amg::Vector3D &	dir,
		bool	forceDir = `false`,
		Trk::BoundaryCheck	bchk = `false`
	)		const

finaloverridevirtualinherited

fast straight line intersection schema - provides closest intersection and (signed) path length

mathematical motivation:

The calculation will be done in the 3-dim frame of the cylinder, i.e. the symmetry axis of the cylinder is the z-axis, x- and y-axis are perpenticular to the the z-axis. In this frame the cylinder is centered around the origin. Therefore the two points describing the line have to be first recalculated into the new frame. Suppose, this is done, the intersection is straight forward:
may \(p_{1}=(p_{1x}, p_{1y}, p_{1z}), p_{2}=(p_{2x}, p_{2y}, p_{2z}) \)@_fakenlthe two points describing the 3D-line, then the line in the \(x-y\)plane can be written as \(y=kx+d\), where \(k =\frac{p_{2y}-p_{1y}}{p_{2x}-p_{1x}}\)such as \(d=\frac{p_{2x}p_{1y}-p_{1x}p_{2y}}{p_{2x}-p_{1x}},\)
and intersects with the corresponding circle \(x^{2}+y^{2} = R^{2}. \)
The solutions can then be found by a simple quadratic equation and reinsertion into the line equation.

Implements Trk::Surface.

Definition at line 277 of file CylinderSurface.cxx.

 {
   bool needsTransform = m_transforms || m_associatedDetElement;
   // create the hep points
   Amg::Vector3D point1 = pos;
   Amg::Vector3D direction = dir;
   if (needsTransform) {
     Amg::Transform3D invTrans = inverseTransformHelper();
     point1 = invTrans * pos;
     direction = invTrans.linear() * dir;
   }
   // the bounds radius
   double R = bounds().r();
   double t1 = 0.;
   double t2 = 0.;
   if (direction.x()) {
     // get line and circle constants
     double idirx = 1. / direction.x();
     double k = direction.y() * idirx;
     double d = point1.y() - point1.x() * k;
     // and solve the qaudratic equation
     Trk::RealQuadraticEquation pquad(1 + k * k, 2 * k * d, d * d - R * R);
     if (pquad.solutions != Trk::none) {
       // the solutions in the 3D frame of the cylinder
       t1 = (pquad.first - point1.x()) * idirx;
       t2 = (pquad.second - point1.x()) * idirx;
     } else // bail out if no solution exists
       return {pos, 0., false};
   } else if (direction.y()) {
     // x value is the one of point1
     // x^2 + y^2 = R^2
     // y = sqrt(R^2-x^2)
     double x = point1.x();
     double r2mx2 = R * R - x * x;
     // bail out if no solution
     if (r2mx2 < 0.)
       return {pos, 0., false};
     double y = sqrt(r2mx2);
     // assign parameters and solutions
     double idiry = 1. / direction.y();
     t1 = (y - point1.y()) * idiry;
     t2 = (-y - point1.y()) * idiry;
   } else {
     return {pos, 0., false};
   }
   Amg::Vector3D sol1raw(point1 + t1 * direction);
   Amg::Vector3D sol2raw(point1 + t2 * direction);
   // now reorder and return
   Amg::Vector3D solution(0, 0, 0);
   double path = 0.;
  
   // first check the validity of the direction
   bool isValid = true;
  
   // both solutions are of same sign, take the smaller, but flag as false if not
   // forward
   if (t1 * t2 > 0 || !forceDir) {
     // asign validity
     isValid = forceDir ? (t1 > 0.) : true;
     // assign the right solution
     if (t1 * t1 < t2 * t2) {
       solution = sol1raw;
       path = t1;
     } else {
       solution = sol2raw;
       path = t2;
     }
   } else {
     if (t1 > 0.) {
       solution = sol1raw;
       path = t1;
     } else {
       solution = sol2raw;
       path = t2;
     }
   }
   // the solution is still in the local 3D frame, direct check
   isValid =
     bchk ? (isValid && m_bounds->inside3D(solution,
                                           Trk::Surface::s_onSurfaceTolerance,
                                           Trk::Surface::s_onSurfaceTolerance))
          : isValid;
  
   // now return
   return needsTransform ? Intersection(transform() * solution, path, isValid)
                         : Intersection(solution, path, isValid);
 }

◆ straightLineIntersection() [2/2]

template<class T >

Intersection Trk::Surface::straightLineIntersection	(	const T &	pars,
		bool	forceDir = `false`,
		const Trk::BoundaryCheck &	bchk = `false`
	)		const

inlineinherited

fst straight line intersection schema - templated for charged and neutral parameters

Definition at line 351 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

   {
     return straightLineIntersection(
       pars.position(), pars.momentum().unit(), forceDir, bchk);
   }

◆ subSurface() [1/2]

const Trk::Surface * Trk::Layer::subSurface ( const Amg::Vector2D & lp ) const

inherited

If no subSurface array is defined or no subSurface can be found to the given Amg::Vector2D, it would return 0.

Definition at line 126 of file Layer.cxx.

                                                                   {
   if (m_surfaceArray) return m_surfaceArray->object(lp);
   return nullptr;
 }

◆ subSurface() [2/2]

const Trk::Surface * Trk::Layer::subSurface ( const Amg::Vector3D & gp ) const

inherited

If no subSurface array is defined or no subSurface can be found to the given Amg::Vector3D, it would return 0.

Definition at line 121 of file Layer.cxx.

                                                                   {
   if (m_surfaceArray) return m_surfaceArray->object(gp);
   return nullptr;
 }

◆ subSurfaceReference()

const Trk::Surface * Trk::Layer::subSurfaceReference ( unsigned int idx = 0 ) const

inherited

Return a reference sub surface of the layer, usually the first one in the array.

if an index is given (and the index is in range), this surface is returned

Definition at line 131 of file Layer.cxx.

                                                                       {
   // the reference surface
   const Trk::Surface* referenceSurface = nullptr;
   if (m_surfaceArray) {
     // get a reference surface
     Trk::BinnedArraySpan<Trk::Surface const * const> surfaces =  std::as_const(*m_surfaceArray).arrayObjects();
     // get a reference surface
     unsigned int rfSurfaces = surfaces.size();
     if (idx && idx < rfSurfaces) return surfaces[idx];
     // get the first one which is non zero
     for (unsigned int rsf = 0; rsf < rfSurfaces; ++rsf) {
       referenceSurface = surfaces[rsf];
       if (referenceSurface) break;
     }
   }
   return referenceSurface;
 }

◆ surfaceArray() [1/2]

SurfaceArray* Trk::Layer::surfaceArray ( )

inherited

Return the entire SurfaceArray, returns nullptr if no SurfaceArray.

◆ surfaceArray() [2/2]

const SurfaceArray* Trk::Layer::surfaceArray ( ) const

inherited

Return the entire SurfaceArray, returns nullptr if no SurfaceArray.

◆ surfaceOnApproach()

const Trk::Surface & Trk::CylinderLayer::surfaceOnApproach	(	const Amg::Vector3D &	pos,
		const Amg::Vector3D &	dir,
		PropDirection	pdir,
		const BoundaryCheck &	bcheck,
		bool	resolveSubSurfaces = `0`,
		const ICompatibilityEstimator *	ice = `nullptr`
	)		const

finaloverridevirtual

Surface seen on approach - if not defined differently, it is the surfaceRepresentation()

Reimplemented from Trk::Layer.

Definition at line 268 of file CylinderLayer.cxx.

                                              {
   // resolve the surfaces
   if (m_approachDescriptor && resolveSubSurfaces) {
     // resolve based on straight line intersection
     return approachSurface(pos, double(pDir) * mom.unit(), bcheck);
   }
   return surfaceRepresentation();
 }

◆ surfaceRepresentation() [1/2]

const Trk::CylinderSurface & Trk::CylinderLayer::surfaceRepresentation ( ) const

finaloverridevirtual

Transforms the layer into a Surface representation for extrapolation.

Implements Trk::Layer.

Definition at line 154 of file CylinderLayer.cxx.

 {
   return (*this);
 }

◆ surfaceRepresentation() [2/2]

Trk::CylinderSurface & Trk::CylinderLayer::surfaceRepresentation ( )

finaloverridevirtual

Implements Trk::Layer.

Definition at line 160 of file CylinderLayer.cxx.

 {
   return (*this);
 }

◆ thickness()

double Trk::Layer::thickness ( ) const

inherited

Return the Thickness of the Layer.

◆ transform()

const Amg::Transform3D& Trk::Surface::transform ( ) const

inherited

Returns HepGeom::Transform3D by reference.

◆ type()

constexpr virtual SurfaceType Trk::CylinderSurface::type ( ) const

constexprfinaloverridevirtualinherited

Return the surface type.

Implements Trk::Surface.

◆ uniqueClone()

std::unique_ptr<Surface> Trk::Surface::uniqueClone ( ) const

inherited

NVI method returning unique_ptr clone.

Member Data Documentation

◆ m_approachDescriptor

std::unique_ptr<IApproachDescriptor> Trk::CylinderLayer::m_approachDescriptor

protected

surfaces on approach to the layer

Definition at line 169 of file CylinderLayer.h.

◆ m_associatedDetElement

const TrkDetElementBase* Trk::Surface::m_associatedDetElement {}

protectedinherited

Not owning Pointer to the TrkDetElementBase.

Definition at line 439 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

◆ m_associatedDetElementId

Identifier Trk::Surface::m_associatedDetElementId {}

protectedinherited

Identifier for the TrkDetElementBase.

Definition at line 442 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

◆ m_associatedLayer

const Layer* Trk::Surface::m_associatedLayer {}

protectedinherited

The associated layer Trk::Layer.

layer in which the Surface is embedded (not owning pointed)

Definition at line 448 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

◆ m_binUtility

const BinUtility* Trk::Layer::m_binUtility

protectedinherited

Enclosing TrackingVolume.

Definition at line 301 of file Layer.h.

◆ m_bounds

SharedObject<const CylinderBounds> Trk::CylinderSurface::m_bounds

protectedinherited

The global reference point (== a point on the surface)

Definition at line 290 of file CylinderSurface.h.

◆ m_enclosingDetachedTrackingVolume

const DetachedTrackingVolume* Trk::Layer::m_enclosingDetachedTrackingVolume

protectedinherited

Definition at line 305 of file Layer.h.

◆ m_enclosingTrackingVolume

const TrackingVolume* Trk::Layer::m_enclosingTrackingVolume

protectedinherited

Enclosing DetachedTrackingVolume.

Definition at line 303 of file Layer.h.

◆ m_index

LayerIndex Trk::Layer::m_index

protectedinherited

LayerIndex.

Definition at line 307 of file Layer.h.

◆ m_layerMaterialProperties

std::unique_ptr<LayerMaterialProperties> Trk::Layer::m_layerMaterialProperties

protectedinherited

thickness of the Layer

Definition at line 289 of file Layer.h.

◆ m_layerThickness

double Trk::Layer::m_layerThickness

protectedinherited

descriptor for overlap/next surface (owning ptr)

Definition at line 291 of file Layer.h.

◆ m_layerType

int Trk::Layer::m_layerType

protectedinherited

active passive layer

Definition at line 308 of file Layer.h.

◆ m_materialLayer

std::shared_ptr<Layer> Trk::Surface::m_materialLayer {}

protectedinherited

Possibility to attach a material descrption.

potentially given as the associated material layer (not owning pointer)

Definition at line 453 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

◆ m_nextLayer

const Layer* Trk::Layer::m_nextLayer

protectedinherited

BinUtility for next/previous decision.

Definition at line 299 of file Layer.h.

◆ m_overlapDescriptor

std::unique_ptr<OverlapDescriptor> Trk::Layer::m_overlapDescriptor

protectedinherited

Definition at line 293 of file Layer.h.

◆ m_owner

SurfaceOwner Trk::Surface::m_owner {SurfaceOwner::noOwn}

protectedinherited

enum for surface owner : 0 free surface

Definition at line 455 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

◆ m_previousLayer

const Layer* Trk::Layer::m_previousLayer

protectedinherited

< the previous Layer according to BinGenUtils

next Layer according to BinGenUtils

Definition at line 297 of file Layer.h.

◆ m_ref

double Trk::Layer::m_ref

protectedinherited

reference measure for local coordinate convertors

Definition at line 309 of file Layer.h.

◆ m_referencePoint

CxxUtils::CachedUniquePtr<Amg::Vector3D> Trk::CylinderSurface::m_referencePoint

protectedinherited

The rotational symmetry axis.

Definition at line 292 of file CylinderSurface.h.

◆ m_rotSymmetryAxis

CxxUtils::CachedUniquePtr<Amg::Vector3D> Trk::CylinderSurface::m_rotSymmetryAxis

protectedinherited

Definition at line 294 of file CylinderSurface.h.

◆ m_surfaceArray

std::unique_ptr<SurfaceArray> Trk::Layer::m_surfaceArray

protectedinherited

MaterialPoperties of this layer Surface.

Definition at line 287 of file Layer.h.

◆ m_transforms

std::unique_ptr<Transforms> Trk::Surface::m_transforms {}

protectedinherited

Definition at line 436 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

◆ s_numberOfInstantiations

std::atomic_size_t Trk::ObjectCounter< Trk::Surface >::s_numberOfInstantiations

inlinestaticinherited

Definition at line 22 of file TrkObjectCounter.h.

◆ s_onSurfaceTolerance

constexpr double Trk::Surface::s_onSurfaceTolerance = 10e-5

staticconstexprprotectedinherited

Tolerance for being on Surface.

Definition at line 458 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

◆ staticType

constexpr SurfaceType Trk::CylinderSurface::staticType = SurfaceType::Cylinder

staticconstexprinherited

The surface type static constexpr.

Definition at line 59 of file CylinderSurface.h.

The documentation for this class was generated from the following files:

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Static Protected Attributes

Private Member Functions

Detailed Description

Member Typedef Documentation

◆ ChargedTrackParametersUniquePtr

◆ NeutralTrackParametersUniquePtr

Constructor & Destructor Documentation

◆ CylinderLayer() [1/10]

◆ CylinderLayer() [2/10]

◆ CylinderLayer() [3/10]

◆ CylinderLayer() [4/10]

◆ CylinderLayer() [5/10]

◆ CylinderLayer() [6/10]

◆ CylinderLayer() [7/10]

◆ CylinderLayer() [8/10]

◆ CylinderLayer() [9/10]

◆ CylinderLayer() [10/10]

◆ ~CylinderLayer()

Member Function Documentation

◆ approachSurface()

◆ assignMaterialProperties()

◆ associatedDetectorElement()

◆ associatedDetectorElementIdentifier()

◆ associatedLayer()

◆ associateLayer()

◆ baseSurface()

◆ binUtility()

◆ bounds()

◆ buildApproachDescriptor()

◆ cachedTransform()

◆ center()

◆ clone()

◆ compactify()

◆ compatibleSurfaces() [1/2]

◆ compatibleSurfaces() [2/2]

◆ createUniqueNeutralParameters() [1/2]

◆ createUniqueNeutralParameters() [2/2]

◆ createUniqueParameters() [1/2]

◆ createUniqueParameters() [2/2]

◆ createUniqueTrackParameters() [1/2]

◆ createUniqueTrackParameters() [2/2]

◆ dump() [1/2]

◆ dump() [2/2]

◆ encloseDetachedTrackingVolume()

◆ encloseTrackingVolume()

◆ enclosingDetachedTrackingVolume()

◆ enclosingTrackingVolume()

◆ fullUpdateMaterialProperties()

◆ getCompatibleSurfaces()

◆ getRef()

◆ globalReferencePoint()

◆ globalToLocal() [1/3]

◆ globalToLocal() [2/3]

◆ globalToLocal() [3/3]

◆ hasBounds()

◆ hasSubStructure()

◆ insideBounds()

◆ insideBoundsCheck()

◆ inverseTransformHelper()

◆ inverseTransformMultHelper()

◆ isActive()

◆ isFree()

◆ isOnLayer()

◆ isOnSurface()

◆ layerIndex()

◆ layerMaterialProperties()

◆ layerType()

◆ localParametersToPosition()

◆ localToGlobal() [1/5]

◆ localToGlobal() [2/5]

◆ localToGlobal() [3/5]

◆ localToGlobal() [4/5]

◆ localToGlobal() [5/5]

◆ materialLayer() [1/2]