#include <ActsDetectorElement.h>

Inheritance diagram for ActsDetectorElement:

Collaboration diagram for ActsDetectorElement:

Public Types
using	DetectorType = ActsTrk::DetectorType
using	AlignmentStore = ActsTrk::GeometryContext::AlignmentStore

Public Member Functions
	ActsDetectorElement (const InDetDD::SiDetectorElement &detElem)
	ActsDetectorElement (const Acts::Transform3 &trf, const InDetDD::TRT_BaseElement &detElem, const Identifier &id)
	Constructor for a straw surface.
	ActsDetectorElement (const InDetDD::HGTD_DetectorElement &detElem, const Identifier &id)
	Constructor for an HGTD surface.
virtual	~ActsDetectorElement ()=default
	Destructor.
Identifier	identify () const override final
	Identifier.
DetectorType	detectorType () const override final
	Detector type.
IdentifierHash	identifyHash () const
	Identifier hash.
virtual unsigned int	storeAlignedTransforms (const ActsTrk::DetectorAlignStore &alignStore) const override
	Caches the aligned transformation in the provided store. Returns the number of cached elements.
virtual const Acts::Transform3 &	localToGlobalTransform (const Acts::GeometryContext &gctx) const final override
virtual const Acts::Surface &	surface () const final override
	Return surface associated with this identifier, which should come from the.
virtual bool	isSensitive () const final override
	Returns whether the detector element is sensitive.
virtual Acts::Surface &	surface () final override
	Mutable surface to this detector element.
const Trk::Surface &	atlasSurface () const
	Return a shared pointer on the ATLAS surface associated with this identifier,.
double	thickness () const
	Returns the thickness of the module.
IdentityHelper	identityHelper () const
const Acts::Transform3 &	getDefaultTransform () const
	Returns default transform.
const GeoVDetectorElement *	upstreamDetectorElement () const
	Returns the underllying GeoModel detectorelement that this one is based on.
Amg::Transform3D	localToGlobal (const ActsTrk::DetectorAlignStore *store) const

Private Attributes
IdentifierHash	m_idHash {}
DetectorType	m_type {DetectorType::UnDefined}
ActsTrk::TransformCacheDetEle< ActsDetectorElement >	m_trfCache {0, this}
const GeoVDetectorElement *	m_detElement {nullptr}
	Detector element as variant.
std::shared_ptr< const Acts::SurfaceBounds >	m_bounds {}
	Boundaries of the detector element.
double	m_thickness {0.}
	Thickness of this detector element.
std::shared_ptr< Acts::Surface >	m_surface {}
	Corresponding Surface.
std::unique_ptr< const Amg::Transform3D >	m_trtTrf {}
Identifier	m_explicitIdentifier {}

Detailed Description

Definition at line 44 of file ActsDetectorElement.h.

Member Typedef Documentation

◆ AlignmentStore

using ActsDetectorElement::AlignmentStore = ActsTrk::GeometryContext::AlignmentStore

Definition at line 47 of file ActsDetectorElement.h.

◆ DetectorType

using ActsDetectorElement::DetectorType = ActsTrk::DetectorType

Definition at line 46 of file ActsDetectorElement.h.

Constructor & Destructor Documentation

◆ ActsDetectorElement() [1/3]

ActsDetectorElement::ActsDetectorElement ( const InDetDD::SiDetectorElement & detElem )

Definition at line 53 of file ActsDetectorElement.cxx.

                                                                                :
  GeoVDetectorElement{detElem.getMaterialGeom()},
  m_idHash(detElem.identifyHash()),
  m_type{detElem.isPixel() ? DetectorType::Pixel : DetectorType::Sct},
  m_detElement{&detElem},
  m_explicitIdentifier(detElem.identify())
{
 
 
  auto boundsType = detElem.bounds().type();
 
  m_thickness = detElem.thickness();
 
 
  if (boundsType == Trk::SurfaceBounds::Rectangle) {
 
    const InDetDD::SiDetectorDesign &design = detElem.design();
    double hlX = design.width() / 2. * length_unit;
    double hlY = design.length() / 2. * length_unit;
 
    auto rectangleBounds = std::make_shared<const Acts::RectangleBounds>(hlX, hlY);
 
    m_bounds = rectangleBounds;
    m_surface = Acts::Surface::makeShared<Acts::PlaneSurface>(rectangleBounds, *this);
    m_surface->assignThickness(thickness());
 
  } else if (boundsType == Trk::SurfaceBounds::Trapezoid) {
 
    const InDetDD::SiDetectorDesign &design = detElem.design();
 
    double minHlX = design.minWidth() / 2. * length_unit;
    double maxHlX = design.maxWidth() / 2. * length_unit;
    double hlY = design.length() / 2. * length_unit;
 
    auto trapezoidBounds =
        std::make_shared<const Acts::TrapezoidBounds>(minHlX, maxHlX, hlY);
 
    m_bounds = trapezoidBounds;
 
    m_surface = Acts::Surface::makeShared<Acts::PlaneSurface>(trapezoidBounds, *this);
    m_surface->assignThickness(thickness());
 
 
  } else if (boundsType == Trk::SurfaceBounds::Annulus) {
 
    const InDetDD::SiDetectorDesign &design = detElem.design();
    const auto *annulus = dynamic_cast<const InDetDD::StripStereoAnnulusDesign *>(&design);
    if (annulus == nullptr) {
      throw std::domain_error("ActsDetectorElement got inconsistent surface");
    }
 
    double phi = annulus->phiWidth();
    double phiS = annulus->stereo();
    double R = annulus->waferCentreR();
    double maxR = annulus->maxR();
    double minR = annulus->minR();
 
    // phiAvg is the bounds-internal local rotation. We don't want one
    double phiAvg = 0; 
    // phi is the total opening angle, set up symmetric phi bounds
    double phiMax = phi / 2.;
    double phiMin = -phiMax;
 
 
    Amg::Vector2D originStripXYRotated(R * (1 - std::cos(phiS)),
                                       R * std::sin(-phiS));
 
    auto annulusBounds = std::make_shared<Acts::AnnulusBounds>(
        minR, maxR, phiMin, phiMax, originStripXYRotated, phiAvg);
    m_bounds = annulusBounds;
 
    m_surface = Acts::Surface::makeShared<Acts::DiscSurface>(annulusBounds, *this);
    m_surface->assignThickness(thickness());
 
  } else {
    std::cout << boundsType << std::endl;
    throw std::domain_error("ActsDetectorElement does not support this surface type");
  }
}

◆ ActsDetectorElement() [2/3]

ActsDetectorElement::ActsDetectorElement	(	const Acts::Transform3 &	trf,
		const InDetDD::TRT_BaseElement &	detElem,
		const Identifier &	id )

Constructor for a straw surface.

Parameters

transform Transform to the straw system

Definition at line 133 of file ActsDetectorElement.cxx.

                                                               :
    GeoVDetectorElement{detElem.getMaterialGeom()},
    m_idHash(detElem.identifyHash()),
    m_type{DetectorType::Trt},
    m_detElement{&detElem}, 
    m_trtTrf{std::make_unique<Amg::Transform3D>(trf)},
    m_explicitIdentifier(id)
{
 
 
  // we know this is a straw
  double length = detElem.strawLength() * 0.5 * length_unit;
 
  // we need to find the radius
  auto ecElem = dynamic_cast<const InDetDD::TRT_EndcapElement *>(&detElem);
  auto brlElem = dynamic_cast<const InDetDD::TRT_BarrelElement *>(&detElem);
  double innerTubeRadius{0.};
  if (ecElem) {
    innerTubeRadius = ecElem->getDescriptor()->innerTubeRadius() * length_unit;
  } else {
    if (brlElem) {
      innerTubeRadius =
          brlElem->getDescriptor()->innerTubeRadius() * length_unit;
    } else {
      THROW_EXCEPTION("Cannot get tube radius for element in ActsDetectorElement c'tor");
    }
  }
 
  auto lineBounds =
      std::make_shared<const Acts::LineBounds>(innerTubeRadius, length);
  m_bounds = lineBounds;
 
  m_surface = Acts::Surface::makeShared<Acts::StrawSurface>(lineBounds, *this);
  m_surface->assignThickness(thickness());
}

◆ ActsDetectorElement() [3/3]

ActsDetectorElement::ActsDetectorElement	(	const InDetDD::HGTD_DetectorElement &	detElem,
		const Identifier &	id )

Constructor for an HGTD surface.

Definition at line 171 of file ActsDetectorElement.cxx.

                                                                                                         :
    GeoVDetectorElement{detElem.getMaterialGeom()},
    m_idHash(detElem.identifyHash()),
    m_type{DetectorType::Hgtd},
    m_detElement{&detElem}, 
    m_thickness{detElem.thickness()}, 
    m_explicitIdentifier{id}
{
 
  auto boundsType = detElem.bounds().type();
 
  if (boundsType == Trk::SurfaceBounds::Rectangle) {
 
    const InDetDD::HGTD_ModuleDesign &design = detElem.design();
    double hlX = design.width() / 2. * length_unit;
    double hlY = design.length() / 2. * length_unit;
 
    auto rectangleBounds =
        std::make_shared<const Acts::RectangleBounds>(hlX, hlY);
 
    m_bounds = rectangleBounds;
 
    m_surface = Acts::Surface::makeShared<Acts::PlaneSurface>(rectangleBounds, *this);
    m_surface->assignThickness(thickness());
        
  } else {
    throw std::domain_error(
        "ActsDetectorElement: the surface type of HGTD is not does not Rectangle, it is wrong");
  }
}

◆ ~ActsDetectorElement()

virtual ActsDetectorElement::~ActsDetectorElement ( )

virtualdefault

Destructor.

Member Function Documentation

◆ atlasSurface()

const Trk::Surface & ActsDetectorElement::atlasSurface ( ) const

Return a shared pointer on the ATLAS surface associated with this identifier,.

Definition at line 274 of file ActsDetectorElement.cxx.

                                                          {
  if (const auto *detElem =
          dynamic_cast<const InDetDD::SiDetectorElement *>(m_detElement);
      detElem != nullptr) {
    return detElem->surface();
  } else {
    throw std::domain_error("Cannot get surface for TRT element");
  }
}

◆ detectorType()

DetectorType ActsDetectorElement::detectorType ( ) const

finaloverridevirtual

Detector type.

Implements ActsTrk::IDetectorElementBase.

Definition at line 294 of file ActsDetectorElement.cxx.

294{ return m_type; }

◆ getDefaultTransform()

const Acts::Transform3 & ActsDetectorElement::getDefaultTransform ( ) const

Returns default transform.

For TRT this is static and set in constructor. For silicon detectors it is calulated from GM, and stored. Thus the method is not const. The store is mutexed.

Definition at line 262 of file ActsDetectorElement.cxx.

                                                                      {
    return m_trfCache.getTransform(nullptr);
}

◆ identify()

Identifier ActsDetectorElement::identify ( ) const

finaloverridevirtual

Identifier.

Implements ActsTrk::IDetectorElementBase.

Definition at line 286 of file ActsDetectorElement.cxx.

                                               {
   return m_explicitIdentifier;
}

◆ identifyHash()

IdentifierHash ActsDetectorElement::identifyHash ( ) const

inline

Identifier hash.

Definition at line 75 of file ActsDetectorElement.h.

75{ return m_idHash; }

◆ identityHelper()

IdentityHelper ActsDetectorElement::identityHelper ( ) const

Definition at line 244 of file ActsDetectorElement.cxx.

                                                         {
  if (detectorType() == DetectorType::Pixel || detectorType() == DetectorType::Sct) {
        return IdentityHelper(static_cast<const InDetDD::SiDetectorElement *>(m_detElement));
  } else {
    throw std::domain_error("Cannot get IdentityHelper for TRT element");
  }
}

◆ isSensitive()

virtual bool ActsDetectorElement::isSensitive ( ) const

inlinefinaloverridevirtual

Returns whether the detector element is sensitive.

Definition at line 85 of file ActsDetectorElement.h.

85{ return true; }

◆ localToGlobal()

Amg::Transform3D ActsDetectorElement::localToGlobal ( const ActsTrk::DetectorAlignStore * store ) const

Pixel or Sct

Definition at line 202 of file ActsDetectorElement.cxx.

                                                                                              {
   
    GeoAlignmentStore* geoModelStore = store ? store->geoModelAlignment.get() : nullptr;
    Amg::Transform3D l2g{Amg::Transform3D::Identity()};
   switch (m_type) {
      case DetectorType::Hgtd:{
         l2g= m_detElement->getMaterialGeom()->getAbsoluteTransform(geoModelStore);
         break;
      } case DetectorType::Trt: {
         l2g = (*m_trtTrf);
         break;
      }
      default: {
          const auto& detElem = static_cast<const InDetDD::SiDetectorElement&>(*m_detElement);
          const InDetDD::SiDetectorDesign&design = detElem.design();
          const Trk::SurfaceBounds::BoundsType boundsType = detElem.bounds().type();
 
          // extra shift for split row modules
          Amg::Transform3D extraTransform{Amg::CLHEPTransformToEigen(detElem.recoToHitTransform())};    
          if (boundsType == Trk::SurfaceBounds::Rectangle &&
              typeid(design) == typeid(InDetDD::StripBoxDesign) ) {            
              extraTransform = design.moduleShift() * extraTransform;
          } else if (boundsType == Trk::SurfaceBounds::Annulus) {
              // need to rotate pi/2 to reproduce ABXY orientation, phiS so that phi=0
              // is center and symmetric
              const double phiShift = M_PI_2 - static_cast<const InDetDD::StripStereoAnnulusDesign&>(design).stereo();
 
              const Amg::Vector2D origin2D = static_cast<const Acts::AnnulusBounds&>(m_surface->bounds()).moduleOrigin();
              const Amg::Translation3D transl{origin2D.x(), origin2D.y(), 0};
              const Amg::Transform3D originTrf{transl * Amg::getRotateZ3D(-phiShift)};
              extraTransform = extraTransform * originTrf.inverse();
          }
          l2g = m_detElement->getMaterialGeom()->getAbsoluteTransform(geoModelStore) * extraTransform;
      }
   };
   // need to make sure translation has correct units
   l2g.translation() *= 1.0 / CLHEP::mm * length_unit;
 
   return l2g;
 
}

◆ localToGlobalTransform()

const Acts::Transform3 & ActsDetectorElement::localToGlobalTransform ( const Acts::GeometryContext & gctx ) const

finaloverridevirtual

Definition at line 252 of file ActsDetectorElement.cxx.

                                                                                                          {
    return m_trfCache.localToGlobalTransform(anygctx);
}

◆ storeAlignedTransforms()

unsigned int ActsDetectorElement::storeAlignedTransforms ( const ActsTrk::DetectorAlignStore & store ) const

overridevirtual

Caches the aligned transformation in the provided store. Returns the number of cached elements.

Implements ActsTrk::IDetectorElement.

Definition at line 256 of file ActsDetectorElement.cxx.

                                                                                                     {
    if (store.detType != detectorType()) return 0;
    m_trfCache.getTransform(&store);
    return 1;
}

◆ surface() [1/2]

const Acts::Surface & ActsDetectorElement::surface ( ) const

finaloverridevirtual

Return surface associated with this identifier, which should come from the.

Definition at line 266 of file ActsDetectorElement.cxx.

                                                      {
  return (*m_surface);
}

◆ surface() [2/2]

Acts::Surface & ActsDetectorElement::surface ( )

finaloverridevirtual

Mutable surface to this detector element.

Definition at line 270 of file ActsDetectorElement.cxx.

                                          {
  return (*m_surface);
}

◆ thickness()

double ActsDetectorElement::thickness ( ) const

Returns the thickness of the module.

Definition at line 284 of file ActsDetectorElement.cxx.

284{ return m_thickness; }

◆ upstreamDetectorElement()

const GeoVDetectorElement * ActsDetectorElement::upstreamDetectorElement ( ) const

Returns the underllying GeoModel detectorelement that this one is based on.

Definition at line 291 of file ActsDetectorElement.cxx.

                                                   {
  return m_detElement;
}

Member Data Documentation

◆ m_bounds

std::shared_ptr<const Acts::SurfaceBounds> ActsDetectorElement::m_bounds {}

private

Boundaries of the detector element.

Definition at line 116 of file ActsDetectorElement.h.

116{};

◆ m_detElement

const GeoVDetectorElement* ActsDetectorElement::m_detElement {nullptr}

private

Detector element as variant.

Definition at line 114 of file ActsDetectorElement.h.

114{nullptr};

◆ m_explicitIdentifier

Identifier ActsDetectorElement::m_explicitIdentifier {}

private

Definition at line 124 of file ActsDetectorElement.h.

124{};

◆ m_idHash

IdentifierHash ActsDetectorElement::m_idHash {}

private

Definition at line 110 of file ActsDetectorElement.h.

110{};

◆ m_surface

std::shared_ptr<Acts::Surface> ActsDetectorElement::m_surface {}

private

Corresponding Surface.

Definition at line 120 of file ActsDetectorElement.h.

120{};

◆ m_thickness

double ActsDetectorElement::m_thickness {0.}

private

Thickness of this detector element.

Definition at line 118 of file ActsDetectorElement.h.

118{0.};

◆ m_trfCache

ActsTrk::TransformCacheDetEle<ActsDetectorElement> ActsDetectorElement::m_trfCache {0, this}

private

Definition at line 112 of file ActsDetectorElement.h.

112{0, this};

◆ m_trtTrf

std::unique_ptr<const Amg::Transform3D> ActsDetectorElement::m_trtTrf {}

private

Definition at line 122 of file ActsDetectorElement.h.

122{};

◆ m_type

DetectorType ActsDetectorElement::m_type {DetectorType::UnDefined}

private

Definition at line 111 of file ActsDetectorElement.h.

111{DetectorType::UnDefined};

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

Public Types

Public Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ AlignmentStore

◆ DetectorType

Constructor & Destructor Documentation

◆ ActsDetectorElement() [1/3]

◆ ActsDetectorElement() [2/3]

◆ ActsDetectorElement() [3/3]

◆ ~ActsDetectorElement()

Member Function Documentation

◆ atlasSurface()

◆ detectorType()

◆ getDefaultTransform()

◆ identify()

◆ identifyHash()

◆ identityHelper()

◆ isSensitive()

◆ localToGlobal()

◆ localToGlobalTransform()

◆ storeAlignedTransforms()

◆ surface() [1/2]

◆ surface() [2/2]

◆ thickness()

◆ upstreamDetectorElement()

Member Data Documentation

◆ m_bounds

◆ m_detElement

◆ m_explicitIdentifier

◆ m_idHash

◆ m_surface

◆ m_thickness

◆ m_trfCache

◆ m_trtTrf

◆ m_type