#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 (ActsTrk::DetectorAlignStore &alignStore) const override
	Caches the aligned transformation in the provided store.
const Acts::Transform3 &	localToGlobalTransform (const Acts::GeometryContext &tgContext) const
const Acts::Transform3 &	localToGlobalTransform (const ActsTrk::GeometryContext &gctx) const override
	Returns the aligned local to global transform from the passed geometry context.
virtual const Acts::Transform3 &	localToGlobalTransform (const ActsTrk::DetectorAlignStore *store) const final override
	Returns the reference to the aligned local to global transform from the DetectorAlignStore actually holding the transform of the Detector element.
const Acts::Surface &	surface () const
	Return surface associated with this identifier, which should come from the.
Acts::Surface &	surface ()
	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.
Acts::Transform3	localToGlobal (const ActsTrk::DetectorAlignStore *store) const

Private Attributes
IdentifierHash	m_idHash {}
DetectorType	m_type {DetectorType::UnDefined}
ActsTrk::ReadoutSurfacePositioning< ActsDetectorElement >	m_trfCache {0, this}
std::shared_ptr< ActsTrk::SurfacePlacement >	m_surfHolder {}
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::unique_ptr< const Acts::Transform3 >	m_trtTrf {}
Identifier	m_explicitIdentifier {}

Detailed Description

Definition at line 58 of file ActsDetectorElement.h.

Member Typedef Documentation

◆ AlignmentStore

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

Definition at line 61 of file ActsDetectorElement.h.

◆ DetectorType

using ActsDetectorElement::DetectorType = ActsTrk::DetectorType

Definition at line 60 of file ActsDetectorElement.h.

Constructor & Destructor Documentation

◆ ActsDetectorElement() [1/3]

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

Definition at line 65 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<Acts::RectangleBounds>(hlX, hlY);
 
    m_bounds = rectangleBounds;
    m_surfHolder = ActsTrk::SurfacePlacement::makeShared<Acts::PlaneSurface>(m_trfCache, std::move(rectangleBounds));
    m_surfHolder->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<Acts::TrapezoidBounds>(minHlX, maxHlX, hlY);
 
    m_bounds = trapezoidBounds;
 
    m_surfHolder = ActsTrk::SurfacePlacement::makeShared<Acts::PlaneSurface>(m_trfCache, std::move(trapezoidBounds));
    m_surfHolder->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_surfHolder = ActsTrk::SurfacePlacement::makeShared<Acts::DiscSurface>(m_trfCache, std::move(annulusBounds));
    m_surfHolder->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 145 of file ActsDetectorElement.cxx.

                                                               :
    GeoVDetectorElement{detElem.getMaterialGeom()},
    m_idHash(detElem.identifyHash()),
    m_type{DetectorType::Trt},
    m_detElement{&detElem}, 
    m_trtTrf{std::make_unique<Acts::Transform3>(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<Acts::LineBounds>(innerTubeRadius, length);
  m_bounds = lineBounds;
 
  m_surfHolder = ActsTrk::SurfacePlacement::makeShared<Acts::StrawSurface>(m_trfCache, std::move(lineBounds));
  m_surfHolder->surface().assignThickness(thickness());
}

◆ ActsDetectorElement() [3/3]

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

Constructor for an HGTD surface.

Definition at line 182 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<Acts::RectangleBounds>(hlX, hlY);
 
    m_bounds = rectangleBounds;
 
    m_surfHolder = ActsTrk::SurfacePlacement::makeShared<Acts::PlaneSurface>(m_trfCache, std::move(rectangleBounds));
    m_surfHolder->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 290 of file ActsDetectorElement.cxx.

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

◆ detectorType()

DetectorType ActsDetectorElement::detectorType ( ) const

finaloverridevirtual

Detector type.

Implements ActsTrk::IDetectorElement.

Definition at line 310 of file ActsDetectorElement.cxx.

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

                                                                      {
    return m_trfCache.getTransform(nullptr);
}

◆ identify()

Identifier ActsDetectorElement::identify ( ) const

finaloverridevirtual

Identifier.

Implements ActsTrk::IDetectorElement.

Definition at line 302 of file ActsDetectorElement.cxx.

                                               {
   return m_explicitIdentifier;
}

◆ identifyHash()

IdentifierHash ActsDetectorElement::identifyHash ( ) const

inline

Identifier hash.

Definition at line 89 of file ActsDetectorElement.h.

89{ return m_idHash; }

◆ identityHelper()

IdentityHelper ActsDetectorElement::identityHelper ( ) const

Definition at line 253 of file ActsDetectorElement.cxx.

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

◆ localToGlobal()

Acts::Transform3 ActsDetectorElement::localToGlobal ( const ActsTrk::DetectorAlignStore * store ) const

Pixel or Sct

Definition at line 212 of file ActsDetectorElement.cxx.

                                                                                              {
   
    GeoAlignmentStore* geoModelStore = store ? store->geoModelAlignment.get() : nullptr;
    Acts::Transform3 l2g{Acts::Transform3::Identity()};
   switch (m_type) {
      case DetectorType::Hgtd:{
         l2g= Amg::toIsometry3D(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
          Acts::Transform3 extraTransform{Amg::CLHEPTransformToEigen(detElem.recoToHitTransform())};    
          if (boundsType == Trk::SurfaceBounds::Rectangle &&
              typeid(design) == typeid(InDetDD::StripBoxDesign) ) {            
              extraTransform = Amg::toIsometry3D(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&>(surface().bounds()).moduleOrigin();
              const Acts::Transform3 originTrf{Amg::getTranslate3D(origin2D.x(), origin2D.y(), 0) * Amg::getRotateZ3D(-phiShift)};
              extraTransform = extraTransform * originTrf.inverse();
          }
          l2g = Amg::toIsometry3D(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() [1/3]

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

Definition at line 267 of file ActsDetectorElement.cxx.

                                                                                                          {
    return m_trfCache.getTransform(anygctx);
}

◆ localToGlobalTransform() [2/3]

const Acts::Transform3 & ActsDetectorElement::localToGlobalTransform ( const ActsTrk::DetectorAlignStore * store ) const

finaloverridevirtual

Returns the reference to the aligned local to global transform from the DetectorAlignStore actually holding the transform of the Detector element.

Parameters

store Pointer to the alignment store. If not provided the nominal transform is returned

Implements ActsTrk::IDetectorElement.

Definition at line 270 of file ActsDetectorElement.cxx.

                                                                                                              {
    return m_trfCache.getTransform(store);
}

◆ localToGlobalTransform() [3/3]

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

overridevirtual

Returns the aligned local to global transform from the passed geometry context.

Parameters

gctx	Reference to the ATLAS geometry context

Implements ActsTrk::IDetectorElement.

Definition at line 264 of file ActsDetectorElement.cxx.

                                                                                                          {
    return m_trfCache.getTransform(gctx);
}

◆ storeAlignedTransforms()

unsigned int ActsDetectorElement::storeAlignedTransforms ( 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 274 of file ActsDetectorElement.cxx.

                                                                                               {
    return m_trfCache.storeTransform(store);
}

◆ surface() [1/2]

Acts::Surface & ActsDetectorElement::surface ( )

Mutable surface to this detector element.

Definition at line 286 of file ActsDetectorElement.cxx.

                                          {
  return m_surfHolder->surface();
}

◆ surface() [2/2]

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

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

Definition at line 282 of file ActsDetectorElement.cxx.

                                                      {
  return m_surfHolder->surface();
}

◆ thickness()

double ActsDetectorElement::thickness ( ) const

Returns the thickness of the module.

Definition at line 300 of file ActsDetectorElement.cxx.

300{ return m_thickness; }

◆ upstreamDetectorElement()

const GeoVDetectorElement * ActsDetectorElement::upstreamDetectorElement ( ) const

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

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

133{};

◆ m_detElement

const GeoVDetectorElement* ActsDetectorElement::m_detElement {nullptr}

private

Detector element as variant.

Definition at line 131 of file ActsDetectorElement.h.

131{nullptr};

◆ m_explicitIdentifier

Identifier ActsDetectorElement::m_explicitIdentifier {}

private

Definition at line 139 of file ActsDetectorElement.h.

139{};

◆ m_idHash

IdentifierHash ActsDetectorElement::m_idHash {}

private

Definition at line 125 of file ActsDetectorElement.h.

125{};

◆ m_surfHolder

std::shared_ptr<ActsTrk::SurfacePlacement> ActsDetectorElement::m_surfHolder {}

private

Definition at line 129 of file ActsDetectorElement.h.

129{};

◆ m_thickness

double ActsDetectorElement::m_thickness {0.}

private

Thickness of this detector element.

Definition at line 135 of file ActsDetectorElement.h.

135{0.};

◆ m_trfCache

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

private

Definition at line 127 of file ActsDetectorElement.h.

127{0, this};

◆ m_trtTrf

std::unique_ptr<const Acts::Transform3> ActsDetectorElement::m_trtTrf {}

private

Definition at line 137 of file ActsDetectorElement.h.

137{};

◆ m_type

DetectorType ActsDetectorElement::m_type {DetectorType::UnDefined}

private

Definition at line 126 of file ActsDetectorElement.h.

126{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()

◆ localToGlobal()

◆ localToGlobalTransform() [1/3]

◆ localToGlobalTransform() [2/3]

◆ localToGlobalTransform() [3/3]

◆ storeAlignedTransforms()

◆ surface() [1/2]

◆ surface() [2/2]

◆ thickness()

◆ upstreamDetectorElement()

Member Data Documentation

◆ m_bounds

◆ m_detElement

◆ m_explicitIdentifier

◆ m_idHash

◆ m_surfHolder

◆ m_thickness

◆ m_trfCache

◆ m_trtTrf

◆ m_type