AnalyticalDerivCalcTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRK_ANALYTICALDERIVCALCTOOL_H
#define TRK_ANALYTICALDERIVCALCTOOL_H
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "GaudiKernel/ToolHandle.h"
#include "EventPrimitives/EventPrimitives.h"
#include "GeoPrimitives/GeoPrimitives.h"
 
#include "TrkAlignInterfaces/IDerivCalcTool.h"
#include "TrkAlignInterfaces/IAlignModuleTool.h"
 
#include "TrkAlignEvent/AlignResidualType.h"
 
#include <vector>
 
 
class AtlasDetectorID;
 
namespace Trk
{
  class AlignModule;
  class AlignTrack;
  class MeasurementTypeID;
 
  class AnalyticalDerivCalcTool : virtual public IDerivCalcTool, public AthAlgTool
  {
 
  public:
    AnalyticalDerivCalcTool(const std::string & type, const std::string & name, const IInterface * parent);
 
    StatusCode initialize() override;
    StatusCode finalize() override;
 
    bool setDerivatives(AlignTrack * alignTrack) override;
 
    void showStatistics() override {}
 
    bool setResidualCovMatrix(AlignTrack * alignTrack) const override;
 
  private:
 
    PublicToolHandle<IAlignModuleTool> m_alignModuleTool{
      this, "AlignModuleTool", "InDet::InDetAlignModuleTool/InDetAlignModuleTool"};
 
    const AtlasDetectorID * m_idHelper = nullptr;
    MeasurementTypeID * m_measTypeIdHelper = nullptr;
 
    bool getMeasErrorMatrix(const AlignTrack * alignTrack, Amg::MatrixX & V) const;
 
    bool getTrkParamCovMatrix(const AlignTrack * alignTrack, Amg::MatrixX & HCH) const;
 
    bool checkValidity(const Amg::MatrixX & R) const;
 
    std::vector<Amg::VectorX> getDerivatives(AlignTrack * alignTrack,  const AlignModule * module);
 
    void checkResidualType(const AlignTrack * alignTrack);
 
    std::vector<std::pair<const AlignModule *, std::vector<Amg::VectorX> > > m_derivatives;
 
    BooleanProperty m_useLocalSetting{this, "UseLocalSetting", false,
      "use local setup for the covariance matrix of the track"};
 
    // Use constant errors for each sub-detector such that it's equivalent
    // to minimize residual distance instead of minizing residual pull.
    // This is only applied if m_useLocalSetting==true
    BooleanProperty m_useIntrinsicPixelErrors{
      this, "UseIntrinsicPixelError", false, "use intrinsic errors for Pixel"};
    BooleanProperty m_useIntrinsicSCTErrors{
      this, "UseIntrinsicSCTError", false, "use intrinsic errors for SCT"};
    BooleanProperty m_useIntrinsicTRTErrors{
      this, "UseIntrinsicTRTError", false, "use intrinsic errors for TRT"};
 
    int m_residualType = Trk::AlignResidualType::HitOnly; 
    bool m_residualTypeSet = false;   
 
    BooleanProperty m_storeDerivatives{this, "StoreDerivatives", false,
      "store derivatives dr/da on AlignTSOS to be filled into ntuple"};
 
  }; // end class
 
} // end namespace
 
#endif // TRK_ANALYTICALDERIVCALCTOOL_H