|
ATLAS Offline Software
|
Properties: VertexContainer - name of the vertex container EventInfo - name of EventInfo container AssociatedTracks - name for attribute holding the list of associated tracks TVATool - tool to do track-vertex association TrkSelTool - tool to select tracks (none ==> no selection)
More...
#include <JetQGTaggerVariableTool.h>
|
| JetQGTaggerVariableTool (const std::string &name) |
|
StatusCode | initialize () override |
| Dummy implementation of the initialisation function. More...
|
|
virtual StatusCode | decorate (const xAOD::JetContainer &jetCont) const override |
| Decorate a jet collection without otherwise modifying it. More...
|
|
const xAOD::Vertex * | findHSVertex (const xAOD::VertexContainer *&) 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 |
|
virtual StatusCode | modify (xAOD::JetContainer &jets) const override final |
| Concrete implementation of the function inherited from IJetModifier. More...
|
|
|
Gaudi::Property< std::string > | m_jetContainerName {this,"JetContainer", "", "SG key for the input jet container"} |
|
SG::ReadHandleKey< xAOD::VertexContainer > | m_vertexContainer_key {this, "VertexContainer", "PrimaryVertices", "SG key for input vertex container"} |
|
SG::ReadHandleKey< xAOD::EventInfo > | m_eventInfo_key {this, "EventInfo", "EventInfo", "SG key for input EventInfo"} |
|
ToolHandle< InDet::IInDetTrackSelectionTool > | m_trkSelectionTool {this, "TrackSelector", "", "Track selector"} |
|
SG::ReadHandleKey< jet::TrackVertexAssociation > | m_tva_key {this, "TrackVertexAssociation", "", "Track vertex association key"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_nTrkKey {this, "NTrksDecorName", "DFCommonJets_QGTagger_NTracks", "SG key for output NTracks decoration"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_trkWidthKey {this, "TracksWidthDecorName", "DFCommonJets_QGTagger_TracksWidth", "SG key for output TracksWidth decoration"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_trkC1Key {this, "TracksC1DecorName", "DFCommonJets_QGTagger_TracksC1", "SG key for output TracksC1 decoration"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_nChargedKey {this, "nChargedDecorName", "DFCommonJets_QGTagger_truthjet_nCharged", "SG key for output truthjet_nCharged decoration"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_truthPtKey {this, "truthPtDecorName", "DFCommonJets_QGTagger_truthjet_pt", "SG key for output truthjet_pt decoration"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_truthEtaKey {this, "truthEtaDecorName", "DFCommonJets_QGTagger_truthjet_eta", "SG key for output truthjet_eta decoration"} |
|
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 |
|
Properties: VertexContainer - name of the vertex container EventInfo - name of EventInfo container AssociatedTracks - name for attribute holding the list of associated tracks TVATool - tool to do track-vertex association TrkSelTool - tool to select tracks (none ==> no selection)
Definition at line 41 of file JetQGTaggerVariableTool.h.
◆ StoreGateSvc_t
◆ JetQGTaggerVariableTool()
JetQGTaggerVariableTool::JetQGTaggerVariableTool |
( |
const std::string & |
name | ) |
|
◆ 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);
253 return PBASE::declareProperty(
name,hndl,
doc);
◆ 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.
229 return PBASE::declareProperty(
name,hndl,
doc);
◆ 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.
338 return PBASE::declareProperty(
name, property,
doc);
◆ 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]
◆ decorate()
Decorate a jet collection without otherwise modifying it.
Implements IJetDecorator.
Definition at line 61 of file JetQGTaggerVariableTool.cxx.
73 if (!vertexContainer.isValid()){
75 return StatusCode::FAILURE;
78 const auto *vertices = vertexContainer.cptr();
82 if (vertices->empty() ) {
86 nTrkHandle(*
jet) = -1;
87 trkWidthHandle(*
jet) = -1.;
88 trkC1Handle(*
jet) = -1.;
89 nChargedHandle(*
jet) = -1;
90 truthPtHandle(*
jet) = -1.;
91 truthEtaHandle(*
jet) = -1.;
94 return StatusCode::SUCCESS;
102 if (!handle_tva.isValid()){
104 return StatusCode::FAILURE;
107 const auto *tva = handle_tva.cptr();
112 int nTracksCount = 0;
113 float TracksWidth = 0., SumTracks_pTs = 0., TracksC1 = 0.,
beta = 0.2;
116 std::vector<const xAOD::IParticle*> jettracks;
120 std::vector<bool> IsGoodTrack;
121 TLorentzVector tracki_TLV, trackj_TLV;
122 TLorentzVector jet_TLV =
jet -> p4();
123 for (
size_t i = 0;
i < jettracks.size();
i++) {
126 IsGoodTrack.push_back(
false);
139 (HSvertex == tva->associatedVertex(trk) || (HSvertex != tva->associatedVertex(trk) && std::abs((trk->
z0()+trk->
vz()-HSvertex->
z())*
sin(trk->
theta()))<3.))
142 IsGoodTrack.push_back(
accept);
147 tracki_TLV = trk -> p4();
148 double DR_tracki_jet = tracki_TLV.DeltaR(jet_TLV);
149 TracksWidth += trk ->
pt() * DR_tracki_jet;
150 SumTracks_pTs += trk ->
pt();
154 if(SumTracks_pTs>0.) TracksWidth = TracksWidth / SumTracks_pTs;
155 else TracksWidth = -1.;
158 for(
size_t i = 0;
i < jettracks.size();
i++) {
160 if( !( IsGoodTrack.at(
i) ) )
continue;
162 for(
size_t j =
i+1; j < jettracks.size(); j++) {
164 if( !( IsGoodTrack.at(j) ) )
continue;
166 tracki_TLV = trki -> p4();
167 trackj_TLV = trkj -> p4();
168 double DR_tracki_trackj = tracki_TLV.DeltaR(trackj_TLV);
169 TracksC1 += trki ->
pt() * trkj ->
pt() *
pow( DR_tracki_trackj,
beta) ;
174 if(SumTracks_pTs>0.) TracksC1 = TracksC1 / (
pow(SumTracks_pTs, 2.) );
180 if (!eventInfoContainer.isValid()){
182 return StatusCode::FAILURE;
185 const auto *eventInfo = eventInfoContainer.cptr();
189 float truthjet_pt = -999.0;
190 float truthjet_eta = -999.0;
195 if(GhostTruthAssociationLinkAcc.isAvailable(*
jet) ){
196 ATH_MSG_DEBUG(
"Accessing GhostTruthAssociationLink: is available");
199 ATH_MSG_DEBUG(
"Accessing GhostTruthAssociationLink: is valid");
211 ATH_MSG_DEBUG(
"Cannot access truth Link: setting weight to 1");
217 truthjet_pt = tjet->
pt();
218 truthjet_eta = tjet->
eta();
228 if( ! (
part->pt()>500.) )
continue;
231 if( !(
part->isCharged()) )
continue;
237 nTrkHandle(*
jet) = nTracksCount;
238 trkWidthHandle(*
jet) = TracksWidth;
239 trkC1Handle(*
jet) = TracksC1;
240 nChargedHandle(*
jet) = tntrk;
241 truthPtHandle(*
jet) = truthjet_pt;
242 truthEtaHandle(*
jet) = truthjet_eta;
246 return StatusCode::SUCCESS;
◆ detStore()
◆ 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
◆ findHSVertex()
Definition at line 251 of file JetQGTaggerVariableTool.cxx.
253 for (
size_t iVertex = 0; iVertex < vertices->
size(); ++iVertex ) {
256 ATH_MSG_VERBOSE(
"JetQGTaggerVariableTool " <<
name() <<
" Found HS vertex at index: "<< iVertex);
257 return vertices->
at(iVertex);
260 ATH_MSG_VERBOSE(
"There is no vertex of type PriVx. Taking default vertex.");
261 return vertices->
at(0);
◆ 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
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
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 JetQGTaggerVariableTool::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 29 of file JetQGTaggerVariableTool.cxx.
33 ATH_MSG_ERROR(
"JetQGTaggerVariableTool needs to have its input jet container configured!");
34 return StatusCode::FAILURE;
56 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.
◆ modify()
|
inlinefinaloverridevirtualinherited |
◆ 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, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
Definition at line 131 of file AsgTool.cxx.
◆ 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.
◆ updateVHKA()
◆ m_detStore
◆ m_eventInfo_key
◆ m_evtStore
◆ m_jetContainerName
Gaudi::Property<std::string> JetQGTaggerVariableTool::m_jetContainerName {this,"JetContainer", "", "SG key for the input jet container"} |
|
private |
◆ m_nChargedKey
◆ m_nTrkKey
◆ m_trkC1Key
◆ m_trkSelectionTool
◆ m_trkWidthKey
◆ m_truthEtaKey
◆ m_truthPtKey
◆ m_tva_key
◆ m_varHandleArraysDeclared
◆ m_vertexContainer_key
◆ m_vhka
The documentation for this class was generated from the following files:
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
StatusCode accept(const xAOD::Muon *mu)
float vz() const
The z origin for the parameters.
float z0() const
Returns the parameter.
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
Helper class to provide constant type-safe access to aux data.
#define ATH_MSG_VERBOSE(x)
const std::string & key() const
Return the StoreGate ID for the referenced object.
Class providing the definition of the 4-vector interface.
@ IS_SIMULATION
true: simulation, false: data
bool isValid() const
Test to see if the link can be dereferenced.
virtual void setOwner(IDataHandleHolder *o)=0
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
bool is_simulation_particle(const T &p)
Method to establish if a particle (or barcode) was created during the simulation (TODO update to be s...
Class describing a truth particle in the MC record.
Handle class for adding a decoration to an object.
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)
float z() const
Returns the z position.
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
virtual double eta() const
The pseudorapidity ( ) of the particle.
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
ElementLink implementation for ROOT usage.
const IParticle * rawConstituent(size_t i) const
Direct access to constituents. WARNING expert use only.
bool isStable(const T &p)
Identify if the particle is stable, i.e. has not decayed.
Class describing a Vertex.
virtual StatusCode decorate(const xAOD::JetContainer &jets) const =0
Decorate a jet collection without otherwise modifying it.
#define ATH_MSG_WARNING(x)
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Class describing a TrackParticle.
const T * at(size_type n) const
Access an element, as an rvalue.
size_t numConstituents() const
Number of constituents in this jets (this is valid even when reading a file where the constituents ha...
virtual double pt() const
The transverse momentum ( ) of the particle.
float theta() const
Returns the parameter, which has range 0 to .
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>