Tracking/TrkEvent/TrkSpacePoint/src/SpacePoint.cxx

Go to the documentation of this file.

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

// SpacePoint.cxx Implementation file for class SpacePoint
// (c) ATLAS Detector software
// Version 1.0 14/10/2003 Veronique Boisvert
// Version 2.0   /12/2005 Martin Siebel
 
#include "TrkSpacePoint/SpacePoint.h"
#include "GaudiKernel/MsgStream.h"
#include "TrkDetElementBase/TrkDetElementBase.h"
#include "TrkEventPrimitives/LocalParameters.h"
#include "TrkPrepRawData/PrepRawData.h"
#include "TrkSurfaces/Surface.h"
#include <ostream>
#include <new>
 
namespace Trk
{
  
  
  // ------------------------------------------------------------------
  
  // Default constructor
  SpacePoint::SpacePoint()
    :
    m_clusList(nullptr, nullptr),
    m_elemIdList(0,0),
    m_position(),
    m_globalCovariance()
  {
  } 
  
  // ------------------------------------------------------------------
  
   
  MsgStream& operator << ( MsgStream& sl, const Trk::SpacePoint& spacePoint)
  { 
    return spacePoint.dump(sl);
  }
  
  // ------------------------------------------------------------------
  
  std::ostream& operator << ( std::ostream& sl, const Trk::SpacePoint& spacePoint)
  {
    return spacePoint.dump(sl);
  }

 
  // We compile this function with optimization, even in debug builds; otherwise,
  // the heavy use of Eigen makes it too slow.  However, from here we may call
  // to out-of-line Eigen code that is linked from other DSOs; in that case,
  // it would not be optimized.  Avoid this by forcing all Eigen code
  // to be inlined here if possible.
  ATH_FLATTEN
  void SpacePoint::setupGlobalFromLocalCovariance() {
    const Amg::MatrixX& lc = this->localCovariance();
 
    Eigen::Matrix3d cov;
 
    cov << lc(0, 0), lc(0, 1), 0., 
           lc(1, 0), lc(1, 1), 0.,
           // cppcheck-suppress constStatement
           0., 0., 0.;
 
    const Amg::RotationMatrix3D& R = associatedSurface().transform().rotation();
    m_globalCovariance = R * cov * R.transpose();
  }
 
  const Surface& SpacePoint::associatedSurface() const {
    assert(m_clusList.first->detectorElement());
    return m_clusList.first->detectorElement()->surface();
  }
 
}  // namespace Trk