DiscOverlapDescriptor.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
// DiscOverlapDescriptor.h, (c) ATLAS Detector software
 
#ifndef INDETTRACKINGGEOMETRY_DISCOVERLAPDESCRIPTOR_H
#define INDETTRACKINGGEOMETRY_DISCOVERLAPDESCRIPTOR_H
 
// Trk
#include "TrkGeometry/OverlapDescriptor.h"
 
// Trk inlcude
#include "TrkDetDescrUtils/BinUtility.h"
#include "TrkDetDescrUtils/BinnedArray1D1D.h"
#include "TrkDetDescrUtils/Intersection.h"
 
// STL include
#include <atomic>
 
#ifndef INDETTRACKINGGEOMETRY_DISCDETADDNEXTPHIETA
#define INDETTRACKINGGEOMETRY_DISCDETADDNEXTPHIETA
 
#define addSurfaceDO(cur, surfaces)                           \
  surfaces.emplace_back(                                      \
      Trk::Intersection(Amg::Vector3D(0., 0., 0.), 0., true), \
      &(cur->surface(cur->identify())))
 
#define addOtherSideDO(cur, surfaces)                         \
  if (cur->otherSide())                                       \
  surfaces.emplace_back(                                      \
      Trk::Intersection(Amg::Vector3D(0., 0., 0.), 0., true), \
      &(cur->otherSide()->surface(cur->otherSide()->identify())))
 
#define addNextInPhiDO(cur, surfaces)       \
  addSurfaceDO(cur->nextInPhi(), surfaces); \
  addOtherSideDO(cur->nextInPhi(), surfaces)
 
#define addPrevInPhiDO(cur, surfaces)       \
  addSurfaceDO(cur->prevInPhi(), surfaces); \
  addOtherSideDO(cur->prevInPhi(), surfaces)
 
#define addNextInEtaDO(cur, surfaces)       \
  addSurfaceDO(cur->nextInEta(), surfaces); \
  addOtherSideDO(cur->nextInEta(), surfaces)
 
#define addPrevInEtaDO(cur, surfaces)       \
  addSurfaceDO(cur->prevInEta(), surfaces); \
  addOtherSideDO(cur->prevInEta(), surfaces)
 
#endif  // INDETTRACKINGGEOMETRY_DISCDETADDNEXTPHIETA
 
class SCT_ID;
class PixelID;
 
namespace Trk {
  class Surface;
}
 
namespace InDetDD {   
     class SiDetectorElement;
}
 
namespace InDet {
    
  class DiscOverlapDescriptor : public Trk::OverlapDescriptor {
  public:
    
    DiscOverlapDescriptor(const Trk::BinnedArray<Trk::Surface>* bin_array = 0,
                          std::vector<Trk::BinUtility*>* singleBinUtils = 0,
                          bool isPixel = false);
    
    virtual ~DiscOverlapDescriptor() {
      if (m_singleBinUtils) {
        std::vector<Trk::BinUtility*>::iterator binIter = m_singleBinUtils->begin();
        for ( ; binIter != m_singleBinUtils->end(); ++binIter) delete *binIter;
      }
      delete m_singleBinUtils;
    }
    
    DiscOverlapDescriptor(const DiscOverlapDescriptor &) = delete;
    
    DiscOverlapDescriptor & operator=(const DiscOverlapDescriptor &) = delete;
    
    virtual DiscOverlapDescriptor* clone() const override;
    
    bool reachableSurfaces(std::vector<Trk::SurfaceIntersection>& surfaces, 
                           const Trk::Surface& sf,
                           const Amg::Vector3D& pos,
                           const Amg::Vector3D& dir) const override;
    
  private:
    bool dumpSurfaces(std::vector<Trk::SurfaceIntersection>& surfaces) const;
    
    const Trk::BinnedArray<Trk::Surface>*           m_bin_array;
    std::vector<Trk::BinUtility*>*                  m_singleBinUtils;
    bool                                            m_pixelCase;
  };
  
  
  inline DiscOverlapDescriptor* DiscOverlapDescriptor::clone() const { return new DiscOverlapDescriptor(); }     
  
}
 
#endif // INDETTRACKINGGEOMETRY_DISCOVERLAPDESCRIPTOR_H