ATLAS Offline Software
Loading...
Searching...
No Matches
InDet::TRT_DriftCircleOnTrackRecalibrateTool Class Referencefinal

This tool creates TRT_DriftCircleOnTrack objects using a given track hypothesis. More...

#include <TRT_DriftCircleOnTrackRecalibrateTool.h>

Inheritance diagram for InDet::TRT_DriftCircleOnTrackRecalibrateTool:
Collaboration diagram for InDet::TRT_DriftCircleOnTrackRecalibrateTool:

Public Member Functions

 TRT_DriftCircleOnTrackRecalibrateTool (const std::string &, const std::string &, const IInterface *)
virtual ~TRT_DriftCircleOnTrackRecalibrateTool ()
virtual StatusCode initialize () override
virtual StatusCode finalize () override
virtual Trk::RIO_OnTrackcorrect (const Trk::PrepRawData &, const Trk::TrackParameters &, const EventContext &ctx=Gaudi::Hive::currentContext()) const override
 abstract base method for the creation of RIO_OnTrack it takes a RIO (PrepRawData) and the given Track Parameter
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Static Public Member Functions

static const InterfaceID & interfaceID ()
 The AlgTool InterfaceID.

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed.

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

ToolHandle< Trk::IRIO_OnTrackCreatorm_riontrackTube
ToolHandle< ITRT_DriftFunctionToolm_drifttool
SG::ReadCondHandleKey< LuminosityCondDatam_lumiDataKey {this, "LumiDataKey", "", "SG key for luminosity data"}
SG::ReadCondHandleKey< RIO_OnTrackErrorScalingm_trtErrorScalingKey {this,"TRTErrorScalingKey", "/Indet/TrkErrorScalingTRT", "Key for TRT error scaling conditions data."}
BooleanProperty m_useToTCorrection
DoubleProperty m_scalefactor
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default)
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

This tool creates TRT_DriftCircleOnTrack objects using a given track hypothesis.

See doxygen to Trk::RIO_OnTrackCreator for details.

Definition at line 33 of file TRT_DriftCircleOnTrackRecalibrateTool.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ TRT_DriftCircleOnTrackRecalibrateTool()

InDet::TRT_DriftCircleOnTrackRecalibrateTool::TRT_DriftCircleOnTrackRecalibrateTool ( const std::string & ty,
const std::string & na,
const IInterface * pa )

Definition at line 25 of file TRT_DriftCircleOnTrackRecalibrateTool.cxx.

27 : AthAlgTool(ty,na,pa)
28{
29 declareInterface<IRIO_OnTrackCreator>(this);
30}
AthAlgTool()
Default constructor:

◆ ~TRT_DriftCircleOnTrackRecalibrateTool()

InDet::TRT_DriftCircleOnTrackRecalibrateTool::~TRT_DriftCircleOnTrackRecalibrateTool ( )
virtualdefault

Member Function Documentation

◆ correct()

Trk::RIO_OnTrack * InDet::TRT_DriftCircleOnTrackRecalibrateTool::correct ( const Trk::PrepRawData & hit,
const Trk::TrackParameters & trk,
const EventContext & ctx = Gaudi::Hive::currentContext() ) const
overridevirtual

abstract base method for the creation of RIO_OnTrack it takes a RIO (PrepRawData) and the given Track Parameter

Implements Trk::IRIO_OnTrackCreator.

Definition at line 74 of file TRT_DriftCircleOnTrackRecalibrateTool.cxx.

76{
77 const InDet::TRT_DriftCircle* DC = dynamic_cast<const InDet::TRT_DriftCircle*>(&rio);
78 if(!DC) return nullptr;
79
80 Trk::RIO_OnTrack *rot=nullptr;
81
82 //recalculate driftradius using the new calibration constants
83 bool isOK=false;
84 double t0=0.;
85 double rawtime = m_drifttool->rawTime(DC->driftTimeBin());
86
88 rawtime -= m_drifttool->driftTimeToTCorrection(DC->timeOverThreshold(), DC->identify());
89 }
90
91 double driftradius = m_drifttool->driftRadius(rawtime,DC->identify(),t0,isOK);
92 double error = m_drifttool->errorOfDriftRadius(rawtime-t0,DC->identify());
93
94 double distance=fabs(fabs(TP.parameters()[0])-driftradius);
95
96
97 msg(MSG::VERBOSE)<<"Old radius: "<< DC->localPosition()[Trk::driftRadius]<<" -> new radius: "<<driftradius<<endmsg;
98
99 if(distance>error*m_scalefactor){
100
101 //tube hit
102 //
103 rot=m_riontrackTube->correct(rio,TP,ctx);
104 }
105 else {
106 // precision hit
107 //
108 const InDetDD::TRT_BaseElement* pE = DC->detectorElement(); if(!pE) return nullptr;
109 IdentifierHash iH = pE->identifyHash();
110 const Amg::VectorX& P = TP.parameters();
111 double se = sin(P[3]);
112 Amg::Vector3D dir(cos(P[2])*se,sin(P[2])*se,cos(P[3]));
113
114 double predictedLocZ=P[1];
115
116 double sign = (TP.parameters()[0] < 0.) ? -1. : 1.;
117
118 // only create objects within the active area
119 //
120 double sl = pE->strawLength()*.5;
121 if (predictedLocZ > sl) predictedLocZ = sl;
122 else if(predictedLocZ <-sl) predictedLocZ =-sl;
123
124
126 Amg::MatrixX cov(DC->localCovariance()*escale2);
127
128 if (!m_trtErrorScalingKey.key().empty()) {
129 double mu=0.;
130 if (!m_lumiDataKey.empty()) {
131 SG::ReadCondHandle<LuminosityCondData> lumiData (m_lumiDataKey,ctx);
132 mu = lumiData->lbAverageInteractionsPerCrossing();
133 }
134
135 bool endcap = false;
136 if(dynamic_cast<const InDetDD::TRT_EndcapElement*>(pE)) endcap = true;
137 SG::ReadCondHandle<RIO_OnTrackErrorScaling> error_scaling( m_trtErrorScalingKey,ctx);
139 ->getScaledCovariance(std::move(cov), endcap, mu);
140 }
141
143 Trk::LocalParameters lp(radius);
144
145 rot=new InDet::TRT_DriftCircleOnTrack(DC,std::move(lp),std::move(cov),iH,predictedLocZ,dir,Trk::DECIDED);
146 }
147 return rot;
148}
#define endmsg
static Double_t P(Double_t *tt, Double_t *par)
static Double_t t0
int sign(int a)
MsgStream & msg() const
virtual const double & strawLength() const =0
Active straw length.
virtual IdentifierHash identifyHash() const override final
identifier hash
SG::ReadCondHandleKey< LuminosityCondData > m_lumiDataKey
SG::ReadCondHandleKey< RIO_OnTrackErrorScaling > m_trtErrorScalingKey
int driftTimeBin() const
returns the leading edge bin defined as in TRT_LoLumRawData to be the first 0-1 transition
double timeOverThreshold() const
returns Time over threshold in ns
virtual const InDetDD::TRT_BaseElement * detectorElement() const override final
return the detector element corresponding to this PRD
const Amg::Vector2D & localPosition() const
return the local position reference
Identifier identify() const
return the identifier
const Amg::MatrixX & localCovariance() const
return const ref to the error matrix
Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic > MatrixX
Dynamic Matrix - dynamic allocation.
float distance(const Amg::Vector3D &p1, const Amg::Vector3D &p2)
calculates the distance between two point in 3D space
Eigen::Matrix< double, 3, 1 > Vector3D
Eigen::Matrix< double, Eigen::Dynamic, 1 > VectorX
Dynamic Vector - dynamic allocation.
@ DECIDED
sign of drift radius has been determined
@ driftRadius
trt, straws
Definition ParamDefs.h:53
@ locX
Definition ParamDefs.h:37
std::pair< double, ParamDefs > DefinedParameter
Typedef to of a std::pair<double, ParamDefs> to identify a passed-through double as a specific type o...
const T_res * ErrorScalingCast(const T_src *src)

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl,
const SG::VarHandleKeyType &  )
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ finalize()

StatusCode InDet::TRT_DriftCircleOnTrackRecalibrateTool::finalize ( )
overridevirtual

Definition at line 65 of file TRT_DriftCircleOnTrackRecalibrateTool.cxx.

66{
67 StatusCode sc = AlgTool::finalize(); return sc;
68}
static Double_t sc
::StatusCode StatusCode
StatusCode definition for legacy code.

◆ initialize()

StatusCode InDet::TRT_DriftCircleOnTrackRecalibrateTool::initialize ( )
overridevirtual

Definition at line 42 of file TRT_DriftCircleOnTrackRecalibrateTool.cxx.

43{
44 StatusCode sc = AlgTool::initialize();
45
46 if(m_drifttool.retrieve().isFailure()){
47 msg(MSG::FATAL) << "Failed to retrieve tool " << m_drifttool << endmsg;
48 return StatusCode::FAILURE;
49 }
50
51 if(m_riontrackTube.retrieve().isFailure()){
52 msg(MSG::FATAL) << "Failed to retrieve tool " << m_riontrackTube << endmsg;
53 return StatusCode::FAILURE;
54 }
55
56 ATH_CHECK(m_trtErrorScalingKey.initialize(!m_trtErrorScalingKey.key().empty()));
57 ATH_CHECK(m_lumiDataKey.initialize ( !m_lumiDataKey.key().empty() && !m_trtErrorScalingKey.key().empty()));
58 return sc;
59}
#define ATH_CHECK
Evaluate an expression and check for errors.

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

const InterfaceID & Trk::IRIO_OnTrackCreator::interfaceID ( )
inlinestaticinherited

The AlgTool InterfaceID.

Definition at line 42 of file IRIO_OnTrackCreator.h.

◆ msg()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24 {
25 return this->msgStream();
26 }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce(T &h)

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )
overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }
std::vector< SG::VarHandleKeyArray * > m_vhka

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_drifttool

ToolHandle<ITRT_DriftFunctionTool> InDet::TRT_DriftCircleOnTrackRecalibrateTool::m_drifttool
private
Initial value:
{this, "DriftFunctionTool",
"TRT_DriftFunctionTool"}

Definition at line 62 of file TRT_DriftCircleOnTrackRecalibrateTool.h.

62 {this, "DriftFunctionTool",
63 "TRT_DriftFunctionTool"};

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_lumiDataKey

SG::ReadCondHandleKey<LuminosityCondData> InDet::TRT_DriftCircleOnTrackRecalibrateTool::m_lumiDataKey {this, "LumiDataKey", "", "SG key for luminosity data"}
private

Definition at line 65 of file TRT_DriftCircleOnTrackRecalibrateTool.h.

66{this, "LumiDataKey", "", "SG key for luminosity data"};

◆ m_riontrackTube

ToolHandle<Trk::IRIO_OnTrackCreator> InDet::TRT_DriftCircleOnTrackRecalibrateTool::m_riontrackTube
private
Initial value:
{this, "RIOonTrackToolTube",
"InDet::TRT_DriftCircleOnTrackNoDriftTimeTool/TRT_DriftCircleOnTrackNoDriftTimeTool"}

Definition at line 60 of file TRT_DriftCircleOnTrackRecalibrateTool.h.

60 {this, "RIOonTrackToolTube",
61 "InDet::TRT_DriftCircleOnTrackNoDriftTimeTool/TRT_DriftCircleOnTrackNoDriftTimeTool"};

◆ m_scalefactor

DoubleProperty InDet::TRT_DriftCircleOnTrackRecalibrateTool::m_scalefactor
private
Initial value:
{this, "ScaleHitUncertainty", 2.,
"scale factor for hit uncertainty"}

Definition at line 74 of file TRT_DriftCircleOnTrackRecalibrateTool.h.

74 {this, "ScaleHitUncertainty", 2.,
75 "scale factor for hit uncertainty"};

◆ m_trtErrorScalingKey

SG::ReadCondHandleKey<RIO_OnTrackErrorScaling> InDet::TRT_DriftCircleOnTrackRecalibrateTool::m_trtErrorScalingKey {this,"TRTErrorScalingKey", "/Indet/TrkErrorScalingTRT", "Key for TRT error scaling conditions data."}
private

Definition at line 69 of file TRT_DriftCircleOnTrackRecalibrateTool.h.

70{this,"TRTErrorScalingKey", "/Indet/TrkErrorScalingTRT", "Key for TRT error scaling conditions data."};

◆ m_useToTCorrection

BooleanProperty InDet::TRT_DriftCircleOnTrackRecalibrateTool::m_useToTCorrection
private
Initial value:
{this, "useDriftTimeToTCorrection", false,
"Shall the Time over Threshold correction be used?"}

Definition at line 72 of file TRT_DriftCircleOnTrackRecalibrateTool.h.

72 {this, "useDriftTimeToTCorrection", false,
73 "Shall the Time over Threshold correction be used?"};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.


The documentation for this class was generated from the following files: