TRT_DriftFunctionTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
 
// TRT_DriftFunctionTool.h
//   Header file for class TRT_DriftFunctionTool
// (c) ATLAS Detector software
// AlgTool used to go from drift time to drift distance
 
#ifndef TRT_DRIFTFUNCTIONTOOL
#define TRT_DRIFTFUNCTIONTOOL
 
#include <string>
#include "AthenaBaseComps/AthAlgTool.h"
#include "TRT_DriftFunctionTool/ITRT_DriftFunctionTool.h"
#include "TRT_ConditionsServices/ITRT_CalDbTool.h"
class TRT_ID;
 
 
#include "TRT_ReadoutGeometry/TRT_DetectorManager.h"
#include "GaudiKernel/ToolHandle.h"
class TRT_DriftFunctionTool: public extends<AthAlgTool, ITRT_DriftFunctionTool>{
 
public:
  TRT_DriftFunctionTool(const std::string& type,
            const std::string& name,
            const IInterface* parent);
 
  virtual ~TRT_DriftFunctionTool();
 
  virtual StatusCode initialize() override;
 
  virtual StatusCode finalize() override;
 
  virtual bool isValidTime(double drifttime) const override;
 
  virtual bool isTestBeamData() const override;
 
  virtual double rawTime(int tdcvalue) const override;
 
 
  virtual double driftRadius(double rawtime, Identifier id, double& t0, bool& isOK, unsigned int word=0) const override;
 
  virtual double driftRadius(double drifttime) const override;
 
  virtual double approxDriftTime(double driftradius) const override;
 
  virtual double errorOfDriftRadius(double drifttime, Identifier id, float mu = -10, unsigned int word=0) const override;
 
  virtual double driftTimeToTCorrection(double tot, Identifier id, bool isArgonStraw=false) const override;
  
  virtual double driftTimeHTCorrection(Identifier id, bool isArgonStraw=false) const override;
  
  void setupRtRelation();
 
 
private:
  
  ToolHandle< ITRT_CalDbTool >   m_TRTCalDbTool;
  ToolHandle< ITRT_CalDbTool >   m_TRTCalDbTool2;
 
  const InDetDD::TRT_DetectorManager* m_manager{};
  const TRT_ID* m_trtid{};
 
 
  double m_drifttimeperbin;            
  double m_error;                      
  
  enum ETimeBins { MaxTimeBin = 50 } ; 
  double m_radius[MaxTimeBin]{};         
  double m_errors[MaxTimeBin]{};         
 
  bool m_ismc;                         
  bool m_isoverlay;                    
  bool m_istestbeam;                   
 
  bool m_dummy;                        
 
  double m_err_fudge;                  
 
  bool  m_allow_digi_version_override; 
  int m_forced_digiversion;            
 
  bool m_override_simcal;              
  bool m_force_universal_errors;       
  double m_uni_error;                  
 
  std::string m_inputfile;             
  std::string m_key;                   
  std::string m_trt_mgr_location;      
  double m_t0_barrel[3]{};               
  double m_t0_endcap[18]{};              
  double m_t0_shift;                   
  double m_ht_correction_barrel_Xe;    
  double m_ht_correction_endcap_Xe;    
  double m_ht_correction_barrel_Ar;    
  double m_ht_correction_endcap_Ar;    
  std::vector<double> m_tot_corrections_barrel_Xe; 
  std::vector<double> m_tot_corrections_endcap_Xe; 
  std::vector<double> m_tot_corrections_barrel_Ar; 
  std::vector<double> m_tot_corrections_endcap_Ar; 
 
  static const size_t s_size_default = 19;
  static constexpr double s_radius_default[s_size_default] = {
    0.   , 0.   , 0.1  , 0.262, 0.466,
    0.607, 0.796, 0.931, 1.065, 1.212,
    1.326, 1.466, 1.585, 1.689, 1.809,
    1.880, 1.940, 1.950, 1.955 };
  static constexpr double s_errors_default[s_size_default] = {
    0.15, 0.15, 0.20, 0.23, 0.21,
    0.18, 0.17, 0.16, 0.15, 0.15,
    0.14, 0.13, 0.12, 0.11, 0.11,
    0.11, 0.13, 0.20, 0.20 };
 
  static const size_t s_size_Comm = 13;
  static constexpr double s_radius_Comm[s_size_Comm] = {
    0.      , 0.     , 0.     , 0.252054, 0.488319,
    0.751514, 1.00173, 1.21851, 1.40886 , 1.68368 ,
    1.85363 , 1.91764, 1.94114 };
  static constexpr double s_errors_Comm[s_size_Comm] = {
    0.10440061, 0.1510298, 0.26130742, 0.260436, 0.246961,
    0.226037,   0.18272  , 0.195482  , 0.213817, 0.157627,
    0.0922559,  0.0463124, 0.0480864 };
 
};
 
inline bool TRT_DriftFunctionTool::isValidTime(double drifttime) const
{ return (drifttime>-10. && drifttime<75.); }
 
inline bool TRT_DriftFunctionTool::isTestBeamData() const
{ return m_istestbeam; }
 
inline double TRT_DriftFunctionTool::rawTime(int tdcvalue) const
{
  double time = (tdcvalue+0.5)*m_drifttimeperbin ;
  return m_istestbeam ? 0.5*time : time;
}
 
 
#endif  // TRT_DRIFTFUNCTIONTOOL_H