SubtractedCylinderLayer.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/

// SubtractedCylinderLayer.cxx, (c) ATLAS Detector software
 
// Trk
#include "TrkGeometry/SubtractedCylinderLayer.h"
 
#include "TrkGeometry/LayerMaterialProperties.h"
#include "TrkGeometry/MaterialProperties.h"
#include "TrkParameters/TrackParameters.h"
// CLHEP
#include "GeoPrimitives/GeoPrimitives.h"
 
Trk::SubtractedCylinderLayer::SubtractedCylinderLayer(
    const Trk::SubtractedCylinderSurface* subCyl,
    const Trk::LayerMaterialProperties& laymatprop, double thickness,
    std::unique_ptr<Trk::OverlapDescriptor> olap, int laytyp)
    : SubtractedCylinderSurface(*subCyl),
      Layer(laymatprop, thickness, std::move(olap), laytyp) {}
 
Trk::SubtractedCylinderLayer::SubtractedCylinderLayer(
    const Trk::SubtractedCylinderLayer& clay) = default;
 
Trk::SubtractedCylinderLayer::SubtractedCylinderLayer(
    const Trk::SubtractedCylinderLayer& clay, const Amg::Transform3D& transf)
    : SubtractedCylinderSurface(clay, transf), Layer(clay) {}
 
Trk::SubtractedCylinderLayer& Trk::SubtractedCylinderLayer::operator=(
    const SubtractedCylinderLayer& clay) {
  if (this != &clay) {
    // call the assignments of the base classes
    Trk::SubtractedCylinderSurface::operator=(clay);
    Trk::Layer::operator=(clay);
  }
  return (*this);
}
 
const Trk::SubtractedCylinderSurface&
Trk::SubtractedCylinderLayer::surfaceRepresentation() const
{
  return (*this);
}
 
Trk::SubtractedCylinderSurface&
Trk::SubtractedCylinderLayer::surfaceRepresentation()
{
  return (*this);
}
 
double Trk::SubtractedCylinderLayer::preUpdateMaterialFactor(
    const Trk::TrackParameters& parm, Trk::PropDirection dir) const {
  if (!Trk::Layer::m_layerMaterialProperties) return 0.;
  // calculate the direction to the normal
  const Amg::Vector3D& parmPos = parm.position();
  Amg::Vector3D pastStep(parmPos + dir * parm.momentum().normalized());
  if (pastStep.perp() > parm.position().perp())
    return Trk::Layer::m_layerMaterialProperties->alongPreFactor();
  return Trk::Layer::m_layerMaterialProperties->oppositePreFactor();
}
 
double Trk::SubtractedCylinderLayer::postUpdateMaterialFactor(
    const Trk::TrackParameters& parm, Trk::PropDirection dir) const {
  if (!Trk::Layer::m_layerMaterialProperties) return 0;
  const Amg::Vector3D& parmPos = parm.position();
  Amg::Vector3D pastStep(parmPos + dir * parm.momentum().normalized());
  if (pastStep.perp() > parm.position().perp())
    return Trk::Layer::m_layerMaterialProperties->alongPostFactor();
  return Trk::Layer::m_layerMaterialProperties->oppositePostFactor();
}
 
void Trk::SubtractedCylinderLayer::moveLayer(Amg::Transform3D& shift) {
  m_transforms =
      std::make_unique<Transforms>(shift * (m_transforms->transform));
}