ATLAS Offline Software
Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
TBExtrapolTrackToCaloTool Class Reference

This class takes a given TrkTrack, extrapolates it to the Calo and returns the variables needed to compare with the cluster position. More...

#include <TBExtrapolTrackToCaloTool.h>

Inheritance diagram for TBExtrapolTrackToCaloTool:
Collaboration diagram for TBExtrapolTrackToCaloTool:

Public Member Functions

 TBExtrapolTrackToCaloTool (const std::string &type, const std::string &name, const IInterface *parent)
 
virtual ~TBExtrapolTrackToCaloTool ()
 
virtual StatusCode initialize () override
 
virtual StatusCode finalize () override
 
void CaloVariables (const Trk::TrackParameters *parm, double &etaCaloLocal, double &phiCaloLocal)
 now return the variables needed to compare with CaloClusters: More...
 
void CaloLocalPoint (const Trk::TrackParameters *parm, Amg::Vector3D *pt_ctb, Amg::Vector3D *pt_local)
 
bool TrackSeenByCalo (const Trk::Track *trk, const CaloCell_ID::CaloSample sample, const double offset, Amg::Vector3D *pt_ctb, Amg::Vector3D *pt_local)
 The "do-it-all" method which combines the 3 steps. More...
 
bool TrackSeenByCalo (const Trk::TrackParameters *parm, const CaloCell_ID::CaloSample sample, const double offset, Amg::Vector3D *pt_ctb, Amg::Vector3D *pt_local)
 
bool TrackSeenByCalo (const Trk::Track *trk, const bool barrel, const CaloCell_ID::SUBCALO subcalo, const int sampling_or_module, const double offset, Amg::Vector3D *pt_ctb, Amg::Vector3D *pt_local)
 
bool TrackSeenByCalo (const Trk::TrackParameters *parm, const bool barrel, const CaloCell_ID::SUBCALO subcalo, const int sampling_or_module, const double offset, Amg::Vector3D *pt_ctb, Amg::Vector3D *pt_local)
 
bool TrackSeenByCalo (const Trk::Track *trk, const CaloCell_ID::CaloSample sample, const double offset, Amg::Vector3D *pt_ctb, Amg::Vector3D *pt_local, double &trketa_atcalo, double &trkphi_atcalo)
 
const Trk::TrackParameters * extrapolate (const Trk::Track &, const CaloCell_ID::CaloSample, const double, const Trk::PropDirection=Trk::alongMomentum, const Trk::ParticleHypothesis=Trk::undefined) const
 
const Trk::TrackParameters * extrapolate (const Trk::TrackParticleBase &, const CaloCell_ID::CaloSample, const double, const Trk::PropDirection=Trk::alongMomentum, const Trk::ParticleHypothesis=Trk::undefined) const
 
const Trk::TrackParameters * extrapolate (const Trk::TrackParameters &, const CaloCell_ID::CaloSample, const double, Trk::ParticleHypothesis, const Trk::PropDirection=Trk::alongMomentum) const
 
const Trk::NeutralParameters * extrapolate (const Trk::NeutralParameters &, const CaloCell_ID::CaloSample, const double, Trk::ParticleHypothesis, const Trk::PropDirection=Trk::alongMomentum) const
 
virtual const Trk::TrackParameters * extrapolate (const xAOD::TrackParticle &, const CaloCell_ID::CaloSample, const double, const Trk::PropDirection=Trk::alongMomentum, const Trk::ParticleHypothesis=Trk::undefined) const
 
virtual const Trk::NeutralParameters * extrapolate (const xAOD::NeutralParticle &, const CaloCell_ID::CaloSample, const double, const Trk::PropDirection=Trk::alongMomentum) const
 get Trk::NeutralParameters at the calorimeter sample from xAOD::NeutralParticle More...
 
virtual const DataVector< const Trk::TrackParameters > * getParametersInCalo (const Trk::TrackParameters &, Trk::ParticleHypothesis, const Trk::PropDirection=Trk::alongMomentum) const
 
virtual const DataVector< const Trk::NeutralParameters > * getParametersInCalo (const Trk::NeutralParameters &, Trk::ParticleHypothesis, const Trk::PropDirection=Trk::alongMomentum) const
 
virtual const DataVector< const Trk::TrackParameters > * getParametersInCalo (const Trk::TrackParticleBase &, const Trk::ParticleHypothesis=Trk::undefined, const Trk::PropDirection=Trk::alongMomentum) const
 
virtual const DataVector< const Trk::TrackParameters > * getParametersInCalo (const Trk::Track &, const Trk::ParticleHypothesis=Trk::undefined, const Trk::PropDirection=Trk::alongMomentum) const
 
virtual const DataVector< const Trk::TrackParameters > * getParametersInCalo (const xAOD::TrackParticle &, const Trk::ParticleHypothesis=Trk::undefined, const Trk::PropDirection=Trk::alongMomentum) const
 
ImpactInCalo * getImpactInCalo (const Trk::Track &, const Trk::ParticleHypothesis=Trk::undefined, const Trk::PropDirection=Trk::alongMomentum) const
 
ImpactInCalo * getImpactInCalo (const Trk::TrackParticleBase &, const Trk::ParticleHypothesis=Trk::undefined, const Trk::PropDirection=Trk::alongMomentum) const
 
virtual ImpactInCalo * getImpactInCalo (const Trk::TrackParameters &, Trk::ParticleHypothesis, const Trk::PropDirection=Trk::alongMomentum) const
 
virtual ImpactInCalo * getImpactInCalo (const Trk::NeutralParameters &, Trk::ParticleHypothesis, const Trk::PropDirection=Trk::alongMomentum) const
 
CaloDepthTool * getCaloDepth ()
 access to the private tool used to define the extrapolation depth, needed to play with delta-eta More...
 
const Trk::TrackParameters * extrapolate (const Trk::TrackParameters *parm, const Trk::Surface *surf)
 to be used on AOD's or to extrapolate from perigee More...
 
const Trk::TrackParameters * extrapolate (const Trk::Track *trk, const Trk::Surface *surf)
 will use the last measured hit, will not work on AOD's More...
 
virtual Amg::Vector3D getMomentumAtVertex (const xAOD::Vertex &vertex, bool) const
 
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysInitialize () override
 Perform system initialization for an algorithm. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 

Static Public Member Functions

static const InterfaceID & interfaceID ()
 

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
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. More...
 

Private Types

typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
 

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

SG::ReadCondHandleKey< CaloDetDescrManager > m_caloMgrKey
 
const CaloCell_ID * m_calo_id {nullptr}
 
ICaloCoordinateTool * m_calo_tb_coord {nullptr}
 
CaloPhiRange m_range
 
ToolHandle< ICaloSurfaceBuilder > m_calosurf
 
ToolHandle< Trk::IExtrapolator > m_extrapolator
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

This class takes a given TrkTrack, extrapolates it to the Calo and returns the variables needed to compare with the cluster position.

It heavily uses the Tracking PreconfiguredExtrapolator class.

The client is asked to select explicitely the CaloCell_ID::CaloSample where he wants to extrapolate. The 4 steps are combined :

The extrapolarion is performed WITH boundary checks.

INTERFACE :

Calorimeter : the enum CaloCell_ID::CaloSample is strongly recommended, the other one being here for backward compatibility.

Trk::Trk is recommended while running the full reconstruction or on ESD. If so, the extrapolation starts from the last measured point.

Trk::Parameter is mandatory while running on AOD's and can be used if you want to start the extrapolation from the perigee.

BEHIND THE SCEENE : the tools of the package CaloTrackingGeometry/CaloSurfaceBuilder are used to create the proper geometry elements ( Trk::Surface, Trk::Layer, Trk::Volume )

Definition at line 82 of file TBExtrapolTrackToCaloTool.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

◆ TBExtrapolTrackToCaloTool()

TBExtrapolTrackToCaloTool::TBExtrapolTrackToCaloTool ( const std::string &  type,
const std::string &  name,
const IInterface *  parent 
)

Definition at line 39 of file TBExtrapolTrackToCaloTool.cxx.

41  :
42  AthAlgTool(type, name, parent),
43  m_calosurf("CaloSurfaceBuilder"),
44  m_extrapolator("Trk::Extrapolator/AtlasExtrapolator")
45 {
46 
47  declareProperty ("CaloSurfaceBuilder",m_calosurf,"Extrapolation ToolHandle");
48  declareProperty ("Extrapolator",m_extrapolator,"CaloSurfaceBuilder ToolHandle");
49 
50 }

◆ ~TBExtrapolTrackToCaloTool()

TBExtrapolTrackToCaloTool::~TBExtrapolTrackToCaloTool ( )
virtual

Definition at line 52 of file TBExtrapolTrackToCaloTool.cxx.

53 {}

Member Function Documentation

◆ CaloLocalPoint()

void TBExtrapolTrackToCaloTool::CaloLocalPoint ( const Trk::TrackParameters *  parm,
Amg::Vector3D *  pt_ctb,
Amg::Vector3D *  pt_local 
)

Definition at line 165 of file TBExtrapolTrackToCaloTool.cxx.

168 {
169 
170  if (!parm) {
171  ATH_MSG_WARNING ("CaloLocalPoint TrkParam called with null Trk::TrackParameters*" );
172  return;
173  }
174 
175  if (!pt_local) {
176  ATH_MSG_WARNING ("CaloLocalPoint TrkParam called with null Amg::Vector3D* pt_local" );
177  return;
178  }
179 
180  if (!pt_ctb) {
181  ATH_MSG_WARNING ("CaloLocalPoint TrkParam called with null Amg::Vector3D* pt_ctb" );
182  return;
183  }
184 
185 
186  *pt_ctb = parm->position();
187 
188  ATH_MSG_DEBUG ("Impact point in ctb coord : x,y,z= "
189  << pt_ctb->x() << " "
190  << pt_ctb->y() << " " << pt_ctb->z() << " R="
191  << std::sqrt( pt_ctb->x()*pt_ctb->x() + pt_ctb->y()*pt_ctb->y()
192  + pt_ctb->z()*pt_ctb->z())
193  << " eta=" << pt_ctb->eta() << " phi=" << pt_ctb->phi() );
194 
195  m_calo_tb_coord->ctb_to_local(*pt_ctb, *pt_local);
196 
197  ATH_MSG_DEBUG ("Impact point in local coord : x,y,z= "
198  << pt_local->x() << " "
199  << pt_local->y() << " " << pt_local->z() << " R="
200  << std::sqrt( pt_local->x()*pt_local->x()
201  + pt_local->y()*pt_local->y()
202  + pt_local->z()*pt_local->z())
203  << " eta=" << pt_local->eta() << " phi=" << pt_local->phi());
204 }

◆ CaloVariables()

void TBExtrapolTrackToCaloTool::CaloVariables ( const Trk::TrackParameters *  parm,
double &  etaCaloLocal,
double &  phiCaloLocal 
)

now return the variables needed to compare with CaloClusters:

Definition at line 122 of file TBExtrapolTrackToCaloTool.cxx.

124 {
125 
126 
127  Amg::Vector3D pt_ctb,pt_local;
128  etaCaloLocal = 0.;
129  phiCaloLocal = 0.;
130 
131  if (!parm) {
132  ATH_MSG_WARNING ("CaloVariables called with null Trk::TrackParameters*" ) ;
133  return;
134  }
135 
136  pt_ctb = parm->position();
137 
138  m_calo_tb_coord->ctb_to_local(pt_ctb, pt_local);
139 
140  if ( msgLvl(MSG::DEBUG) ){
141 
142  ATH_MSG_DEBUG ("Impact point in ctb coord : x,y,z= "
143  << pt_ctb.x() << " "
144  << pt_ctb.y() << " " << pt_ctb.z() << " R="
145  << std::sqrt( pt_ctb.x()*pt_ctb.x() + pt_ctb.y()*pt_ctb.y()
146  + pt_ctb.z()*pt_ctb.z())
147  << " eta=" << pt_ctb.eta() << " phi=" << pt_ctb.phi() );
148 
149 
150 
151  ATH_MSG_DEBUG ("Impact point in local coord : x,y,z= "
152  << pt_local.x() << " "
153  << pt_local.y() << " " << pt_local.z() << " R="
154  << std::sqrt( pt_local.x()*pt_local.x() + pt_local.y()*pt_local.y()
155  + pt_local.z()*pt_local.z())
156  << " eta=" << pt_local.eta() << " phi=" << pt_local.phi());
157  }
158 
159  etaCaloLocal = pt_local.eta();
160  phiCaloLocal = pt_local.phi();
161 
162 }

◆ declareGaudiProperty() [1/4]

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

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

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

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

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

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

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

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

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  t,
const SG::NotHandleType &   
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleBase &  hndl,
const std::string &  doc,
const SG::VarHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKey &  hndl,
const std::string &  doc,
const SG::VarHandleKeyType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray &  hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType &   
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

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

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
147  return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
148  }

◆ 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.

95 { return m_detStore; }

◆ evtStore() [1/2]

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.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ 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

◆ extrapolate() [1/8]

const Trk::NeutralParameters* TBExtrapolTrackToCaloTool::extrapolate ( const Trk::NeutralParameters &  ,
const CaloCell_ID::CaloSample  ,
const double  ,
Trk::ParticleHypothesis  ,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inline

Definition at line 162 of file TBExtrapolTrackToCaloTool.h.

166  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ extrapolate() [2/8]

const Trk::TrackParameters* TBExtrapolTrackToCaloTool::extrapolate ( const Trk::Track &  ,
const CaloCell_ID::CaloSample  ,
const double  ,
const Trk::PropDirection  = Trk::alongMomentum,
const Trk::ParticleHypothesis  = Trk::undefined 
) const
inline

Definition at line 144 of file TBExtrapolTrackToCaloTool.h.

148  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ extrapolate() [3/8]

const Trk::TrackParameters * TBExtrapolTrackToCaloTool::extrapolate ( const Trk::Track *  trk,
const Trk::Surface *  surf 
)

will use the last measured hit, will not work on AOD's

Definition at line 460 of file TBExtrapolTrackToCaloTool.cxx.

462 {
463 
464  const Trk::TrackParameters* param=0;
465  if (!trk) {
466  ATH_MSG_WARNING ("extrapolate Trk called with null Trk::Track*");
467  return param;
468  }
469 
470  if (!surf) {
471  ATH_MSG_WARNING ("extrapolate Trk called with null Trk::Surface*" );
472  return param;
473  }
474 
475 
476 
477 
478  ATH_MSG_DEBUG ("Trying to propagate to Surface ... " << (*surf) );
479 
480  param = m_extrapolator->extrapolateTrack(
481  Gaudi::Hive::currentContext(),
482  *trk, *surf,
483  Trk::alongMomentum,
484  true, Trk::nonInteracting).release();
485 
486  // the other way to do it:
487  //
488  // const TrackParameters* extrapolate(const Trk::Propagator& prop,
489  // const Trk::Track& trk,
490  // const Trk::Surface& sf,
491  // Trk::PropDirectiondir=anyDirection,Trk::BoundaryCheck bcheck = true)
492  // const;
493 
494  if (param)
495  ATH_MSG_DEBUG ("Propagation successful " );
496  else
497  ATH_MSG_DEBUG ("Propagation failed ");
498 
499  return param;
500 }

◆ extrapolate() [4/8]

const Trk::TrackParameters* TBExtrapolTrackToCaloTool::extrapolate ( const Trk::TrackParameters &  ,
const CaloCell_ID::CaloSample  ,
const double  ,
Trk::ParticleHypothesis  ,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inline

Definition at line 156 of file TBExtrapolTrackToCaloTool.h.

160  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ extrapolate() [5/8]

const Trk::TrackParameters * TBExtrapolTrackToCaloTool::extrapolate ( const Trk::TrackParameters *  parm,
const Trk::Surface *  surf 
)

to be used on AOD's or to extrapolate from perigee

Definition at line 421 of file TBExtrapolTrackToCaloTool.cxx.

423 {
424 
425 
426  const Trk::TrackParameters* param=0;
427 
428  if (!parm) {
429  ATH_MSG_WARNING ("extrapolate TrkParam called with null Trk::TrackParameters*" ) ;
430  return param;
431  }
432 
433  if (!surf) {
434  ATH_MSG_WARNING ("extrapolate TrkParam called with null Trk::Surface*" << endmsg );
435  return param;
436  }
437 
438 
439 
440 
441  ATH_MSG_DEBUG ( "Trying to propagate TrackParameters to Surface ... "
442  << (*surf));
443 
444  // for the moment use nonInteracting to avoid the navigation
445  param = m_extrapolator->extrapolate(
446  Gaudi::Hive::currentContext(),
447  *parm, *surf,
448  Trk::alongMomentum,
449  true, Trk::nonInteracting).release();
450 
451  if (param)
452  ATH_MSG_DEBUG ("Propagation successful ");
453  else
454  ATH_MSG_DEBUG ("Propagation failed ");
455 
456  return param;
457 }

◆ extrapolate() [6/8]

const Trk::TrackParameters* TBExtrapolTrackToCaloTool::extrapolate ( const Trk::TrackParticleBase &  ,
const CaloCell_ID::CaloSample  ,
const double  ,
const Trk::PropDirection  = Trk::alongMomentum,
const Trk::ParticleHypothesis  = Trk::undefined 
) const
inline

Definition at line 150 of file TBExtrapolTrackToCaloTool.h.

154  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ extrapolate() [7/8]

virtual const Trk::NeutralParameters* TBExtrapolTrackToCaloTool::extrapolate ( const xAOD::NeutralParticle &  ,
const CaloCell_ID::CaloSample  ,
const double  ,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inlinevirtual

get Trk::NeutralParameters at the calorimeter sample from xAOD::NeutralParticle

Definition at line 176 of file TBExtrapolTrackToCaloTool.h.

180  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ extrapolate() [8/8]

virtual const Trk::TrackParameters* TBExtrapolTrackToCaloTool::extrapolate ( const xAOD::TrackParticle &  ,
const CaloCell_ID::CaloSample  ,
const double  ,
const Trk::PropDirection  = Trk::alongMomentum,
const Trk::ParticleHypothesis  = Trk::undefined 
) const
inlinevirtual

Definition at line 168 of file TBExtrapolTrackToCaloTool.h.

173  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ finalize()

StatusCode TBExtrapolTrackToCaloTool::finalize ( )
overridevirtual

Definition at line 115 of file TBExtrapolTrackToCaloTool.cxx.

116 {
117  StatusCode sc = StatusCode::SUCCESS;
118  return sc;
119 }

◆ getCaloDepth()

CaloDepthTool * TBExtrapolTrackToCaloTool::getCaloDepth ( )

access to the private tool used to define the extrapolation depth, needed to play with delta-eta

Definition at line 412 of file TBExtrapolTrackToCaloTool.cxx.

413 {
414  if (m_calosurf==0) return nullptr;
415 
416  // getCaloDepth from the surface builder
417  return m_calosurf->getCaloDepth();
418 }

◆ getImpactInCalo() [1/4]

virtual ImpactInCalo* TBExtrapolTrackToCaloTool::getImpactInCalo ( const Trk::NeutralParameters &  ,
Trk::ParticleHypothesis  ,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inlinevirtual

Definition at line 215 of file TBExtrapolTrackToCaloTool.h.

217  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ getImpactInCalo() [2/4]

ImpactInCalo* TBExtrapolTrackToCaloTool::getImpactInCalo ( const Trk::Track &  ,
const Trk::ParticleHypothesis  = Trk::undefined,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inline

Definition at line 203 of file TBExtrapolTrackToCaloTool.h.

205  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ getImpactInCalo() [3/4]

virtual ImpactInCalo* TBExtrapolTrackToCaloTool::getImpactInCalo ( const Trk::TrackParameters &  ,
Trk::ParticleHypothesis  ,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inlinevirtual

Definition at line 211 of file TBExtrapolTrackToCaloTool.h.

213  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ getImpactInCalo() [4/4]

ImpactInCalo* TBExtrapolTrackToCaloTool::getImpactInCalo ( const Trk::TrackParticleBase &  ,
const Trk::ParticleHypothesis  = Trk::undefined,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inline

Definition at line 207 of file TBExtrapolTrackToCaloTool.h.

209  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ getMomentumAtVertex()

Amg::Vector3D TBExtrapolTrackToCaloTool::getMomentumAtVertex ( const xAOD::Vertex &  vertex,
bool   
) const
virtual

Definition at line 503 of file TBExtrapolTrackToCaloTool.cxx.

504 {
505  const EventContext& ctx = Gaudi::Hive::currentContext();
506  Amg::Vector3D momentum(0., 0., 0.);
507  if (vertex.vxTrackAtVertexAvailable())
508  {
509  // Use the parameters at the vertex
510  // (the tracks should be parallel but we will do the sum anyway)
511  for (const auto& trkAtVertex : vertex.vxTrackAtVertex())
512  {
513  const Trk::TrackParameters* paramAtVertex = trkAtVertex.perigeeAtVertex();
514  if (!paramAtVertex)
515  ATH_MSG_WARNING("VxTrackAtVertex does not have perigee at vertex");
516  else
517  momentum += paramAtVertex->momentum();
518  }
519  }
520  else if (vertex.nTrackParticles() == 1)
521  {
522  // Use the first measurement
523  const xAOD::TrackParticle *tp = vertex.trackParticle(0);
524  unsigned int index(0);
525  if (!tp || !tp->indexOfParameterAtPosition(index, xAOD::FirstMeasurement))
526  {
527  ATH_MSG_WARNING("No TrackParticle or no have first measurement");
528  }
529  else
530  momentum += tp->curvilinearParameters(index).momentum();
531  // OR last 3 values of trackParameters(index)
532  }
533  else
534  {
535  // Extrapolate track particles to vertex
536  // (the tracks should be parallel but we will do the sum anyway)
537  const Trk::PerigeeSurface *surface = new Trk::PerigeeSurface(vertex.position());
538  for (unsigned int i = 0; i < vertex.nTrackParticles(); ++i)
539  {
540  const xAOD::TrackParticle* tp = vertex.trackParticle( i );
541  if (!tp)
542  {
543  ATH_MSG_WARNING("NULL pointer to TrackParticle in vertex");
544  continue;
545  }
546  const Trk::TrackParameters* params = m_extrapolator->extrapolate(ctx,
547  tp->perigeeParameters(),
548  *surface,
549  Trk::alongMomentum).release();
550  if (!params)
551  ATH_MSG_DEBUG("Extrapolation to vertex (perigee) failed");
552  else
553  {
554  momentum += params->momentum();
555  delete params;
556  }
557  }
558  delete surface;
559  }
560 
561  return momentum;
562 }

◆ getParametersInCalo() [1/5]

virtual const DataVector< const Trk::NeutralParameters >* TBExtrapolTrackToCaloTool::getParametersInCalo ( const Trk::NeutralParameters &  ,
Trk::ParticleHypothesis  ,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inlinevirtual

Definition at line 186 of file TBExtrapolTrackToCaloTool.h.

188  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ getParametersInCalo() [2/5]

virtual const DataVector< const Trk::TrackParameters >* TBExtrapolTrackToCaloTool::getParametersInCalo ( const Trk::Track &  ,
const Trk::ParticleHypothesis  = Trk::undefined,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inlinevirtual

Definition at line 194 of file TBExtrapolTrackToCaloTool.h.

196  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ getParametersInCalo() [3/5]

virtual const DataVector< const Trk::TrackParameters >* TBExtrapolTrackToCaloTool::getParametersInCalo ( const Trk::TrackParameters &  ,
Trk::ParticleHypothesis  ,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inlinevirtual

Definition at line 182 of file TBExtrapolTrackToCaloTool.h.

184  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ getParametersInCalo() [4/5]

virtual const DataVector< const Trk::TrackParameters >* TBExtrapolTrackToCaloTool::getParametersInCalo ( const Trk::TrackParticleBase &  ,
const Trk::ParticleHypothesis  = Trk::undefined,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inlinevirtual

Definition at line 190 of file TBExtrapolTrackToCaloTool.h.

192  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ getParametersInCalo() [5/5]

virtual const DataVector< const Trk::TrackParameters >* TBExtrapolTrackToCaloTool::getParametersInCalo ( const xAOD::TrackParticle &  ,
const Trk::ParticleHypothesis  = Trk::undefined,
const Trk::PropDirection  = Trk::alongMomentum 
) const
inlinevirtual

Definition at line 198 of file TBExtrapolTrackToCaloTool.h.

200  {std::cout<<"Not Implemented !!!!"<<std::endl; return 0;}

◆ initialize()

StatusCode TBExtrapolTrackToCaloTool::initialize ( )
overridevirtual

Definition at line 56 of file TBExtrapolTrackToCaloTool.cxx.

57 {
58 
59  // at this point the joboption should have been read
60  ATH_CHECK(m_caloMgrKey.initialize());
61  ATH_CHECK(detStore()->retrieve(m_calo_id,"CaloCell_ID"));
62 
63 
64 
65  ATH_MSG_DEBUG ( " getting started ");
66 
67 
68 
69  // Get TrkExtrapolator from ToolService
70  StatusCode s=m_extrapolator.retrieve();
71 
72  if (s.isFailure())
73  {
74  ATH_MSG_FATAL ( "Could not find Tool "
75  << m_extrapolator.name() << " with parent " << m_extrapolator.parentName() );
76  return s;
77  } else
78  {
79  ATH_MSG_INFO ( "Successfully retrieved tool " << m_extrapolator.typeAndName() << " with parent " << m_extrapolator.parentName() ) ;
80  }
81 
82 
83  IAlgTool* tool;
84  s = toolSvc()->retrieveTool("TBCaloCoordinate",tool);
85  if (s.isFailure())
86  {
87  ATH_MSG_FATAL ("Could not find Tool TBCaloCoordinate" );
88 
89  return s;
90  }
91  m_calo_tb_coord = dynamic_cast<ICaloCoordinateTool*>(tool);
92 
93 
94 
95  // Get the CaloSurface tool and sets the depthtool it should use:
96  s=m_calosurf.retrieve();
97 
98  if (s.isFailure())
99  {
100  ATH_MSG_FATAL ("Could not find Tool CaloSurfaceBuilder"
101  << m_calosurf.typeAndName());
102  return s;
103  } else {
104  ATH_MSG_INFO ("Successfully retrieved Tool CaloSurfaceBuilder "
105  << m_calosurf.typeAndName());
106 
107  }
108 
109 
110  StatusCode sc = StatusCode::SUCCESS;
111  return sc;
112 }

◆ 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()

static const InterfaceID& TBExtrapolTrackToCaloTool::interfaceID ( )
inlinestatic

Definition at line 87 of file TBExtrapolTrackToCaloTool.h.

87 { return IID_TBExtrapolTrackToCaloTool; }

◆ msg() [1/2]

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ 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();
383  PBASE::renounce (h);
384  }

◆ 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  {
365  handlesArray.renounce();
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 DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ 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.

◆ TrackSeenByCalo() [1/5]

bool TBExtrapolTrackToCaloTool::TrackSeenByCalo ( const Trk::Track *  trk,
const bool  barrel,
const CaloCell_ID::SUBCALO  subcalo,
const int  sampling_or_module,
const double  offset,
Amg::Vector3D *  pt_ctb,
Amg::Vector3D *  pt_local 
)

Definition at line 312 of file TBExtrapolTrackToCaloTool.cxx.

319 {
320 
321 
322  CaloCell_ID::CaloSample sample;
323  CaloDetDescrManager::build_sample ( subcalo, barrel ,sampling_or_module, sample);
324  return TrackSeenByCalo(trk,sample,offset,pt_ctb,pt_local);
325 }

◆ TrackSeenByCalo() [2/5]

bool TBExtrapolTrackToCaloTool::TrackSeenByCalo ( const Trk::Track *  trk,
const CaloCell_ID::CaloSample  sample,
const double  offset,
Amg::Vector3D *  pt_ctb,
Amg::Vector3D *  pt_local 
)

The "do-it-all" method which combines the 3 steps.

Definition at line 207 of file TBExtrapolTrackToCaloTool.cxx.

212 {
213 
214  if (!trk) {
215  ATH_MSG_WARNING ("TrackSeenByCalo Trk called with null Trk::TrackParameters*" );
216  return false;
217  }
218 
219  if (!pt_local) {
220  ATH_MSG_WARNING ("TrackSeenByCalo Trk called with null Amg::Vector3D* pt_local");
221 
222  return false;
223  }
224 
225  if (!pt_ctb) {
226  ATH_MSG_WARNING ("TrackSeenByCalo Trk called with null Amg::Vector3D* pt_ctb");
227 
228  return false;
229  }
230 
231 
232  Trk::Surface* surf = 0;
233  bool success = false;
234  double trketa = 0.;
235  SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{ m_caloMgrKey };
236  const CaloDetDescrManager* caloDDMgr = *caloMgrHandle;
237  // Take eta of the last measured hit as best guess and create surface :
238  const DataVector<const Trk::TrackParameters>* paramvec =
239  trk->trackParameters();
240  if (paramvec) {
241  for (const Trk::TrackParameters* params : *paramvec)
242  trketa = params->eta();
243  surf = m_calosurf->CreateUserSurface(sample, offset, trketa, caloDDMgr);
244  }
245  if (!surf)
246  return success;
247 
248  //std::cout << " here we go : " << std::endl;
249  const Trk::TrackParameters* param = extrapolate(trk,surf);
250 
251 
252 
253  if(param) {
254  CaloLocalPoint (param,pt_ctb, pt_local);
255  success = true;
256  delete param;
257  }
258  delete surf;
259 
260  return success;
261 }

◆ TrackSeenByCalo() [3/5]

bool TBExtrapolTrackToCaloTool::TrackSeenByCalo ( const Trk::Track *  trk,
const CaloCell_ID::CaloSample  sample,
const double  offset,
Amg::Vector3D *  pt_ctb,
Amg::Vector3D *  pt_local,
double &  trketa_atcalo,
double &  trkphi_atcalo 
)

Definition at line 342 of file TBExtrapolTrackToCaloTool.cxx.

348 {
349 
350 
351  if (!trk) {
352  ATH_MSG_WARNING ("TrackSeenByCalo Trk called with null Trk::TrackParameters*" );
353  return false;
354  }
355 
356  if (!pt_local) {
357  ATH_MSG_WARNING ("TrackSeenByCalo Trk called with null Amg::Vector3D* pt_local" ) ;
358  return false;
359  }
360 
361  if (!pt_ctb) {
362  ATH_MSG_WARNING ("TrackSeenByCalo Trk called with null Amg::Vector3D* pt_ctb" );
363  return false;
364  }
365 
366  Trk::Surface* surf = 0;
367  bool success = false;
368  double trketa = 0.;
369 
370  // Take eta of the last measured hit as best guess and create surface :
371  const DataVector<const Trk::TrackParameters>* paramvec =
372  trk->trackParameters();
373  if (paramvec) {
374  SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{ m_caloMgrKey };
375  if ((*caloMgrHandle)->lar_geometry() == "H8") {
376  trketa = m_calo_tb_coord->beam_local_eta();
377  } else {
378  for (const Trk::TrackParameters* params : *paramvec)
379  trketa = params->eta();
380  }
381 
382  surf = m_calosurf->CreateUserSurface(sample, offset, trketa, (*caloMgrHandle));
383  }
384  if (!surf)
385  return success;
386 
387  const Trk::TrackParameters* param = extrapolate(trk,surf);
388 
389 
390 
391  if(param) {
392  CaloLocalPoint (param,pt_ctb, pt_local);
393  trketa_atcalo = param->eta();
394  if ( sin(param->parameters()[Trk::phi]) > 0.)
395  trkphi_atcalo = std::acos(cos(param->parameters()[Trk::phi]));
396  else
397  trkphi_atcalo = -1. * std::acos(cos(param->parameters()[Trk::phi]));
398  trkphi_atcalo = m_range.fix(trkphi_atcalo);
399  success = true;
400  delete param;
401  }
402  else {
403  trketa_atcalo = 999999.;
404  trkphi_atcalo = 999999.;
405  }
406  delete surf;
407 
408  return success;
409 }

◆ TrackSeenByCalo() [4/5]

bool TBExtrapolTrackToCaloTool::TrackSeenByCalo ( const Trk::TrackParameters *  parm,
const bool  barrel,
const CaloCell_ID::SUBCALO  subcalo,
const int  sampling_or_module,
const double  offset,
Amg::Vector3D *  pt_ctb,
Amg::Vector3D *  pt_local 
)

Definition at line 328 of file TBExtrapolTrackToCaloTool.cxx.

335 {
336  CaloCell_ID::CaloSample sample;
337  CaloDetDescrManager::build_sample ( subcalo, barrel ,sampling_or_module, sample);
338  return TrackSeenByCalo(parm,sample,offset,pt_ctb,pt_local);
339 }

◆ TrackSeenByCalo() [5/5]

bool TBExtrapolTrackToCaloTool::TrackSeenByCalo ( const Trk::TrackParameters *  parm,
const CaloCell_ID::CaloSample  sample,
const double  offset,
Amg::Vector3D *  pt_ctb,
Amg::Vector3D *  pt_local 
)

Definition at line 264 of file TBExtrapolTrackToCaloTool.cxx.

269 {
270  bool success = false;
271 
272 
273  if (!parm) {
274  ATH_MSG_WARNING ("TrackSeenByCalo TrkParam called with null Trk::TrackParameters*" );
275  return false;
276  }
277 
278  if (!pt_local) {
279 
280  ATH_MSG_WARNING("TrackSeenByCalo TrkParam called with null Amg::Vector3D* pt_local" );
281  return false;
282  }
283 
284  if (!pt_ctb) {
285  ATH_MSG_WARNING ("TrackSeenByCalo TrkParam called with null Amg::Vector3D* pt_ctb" );
286  return false;
287  }
288 
289 
290  // Take eta as best guess of track direction and create surface :
291  double trketa = parm->eta();
292  SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{ m_caloMgrKey };
293  const CaloDetDescrManager* caloDDMgr = *caloMgrHandle;
294  Trk::Surface* surf = m_calosurf->CreateUserSurface (sample,offset,trketa,caloDDMgr);
295  if (!surf) return success;
296 
297  const Trk::TrackParameters* resparam = extrapolate(parm,surf);
298 
299 
300 
301  if(resparam) {
302  CaloLocalPoint (resparam,pt_ctb, pt_local);
303  success = true;
304  delete resparam;
305  }
306  delete surf;
307 
308  return success;
309 }

◆ 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) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ m_calo_id

const CaloCell_ID* TBExtrapolTrackToCaloTool::m_calo_id {nullptr}
private

Definition at line 248 of file TBExtrapolTrackToCaloTool.h.

◆ m_calo_tb_coord

ICaloCoordinateTool* TBExtrapolTrackToCaloTool::m_calo_tb_coord {nullptr}
private

Definition at line 249 of file TBExtrapolTrackToCaloTool.h.

◆ m_caloMgrKey

SG::ReadCondHandleKey<CaloDetDescrManager> TBExtrapolTrackToCaloTool::m_caloMgrKey
private
Initial value:
{ this
, "CaloDetDescrManager"
, "CaloDetDescrManager"
, "SG Key for CaloDetDescrManager in the Condition Store" }

Definition at line 243 of file TBExtrapolTrackToCaloTool.h.

◆ m_calosurf

ToolHandle<ICaloSurfaceBuilder> TBExtrapolTrackToCaloTool::m_calosurf
private

Definition at line 254 of file TBExtrapolTrackToCaloTool.h.

◆ 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_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_extrapolator

ToolHandle< Trk::IExtrapolator > TBExtrapolTrackToCaloTool::m_extrapolator
private

Definition at line 258 of file TBExtrapolTrackToCaloTool.h.

◆ m_range

CaloPhiRange TBExtrapolTrackToCaloTool::m_range
private

Definition at line 251 of file TBExtrapolTrackToCaloTool.h.

◆ 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:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
ICaloCoordinateTool::beam_local_eta
virtual double beam_local_eta()=0
the most common use-case : the H8 beam goes along the x axis -> that are the local calorimeters eta a...
TBExtrapolTrackToCaloTool::m_calo_tb_coord
ICaloCoordinateTool * m_calo_tb_coord
Definition: TBExtrapolTrackToCaloTool.h:249
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
python.SystemOfUnits.s
int s
Definition: SystemOfUnits.py:131
python.PerfMonSerializer.p
def p
Definition: PerfMonSerializer.py:743
TBExtrapolTrackToCaloTool::m_caloMgrKey
SG::ReadCondHandleKey< CaloDetDescrManager > m_caloMgrKey
Definition: TBExtrapolTrackToCaloTool.h:243
TBExtrapolTrackToCaloTool::m_calo_id
const CaloCell_ID * m_calo_id
Definition: TBExtrapolTrackToCaloTool.h:248
SG::ReadCondHandle
Definition: ReadCondHandle.h:44
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
Trk::PerigeeSurface
Definition: PerigeeSurface.h:43
Trk::ParametersBase::position
const Amg::Vector3D & position() const
Access method for the position.
index
Definition: index.py:1
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
ParticleTest.tp
tp
Definition: ParticleTest.py:25
AthCommonMsg< AlgTool >::msgLvl
bool msgLvl(const MSG::Level lvl) const
Definition: AthCommonMsg.h:30
Trk::alongMomentum
@ alongMomentum
Definition: PropDirection.h:20
TBExtrapolTrackToCaloTool::extrapolate
const Trk::TrackParameters * extrapolate(const Trk::Track &, const CaloCell_ID::CaloSample, const double, const Trk::PropDirection=Trk::alongMomentum, const Trk::ParticleHypothesis=Trk::undefined) const
Definition: TBExtrapolTrackToCaloTool.h:144
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
TBExtrapolTrackToCaloTool::TrackSeenByCalo
bool TrackSeenByCalo(const Trk::Track *trk, const CaloCell_ID::CaloSample sample, const double offset, Amg::Vector3D *pt_ctb, Amg::Vector3D *pt_local)
The "do-it-all" method which combines the 3 steps.
Definition: TBExtrapolTrackToCaloTool.cxx:207
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:95
CaloDetDescrManager_Base::build_sample
static void build_sample(CaloCell_ID::SUBCALO subCalo, bool barrel, int sampling_or_module, CaloCell_ID::CaloSample &sample)
translate between the 2 ways to label a sub-detector:
Definition: CaloDetDescrManager.cxx:1425
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
TBExtrapolTrackToCaloTool::m_calosurf
ToolHandle< ICaloSurfaceBuilder > m_calosurf
Definition: TBExtrapolTrackToCaloTool.h:254
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
ParticleGun_EoverP_Config.momentum
momentum
Definition: ParticleGun_EoverP_Config.py:63
FullCPAlgorithmsTest_eljob.sample
sample
Definition: FullCPAlgorithmsTest_eljob.py:100
lumiFormat.i
int i
Definition: lumiFormat.py:92
CaloSampling::CaloSample
CaloSample
Definition: Calorimeter/CaloGeoHelpers/CaloGeoHelpers/CaloSampling.h:22
xAOD::FirstMeasurement
@ FirstMeasurement
Parameter defined at the position of the 1st measurement.
Definition: TrackingPrimitives.h:213
endmsg
#define endmsg
Definition: AnalysisConfig_Ntuple.cxx:63
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
test_pyathena.parent
parent
Definition: test_pyathena.py:15
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
Trk::ParametersBase
Definition: ParametersBase.h:55
DataVector
Derived DataVector<T>.
Definition: DataVector.h:581
CaloPhiRange::fix
static double fix(double phi)
Definition: CaloPhiRange.cxx:14
ICaloCoordinateTool
This (clean) interface is implemented in the (rather dirty) ICaloCoordinateTool class,...
Definition: ICaloCoordinateTool.h:65
AthAlgTool::AthAlgTool
AthAlgTool()
Default constructor:
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
Trk::Track::trackParameters
const DataVector< const TrackParameters > * trackParameters() const
Return a pointer to a vector of TrackParameters.
Definition: Tracking/TrkEvent/TrkTrack/src/Track.cxx:97
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:192
Trk::nonInteracting
@ nonInteracting
Definition: ParticleHypothesis.h:25
AtlCoolConsole.tool
tool
Definition: AtlCoolConsole.py:453
SG::CondHandleKey::initialize
StatusCode initialize(bool used=true)
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
TBExtrapolTrackToCaloTool::m_range
CaloPhiRange m_range
Definition: TBExtrapolTrackToCaloTool.h:251
Trk::ParametersBase::momentum
const Amg::Vector3D & momentum() const
Access method for the momentum.
Trk::vertex
@ vertex
Definition: MeasurementType.h:21
DeMoScan.index
string index
Definition: DeMoScan.py:362
a
TList * a
Definition: liststreamerinfos.cxx:10
h
CaloDetDescrManager
This class provides the client interface for accessing the detector description information common to...
Definition: CaloDetDescrManager.h:473
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
python.CaloScaleNoiseConfig.type
type
Definition: CaloScaleNoiseConfig.py:78
DEBUG
#define DEBUG
Definition: page_access.h:11
DetectorZone::barrel
@ barrel
convertTimingResiduals.offset
offset
Definition: convertTimingResiduals.py:71
TBExtrapolTrackToCaloTool::CaloLocalPoint
void CaloLocalPoint(const Trk::TrackParameters *parm, Amg::Vector3D *pt_ctb, Amg::Vector3D *pt_local)
Definition: TBExtrapolTrackToCaloTool.cxx:165
ICaloCoordinateTool::ctb_to_local
virtual void ctb_to_local(Amg::Vector3D &pt_ctb, Amg::Vector3D &pt_local)=0
General use method: basic translation between the 2 cartesian coordinate systems:
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:616
TBExtrapolTrackToCaloTool::m_extrapolator
ToolHandle< Trk::IExtrapolator > m_extrapolator
Definition: TBExtrapolTrackToCaloTool.h:258
declareProperty
#define declareProperty(n, p, h)
Definition: BaseFakeBkgTool.cxx:15
Trk::ParametersBase::eta
double eta() const
Access method for pseudorapidity - from momentum.
Trk::phi
@ phi
Definition: ParamDefs.h:81
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:790
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
PowhegControl_ttFCNC_NLO.params
params
Definition: PowhegControl_ttFCNC_NLO.py:226
Trk::Surface
Definition: Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h:75
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
fitman.k
k
Definition: fitman.py:528
Generated on Fri Jun 28 2024 22:00:48 for ATLAS Offline Software by   doxygen 1.8.18