TRTTimeCorrection.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRT_DIGITIZATION_TRTTIMECORRECTION_H
#define TRT_DIGITIZATION_TRTTIMECORRECTION_H
 
#include "AthenaBaseComps/AthMessaging.h"
#include "Identifier/Identifier.h"
#include "TRT_ConditionsServices/ITRT_CalDbTool.h"
 
#include <vector>
 
namespace InDetDD {
  class TRT_DetectorManager;
  class TRT_DetElementContainer;
}
 
class TRTDigSettings;
class TRT_ID;
 
#include "GeoPrimitives/GeoPrimitives.h"

class TRTTimeCorrection : public AthMessaging {
 
public:
 
  TRTTimeCorrection(const TRTDigSettings* digset,
                    const InDetDD::TRT_DetectorManager* detmgr,
                    const TRT_ID*,
                    const ITRT_CalDbTool*);
 
  ~TRTTimeCorrection();

  double TimeShift(int strawID,
                   const InDetDD::TRT_DetElementContainer* detElements);

  void PropagationTime(int strawID,
                       double meanZ,
                       double& propagationTime1,
                       double& propagationTime2);
 
private:
 
  Identifier getIdentifier ( int hitID,
                             bool& statusok) ; //Note: Changed from const due to message service hick ups
 
 
  void Initialize();

  double calculateTimeShift_Barrel( unsigned int iPhi,
                                    unsigned int iRing,
                                    unsigned int iLayer,
                                    unsigned int iStraw,
                                    int strawID,
                                    const InDetDD::TRT_DetElementContainer* detElements) ; //Note: Changed from const due to message service hick ups

  double calculateTimeShift_EndCap( unsigned int iPhi,
                                    unsigned int iWheel,
                                    unsigned int iLayer,
                                    int strawID,
                                    const InDetDD::TRT_DetElementContainer* detElements) ; //Note: Changed from const due to message service hick ups

  double calculateTimeShiftFromStrawEnds( const Amg::Vector3D& strawend1_globalcoord,
                                          const Amg::Vector3D& strawend2_globalcoord,
                                          int strawID) ;
  void calculateSignalDists_Barrel(unsigned int iRing,
                   unsigned int iLayer,
                                   double& direct_dist,
                                   double& reflect_dist ) const;
  void calculateSignalDists_EndCap(unsigned int iWheel,
                                   double& direct_dist,
                                   double& reflect_dist ) const;
 
  const TRTDigSettings* m_settings;
  const InDetDD::TRT_DetectorManager* m_detmgr;
protected:
  const TRT_ID* m_trt_id;
private:
  std::vector< std::vector<  std::vector<double> > >
  m_timeShiftForEndCapPlanes;
  std::vector< std::vector< std::vector< std::vector<double> > > >
  m_timeShiftForBarrelStraws;

  std::vector<double> m_directDistsForEndCapWheels;
  std::vector<double> m_reflectedDistsForEndCapWheels;
  std::vector< std::vector<double> > m_directDistsForBarrelLayers;
  std::vector< std::vector<double> > m_reflectedDistsForBarrelLayers;
 
  //Substitute these with use of hitidhelper?:
  const unsigned int m_subdetectorMask;
  const unsigned int m_right5Bits;
  const unsigned int m_shift5Bits;
  const unsigned int m_shift10Bits;
  const unsigned int m_shift15Bits;

  const double m_notInitVal;
 
  //Digitization parameters
  double m_signalPropagationSpeed = 0.0;
  double m_lengthDeadRegion = 0.0;
  double m_maxVertexDisplacement = 0.0;
  bool m_timeShiftPhiSectSymmetry = false;
  bool m_getT0FromData = false;
 
  const ITRT_CalDbTool* m_trtcaldbtool;
};
 
#endif