ActsDetectorElement Class Reference

#include <ActsDetectorElement.h>

Inheritance diagram for ActsDetectorElement:

Collaboration diagram for ActsDetectorElement:

Public Types
using	DetectorType = ActsTrk::DetectorType
using	AlignmentStore = ActsGeometryContext::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. More...
	ActsDetectorElement (const InDetDD::HGTD_DetectorElement &detElem, const Identifier &id)
	Constructor for an HGTD surface. More...
virtual	~ActsDetectorElement ()=default
	Destructor. More...
Identifier	identify () const override final
	Identifier. More...
DetectorType	detectorType () const override final
	Detector type. More...
virtual unsigned int	storeAlignedTransforms (const ActsTrk::DetectorAlignStore &alignStore) const override
	Caches the aligned transformation in the provided store. Returns the number of cached elements. More...
virtual const Acts::Transform3 &	transform (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. More...
virtual Acts::Surface &	surface () final override
	Mutable surface to this detector element. More...
const Trk::Surface &	atlasSurface () const
	Return a shared pointer on the ATLAS surface associated with this identifier,. More...
virtual double	thickness () const final override
	Returns the thickness of the module. More...
IdentityHelper	identityHelper () const
const Acts::Transform3 &	getDefaultTransform () const
	Returns default transform. More...
const GeoVDetectorElement *	upstreamDetectorElement () const
	Returns the underllying GeoModel detectorelement that this one is based on. More...
Amg::Transform3D	transform (const ActsTrk::DetectorAlignStore *store) const

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

Detailed Description

Definition at line 42 of file ActsDetectorElement.h.

Member Typedef Documentation

AlignmentStore

using ActsDetectorElement::AlignmentStore = ActsGeometryContext::AlignmentStore

Definition at line 45 of file ActsDetectorElement.h.

DetectorType

using ActsDetectorElement::DetectorType = ActsTrk::DetectorType

Definition at line 44 of file ActsDetectorElement.h.

Constructor & Destructor Documentation

ActsDetectorElement() [1/3]

ActsDetectorElement::ActsDetectorElement

(

const InDetDD::SiDetectorElement &

detElem

)

Definition at line 51 of file ActsDetectorElement.cxx.

                                                                                :
  GeoVDetectorElement{detElem.getMaterialGeom()},
  m_type{detElem.isPixel() ? DetectorType::Pixel : DetectorType::Sct},
  m_detElement{&detElem} {
 
 
  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);
 
  } 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);
 
 
  } 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);
 
  } 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 125 of file ActsDetectorElement.cxx.

                                                               :
    GeoVDetectorElement{detElem.getMaterialGeom()},
    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);
}

ActsDetectorElement() [3/3]

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

Constructor for an HGTD surface.

Definition at line 159 of file ActsDetectorElement.cxx.

                                                                                                         :
    GeoVDetectorElement{detElem.getMaterialGeom()},
    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);
        
  } 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 259 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::IDetectorElement.

Definition at line 291 of file ActsDetectorElement.cxx.

{ 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 247 of file ActsDetectorElement.cxx.

                                                                      {
    return m_trfCache.getTransform(nullptr);
}

identify()

Identifier ActsDetectorElement::identify ( ) const

finaloverridevirtual

Identifier.

Implements ActsTrk::IDetectorElement.

Definition at line 271 of file ActsDetectorElement.cxx.

                                               {
  if (const auto *detElem =
          dynamic_cast<const InDetDD::SiDetectorElement *>(m_detElement);
      detElem != nullptr) {
    return detElem->identify();
  } else if (dynamic_cast<const InDetDD::TRT_BaseElement *>(m_detElement) !=
             nullptr) {
    return m_explicitIdentifier;
  } else if (dynamic_cast<const InDetDD::HGTD_DetectorElement *>(m_detElement) !=
             nullptr) {
    return m_explicitIdentifier;
  } else {
    THROW_EXCEPTION("Unknown detector element type");
  }
}

identityHelper()

IdentityHelper ActsDetectorElement::identityHelper

(

)

const

Definition at line 229 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");
  }
}

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 241 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 251 of file ActsDetectorElement.cxx.

                                                      {
  return (*m_surface);
}

surface() [2/2]

Acts::Surface & ActsDetectorElement::surface ( )

finaloverridevirtual

Mutable surface to this detector element.

Definition at line 255 of file ActsDetectorElement.cxx.

                                          {
  return (*m_surface);
}

thickness()

double ActsDetectorElement::thickness ( ) const

finaloverridevirtual

Returns the thickness of the module.

Definition at line 269 of file ActsDetectorElement.cxx.

{ return m_thickness; }

transform() [1/2]

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

finaloverridevirtual

Definition at line 237 of file ActsDetectorElement.cxx.

                                                                                             {
    return m_trfCache.transform(anygctx);
}

transform() [2/2]

Amg::Transform3D ActsDetectorElement::transform

(

const ActsTrk::DetectorAlignStore *

store

)

const

Pixel or Sct

Definition at line 187 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;
 
}

upstreamDetectorElement()

const GeoVDetectorElement * ActsDetectorElement::upstreamDetectorElement

(

)

const

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

Definition at line 288 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 109 of file ActsDetectorElement.h.

m_detElement

const GeoVDetectorElement* ActsDetectorElement::m_detElement {nullptr}

private

Detector element as variant.

Definition at line 107 of file ActsDetectorElement.h.

m_explicitIdentifier

Identifier ActsDetectorElement::m_explicitIdentifier {}

private

Definition at line 118 of file ActsDetectorElement.h.

m_surface

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

private

Corresponding Surface.

Definition at line 113 of file ActsDetectorElement.h.

m_surfaces

std::vector<std::shared_ptr<const Acts::Surface> > ActsDetectorElement::m_surfaces {}

private

Definition at line 114 of file ActsDetectorElement.h.

m_thickness

double ActsDetectorElement::m_thickness {0.}

private

Thickness of this detector element.

Definition at line 111 of file ActsDetectorElement.h.

m_trfCache

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

private

Definition at line 105 of file ActsDetectorElement.h.

m_trtTrf

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

private

Definition at line 116 of file ActsDetectorElement.h.

m_type

DetectorType ActsDetectorElement::m_type {DetectorType::UnDefined}

private

Definition at line 103 of file ActsDetectorElement.h.

---

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

• ActsDetectorElement.h
• ActsDetectorElement.cxx