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"
 
#include "CLHEP/Units/SystemOfUnits.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 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 {this, "TRTCalDbTool",  "TRT_CalDbTool"};
  ToolHandle< ITRT_CalDbTool > m_TRTCalDbTool2{this, "TRTCalDbTool2", ""};
 
  const InDetDD::TRT_DetectorManager* m_manager{};
  const TRT_ID* m_trtid{};
 
  static constexpr double s_drifttimeperbin = 3.125 * CLHEP::ns;
  
  enum ETimeBins { MaxTimeBin = 50 } ;   
  double m_radius[MaxTimeBin]{};         
  double m_errors[MaxTimeBin]{};         
 
  BooleanProperty m_ismc{this, "IsMC", true};
  BooleanProperty m_isoverlay{this, "IsOverlay", false};
  BooleanProperty m_dummy{this, "DummyMode", false,
    "flag for ignoring drift time info"};
 
  bool m_istestbeam = false;
 
  DoubleProperty m_err_fudge{this, "ErrorFudgeFactor", 1.0,
    "fudge_factor for error scaling"};
 
  BooleanProperty m_allow_digi_version_override
    {this, "AllowDigiVersionOverride", false}; 
  IntegerProperty m_forced_digiversion
    {this, "ForcedDigiVersion", 11};           
 
  BooleanProperty m_force_universal_errors
    {this, "ForceUniversalErrors", false}; 
  DoubleProperty m_uni_error{this, "UniversalError", 0.136}; 
 
  std::string m_key = ""; 
  StringProperty m_trt_mgr_location
    {this, "TrtDescrManageLocation", "TRT", "Name of TRT detector manager"};
 
  double m_t0_barrel[3]{};               
  double m_t0_endcap[18]{};              
  double m_t0_shift;                   
 
  DoubleProperty m_ht_correction_barrel_Xe{this, "HTCorrectionBarrelXe", 0.,
    "HT correction for Xe straws in barrel"};
  DoubleProperty m_ht_correction_endcap_Xe{this, "HTCorrectionEndcapXe", 0.,
    "HT correction for Xe straws in endcaps"};
  DoubleProperty m_ht_correction_barrel_Ar{this, "HTCorrectionBarrelAr", 0.,
    "HT correction for Ar straws in barrel"};
  DoubleProperty m_ht_correction_endcap_Ar{this, "HTCorrectionEndcapAr", 0.,
    "HT correction for Ar straws in endcaps"};
 
  DoubleArrayProperty m_tot_corrections_barrel_Xe
    {this, "ToTCorrectionsBarrelXe", {},
     "ToT corrections for 20 ToT bins in Xe barrel straws"};
  DoubleArrayProperty m_tot_corrections_endcap_Xe
    {this, "ToTCorrectionsEndcapXe", {},
     "ToT corrections for 20 ToT bins in Xe endcap straws"};
  DoubleArrayProperty m_tot_corrections_barrel_Ar
    {this, "ToTCorrectionsBarrelAr", {},
     "ToT corrections for 20 ToT bins in Ar barrel straws"};
  DoubleArrayProperty m_tot_corrections_endcap_Ar
    {this, "ToTCorrectionsEndcapAr", {},
     "ToT corrections for 20 ToT bins in Ar endcap straws"};
 
  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 double TRT_DriftFunctionTool::rawTime(int tdcvalue) const
{
  double time = (tdcvalue+0.5)*s_drifttimeperbin ;
  return m_istestbeam ? 0.5*time : time;
}
 
 
#endif  // TRT_DRIFTFUNCTIONTOOL_H