ATLAS Offline Software
Loading...
Searching...
No Matches
Trk::SubtractedDiscSurface Class Referencefinal

Class for a planar subtracted/shared surface in the ATLAS detector. More...

#include <SubtractedDiscSurface.h>

Inheritance diagram for Trk::SubtractedDiscSurface:
Collaboration diagram for Trk::SubtractedDiscSurface:

Public Types

using ChargedTrackParametersUniquePtr
 Unique ptr types.
using NeutralTrackParametersUniquePtr

Public Member Functions

 SubtractedDiscSurface ()=default
 Defaults.
 SubtractedDiscSurface (SubtractedDiscSurface &&)=default
SubtractedDiscSurfaceoperator= (SubtractedDiscSurface &&)=default
virtual ~SubtractedDiscSurface ()=default
 SubtractedDiscSurface (const DiscSurface &ps, std::shared_ptr< const AreaExcluder > vol, bool shared)
 Constructor.
 SubtractedDiscSurface (const SubtractedDiscSurface &psf)
 Copy Constructor.
 SubtractedDiscSurface (const SubtractedDiscSurface &psf, const Amg::Transform3D &shift)
 Copy Constructor.
SubtractedDiscSurfaceoperator= (const SubtractedDiscSurface &psf)
 Assignment operator.
virtual bool operator== (const Surface &sf) const override
 Equality operator.
bool shared () const
 This method indicates the subtraction mode.
virtual bool insideBounds (const Amg::Vector2D &locpos, double tol1=0., double tol2=0.) const override final
 This method calls the inside() method of the Bounds.
const AreaExcludersubtractedVolume () const
 This method allows access to the subtracted part.
virtual std::string name () const override final
 Return properly formatted class name for screen output.
bool operator== (const DiscSurface &cf) const
virtual DiscSurfaceclone () const override
 Virtual constructor.
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.
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.
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.
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.
template<int DIM, class T>
std::unique_ptr< ParametersT< DIM, T, DiscSurface > > 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.
template<int DIM, class T>
std::unique_ptr< ParametersT< DIM, T, DiscSurface > > 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.
virtual constexpr SurfaceType type () const override final
 Return the surface type.
virtual const Amg::Vector3DglobalReferencePoint () const override final
 Returns a global reference point: For the Disc this is \( (R*cos(\phi), R*sin(\phi),0)*transform() \) Where \( r, \phi \) denote the r(), averagePhi() of the Bounds.
const SurfaceBoundsbounds () const override final
 This method returns the bounds by reference.
virtual bool insideBoundsCheck (const Amg::Vector2D &locpos, const BoundaryCheck &bchk) const override final
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.
virtual Amg::Vector2D localParametersToPosition (const LocalParameters &locpars) const override final
 Specialized for DiscSurface : LocalParameters to Vector2D.
virtual void localToGlobal (const Amg::Vector2D &locp, const Amg::Vector3D &mom, Amg::Vector3D &glob) const override
 Specialized for DiscSurface: LocalToGlobal method without dynamic memory allocation.
Amg::Vector3D localToGlobal (const Amg::Vector2D &locpos) const
 This method returns the GlobalPosition from a LocalPosition uses the per surface localToGlobal.
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.
Amg::Vector3D localToGlobal (const LocalParameters &locpars) const
 This method returns the GlobalPosition from LocalParameters The LocalParameters can be outside Surface bounds.
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.
virtual bool globalToLocal (const Amg::Vector3D &glob, const Amg::Vector3D &mom, Amg::Vector2D &loc) const override
 Specialized for DiscSurface: GlobalToLocal method without dynamic memory allocation - boolean checks if on surface.
std::optional< Amg::Vector2DglobalToLocal (const Amg::Vector3D &glopos, double tol=0.) const
 This method returns the LocalPosition from a provided GlobalPosition.
std::optional< Amg::Vector2DglobalToLocal (const Amg::Vector3D &glopos, const Amg::Vector3D &glomom) const
 This method returns the LocalPosition from a provided GlobalPosition.
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 - standard: provides closest intersection and (signed) path length forceDir is to provide the closest forward solution
template<typename 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
virtual DistanceSolution straightLineDistanceEstimate (const Amg::Vector3D &pos, const Amg::Vector3D &dir) const override
 fast straight line distance evaluation to Surface
virtual DistanceSolution straightLineDistanceEstimate (const Amg::Vector3D &pos, const Amg::Vector3D &dir, bool Bound) const override
 fast straight line distance evaluation to Surface - with bound option
bool operator!= (const Surface &sf) const
 Non-equality operator.
std::unique_ptr< SurfaceuniqueClone () const
 NVI method returning unique_ptr clone.
const Amg::Transform3DcachedTransform () const
 Return the cached transformation directly.
const Amg::Transform3Dtransform () const
 Returns HepGeom::Transform3D by reference.
const Amg::Vector3Dcenter () const
 Returns the center position of the Surface.
virtual const Amg::Vector3Dnormal () const
 Returns the normal vector of the Surface (i.e.
virtual Amg::Vector3D normal (const Amg::Vector2D &lp) const
 Returns a normal vector at a specific local position.
const TrkDetElementBaseassociatedDetectorElement () const
 return associated Detector Element
Identifier associatedDetectorElementIdentifier () const
 return Identifier of the associated Detector Element
const Trk::LayerassociatedLayer () const
 return the associated Layer
const Trk::MaterialLayermaterialLayer () const
 return the material Layer
Trk::MaterialLayermaterialLayer ()
virtual const Trk::SurfacebaseSurface () const
 return the base surface (simplified for persistification)
std::optional< Amg::Vector2DpositionOnSurface (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).
template<class T>
bool onSurface (const T &parameters, const BoundaryCheck &bchk=BoundaryCheck(true)) const
 The templated Parameters OnSurface method - checks on surface pointer first.
virtual double pathCorrection (const Amg::Vector3D &pos, const Amg::Vector3D &mom) const
 the pathCorrection for derived classes with thickness - it reflects if the direction projection is positive or negative
virtual Amg::RotationMatrix3D measurementFrame (const Amg::Vector3D &glopos, const Amg::Vector3D &glomom) const
 Return the measurement frame - this is needed for alignment, in particular for StraightLine and Perigee Surface.
bool isFree () const
 Returns 'true' if this surface is 'free', i.e.
bool isActive () const
 Return 'true' if this surface is owned by the detector element.
void setTransform (const Amg::Transform3D &trans)
 Set the transform updates center and normal.
void setOwner (SurfaceOwner x)
 set Ownership
SurfaceOwner owner () const
 return ownership
void setMaterialLayer (std::shared_ptr< Trk::MaterialLayer > mlay)
 set material layer
virtual MsgStream & dump (MsgStream &sl) const
 Output Method for MsgStream, to be overloaded by child classes.
virtual std::ostream & dump (std::ostream &sl) const
 Output Method for std::ostream, to be overloaded by child classes.
void associateLayer (const Layer &lay)
 method to associate a Trk::Layer.

Static Public Member Functions

static std::size_t numberOfInstantiations ()

Static Public Attributes

static constexpr SurfaceType staticType = SurfaceType::Disc
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.
Amg::Vector3D inverseTransformMultHelper (const Amg::Vector3D &glopos) const

Protected Attributes

std::shared_ptr< const AreaExcluderm_subtrVol {nullptr}
bool m_shared {true}
std::shared_ptr< const SurfaceBoundsm_bounds
 reference Point on the Surface
CxxUtils::CachedUniquePtr< Amg::Vector3Dm_referencePoint
 static member for boundless approach
std::unique_ptr< Transformsm_transforms {}
 Unique Pointer to the Transforms struct.
const TrkDetElementBasem_associatedDetElement {}
 Not owning Pointer to the Detector Element.
Identifier m_associatedDetElementId {}
 Identifier to the Detector Element.
const Layerm_associatedLayer {}
 The associated layer Trk::Layer This is the layer in which the Surface is embedded.
std::shared_ptr< MaterialLayerm_materialLayer {}
 Possibility to attach a material description to the surface.
SurfaceOwner m_owner {SurfaceOwner::noOwn}
 enum for surface owner : default free surface

Static Protected Attributes

static const NoBounds s_boundless
static constexpr double s_onSurfaceTolerance = 10e-5
 Tolerance for being on Surface.

Detailed Description

Class for a planar subtracted/shared surface in the ATLAS detector.

It owns its surface bounds and subtracted volume.

Author
Sarka.nosp@m..Tod.nosp@m.orova.nosp@m.@cer.nosp@m.n.ch

Definition at line 32 of file SubtractedDiscSurface.h.

Member Typedef Documentation

◆ ChargedTrackParametersUniquePtr

Initial value:
std::unique_ptr<ParametersBase<5, Trk::Charged>>

Unique ptr types.

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

◆ NeutralTrackParametersUniquePtr

Initial value:
std::unique_ptr<ParametersBase<5, Trk::Neutral>>

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

Constructor & Destructor Documentation

◆ SubtractedDiscSurface() [1/5]

Trk::SubtractedDiscSurface::SubtractedDiscSurface ( )
default

Defaults.

Copies need care due to unique_ptr

◆ SubtractedDiscSurface() [2/5]

Trk::SubtractedDiscSurface::SubtractedDiscSurface ( SubtractedDiscSurface && )
default

◆ ~SubtractedDiscSurface()

virtual Trk::SubtractedDiscSurface::~SubtractedDiscSurface ( )
virtualdefault

◆ SubtractedDiscSurface() [3/5]

Trk::SubtractedDiscSurface::SubtractedDiscSurface ( const DiscSurface & ps,
std::shared_ptr< const AreaExcluder > vol,
bool shared )

Constructor.

Definition at line 19 of file SubtractedDiscSurface.cxx.

23 : Trk::DiscSurface(ps), m_subtrVol(std::move(vol)), m_shared(shared) {}
std::shared_ptr< const AreaExcluder > m_subtrVol
bool shared() const
This method indicates the subtraction mode.

◆ SubtractedDiscSurface() [4/5]

Trk::SubtractedDiscSurface::SubtractedDiscSurface ( const SubtractedDiscSurface & psf)
default

Copy Constructor.

◆ SubtractedDiscSurface() [5/5]

Trk::SubtractedDiscSurface::SubtractedDiscSurface ( const SubtractedDiscSurface & psf,
const Amg::Transform3D & shift )

Copy Constructor.

Definition at line 29 of file SubtractedDiscSurface.cxx.

30 : Trk::DiscSurface(psf, shift)
31 , m_subtrVol{psf.m_subtrVol}
32 , m_shared(psf.m_shared)
33{}

Member Function Documentation

◆ 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 a Trk::Layer.

We do not take ownership as this is typically already owned by the Geometry

◆ baseSurface()

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

return the base surface (simplified for persistification)

◆ bounds()

const SurfaceBounds & Trk::DiscSurface::bounds ( ) const
finaloverridevirtualinherited

This method returns the bounds by reference.

Implements Trk::Surface.

◆ 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 DiscSurface * Trk::DiscSurface::clone ( ) const
overridevirtualinherited

Virtual constructor.

Implements Trk::Surface.

Reimplemented in Trk::SlidingDiscSurface.

◆ createUniqueNeutralParameters() [1/2]

Trk::Surface::NeutralTrackParametersUniquePtr Trk::DiscSurface::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 201 of file DiscSurface.cxx.

203 {
204 return std::make_unique<ParametersT<5, Neutral, DiscSurface>>(
205 position, momentum, charge, *this, std::move(cov));
206}
double charge(const T &p)
Definition AtlasPID.h:997

◆ createUniqueNeutralParameters() [2/2]

Trk::Surface::NeutralTrackParametersUniquePtr Trk::DiscSurface::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 191 of file DiscSurface.cxx.

193 {
194 return std::make_unique<ParametersT<5, Neutral, DiscSurface>>(
195 l1, l2, phi, theta, qop, *this, std::move(cov));
196}
@ theta
Definition ParamDefs.h:66
@ phi
Definition ParamDefs.h:75

◆ createUniqueParameters() [1/2]

template<int DIM, class T>
std::unique_ptr< ParametersT< DIM, T, DiscSurface > > Trk::DiscSurface::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, DiscSurface > > Trk::DiscSurface::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::DiscSurface::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 181 of file DiscSurface.cxx.

183 {
184 return std::make_unique<ParametersT<5, Charged, DiscSurface>>(
185 position, momentum, charge, *this, std::move(cov));
186}

◆ createUniqueTrackParameters() [2/2]

Trk::Surface::ChargedTrackParametersUniquePtr Trk::DiscSurface::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 172 of file DiscSurface.cxx.

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

◆ dump() [1/2]

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

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

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

Definition at line 157 of file Surface.cxx.

158{
159 sl << std::setiosflags(std::ios::fixed);
160 sl << std::setprecision(4);
161 sl << name() << std::endl;
162 if (associatedDetectorElement()!=nullptr){
163 sl<<" Detector Type = "<<associatedDetectorElement()->detectorTypeString()<<std::endl;
164 }
165 sl << " Center position (x, y, z) = (" << center().x() << ", " << center().y() << ", " << center().z() << ")"
166 << std::endl;
168 Amg::Vector3D rotX(rot.col(0));
169 Amg::Vector3D rotY(rot.col(1));
170 Amg::Vector3D rotZ(rot.col(2));
171 sl << std::setprecision(6);
172 sl << " Rotation: colX = (" << rotX(0) << ", " << rotX(1) << ", " << rotX(2) << ")" << std::endl;
173 sl << " colY = (" << rotY(0) << ", " << rotY(1) << ", " << rotY(2) << ")" << std::endl;
174 sl << " colZ = (" << rotZ(0) << ", " << rotZ(1) << ", " << rotZ(2) << ")" << std::endl;
175 sl << " Bounds : " << bounds();
176 if (!checkTransform(*this)) {
177 sl << std::endl << " NOT a strict rotation matrix." << std::endl;
178 }
179 sl << std::setprecision(-1);
180 return sl;
181}
const TrkDetElementBase * associatedDetectorElement() const
return associated Detector Element
const Amg::Transform3D & transform() const
Returns HepGeom::Transform3D by reference.
virtual const SurfaceBounds & bounds() const =0
Surface Bounds method.
const Amg::Vector3D & center() const
Returns the center position of the Surface.
virtual std::string name() const =0
Return properly formatted class name.
std::string detectorTypeString() const
Returns a string of the Detector element type.
Eigen::Matrix< double, 3, 3 > RotationMatrix3D
Eigen::Matrix< double, 3, 1 > Vector3D

◆ 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::DetElementSurface, and Trk::PerigeeSurface.

Definition at line 185 of file Surface.cxx.

186{
187 sl << std::setiosflags(std::ios::fixed);
188 sl << std::setprecision(4);
189 sl << name() << std::endl;
190 if (associatedDetectorElement()!=nullptr){
191 sl<<" Detector Type = "<<associatedDetectorElement()->detectorTypeString()<<std::endl;
192 }
193 sl << " Center position (x, y, z) = (" << center().x() << ", " << center().y() << ", " << center().z() << ")"
194 << std::endl;
196 Amg::Vector3D rotX(rot.col(0));
197 Amg::Vector3D rotY(rot.col(1));
198 Amg::Vector3D rotZ(rot.col(2));
199 sl << std::setprecision(6);
200 sl << " Rotation: colX = (" << rotX(0) << ", " << rotX(1) << ", " << rotX(2) << ")" << std::endl;
201 sl << " colY = (" << rotY(0) << ", " << rotY(1) << ", " << rotY(2) << ")" << std::endl;
202 sl << " colZ = (" << rotZ(0) << ", " << rotZ(1) << ", " << rotZ(2) << ")" << std::endl;
203 sl << " Bounds : " << bounds();
204 if (!checkTransform(*this)) {
205 sl << std::endl << " NOT a strict rotation matrix." << std::endl;
206 }
207 sl << std::setprecision(-1);
208 return sl;
209}

◆ globalReferencePoint()

const Amg::Vector3D & Trk::DiscSurface::globalReferencePoint ( ) const
finaloverridevirtualinherited

Returns a global reference point: For the Disc this is \( (R*cos(\phi), R*sin(\phi),0)*transform() \) Where \( r, \phi \) denote the r(), averagePhi() of the Bounds.

Reimplemented from Trk::Surface.

Definition at line 209 of file DiscSurface.cxx.

210{
211 if (!m_referencePoint) {
212 const Trk::DiscBounds* dbo =
213 dynamic_cast<const Trk::DiscBounds*>(&(bounds()));
214 if (dbo) {
215 double rMedium = bounds().r();
216 double phi = dbo->averagePhi();
217 Amg::Vector3D gp(rMedium * cos(phi), rMedium * sin(phi), 0.);
218 m_referencePoint.set(std::make_unique<Amg::Vector3D>(transform() * gp));
219 } else {
220 const Trk::DiscTrapezoidalBounds* dtbo =
221 dynamic_cast<const Trk::DiscTrapezoidalBounds*>(&(bounds()));
222 // double rMedium = dtbo ? bounds().r() : dtbo->rCenter() ; //nonsense, or
223 // logic inverted?
224 double rMedium = bounds().r();
225 double phi = dtbo ? dtbo->averagePhi() : 0.;
226 Amg::Vector3D gp(rMedium * cos(phi), rMedium * sin(phi), 0.);
227 m_referencePoint.set(std::make_unique<Amg::Vector3D>(transform() * gp));
228 }
229 }
230 return (*m_referencePoint);
231}
double averagePhi() const
This method returns the average phi.
CxxUtils::CachedUniquePtr< Amg::Vector3D > m_referencePoint
static member for boundless approach
const SurfaceBounds & bounds() const override final
This method returns the bounds by reference.
double averagePhi() const
This method returns the average phi.
virtual double r() const =0
Interface method for the maximal extension or the radius.

◆ globalToLocal() [1/3]

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

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

local<->global transformation in case of polar local coordinates

Implements Trk::Surface.

Reimplemented in Trk::SlidingDiscSurface.

Definition at line 250 of file DiscSurface.cxx.

253{
254 Amg::Vector3D loc3Dframe = inverseTransformMultHelper(glopos);
255 locpos = Amg::Vector2D(loc3Dframe.perp(), loc3Dframe.phi());
256 return (std::fabs(loc3Dframe.z()) <= s_onSurfaceTolerance);
257}
static constexpr double s_onSurfaceTolerance
Tolerance for being on Surface.
Amg::Vector3D inverseTransformMultHelper(const Amg::Vector3D &glopos) const
Eigen::Matrix< double, 2, 1 > Vector2D

◆ 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

◆ insideBounds()

bool Trk::SubtractedDiscSurface::insideBounds ( const Amg::Vector2D & locpos,
double tol1 = 0.,
double tol2 = 0. ) const
inlinefinaloverridevirtual

This method calls the inside() method of the Bounds.

Reimplemented from Trk::DiscSurface.

Definition at line 80 of file SubtractedDiscSurface.h.

83{
84 // no subtracted Volume exists
85 if (!m_subtrVol.get()) {
86 return (this->bounds().inside(locpos, tol1, tol2));
87 }
88 // subtracted Volume exists, needs to be checked
89 double rPos = locpos[Trk::locR];
90 double phiPos = locpos[Trk::locPhi];
91 Amg::Vector3D gp(rPos * cos(phiPos), rPos * sin(phiPos), 0.);
92 if (m_shared) {
93 return (this->bounds().inside(locpos, tol1, tol2) &&
94 m_subtrVol.get()->inside(gp, 0.));
95 }
96 bool in(this->bounds().inside(locpos, tol1, tol2) &&
97 !m_subtrVol.get()->inside(gp, 0.));
98
99 return in;
100}
@ locR
Definition ParamDefs.h:44
@ locPhi
local polar
Definition ParamDefs.h:45

◆ insideBoundsCheck()

virtual bool Trk::DiscSurface::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

◆ isOnSurface()

bool Trk::DiscSurface::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::SlidingDiscSurface.

Definition at line 287 of file DiscSurface.cxx.

291{
292 Amg::Vector3D loc3Dframe = inverseTransformMultHelper(glopo);
293 if (std::abs(loc3Dframe.z()) > (s_onSurfaceTolerance + tol1)) {
294 return false;
295 }
296 return (bchk
297 ? bounds().inside(
298 Amg::Vector2D(loc3Dframe.perp(), loc3Dframe.phi()), tol1, tol2)
299 : true);
300}

◆ localParametersToPosition()

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

Specialized for DiscSurface : LocalParameters to Vector2D.

Reimplemented from Trk::Surface.

◆ localToGlobal() [1/5]

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

Specialized for DiscSurface: LocalToGlobal method without dynamic memory allocation.

Implements Trk::Surface.

Reimplemented in Trk::SlidingDiscSurface.

Definition at line 236 of file DiscSurface.cxx.

239{
240 // create the position in the local 3d frame
241 Amg::Vector3D loc3Dframe(locpos[Trk::locR] * cos(locpos[Trk::locPhi]),
242 locpos[Trk::locR] * sin(locpos[Trk::locPhi]),
243 0.);
244 // transport it to the globalframe
245 glopos = transform() * loc3Dframe;
246}

◆ 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::MaterialLayer * Trk::Surface::materialLayer ( )
inherited

◆ materialLayer() [2/2]

const Trk::MaterialLayer * Trk::Surface::materialLayer ( ) const
inherited

return the material Layer

◆ measurementFrame()

Amg::RotationMatrix3D Surface::measurementFrame ( const Amg::Vector3D & glopos,
const Amg::Vector3D & glomom ) const
virtualinherited

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

Reimplemented in GXF::CylinderSurface, Trk::ConeSurface, Trk::CylinderSurface, Trk::PerigeeSurface, and Trk::StraightLineSurface.

Definition at line 135 of file Surface.cxx.

136{
137 return transform().linear();
138}

◆ name()

virtual std::string Trk::SubtractedDiscSurface::name ( ) const
inlinefinaloverridevirtual

Return properly formatted class name for screen output.

Reimplemented from Trk::DiscSurface.

Definition at line 69 of file SubtractedDiscSurface.h.

70 {
71 return "Trk::SubtractedDiscSurface";
72 }

◆ 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 GXF::CylinderSurface, Trk::ConeSurface, Trk::CylinderSurface, and Trk::PerigeeSurface.

◆ normal() [2/2]

virtual Amg::Vector3D Trk::Surface::normal ( const Amg::Vector2D & lp) const
virtualinherited

Returns a normal vector at a specific local position.

Reimplemented in GXF::CylinderSurface, Trk::ConeSurface, Trk::CylinderSurface, and Trk::PerigeeSurface.

◆ numberOfInstantiations()

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

Definition at line 25 of file TrkObjectCounter.h.

26 {
27#ifndef NDEBUG
28 return s_numberOfInstantiations.load();
29#endif
30 return 0;
31 }
Helper to enable counting number of instantiations in debug builds.

◆ 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=() [1/2]

Trk::SubtractedDiscSurface & Trk::SubtractedDiscSurface::operator= ( const SubtractedDiscSurface & psf)
default

Assignment operator.

◆ operator=() [2/2]

SubtractedDiscSurface & Trk::SubtractedDiscSurface::operator= ( SubtractedDiscSurface && )
default

◆ operator==() [1/2]

bool Trk::DiscSurface::operator== ( const DiscSurface & cf) const
inherited

◆ operator==() [2/2]

bool Trk::SubtractedDiscSurface::operator== ( const Surface & sf) const
overridevirtual

Equality operator.

Reimplemented from Trk::DiscSurface.

Definition at line 40 of file SubtractedDiscSurface.cxx.

41{
42 // first check the type not to compare apples with oranges
43 const Trk::SubtractedDiscSurface* sdsf = dynamic_cast<const Trk::SubtractedDiscSurface*>(&sf);
44 if (!sdsf)
45 return false;
46 bool surfaceEqual = Trk::DiscSurface::operator==(sf);
47 bool sharedEqual = (surfaceEqual) ? (shared() == sdsf->shared()) : false;
48 return sharedEqual;
49}
virtual bool operator==(const Surface &sf) const override
Equality operator.

◆ owner()

SurfaceOwner Trk::Surface::owner ( ) const
inherited

return ownership

◆ pathCorrection()

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

the pathCorrection for derived classes with thickness - it reflects if the direction projection is positive or negative

Reimplemented in GXF::CylinderSurface, Trk::ConeSurface, Trk::CylinderSurface, Trk::PerigeeSurface, and Trk::StraightLineSurface.

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

110{
111 std::optional<Amg::Vector2D> posOnSurface = globalToLocal(glopo, tol1);
112 if (!bchk){
113 return posOnSurface;
114 }
115 if (posOnSurface && insideBounds(*posOnSurface, tol1, tol2)){
116 return posOnSurface;
117 }
118 return std::nullopt;
119}
virtual bool globalToLocal(const Amg::Vector3D &glob, const Amg::Vector3D &mom, Amg::Vector2D &loc) const =0
Specified by each surface type: GlobalToLocal method without dynamic memory allocation - boolean chec...
virtual bool insideBounds(const Amg::Vector2D &locpos, double tol1=0., double tol2=0.) const =0
virtual methods to be overwritten by the inherited surfaces

◆ setMaterialLayer()

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

set material layer

◆ setOwner()

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

set Ownership

◆ setTransform()

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

Set the transform updates center and normal.

◆ shared()

bool Trk::SubtractedDiscSurface::shared ( ) const
inline

This method indicates the subtraction mode.

Definition at line 103 of file SubtractedDiscSurface.h.

104{
105 return m_shared;
106}

◆ straightLineDistanceEstimate() [1/2]

Trk::DistanceSolution Trk::DiscSurface::straightLineDistanceEstimate ( const Amg::Vector3D & pos,
const Amg::Vector3D & dir ) const
overridevirtualinherited

fast straight line distance evaluation to Surface

distance to surface

Implements Trk::Surface.

Reimplemented in Trk::SlidingDiscSurface.

Definition at line 304 of file DiscSurface.cxx.

306{
307 double tol = 0.001;
308
309 const Amg::Vector3D& C = center();
310 const Amg::Vector3D& N = normal();
311
312 double S = C.dot(N);
313 double b = S < 0. ? -1 : 1;
314 double d = (pos - C).dot(N); // distance to surface
315
316 double A = b * dir.dot(N);
317 if (A == 0.) { // direction parallel to surface
318 if (fabs(d) < tol) {
319 return {1, 0., true, 0.};
320 }
321 return {0, d, true, 0.};
322 }
323
324 double D = b * (S - (pos.dot(N))) / A;
325 return {1, d, true, D};
326}
virtual const Amg::Vector3D & normal() const
Returns the normal vector of the Surface (i.e.
struct color C
dot(G, fn, nodesToHighlight=[])
Definition dot.py:5

◆ straightLineDistanceEstimate() [2/2]

Trk::DistanceSolution Trk::DiscSurface::straightLineDistanceEstimate ( const Amg::Vector3D & pos,
const Amg::Vector3D & dir,
bool Bound ) const
overridevirtualinherited

fast straight line distance evaluation to Surface - with bound option

Implements Trk::Surface.

Reimplemented in Trk::SlidingDiscSurface.

Definition at line 329 of file DiscSurface.cxx.

332{
333 const Amg::Transform3D& T = transform();
334 const double Az[3] = { T(0, 2), T(1, 2), T(2, 2) };
335
336 // Transformation to cylinder system coordinates
337 //
338 const double dx = pos[0] - T(0, 3);
339 const double dy = pos[1] - T(1, 3);
340 const double dz = pos[2] - T(2, 3);
341 const double z = dx * Az[0] + dy * Az[1] + dz * Az[2];
342 const double az = dir[0] * Az[0] + dir[1] * Az[1] + dir[2] * Az[2];
343
344 // Step to surface
345 //
346 int ns = 0;
347 double s = 0.;
348 if (az != 0.) {
349 s = -z / az;
350 ns = 1;
351 }
352 double dist = std::abs(z);
353 if (!bound) {
354 return {ns, dist, true, s};
355 }
356
357 // Min distance to surface
358 //
359 const double x = dx * T(0, 0) + dy * T(1, 0) + dz * T(2, 0);
360 const double y = dx * T(0, 1) + dy * T(1, 1) + dz * T(2, 1);
361
362 Amg::Vector2D lp(sqrt(x * x + y * y), atan2(y, x));
363
364 double d = bounds().minDistance(lp);
365 if (d > 0.) {
366 dist = std::sqrt(dist * dist + d * d);
367 }
368
369 return {ns, dist, true, s};
370}
virtual double minDistance(const Amg::Vector2D &pos) const =0
Minimal distance to boundary ( > 0 if outside and <=0 if inside)
Eigen::Affine3d Transform3D
unsigned long long T
@ x
Definition ParamDefs.h:55
@ z
global position (cartesian)
Definition ParamDefs.h:57
@ y
Definition ParamDefs.h:56

◆ straightLineIntersection() [1/2]

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

fast straight line intersection schema - standard: provides closest intersection and (signed) path length forceDir is to provide the closest forward solution

mathematical motivation:

the equation of the plane is given by:
\( \vec n \cdot \vec x = \vec n \cdot \vec p,\)
where \( \vec n = (n_{x}, n_{y}, n_{z})\) denotes the normal vector of the plane, \( \vec p = (p_{x}, p_{y}, p_{z})\) one specific point on the plane and \( \vec x = (x,y,z) \) all possible points on the plane.
Given a line with:
\( \vec l(u) = \vec l_{1} + u \cdot \vec v \),
the solution for \( u \) can be written: \( u = \frac{\vec n (\vec p - \vec l_{1})}{\vec n \vec v}\)
If the denominator is 0 then the line lies:

  • either in the plane
  • perpenticular to the normal of the plane

Implements Trk::Surface.

Definition at line 260 of file DiscSurface.cxx.

263 {
264 double denom = dir.dot(normal());
265 if (denom) {
266 double u = (normal().dot((center() - pos))) / (denom);
267 Amg::Vector3D intersectPoint(pos + u * dir);
268 // evaluate the intersection in terms of direction
269 bool isValid = forceDir ? (u > 0.) : true;
270 // evaluate (if necessary in terms of boundaries)
271 isValid = bchk ? (isValid && isOnSurface(intersectPoint)) : isValid;
272 // return the result
273 return Trk::Intersection(intersectPoint, u, isValid);
274 }
275 return Trk::Intersection(pos, 0., false);
276}
bool isValid(const T &p)
Av: we implement here an ATLAS-sepcific convention: all particles which are 99xxxxx are fine.
Definition AtlasPID.h:878
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...
@ u
Enums for curvilinear frames.
Definition ParamDefs.h:77

◆ straightLineIntersection() [2/2]

template<typename 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 352 of file Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h.

356 {
358 pars.position(), pars.momentum().unit(), forceDir, bchk);
359 }
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

◆ subtractedVolume()

const AreaExcluder * Trk::SubtractedDiscSurface::subtractedVolume ( ) const
inline

This method allows access to the subtracted part.

Definition at line 109 of file SubtractedDiscSurface.h.

110{
111 return m_subtrVol.get();
112}

◆ transform()

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

Returns HepGeom::Transform3D by reference.

◆ type()

virtual constexpr SurfaceType Trk::DiscSurface::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_associatedDetElement

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

Not owning Pointer to the Detector Element.

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

441{};

◆ m_associatedDetElementId

Identifier Trk::Surface::m_associatedDetElementId {}
protectedinherited

Identifier to the Detector Element.

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

443{};

◆ m_associatedLayer

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

The associated layer Trk::Layer This is the layer in which the Surface is embedded.

so not owning pointer.

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

448{};

◆ m_bounds

std::shared_ptr<const SurfaceBounds> Trk::DiscSurface::m_bounds
protectedinherited

reference Point on the Surface

Definition at line 284 of file DiscSurface.h.

◆ m_materialLayer

std::shared_ptr<MaterialLayer> Trk::Surface::m_materialLayer {}
protectedinherited

Possibility to attach a material description to the surface.

In this case the surface is what holds the MaterialLayer. This is usually done for boundary surfaces so shared ptr.

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

453{};

◆ m_owner

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

enum for surface owner : default free surface

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

◆ m_referencePoint

CxxUtils::CachedUniquePtr<Amg::Vector3D> Trk::DiscSurface::m_referencePoint
protectedinherited

static member for boundless approach

Definition at line 286 of file DiscSurface.h.

◆ m_shared

bool Trk::SubtractedDiscSurface::m_shared {true}
protected

Definition at line 76 of file SubtractedDiscSurface.h.

76{true};

◆ m_subtrVol

std::shared_ptr<const AreaExcluder> Trk::SubtractedDiscSurface::m_subtrVol {nullptr}
protected

Definition at line 75 of file SubtractedDiscSurface.h.

75{nullptr};

◆ m_transforms

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

Unique Pointer to the Transforms struct.

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

439{};

◆ s_boundless

const Trk::NoBounds Trk::DiscSurface::s_boundless
staticprotectedinherited

Definition at line 288 of file DiscSurface.h.

◆ s_numberOfInstantiations

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

Definition at line 22 of file TrkObjectCounter.h.

◆ s_onSurfaceTolerance

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

SurfaceType Trk::DiscSurface::staticType = SurfaceType::Disc
staticconstexprinherited

Definition at line 57 of file DiscSurface.h.


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