ATLAS Offline Software
Loading...
Searching...
No Matches
HIJetConstituentModifierTool Class Reference

#include <HIJetConstituentModifierTool.h>

Inheritance diagram for HIJetConstituentModifierTool:
Collaboration diagram for HIJetConstituentModifierTool:

Public Member Functions

 HIJetConstituentModifierTool (const std::string &myname)
virtual StatusCode initialize () override
 Dummy implementation of the initialisation function.
virtual int modifyJet (xAOD::Jet &jet) const override
 Modify a single jet. This is obsolete and set to be removed.
virtual StatusCode modify (xAOD::JetContainer &jets) const
 Loop over calls to modifyJet.
virtual void print () const
 Print the state of the tool.
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
Additional helper functions, not directly mimicking Athena
template<class T>
const T * getProperty (const std::string &name) const
 Get one of the tool's properties.
const std::string & msg_level_name () const __attribute__((deprecated))
 A deprecated function for getting the message level's name.
const std::string & getName (const void *ptr) const
 Get the name of an object that is / should be in the event store.
SG::sgkey_t getKey (const void *ptr) const
 Get the (hashed) key of an object that is in the event store.

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

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

Private Attributes

SG::ReadHandleKey< xAOD::CaloClusterContainerm_clusterKey { this, "ClusterKey", "ClusterKey", "Name of the input Cluster Container"}
 Name of input cluster container.
ToolHandle< IHISubtractorToolm_subtractorTool { this, "Subtractor", "HIJetSubtractorToolBase", "" }
 handle to IHISubtractorTool that determines the subtracted kinematics for each constituent
Gaudi::Property< bool > m_originCorrection { this, "ApplyOriginCorrection", false, "Apply Origin Correction boolean switch"}
 |brief boolean switch to drive the JetScale settings after constituents are added
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

Definition at line 24 of file HIJetConstituentModifierTool.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

◆ HIJetConstituentModifierTool()

HIJetConstituentModifierTool::HIJetConstituentModifierTool ( const std::string & myname)

Definition at line 15 of file HIJetConstituentModifierTool.cxx.

15 : JetModifierBase(myname)
16{
17}
JetModifierBase(const std::string &myname)
Ctor.

Member Function Documentation

◆ ASG_TOOL_CLASS()

HIJetConstituentModifierTool::ASG_TOOL_CLASS ( HIJetConstituentModifierTool ,
IJetModifier  )
private

◆ 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

◆ getKey()

SG::sgkey_t asg::AsgTool::getKey ( const void * ptr) const
inherited

Get the (hashed) key of an object that is in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::TEvent both provide ways for getting the SG::sgkey_t key for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also
asg::AsgTool::getName
Parameters
ptrThe bare pointer to the object that the event store should know about
Returns
The hashed key of the object in the store. If not found, an invalid (zero) key.

Definition at line 119 of file AsgTool.cxx.

119 {
120
121#ifdef XAOD_STANDALONE
122 // In case we use @c xAOD::TEvent, we have a direct function call
123 // for this.
124 return evtStore()->event()->getKey( ptr );
125#else
126 const SG::DataProxy* proxy = evtStore()->proxy( ptr );
127 return ( proxy == nullptr ? 0 : proxy->sgkey() );
128#endif // XAOD_STANDALONE
129 }
ServiceHandle< StoreGateSvc > & evtStore()

◆ getName()

const std::string & asg::AsgTool::getName ( const void * ptr) const
inherited

Get the name of an object that is / should be in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::TEvent both provide ways for getting the std::string name for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also
asg::AsgTool::getKey
Parameters
ptrThe bare pointer to the object that the event store should know about
Returns
The string name of the object in the store. If not found, an empty string.

Definition at line 106 of file AsgTool.cxx.

106 {
107
108#ifdef XAOD_STANDALONE
109 // In case we use @c xAOD::TEvent, we have a direct function call
110 // for this.
111 return evtStore()->event()->getName( ptr );
112#else
113 const SG::DataProxy* proxy = evtStore()->proxy( ptr );
114 static const std::string dummy = "";
115 return ( proxy == nullptr ? dummy : proxy->name() );
116#endif // XAOD_STANDALONE
117 }

◆ getProperty()

template<class T>
const T * asg::AsgTool::getProperty ( const std::string & name) const
inherited

Get one of the tool's properties.

◆ initialize()

StatusCode HIJetConstituentModifierTool::initialize ( void )
overridevirtual

Dummy implementation of the initialisation function.

It's here to allow the dual-use tools to skip defining an initialisation function. Since many are doing so...

Reimplemented from asg::AsgTool.

Definition at line 19 of file HIJetConstituentModifierTool.cxx.

19 {
20
21 ATH_MSG_INFO("Initializing HIJetConstituentModifierTool w/ Cluster Key " << m_clusterKey.key() );
22 ATH_CHECK( m_clusterKey.initialize() );
23 return StatusCode::SUCCESS;
24
25}
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_INFO(x)
SG::ReadHandleKey< xAOD::CaloClusterContainer > m_clusterKey
Name of input cluster container.

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

◆ modify()

StatusCode JetModifierBase::modify ( xAOD::JetContainer & jets) const
virtualinherited

Loop over calls to modifyJet.

Implements IJetModifier.

Reimplemented in CenterOfMassShapesTool, EnergyCorrelatorGeneralizedTool, EnergyCorrelatorTool, HIJetClusterIndexAssociationTool, HIJetConstituentSubtractionTool, HIJetDiscriminatorTool, HIJetDRAssociationTool, KTSplittingScaleTool, NSubjettinessTool, and QwTool.

Definition at line 19 of file JetModifierBase.cxx.

19 {
20 for ( xAOD::JetContainer::iterator ijet=jets.begin(); ijet!=jets.end(); ++ijet) {
21 modifyJet(**ijet);
22 }
23 return StatusCode::SUCCESS;
24}
DataModel_detail::iterator< DataVector > iterator
Definition DataVector.h:842
virtual int modifyJet(xAOD::Jet &jet) const
Modify a single jet. This is obsolete and set to be removed.

◆ modifyJet()

int HIJetConstituentModifierTool::modifyJet ( xAOD::Jet & jet) const
overridevirtual

Modify a single jet. This is obsolete and set to be removed.

The accessor for the cluster element links

Reimplemented from JetModifierBase.

Definition at line 27 of file HIJetConstituentModifierTool.cxx.

27 {
28
29 float E_min=m_subtractorTool->minEnergyForMoments();
30 const xAOD::JetConstituentVector constituents = jet.getConstituents();
31 std::vector<size_t> cluster_indices;
32 cluster_indices.reserve(constituents.size());
33
34 for (const auto *constituent : constituents)
35 {
36 cluster_indices.push_back(constituent->rawConstituent()->index());
37 }
38
41 constituentAcc( "constituentLinks" );
43 constituentWeightAcc( "constituentWeights" );
44
45 if( constituentAcc.isAvailable(jet) ) constituentAcc( jet ).resize(0);
46 if( constituentWeightAcc.isAvailable(jet) ) constituentWeightAcc( jet ).resize(0);
47
48 //save unsubtracted kinematics as moment if they don't exist already...
49 xAOD::IParticle::FourMom_t unsubtractedP4;
50 unsubtractedP4 = jet.p4();
52
53 xAOD::IParticle::FourMom_t subtractedP4;
54 xAOD::IParticle::FourMom_t subtractedOriginCorrP4;
55 xAOD::JetFourMom_t subtractedJet4vec;
56 xAOD::JetFourMom_t subtractedOriginCorrJet4vec;
57 //Cluster container access (from the e-gamma + cluster subtracted deep copy)
58 SG::ReadHandle<xAOD::CaloClusterContainer> readHandleSubtractedClusters ( m_clusterKey );
59
60 const xAOD::CaloClusterContainer* ccl=readHandleSubtractedClusters.get();
61
62 for(auto index : cluster_indices)
63 {
64 const xAOD::CaloCluster* cl= ccl->at(index);
65
66 jet.addConstituent(cl);
67 xAOD::IParticle::FourMom_t subtractedClusterP4;
69 subtractedP4+=subtractedClusterP4;
70 ATH_MSG_DEBUG("Subracted Cluster #: " << cl->index() << " :: E: " << subtractedClusterP4.E() << " :: Eta: " << subtractedClusterP4.Eta() << " :: Phi: " << subtractedClusterP4.Phi());
71
73 xAOD::IParticle::FourMom_t subtractedOriginCorrClusterP4;
75 subtractedOriginCorrP4+=subtractedOriginCorrClusterP4;
76 ATH_MSG_DEBUG("Subrtracted OC Cluster #: " << cl->index() << " :: E: " << subtractedOriginCorrClusterP4.E() << " :: Eta: " << subtractedOriginCorrClusterP4.Eta() << " :: Phi: " << subtractedOriginCorrClusterP4.Phi());
77 }
78
79 }
80
81 //Check for subtracted Kinematics
82 if(subtractedP4.E()/std::cosh(subtractedP4.Eta()) < E_min)
83 {
84 subtractedP4=unsubtractedP4;
85 subtractedP4*=1e-7;//ghost scale
86 }
87 //Check for subtracted + Origin Correction Kinematics
88 if(subtractedOriginCorrP4.E()/std::cosh(subtractedOriginCorrP4.Eta()) < E_min && m_originCorrection )
89 {
90 subtractedOriginCorrP4=unsubtractedP4;
91 subtractedOriginCorrP4*=1e-7;//ghost scale
92 }
93
94 subtractedJet4vec.SetCoordinates(subtractedP4.Pt(),subtractedP4.Eta(),
95 subtractedP4.Phi(),subtractedP4.M());
96 jet.setJetP4(HIJetRec::subtractedJetState(), subtractedJet4vec);
97
98 if( m_originCorrection ) {
99 subtractedOriginCorrJet4vec.SetCoordinates(subtractedOriginCorrP4.Pt(),subtractedOriginCorrP4.Eta(),
100 subtractedOriginCorrP4.Phi(),subtractedOriginCorrP4.M());
101 jet.setJetP4(HIJetRec::subtractedOriginCorrectedJetState(), subtractedOriginCorrJet4vec);
102 jet.setJetP4(subtractedOriginCorrJet4vec);
104
105 }
106 else {
107 jet.setJetP4(subtractedJet4vec);
109 }
110
111 return 0;
112}
#define ATH_MSG_DEBUG(x)
const T * get(size_type n) const
Access an element, as an rvalue.
const T * at(size_type n) const
Access an element, as an rvalue.
Gaudi::Property< bool > m_originCorrection
|brief boolean switch to drive the JetScale settings after constituents are added
ToolHandle< IHISubtractorTool > m_subtractorTool
handle to IHISubtractorTool that determines the subtracted kinematics for each constituent
SG::Accessor< T, ALLOC > Accessor
Definition AuxElement.h:572
TLorentzVector FourMom_t
Definition of the 4-momentum type.
size_t size() const
number of constituents
void addConstituent(const ElementLink< IParticleContainer > &link, float weight=1.0)
Add a constituent directly in the ElementLink format.
Definition Jet_v1.cxx:111
virtual FourMom_t p4() const
The full 4-momentum of the particle.
Definition Jet_v1.cxx:71
void setJetP4(const JetFourMom_t &p4)
Definition Jet_v1.cxx:171
void setConstituentsSignalState(JetConstitScale t)
Set the state at which constituents were when this jet was found. This function is called by jet buil...
Definition Jet_v1.cxx:141
JetConstituentVector getConstituents() const
Return a vector of consituents. The object behaves like vector<const IParticle*>. See JetConstituentV...
Definition Jet_v1.cxx:147
JetFourMom_t jetP4() const
The full 4-momentum of the particle : internal jet type.
Definition Jet_v1.cxx:76
constexpr const char * unsubtractedJetState()
constexpr xAOD::JetConstitScale subtractedOriginCorrectedConstitState()
constexpr const char * subtractedJetState()
TLorentzVector getClusterP4(const xAOD::CaloCluster *cl, xAOD::CaloCluster::State s)
: getClusterP4 Added during Rel22 migration.
constexpr const char * subtractedOriginCorrectedJetState()
constexpr xAOD::CaloCluster::State subtractedClusterState()
constexpr xAOD::JetConstitScale subtractedConstitState()
constexpr xAOD::CaloCluster::State subtractedOriginCorrectedClusterState()
cl
print [x.__class__ for x in toList(dqregion.getSubRegions()) ]
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
CaloClusterContainer_v1 CaloClusterContainer
Define the latest version of the calorimeter cluster container.
ROOT::Math::LorentzVector< ROOT::Math::PtEtaPhiM4D< double > > JetFourMom_t
Base 4 Momentum type for Jet.
Definition JetTypes.h:17

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg_level_name()

const std::string & asg::AsgTool::msg_level_name ( ) const
inherited

A deprecated function for getting the message level's name.

Instead of using this, weirdly named function, user code should get the string name of the current minimum message level (in case they really need it...), with:

MSG::name( msg().level() )

This function's name doesn't follow the ATLAS coding rules, and as such will be removed in the not too distant future.

Returns
The string name of the current minimum message level that's printed

Definition at line 101 of file AsgTool.cxx.

101 {
102
103 return MSG::name( msg().level() );
104 }
MsgStream & msg() const
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition MsgLevel.cxx:19

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

◆ print()

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

SG::ReadHandleKey< xAOD::CaloClusterContainer > HIJetConstituentModifierTool::m_clusterKey { this, "ClusterKey", "ClusterKey", "Name of the input Cluster Container"}
private

Name of input cluster container.

Definition at line 39 of file HIJetConstituentModifierTool.h.

39{ this, "ClusterKey", "ClusterKey", "Name of the input Cluster Container"};

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

Gaudi::Property< bool > HIJetConstituentModifierTool::m_originCorrection { this, "ApplyOriginCorrection", false, "Apply Origin Correction boolean switch"}
private

|brief boolean switch to drive the JetScale settings after constituents are added

Definition at line 43 of file HIJetConstituentModifierTool.h.

43{ this, "ApplyOriginCorrection", false, "Apply Origin Correction boolean switch"};

◆ m_subtractorTool

ToolHandle<IHISubtractorTool> HIJetConstituentModifierTool::m_subtractorTool { this, "Subtractor", "HIJetSubtractorToolBase", "" }
private

handle to IHISubtractorTool that determines the subtracted kinematics for each constituent

Definition at line 41 of file HIJetConstituentModifierTool.h.

41{ this, "Subtractor", "HIJetSubtractorToolBase", "" };

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