![Logo](../../ATLAS-Logo-Square-Blue-RGB.png) |
ATLAS Offline Software
|
#include <JetQuarkLabel.h>
|
| JetQuarkLabel (const std::string &name) |
|
virtual | ~JetQuarkLabel () |
|
virtual StatusCode | initialize () override |
| Dummy implementation of the initialisation function. More...
|
|
virtual bool | matchJet (const xAOD::Jet &myJet, MatchInfo *info=nullptr) const override |
| AlgTool interface methods. More...
|
|
virtual void | printParameterSettings () const override |
| print parameter settings of the truth match tool More...
|
|
void | EventSelection (short s) |
|
short | EventSelection () const |
|
virtual void | print () const |
| Print the state of the tool. 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 | 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 |
|
Definition at line 38 of file JetQuarkLabel.h.
◆ StoreGateSvc_t
◆ JetQuarkLabel()
Analysis::JetQuarkLabel::JetQuarkLabel |
( |
const std::string & |
name | ) |
|
◆ ~JetQuarkLabel()
Analysis::JetQuarkLabel::~JetQuarkLabel |
( |
| ) |
|
|
virtual |
◆ 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()
◆ EventSelection() [1/2]
short Analysis::JetQuarkLabel::EventSelection |
( |
| ) |
const |
|
inline |
◆ EventSelection() [2/2]
void Analysis::JetQuarkLabel::EventSelection |
( |
short |
s | ) |
|
|
inline |
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ 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
◆ getKey()
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
-
ptr | The 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.
121 #ifdef XAOD_STANDALONE
124 return evtStore()->event()->getKey( ptr );
127 return (
proxy ==
nullptr ? 0 :
proxy->sgkey() );
128 #endif // XAOD_STANDALONE
◆ 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
-
ptr | The 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.
108 #ifdef XAOD_STANDALONE
111 return evtStore()->event()->getName( ptr );
114 static const std::string
dummy =
"";
116 #endif // XAOD_STANDALONE
◆ getProperty()
template<class T >
const T* asg::AsgTool::getProperty |
( |
const std::string & |
name | ) |
const |
|
inherited |
Get one of the tool's properties.
◆ initialize()
StatusCode Analysis::JetQuarkLabel::initialize |
( |
| ) |
|
|
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 42 of file JetQuarkLabel.cxx.
44 return StatusCode::SUCCESS;
◆ 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.
◆ matchJet()
AlgTool interface methods.
Method to truth tag a jet:
Implements Analysis::IJetTruthMatching.
Definition at line 47 of file JetQuarkLabel.cxx.
55 if (!truthEventContainerReadHandle.isValid()){
56 ATH_MSG_DEBUG(
" Invalid ReadHandle for xAOD::TruthEventContainer with key: " << truthEventContainerReadHandle.key());
60 return testJet(myJet,truthEventContainerReadHandle.cptr(),
info);
◆ msg() [1/2]
◆ msg() [2/2]
◆ 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.
◆ 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.
◆ print()
void asg::AsgTool::print |
( |
| ) |
const |
|
virtualinherited |
Print the state of the tool.
Implements asg::IAsgTool.
Reimplemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, top::TopObjectSelection, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
Definition at line 131 of file AsgTool.cxx.
◆ printParameterSettings()
void Analysis::JetQuarkLabel::printParameterSettings |
( |
| ) |
const |
|
overridevirtual |
◆ renounce()
◆ renounceArray()
◆ sysInitialize()
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ testJet()
Definition at line 64 of file JetQuarkLabel.cxx.
69 TLorentzVector jet_hlv = myJet.
p4();
71 ATH_MSG_DEBUG(
"Truth matching for jet " << myJet.
pt() <<
" " << myJet.
eta() <<
" " << myJet.
phi() <<
" using pile-up event selection " <<
m_inTime);
75 int NEventInCollection = truthEventContainer->
size();
77 info->NEventInCollection =NEventInCollection;
82 double deltaRC = 999.;
83 double deltaRB = 999.;
84 double deltaRT = 999.;
88 unsigned int nLab = 0;
94 for (; itEvt != itEvtE; ++itEvt) {
96 std::vector<ElementLink<xAOD::TruthParticleContainer> >::const_iterator ELpitr = GenEvent->
truthParticleLinks().begin();
97 std::vector<ElementLink<xAOD::TruthParticleContainer> >::const_iterator pitrE = GenEvent->
truthParticleLinks().end();
98 for (; ELpitr != pitrE; ++ELpitr) {
100 if(!pitr.
isValid() || !(*pitr)) {
105 int pdg = (*pitr)->pdgId();
108 TLorentzVector part_momentum_lv = (*pitr)->p4();
111 if (std::abs(pdg) == 5 || std::abs(pdg) == 4) {
112 double pt = part_momentum_lv.Pt();
117 <<
" pT= " <<part_momentum_lv.Pt()
118 <<
" eta= " <<part_momentum_lv.Eta()
119 <<
" phi= " <<part_momentum_lv.Phi()
120 <<
" dR= " <<part_momentum_lv.DeltaR(jet_hlv));
122 bool afterFSR =
true;
123 if ((*pitr)->hasDecayVtx()) {
124 std::vector<ElementLink<xAOD::TruthParticleContainer> >::const_iterator firstChild = (*pitr)->decayVtx()->outgoingParticleLinks().begin();
125 std::vector<ElementLink<xAOD::TruthParticleContainer> >::const_iterator endChild = (*pitr)->decayVtx()->outgoingParticleLinks().begin();
126 std::vector<ElementLink<xAOD::TruthParticleContainer> >::const_iterator thisChild = firstChild;
127 for(; thisChild != endChild; ++thisChild){
128 if(!thisChild->isValid() || !(**thisChild)) {
133 if ((**thisChild)->pdgId() == pdg) afterFSR =
false;
142 deltaR=part_momentum_lv.DeltaR(jet_hlv);
148 if (std::abs(pdg) == 15) {
149 double pt = part_momentum_lv.Pt();
152 <<
" pT= " <<part_momentum_lv.Pt()
153 <<
" eta= " <<part_momentum_lv.Eta()
154 <<
" phi= " <<part_momentum_lv.Phi()
155 <<
" dR= " <<part_momentum_lv.DeltaR(jet_hlv) );
156 deltaR=part_momentum_lv.DeltaR(jet_hlv);
163 ATH_MSG_DEBUG(
"Number of events in the EventCollection : " << NEventInCollection <<
" and used for labelling " << nLab);
166 info->distanceToQuarks.insert(std::make_pair(
"B",deltaRB));
167 info->distanceToQuarks.insert(std::make_pair(
"C",deltaRC));
168 info->distanceToQuarks.insert(std::make_pair(
"T",deltaRT));
170 ATH_MSG_VERBOSE(
"DeltaR " << deltaRB <<
" " << deltaRC <<
" " << deltaRT);
174 info->barcode = barcb;
177 ATH_MSG_VERBOSE(
"Jet matched with a b "<<barcb<<
" after FSR, dR: " << deltaRB);
180 if (LabellingParticle == 0) {
181 ATH_MSG_WARNING(
"A b labelled jet without a labelling particle ? Should not exist");
188 info->barcode = barcc;
191 ATH_MSG_VERBOSE(
"Jet matched with a c "<<barcc<<
" after FSR, dR: " << deltaRC);
◆ updateVHKA()
◆ m_deltaRCut
double Analysis::JetQuarkLabel::m_deltaRCut |
|
private |
deltaR cut value of the cone matching (max distance between Jet axis and momentum of truth particel)
Definition at line 62 of file JetQuarkLabel.h.
◆ m_detStore
◆ m_evtStore
◆ m_inTime
short Analysis::JetQuarkLabel::m_inTime |
|
private |
◆ m_mcEventCollection
std::string Analysis::JetQuarkLabel::m_mcEventCollection |
|
private |
◆ m_noDoc
bool Analysis::JetQuarkLabel::m_noDoc |
|
private |
◆ m_ptCut
double Analysis::JetQuarkLabel::m_ptCut |
|
private |
◆ m_truthEventContainerKey
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
Const iterator class for DataVector/DataList.
double m_deltaRCut
deltaR cut value of the cone matching (max distance between Jet axis and momentum of truth particel)
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
bool testJet(const xAOD::Jet &, const xAOD::TruthEventContainer *, MatchInfo *info) const
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
virtual double phi() const
The azimuthal angle ( ) of the particle.
#define ATH_MSG_VERBOSE(x)
bool isValid() const
Test to see if the link can be dereferenced.
virtual void setOwner(IDataHandleHolder *o)=0
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
constexpr int SPECIALSTATUS
Constant that the meaning of which is currently lost, to be recovered...
Class describing a truth particle in the MC record.
SG::ReadHandleKey< xAOD::TruthEventContainer > m_truthEventContainerKey
Class describing a signal truth event in the MC record.
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
virtual double eta() const
The pseudorapidity ( ) of the particle.
virtual void renounce()=0
constexpr int UNDEFINED_ID
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
ElementLink implementation for ROOT usage.
#define GeVtoMeV
Purpose: label jets with b or c quarks.
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
virtual FourMom_t p4() const
The full 4-momentum of the particle.
#define ATH_MSG_WARNING(x)
double m_ptCut
pT cut for partons
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
const TruthParticleLinks_t & truthParticleLinks() const
Get all the truth particles.
virtual double pt() const
The transverse momentum ( ) of the particle.
size_type size() const noexcept
Returns the number of elements in the collection.
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.