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(const int& strawID,
           const InDetDD::TRT_DetElementContainer* detElements);

  void PropagationTime(const int& strawID,
                       const 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( const unsigned int& iPhi,
                                    const unsigned int& iRing,
                                    const unsigned int& iLayer,
                                    const unsigned int& iStraw,
                                    const int strawID,
                    const InDetDD::TRT_DetElementContainer* detElements) ; //Note: Changed from const due to message service hick ups

  double calculateTimeShift_EndCap( const unsigned int& iPhi,
                                    const unsigned int& iWheel,
                                    const unsigned int& iLayer,
                                    const 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,
                                          const int strawID) ;
  void calculateSignalDists_Barrel(const unsigned int& iRing,
                                   const unsigned int& iLayer,
                                   double& direct_dist,
                                   double& reflect_dist ) const;
  void calculateSignalDists_EndCap(const 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