![Logo](../../ATLAS-Logo-Square-Blue-RGB.png) |
ATLAS Offline Software
|
#include <TrackParticleClusterAssociationAlg.h>
|
| TrackParticleClusterAssociationAlg (const std::string &name, ISvcLocator *pSvcLocator) |
|
StatusCode | initialize () override |
|
StatusCode | execute () override |
|
virtual StatusCode | sysInitialize () override |
| Override sysInitialize. More...
|
|
virtual const DataObjIDColl & | extraOutputDeps () const override |
| Return the list of extra output dependencies. More...
|
|
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 | 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 |
|
|
Gaudi::Property< float > | m_dr {this, "DeltaR", 0.1, "max deltaR to match track and cluster" } |
|
Gaudi::Property< bool > | m_useCovariance {this, "UseCovariance", true, "use covariance from TrkParameters to evaluate angular uncertainties"} |
|
Gaudi::Property< double > | m_ptCut {this, "PtCut", 25000., ""} |
|
SG::ReadHandleKey< CaloExtensionCollection > | m_caloExtKey {this, "CaloExtensionName", "", ""} |
|
SG::ReadHandleKey< xAOD::TrackParticleContainer > | m_trackParticleCollectionHandle {this,"TrackParticleContainerName", "InDetTrackParticles", "input tracks" } |
|
SG::ReadHandleKey< xAOD::CaloClusterContainer > | m_caloClusters {this, "CaloClusterLocation", "CaloCalTopoClusters","input calo clusters"} |
|
SG::ReadDecorHandleKey< xAOD::CaloClusterContainer > | m_detectorEtaDecor { this, "DetectorEtaName", "", "Decoration for CaloCluster DetectorEta" } |
|
bool | m_doDetEta |
|
SG::ReadHandleKey< xAOD::VertexContainer > | m_vertexContHandle {this, "VertexContainerName", "", "if empty all tracks will be decorated. if not only those corresponding to the PV[0] will."} |
|
ToolHandle< CP::ITrackVertexAssociationTool > | m_trackvertexassoTool {this, "TrackVertexAssoTool", "" } |
|
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > | m_assocClustersDecor {this, "AssociatedClusterDecorKey" , "", "Decoration key to store associated clusters. IMPORTANT must be consistent with TrackParticleContainer" } |
|
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > | m_caloEntryParsDecor {this, "CaloEntryParsDecorKey" , "", "Decoration name to store trk parameters to calo entry (if non blank). IMPORTANT : must be consistent with TrackParticleContainer" } |
|
DataObjIDColl | m_extendedExtraObjects |
|
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 |
|
Associate clusters to each track based on its extrapolation to the calo.
This alg reads in a CaloExtensionCollection corresponding to a TrackParticleContainer. It then decorates each TrackParticleContainer with a vector<EL<CaloCluster>> containing the clusters close enough (DeltaR property) w.r.t CaloExtension corresponding to the track. Optionnaly tracks are also decorated with a const Trk::TrackParameters* obtain from caloExtension->caloEntryLayerIntersection().
Note: In earlier versions this alg was creating 2 maps a xAODAssociations/TrackParticleClusterAssociationContainer and a Rec::ParticleClusterAssociation and was using ParticleCaloClusterAssociationTool to retrieve the extrapolation and perform the clusters association. All of this was suppress or replaced in order to simplify the procedure.
Definition at line 38 of file TrackParticleClusterAssociationAlg.h.
◆ StoreGateSvc_t
◆ TrackParticleClusterAssociationAlg()
TrackParticleClusterAssociationAlg::TrackParticleClusterAssociationAlg |
( |
const std::string & |
name, |
|
|
ISvcLocator * |
pSvcLocator |
|
) |
| |
◆ associatedClusters()
returns the clusters from allClusters which are close enough to caloExtensions
Definition at line 144 of file TrackParticleClusterAssociationAlg.cxx.
146 std::vector<const xAOD::CaloCluster* >
clusters;
150 ATH_MSG_WARNING(
" NO TrackParameters caloExtension.caloEntryLayerIntersection() ");
154 float eta =
pars->position().eta();
155 float phi =
pars->position().phi();
157 double uncertEta = 0.;
158 double uncertPhi = 0.;
159 if(
pars->covariance()) {
163 double uncertExtrp = uncertEta*uncertEta + uncertPhi*uncertPhi;
177 cl_eta = detEta(*
cl);
184 float dr2Cut = dr2Cut0;
188 double sigmaWidth = sig_acc(*
cl);
189 double uncertClus = 2.*sigmaWidth*sigmaWidth;
190 if(uncertExtrp>uncertClus){
191 ATH_MSG_DEBUG(
"Extrapolation uncertainty larger than cluster width! Returning without association.");
195 dr2Cut = (sigmaWidth+uncertEta)*(sigmaWidth+uncertEta)+(sigmaWidth+uncertPhi)*(sigmaWidth+uncertPhi);
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation 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.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
◆ declareProperty() [2/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation 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.
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation 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.
◆ declareProperty() [5/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation 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.
◆ declareProperty() [6/6]
◆ detStore()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode TrackParticleClusterAssociationAlg::execute |
( |
| ) |
|
|
override |
Definition at line 46 of file TrackParticleClusterAssociationAlg.cxx.
74 double sigmaWidth = 0.0;
75 if(cent > 0) sigmaWidth =
atan(sqrt(
rad)/cent)*cosh(cl_eta);
76 sig_dec(*
cl) = sigmaWidth;
87 if(!vxCont->empty()) pv0=(*vxCont)[0];
97 unsigned int ntracks = 0;
101 std::vector< ElementLink< xAOD::CaloClusterContainer > > & caloClusterLinks = assoClustDecor(*
tp);
107 ATH_MSG_DEBUG(
" Selected track " <<
tp->index() <<
" pt " <<
tp->pt() <<
" eta " <<
tp->eta() <<
" phi " <<
tp->phi() );
111 if (caloExtension ==
nullptr ) {
112 ATH_MSG_DEBUG(
" Selected track "<<
tp->index() <<
" has no caloExtension ");
117 std::vector<const xAOD::CaloCluster*> assoClusters =
associatedClusters( *caloExtension, *clusterContainer);
121 caloClusterLinks.reserve( assoClusters.size() );
122 for(
const xAOD::CaloCluster* cluster : assoClusters) caloClusterLinks.emplace_back( *clusterContainer,cluster->index() );
133 if (caloExtension ==
nullptr ) trkParamDecor( *
tp ) = nullptr ;
138 ATH_MSG_DEBUG(
" Total number of selected tracks: " << ntracks );
140 return StatusCode::SUCCESS;
◆ extraDeps_update_handler()
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
◆ extraOutputDeps()
const DataObjIDColl & AthAlgorithm::extraOutputDeps |
( |
| ) |
const |
|
overridevirtualinherited |
Return the list of extra output dependencies.
This list is extended to include symlinks implied by inheritance relations.
Definition at line 50 of file AthAlgorithm.cxx.
57 return Algorithm::extraOutputDeps();
◆ initialize()
StatusCode TrackParticleClusterAssociationAlg::initialize |
( |
| ) |
|
|
override |
◆ inputHandles()
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.
◆ msg() [1/2]
◆ msg() [2/2]
◆ msgLvl()
◆ outputHandles()
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()
◆ renounceArray()
◆ sysInitialize()
StatusCode AthAlgorithm::sysInitialize |
( |
| ) |
|
|
overridevirtualinherited |
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ updateVHKA()
◆ m_assocClustersDecor
◆ m_caloClusters
◆ m_caloEntryParsDecor
◆ m_caloExtKey
◆ m_detectorEtaDecor
◆ m_detStore
◆ m_doDetEta
bool TrackParticleClusterAssociationAlg::m_doDetEta |
|
private |
◆ m_dr
Gaudi::Property<float> TrackParticleClusterAssociationAlg::m_dr {this, "DeltaR", 0.1, "max deltaR to match track and cluster" } |
|
private |
◆ m_evtStore
◆ m_extendedExtraObjects
DataObjIDColl AthAlgorithm::m_extendedExtraObjects |
|
privateinherited |
◆ m_ptCut
Gaudi::Property<double> TrackParticleClusterAssociationAlg::m_ptCut {this, "PtCut", 25000., ""} |
|
private |
◆ m_trackParticleCollectionHandle
◆ m_trackvertexassoTool
◆ m_useCovariance
Gaudi::Property<bool> TrackParticleClusterAssociationAlg::m_useCovariance {this, "UseCovariance", true, "use covariance from TrkParameters to evaluate angular uncertainties"} |
|
private |
◆ m_varHandleArraysDeclared
◆ m_vertexContHandle
◆ m_vhka
The documentation for this class was generated from the following files:
@ CENTER_MAG
Cluster Centroid ( )
@ SECOND_R
Second Moment in .
Gaudi::Property< bool > m_useCovariance
SG::ReadHandleKey< xAOD::VertexContainer > m_vertexContHandle
Scalar phi() const
phi method
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_trackParticleCollectionHandle
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_assocClustersDecor
Tracking class to hold the extrapolation from a particle from the ID to the muon system (or the other...
Scalar eta() const
pseudorapidity method
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
TrackParticleContainer_v1 TrackParticleContainer
Definition of the current "TrackParticle container version".
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
const TrackParameters * caloEntryLayerIntersection() const
access to intersection with the calorimeter entry layer return NULL if the intersection failed
Helper class to provide constant type-safe access to aux data.
const std::string & key() const
Return the StoreGate ID for the referenced object.
bool empty() const
Test if the key is blank.
virtual void setOwner(IDataHandleHolder *o)=0
SG::ReadDecorHandleKey< xAOD::CaloClusterContainer > m_detectorEtaDecor
Description of a calorimeter cluster.
double deltaPhi(double phiA, double phiB)
delta Phi in range [-pi,pi[
virtual StatusCode sysInitialize() override
Override sysInitialize.
Handle class for reading a decoration on an object.
Helper class to provide type-safe access to aux data.
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_caloEntryParsDecor
Gaudi::Property< float > m_dr
::StatusCode StatusCode
StatusCode definition for legacy code.
bool dPhi(const xAOD::TauJet &tau, const xAOD::TauTrack &track, double &out)
Handle class for adding a decoration to an object.
SG::ReadHandleKey< CaloExtensionCollection > m_caloExtKey
Gaudi::Property< double > m_ptCut
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
std::vector< const xAOD::CaloCluster * > associatedClusters(const Trk::CaloExtension &caloExtensions, const xAOD::CaloClusterContainer &allClusters)
returns the clusters from allClusters which are close enough to caloExtensions
DataObjIDColl m_extendedExtraObjects
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Class describing a Vertex.
#define ATH_MSG_WARNING(x)
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
AthAlgorithm()
Default constructor:
Class describing a TrackParticle.
ToolHandle< CP::ITrackVertexAssociationTool > m_trackvertexassoTool
bool dEta(const xAOD::TauJet &tau, const xAOD::TauTrack &track, double &out)
cl
print [x.__class__ for x in toList(dqregion.getSubRegions()) ]
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
SG::ReadHandleKey< xAOD::CaloClusterContainer > m_caloClusters