CutoutCylinderVolumeBounds.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2019-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/Definitions/Algebra.hpp"
#include "Acts/Geometry/Volume.hpp"
#include "Acts/Geometry/VolumeBounds.hpp"

#include <array>
#include <iosfwd>
#include <memory>
#include <stdexcept>
#include <vector>

namespace Acts {

class CylinderBounds;
class DiscBounds;

class CutoutCylinderVolumeBounds : public VolumeBounds {
 public:
  enum BoundValues : int {
    eMinR = 0,
    eMedR = 1,
    eMaxR = 2,
    eHalfLengthZ = 3,
    eHalfLengthZcutout = 4,
    eSize
  };

  CutoutCylinderVolumeBounds() = delete;

  CutoutCylinderVolumeBounds(double rmin, double rmed, double rmax, double hlZ,
                             double hlZc) noexcept(false)
      : m_values({rmin, rmed, rmax, hlZ, hlZc}) {
    checkConsistency();
    buildSurfaceBounds();
  }

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

  ~CutoutCylinderVolumeBounds() override = default;

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

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

  bool inside(const Vector3& gpos, double tol = 0) const override;

  OrientedSurfaces orientedSurfaces(
      const Transform3& transform = Transform3::Identity()) const override;

  Volume::BoundingBox boundingBox(const Transform3* trf = nullptr,
                                  const Vector3& 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:
  std::array<double, eSize> m_values;

  // The surface bound objects
  std::shared_ptr<const CylinderBounds> m_innerCylinderBounds{nullptr};
  std::shared_ptr<const CylinderBounds> m_cutoutCylinderBounds{nullptr};
  std::shared_ptr<const CylinderBounds> m_outerCylinderBounds{nullptr};
  std::shared_ptr<const DiscBounds> m_outerDiscBounds{nullptr};
  std::shared_ptr<const DiscBounds> m_innerDiscBounds{nullptr};

  void buildSurfaceBounds();

  void checkConsistency() noexcept(false);
};

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

inline void CutoutCylinderVolumeBounds::checkConsistency() noexcept(false) {
  if (get(eMinR) < 0. or get(eMedR) <= 0. or get(eMaxR) <= 0. or
      get(eMinR) >= get(eMedR) or get(eMinR) >= get(eMaxR) or
      get(eMedR) >= get(eMaxR)) {
    throw std::invalid_argument(
        "CutoutCylinderVolumeBounds: invalid radial input.");
  }
  if (get(eHalfLengthZ) <= 0 or get(eHalfLengthZcutout) <= 0. or
      get(eHalfLengthZcutout) > get(eHalfLengthZ)) {
    throw std::invalid_argument(
        "CutoutCylinderVolumeBounds: invalid longitudinal input.");
  }
}

}  // namespace Acts