CuboidVolumeBounds.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file CuboidVolumeBounds.hpp

// This file is part of the Acts project.
//
// Copyright (C) 2016-2020 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.

#pragma once

#include "Acts/Geometry/Volume.hpp"
#include "Acts/Geometry/VolumeBounds.hpp"
#include "Acts/Utilities/BoundingBox.hpp"
#include "Acts/Utilities/Definitions.hpp"

#include <array>
#include <cmath>
#include <exception>
#include <vector>

namespace Acts {

class RectangleBounds;
class Volume;
class Surface;


class CuboidVolumeBounds : public VolumeBounds {
 public:
  enum BoundValues : int {
    eHalfLengthX = 0,
    eHalfLengthY = 1,
    eHalfLengthZ = 2,
    eSize = 3
  };

  CuboidVolumeBounds() = delete;

  CuboidVolumeBounds(double halex, double haley, double halez) noexcept(false);

  CuboidVolumeBounds(const std::array<double, eSize>& values) noexcept(false)
      : m_values(values) {
    checkConsistency();
  }

  CuboidVolumeBounds(const CuboidVolumeBounds& bobo);

  CuboidVolumeBounds& operator=(const CuboidVolumeBounds& bobo);

  ~CuboidVolumeBounds() override = default;

  VolumeBounds::BoundsType type() const final { return VolumeBounds::eCuboid; }

  std::vector<double> values() const final;

  bool inside(const Vector3D& pos, double tol = 0.) const override;

  SurfacePtrVector decomposeToSurfaces(
      const Transform3D* transformPtr) const override;

  Volume::BoundingBox boundingBox(const Transform3D* trf = nullptr,
                                  const Vector3D& envelope = {0, 0, 0},
                                  const Volume* entity = nullptr) const final;

  std::ostream& toStream(std::ostream& sl) const override;

  double get(BoundValues bValue) const { return m_values[bValue]; }

 private:
  template <class T>
  T& dumpT(T& dt) const;

  std::array<double, eSize> m_values;

  std::shared_ptr<const RectangleBounds> m_xyBounds = nullptr;
  std::shared_ptr<const RectangleBounds> m_yzBounds = nullptr;
  std::shared_ptr<const RectangleBounds> m_zxBounds = nullptr;

  void checkConsistency() noexcept(false);
};

inline bool CuboidVolumeBounds::inside(const Vector3D& pos, double tol) const {
  return (std::abs(pos.x()) <= get(eHalfLengthX) + tol &&
          std::abs(pos.y()) <= get(eHalfLengthY) + tol &&
          std::abs(pos.z()) <= get(eHalfLengthZ) + tol);
}

inline std::vector<double> CuboidVolumeBounds::values() const {
  std::vector<double> valvector;
  valvector.insert(valvector.begin(), m_values.begin(), m_values.end());
  return valvector;
}

inline void CuboidVolumeBounds::checkConsistency() noexcept(false) {
  if (get(eHalfLengthX) <= 0 or get(eHalfLengthY) <= 0 or
      get(eHalfLengthZ) <= 0.) {
    throw std::invalid_argument(
        "CuboidVolumeBounds: invalid input, zero or negative.");
  }
}

template <class T>
T& CuboidVolumeBounds::dumpT(T& dt) const {
  dt << std::setiosflags(std::ios::fixed);
  dt << std::setprecision(5);
  dt << "Acts::CuboidVolumeBounds: (halfLengthX, halfLengthY, halfLengthZ) = ";
  dt << "(" << get(eHalfLengthX) << ", " << get(eHalfLengthY) << ", "
     << get(eHalfLengthZ) << ")";
  return dt;
}
}  // namespace Acts