AnnulusBounds.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 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/Surfaces/DiscBounds.hpp"
#include "Acts/Utilities/Definitions.hpp"
#include "Acts/Utilities/ParameterDefinitions.hpp"
#include "Acts/Utilities/detail/periodic.hpp"

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

namespace Acts {

class AnnulusBounds : public DiscBounds {
 public:
  using Transform2D = Eigen::Affine2d;

  enum BoundValues : int {
    eMinR = 0,
    eMaxR = 1,
    eMinPhiRel = 2,
    eMaxPhiRel = 3,
    eAveragePhi = 4,
    eOriginX = 5,
    eOriginY = 6,
    eSize = 7
  };

  AnnulusBounds() = delete;

  AnnulusBounds(double minR, double maxR, double minPhiRel, double maxPhiRel,
                const Vector2D& moduleOrigin = {0, 0},
                double avgPhi = 0) noexcept(false)
      : AnnulusBounds({minR, maxR, minPhiRel, maxPhiRel, avgPhi,
                       moduleOrigin.x(), moduleOrigin.y()}) {}

  AnnulusBounds(const std::array<double, eSize>& values) noexcept(false);

  AnnulusBounds(const AnnulusBounds& source) = default;

  SurfaceBounds::BoundsType type() const final;

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

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

  virtual double distanceToBoundary(const Vector2D& lposition) const final;

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

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

  double phiMin() const;

  double phiMax() const;

  bool coversFullAzimuth() const final;

  bool insideRadialBounds(double R, double tolerance = 0.) const final;

  double binningValueR() const final;

  double binningValuePhi() const final;

  Vector2D moduleOrigin() const;

  std::vector<Vector2D> corners() const;

  std::vector<Vector2D> vertices(unsigned int lseg) const;

  double rMin() const final;

  double rMax() const final;

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

  // @TODO: Does this need to be in bound values?
  Vector2D m_moduleOrigin;
  Vector2D m_shiftXY;  // == -m_moduleOrigin
  Vector2D m_shiftPC;
  double m_phiAvg;
  Transform2D m_rotationStripPC;
  Transform2D m_translation;

  // Vectors needed for inside checking
  Vector2D m_outLeftStripPC;
  Vector2D m_inLeftStripPC;
  Vector2D m_outRightStripPC;
  Vector2D m_inRightStripPC;

  Vector2D m_outLeftModulePC;
  Vector2D m_inLeftModulePC;
  Vector2D m_outRightModulePC;
  Vector2D m_inRightModulePC;

  Vector2D m_outLeftStripXY;
  Vector2D m_inLeftStripXY;
  Vector2D m_outRightStripXY;
  Vector2D m_inRightStripXY;

  void checkConsistency() noexcept(false);

  virtual bool inside(const Vector2D& lposition, double tolR,
                      double tolPhi) const final;

  Vector2D stripXYToModulePC(const Vector2D& vStripXY) const;

  Vector2D closestOnSegment(const Vector2D& a, const Vector2D& b,
                            const Vector2D& p,
                            const Eigen::Matrix<double, 2, 2>& weight) const;

  double squaredNorm(const Vector2D& v,
                     const Eigen::Matrix<double, 2, 2>& weight) const;
};

inline SurfaceBounds::BoundsType AnnulusBounds::type() const {
  return SurfaceBounds::eAnnulus;
}

inline double AnnulusBounds::rMin() const {
  return get(eMinR);
}

inline double AnnulusBounds::rMax() const {
  return get(eMaxR);
}

inline double AnnulusBounds::phiMin() const {
  return get(eMinPhiRel) + get(eAveragePhi);
}

inline double AnnulusBounds::phiMax() const {
  return get(eMaxPhiRel) + get(eAveragePhi);
}

inline bool AnnulusBounds::coversFullAzimuth() const {
  return (std::abs((get(eMinPhiRel) - get(eMaxPhiRel)) - M_PI) <
          s_onSurfaceTolerance);
}

inline bool AnnulusBounds::insideRadialBounds(double R,
                                              double tolerance) const {
  return ((R + tolerance) > get(eMinR) and (R - tolerance) < get(eMaxR));
}

inline double AnnulusBounds::binningValueR() const {
  return 0.5 * (get(eMinR) + get(eMaxR));
}

inline double AnnulusBounds::binningValuePhi() const {
  return get(eAveragePhi);
}

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

inline void AnnulusBounds::checkConsistency() noexcept(false) {
  if (get(eMinR) < 0. or get(eMaxR) < 0. or get(eMinR) > get(eMaxR) or
      std::abs(get(eMinR) - get(eMaxR)) < s_epsilon) {
    throw std::invalid_argument("AnnulusBounds: invalid radial setup.");
  }
  if (get(eMinPhiRel) != detail::radian_sym(get(eMinPhiRel)) or
      get(eMaxPhiRel) != detail::radian_sym(get(eMaxPhiRel)) or
      get(eMinPhiRel) > get(eMaxPhiRel)) {
    throw std::invalid_argument("AnnulusBounds: invalid phi boundary setup.");
  }
  if (get(eAveragePhi) != detail::radian_sym(get(eAveragePhi))) {
    throw std::invalid_argument("AnnulusBounds: invalid phi positioning.");
  }
}

}  // namespace Acts