BoundaryCheck.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/

// BoundaryCheck.h, (c) ATLAS Detector software
 
#ifndef TRKSURFACES_BOUNDARYCHECK_H
#define TRKSURFACES_BOUNDARYCHECK_H
 
#include "EventPrimitives/EventPrimitives.h"
#include "GeoPrimitives/GeoPrimitives.h"
#include "TrkEventPrimitives/ParamDefs.h"
namespace Trk {

 
// maint struct for comparing the extent of arbitrary convex polygons relative to prefixed axes
struct KDOP
{
  float min;
  // Minimum distance (from origin) along axes
  float max;
  // Maximum distance (from origin) along axes
};
 
// struct needed for FastSinCos method (see below)
struct sincosCache
{
  double sinC;
  double cosC;
};
 
class BoundaryCheck
{
  // saves us a lot of function calls in the EllipseToPoly method
  static constexpr double s_cos22 = 0.923879532511286756128183189396788286822416625863642486115097;
  static constexpr double s_cos45 = 0.707106781186547524400844362104849039284835937688474036588339;
  static constexpr double s_cos67 = 0.382683432365089771728459984030398866761344562485627041433800;
 
public:
  enum BoundaryCheckType
  {
    absolute = 0, 
    chi2corr = 1  
  };
 
 
  bool checkLoc1; 
  bool checkLoc2; 
  int nSigmas;          
  double toleranceLoc1; 
  double toleranceLoc2; 
  BoundaryCheckType bcType;
  AmgSymMatrix(2) lCovariance; 
 

  BoundaryCheck(bool sCheck)
    : checkLoc1(sCheck)
    , checkLoc2(sCheck)
    , nSigmas(-1)
    , toleranceLoc1(0.)
    , toleranceLoc2(0.)
    , bcType(absolute)
    , lCovariance(AmgSymMatrix(2)::Identity())
  {}

  BoundaryCheck(bool chkL1, bool chkL2, double tloc1 = 0., double tloc2 = 0.)
    : checkLoc1(chkL1)
    , checkLoc2(chkL2)
    , nSigmas(-1)
    , toleranceLoc1(tloc1)
    , toleranceLoc2(tloc2)
    , bcType(absolute)
    , lCovariance(AmgSymMatrix(2)::Identity())
  {}

  BoundaryCheck(const AmgSymMatrix(2) & lCov, int nsig = 1, bool chkL1 = true, bool chkL2 = true)
    : checkLoc1(chkL1)
    , checkLoc2(chkL2)
    , nSigmas(nsig)
    , toleranceLoc1(0.)
    , toleranceLoc2(0.)
    , bcType(chi2corr)
    , lCovariance(lCov)
  {}

  operator bool() const { return (checkLoc1 || checkLoc2); }

  void ComputeKDOP(const std::vector<Amg::Vector2D> &v, const std::vector<Amg::Vector2D> &KDOPAxes, 
                   std::vector<KDOP>& kdop) const;
 
  std::vector<Amg::Vector2D> EllipseToPoly(int resolution = 3) const;
 
  bool TestKDOPKDOP(const std::vector<KDOP>& a, const std::vector<KDOP>& b) const;
 
  double FastArcTan(double x) const;
 
  sincosCache FastSinCos(double x) const;
};
 
}
 
#include "TrkSurfaces/BoundaryCheck.icc"
#endif // TRKSURFACES_BOUNDARYCHECK_H