AnnulusBoundsPC.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRKSURFACES_ANNULUSBOUNDSPC_H
#define TRKSURFACES_ANNULUSBOUNDSPC_H
 
#include "TrkSurfaces/SurfaceBounds.h"
#include "GeoPrimitives/GeoPrimitives.h"
 
#include <iosfwd>
 
class BoundaryCheck;
class MsgStream;
 
#ifdef TRKDETDESCR_USEFLOATPRECISON
typedef float TDD_real_t;
#else
typedef double TDD_real_t;
#endif
 
class AnnulusTestAlg;
 
namespace Trk {
 
class AnnulusBounds;

class AnnulusBoundsPC final : public SurfaceBounds {
 
  using Rotation2D = Eigen::Rotation2D<double>;
  using Transform2D = Eigen::Transform<double, 2, Eigen::Affine>;
  using Matrix2D = Eigen::Matrix<double, 2, 2>;
 
  friend AnnulusTestAlg;
 
public:
 
  enum BoundValues {
    bv_minR,
    bv_maxR,
    bv_phiMin,
    bv_phiMax,
    bv_phiAvg,
    bv_originX,
    bv_originY
  };

  AnnulusBoundsPC(double minR,
                  double maxR,
                  double phiMin,
                  double phiMax,
                  Amg::Vector2D moduleOrigin = {0, 0},
                  double phiAvg = 0);

  AnnulusBoundsPC(const AnnulusBoundsPC& annbo) = default;
  AnnulusBoundsPC& operator=(const AnnulusBoundsPC& sbo) = default;
  AnnulusBoundsPC(AnnulusBoundsPC&& annbo) = default;
  AnnulusBoundsPC& operator=(AnnulusBoundsPC&& sbo) = default;
 
 

  static
  std::pair<AnnulusBoundsPC, double> fromCartesian(const AnnulusBounds& annbo);
 
  virtual AnnulusBoundsPC* clone() const override;
  bool operator==(const SurfaceBounds& sb) const override;
 
  SurfaceBounds::BoundsType type() const override;

  bool inside(const Amg::Vector2D& locpo, double tol1 = 0., double tol2 = 0.) const override final;

  bool inside(const Amg::Vector2D& locpo, const BoundaryCheck& bchk) const override final;

  bool insideLoc1(const Amg::Vector2D& locpo, double tol1 = 0.) const override final;

  bool insideLoc2(const Amg::Vector2D& locpo, double tol2 = 0.) const override final;

  double minDistance(const Amg::Vector2D& locpo) const override final;

  double r() const override;

  double rMin() const { return m_minR; }

  double rMax() const { return m_maxR; }

  double phiMin() const { return m_phiMin + m_phiAvg; }

  double phiMax() const { return m_phiMax + m_phiAvg; }

  Amg::Vector2D moduleOrigin() const;

  MsgStream& dump(MsgStream& sl) const override;
  std::ostream& dump(std::ostream& sl) const override;

  std::array<std::pair<double, double>, 4> corners() const;
 
private:
  double m_minR;
  double m_maxR;
  double m_phiMin;
  double m_phiMax;
  // @TODO: Does this need to be in bound values?
  Amg::Vector2D m_moduleOrigin;
  Amg::Vector2D m_shiftXY; // == -m_moduleOrigin
  Amg::Vector2D m_shiftPC;
  double m_phiAvg;
  Transform2D m_rotationStripPC;
  Transform2D m_translation;
 
  // vectors needed for inside checking
  Amg::Vector2D m_outLeftStripPC;
  Amg::Vector2D m_inLeftStripPC;
  Amg::Vector2D m_outRightStripPC;
  Amg::Vector2D m_inRightStripPC;
 
  Amg::Vector2D m_outLeftModulePC;
  Amg::Vector2D m_inLeftModulePC;
  Amg::Vector2D m_outRightModulePC;
  Amg::Vector2D m_inRightModulePC;
 
  Amg::Vector2D m_outLeftStripXY;
  Amg::Vector2D m_inLeftStripXY;
  Amg::Vector2D m_outRightStripXY;
  Amg::Vector2D m_inRightStripXY;
 
  std::vector<TDD_real_t>   m_boundValues;
 
  Amg::Vector2D stripXYToModulePC(const Amg::Vector2D& vStripXY) const;
  static Amg::Vector2D
  closestOnSegment(const Amg::Vector2D& a,
                   const Amg::Vector2D& b,
                   const Amg::Vector2D& p,
                   const Eigen::Matrix<double, 2, 2>& weight) ;
  static double
  squaredNorm(const Amg::Vector2D& v, const Eigen::Matrix<double, 2, 2>& weight) ;
};
 
} // End of Trk Namespace
 
#include "TrkSurfaces/AnnulusBoundsPC.icc"
 
#endif