TrapezoidBounds.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file TrapezoidBounds.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/Definitions/Algebra.hpp"
#include "Acts/Definitions/TrackParametrization.hpp"
#include "Acts/Surfaces/BoundaryCheck.hpp"
#include "Acts/Surfaces/PlanarBounds.hpp"
#include "Acts/Surfaces/RectangleBounds.hpp"
#include "Acts/Surfaces/SurfaceBounds.hpp"

#include <algorithm>
#include <array>
#include <cmath>
#include <iosfwd>
#include <stdexcept>
#include <vector>

namespace Acts {


class TrapezoidBounds : public PlanarBounds {
 public:
  enum BoundValues {
    eHalfLengthXnegY = 0,
    eHalfLengthXposY = 1,
    eHalfLengthY = 2,
    eSize = 3
  };

  TrapezoidBounds() = delete;

  TrapezoidBounds(double halfXnegY, double halfXposY,
                  double halfY) noexcept(false)
      : m_values({halfXnegY, halfXposY, halfY}),
        m_boundingBox(std::max(halfXnegY, halfXposY), halfY) {
    checkConsistency();
  }

  TrapezoidBounds(const std::array<double, eSize>& values) noexcept(false)
      : m_values(values),
        m_boundingBox(
            std::max(values[eHalfLengthXnegY], values[eHalfLengthXposY]),
            values[eHalfLengthY]) {
    checkConsistency();
  }

  ~TrapezoidBounds() override;

  BoundsType type() const final;

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

  bool inside(const Vector2& lposition,
              const BoundaryCheck& bcheck) const final;

  std::vector<Vector2> vertices(unsigned int lseg = 1) const final;

  // Bounding box representation
  const RectangleBounds& boundingBox() const final;

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

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

 private:
  std::array<double, eSize> m_values;
  RectangleBounds m_boundingBox;

  void checkConsistency() noexcept(false);
};

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

inline void TrapezoidBounds::checkConsistency() noexcept(false) {
  if (get(eHalfLengthXnegY) <= 0. or get(eHalfLengthXposY) <= 0.) {
    throw std::invalid_argument("TrapezoidBounds: invalid local x setup");
  }
  if (get(eHalfLengthY) <= 0.) {
    throw std::invalid_argument("TrapezoidBounds: invalid local y setup");
  }
}

}  // namespace Acts