Volume.cxx

Go to the documentation of this file.

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

// Volume.cxx, (c) ATLAS Detector software
 
// Trk
#include "TrkVolumes/Volume.h"
#include "TrkDetDescrUtils/ObjectAccessor.h"
#include "TrkVolumes/VolumeBounds.h"
//CxxUtils
#include "CxxUtils/inline_hints.h"
// Gaudi
#include "GaudiKernel/MsgStream.h"
// STD
#include <iostream>
 
 
// constructor with Amg::Transform3D
Trk::Volume::Volume(std::unique_ptr<Amg::Transform3D> htrans,
                    std::shared_ptr<Trk::VolumeBounds> volbounds)
  : m_transform(std::move(htrans))
  , m_center(m_transform ? m_transform->translation() : Trk::s_origin)
  , m_volumeBounds(std::move(volbounds))
{}
 
// copy constructor
Trk::Volume::Volume(const Trk::Volume& vol)
    : m_transform(vol.m_transform
                      ? std::make_unique<Amg::Transform3D>(*(vol.m_transform))
                      : nullptr),
      m_center(vol.m_center),
      m_volumeBounds(vol.m_volumeBounds) {}
 
// copy constructor with shift
Trk::Volume::Volume(const Trk::Volume& vol, const Amg::Transform3D& shift)
    : m_transform(vol.m_transform ? std::make_unique<Amg::Transform3D>(
                                        shift * (*(vol.m_transform)))
                                  : std::make_unique<Amg::Transform3D>(shift)),
      m_center(shift * vol.m_center),
      m_volumeBounds(vol.m_volumeBounds) {}
 
// assignment operator
Trk::Volume&
Trk::Volume::operator=(const Trk::Volume& vol)
{
  if (this != &vol) {
    m_transform = vol.m_transform
                    ? std::make_unique<Amg::Transform3D>(*vol.m_transform)
                    : nullptr;
    m_center = vol.m_center;
    m_volumeBounds = vol.m_volumeBounds;
  }
  return *this;
}
 
Trk::Volume*
Trk::Volume::clone() const
{
  return new Trk::Volume(*this);
}
 
// We compile this package with optimization, even in debug builds; otherwise,
// the heavy use of Eigen makes it too slow.  However, from here we may call
// to out-of-line Eigen code that is linked from other DSOs; in that case,
// it would not be optimized.  Avoid this by forcing all Eigen code
// to be inlined here if possible.
ATH_FLATTEN
bool
Trk::Volume::inside(const Amg::Vector3D& gp, double tol) const
{
  if (!m_transform) {
    return (volumeBounds()).inside(gp, tol);
  }
  Amg::Vector3D posInVolFrame((transform().inverse()) * gp);
  return (volumeBounds()).inside(posInVolFrame, tol);
}
 
Trk::ObjectAccessor
Trk::Volume::boundarySurfaceAccessor(
  const Amg::Vector3D& gp,
  const Amg::Vector3D& dir,
  bool forceInside) const
{
  if (!m_transform) {
    return Trk::ObjectAccessor(
      volumeBounds().boundarySurfaceAccessor(gp, dir, forceInside));
  }
  return Trk::ObjectAccessor(volumeBounds().boundarySurfaceAccessor(
    transform().inverse() * gp, dir, forceInside));
}

MsgStream&
Trk::operator<<(MsgStream& sl, const Trk::Volume& vol)
{
  sl << "Trk::Volume with VolumeBounds :" << vol.volumeBounds() << endmsg;
  return sl;
}
 
std::ostream&
Trk::operator<<(std::ostream& sl, const Trk::Volume& vol)
{
  sl << "Trk::Volume with VolumeBounds :" << vol.volumeBounds() << std::endl;
  return sl;
}