 |
ATLAS Offline Software
|
Retrieves the Jet
objects
(AntiKt4TopoEMJets etc.)
More...
#include <xAODJetRetriever.h>
|
| xAODJetRetriever (const std::string &type, const std::string &name, const IInterface *parent) |
| Standard Constructor. More...
|
|
virtual StatusCode | retrieve (ToolHandle< IFormatTool > &FormatTool) |
| For each jet collections retrieve basic parameters. More...
|
|
virtual StatusCode | initialize () |
|
const DataMap | getData (const xAOD::JetContainer *, const std::string &jetkey) |
| Puts the variables into a DataMap. More...
|
|
virtual std::string | dataTypeName () const |
| Return the name of the data type that is generated by this retriever. 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, V, H > &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::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &) |
| specialization for handling Gaudi::Property<SG::VarHandleKey> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &) |
| specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &) |
| specialization for handling Gaudi::Property<SG::VarHandleBase> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &) |
| specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
|
|
Retrieves the Jet
objects
(AntiKt4TopoEMJets etc.)
- Properties
- JetCollections
- WriteJetQuality
- TracksName
- BTaggerNames
- CDIPaths
- Retrieved Data
- Usual four-vectors: phi, eta, et etc.
- Associations for clusters and tracks via ElementLink: key/index scheme
Definition at line 36 of file xAODJetRetriever.h.
◆ StoreGateSvc_t
◆ xAODJetRetriever()
JiveXML::xAODJetRetriever::xAODJetRetriever |
( |
const std::string & |
type, |
|
|
const std::string & |
name, |
|
|
const IInterface * |
parent |
|
) |
| |
Standard Constructor.
This is the standard AthAlgTool constructor.
- Parameters
-
type | AlgTool type name |
name | AlgTool instance name |
parent | AlgTools parent owning this tool |
Definition at line 25 of file xAODJetRetriever.cxx.
◆ dataTypeName()
virtual std::string JiveXML::xAODJetRetriever::dataTypeName |
( |
| ) |
const |
|
inlinevirtual |
◆ 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]
◆ 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
◆ getData()
Puts the variables into a DataMap.
Retrieve basic parameters, mainly four-vectors, for each collection.
Also association with clusters and tracks (ElementLink).
Definition at line 84 of file xAODJetRetriever.cxx.
107 DataVect trackContKey; trackContKey.reserve(jetCont->
size());
108 DataVect trackLinkCount; trackLinkCount.reserve(jetCont->
size());
117 for (
const auto jet : *jetCont) {
122 size_t numConstit =
jet->numConstituents();
123 std::vector<std::string> tempCellID;
124 size_t trackcounter = 0;
133 for (
size_t i = 0;
i < numConstit;
i++) {
145 size_t numConstit =
jet->numConstituents();
146 for (
size_t i = 0;
i < numConstit;
i++) {
153 if (thisTrack !=
nullptr) {
165 for (
size_t j = 0; j < numConstit; ++j) {
170 if (
std::find(tempCellID.begin(), tempCellID.end(),
std::to_string(
cc->caloDDE()->identify().get_compact())) != tempCellID.end()) {
173 cellID.emplace_back(
DataType(
cc->caloDDE()->identify().get_compact()));
174 tempCellID.emplace_back(
std::to_string(
cc->caloDDE()->identify().get_compact()));
180 std::vector<const xAOD::TrackParticle*> ghosttracks;
182 if (ghosttracks.empty()) {
185 for (
size_t i = 0;
i < ghosttracks.size();
i++) {
193 ATH_MSG_VERBOSE(
" Associated track: d0 = " << ghosttracks[
i]->
d0() <<
", pt = " << ghosttracks[
i]->
pt());
195 trackcounter = ghosttracks.size();
198 for (
size_t j = 0; j < numConstit; ++j) {
213 trackLinkCount.emplace_back(
DataType(trackcounter));
220 idVec.emplace_back(
DataType( ++
id ));
230 if (jetkey!=
"AntiKt4EMPFlowJets" || (
m_nTaggers==0)){
232 bTagName.emplace_back(
DataType(
"None"));
233 bTagValue.emplace_back(
DataType(0.));
240 ATH_MSG_DEBUG(
"Failed to get btagging weight for tagger " << taggerName);
243 bTagName.emplace_back(
DataType(taggerName));
244 bTagValue.emplace_back(
DataType(btagValue));
258 std::vector<float> jvfread;
263 jvf.emplace_back(
DataType( jvfread[0] ));
271 jvt.emplace_back(
DataType(jvtread));
274 float emfracread = 0;
277 emfrac.emplace_back(
DataType( 0. ));
279 emfrac.emplace_back(
DataType( emfracread ));
291 DataMap[
"bTagName multiple=\""+str_nTaggers+
"\""] = bTagName;
292 DataMap[
"bTagValue multiple=\""+str_nTaggers+
"\""] = bTagValue;
302 if ((trackKey.size()) != 0){
303 double NTracksPerVertex = trackKey.size()*1./jetCont->size();
304 std::string
tag =
"trackIndex multiple=\"" +
DataType(NTracksPerVertex).toString()+
"\"";
306 tag =
"trackKey multiple=\"" +
DataType(NTracksPerVertex).toString()+
"\"";
310 if ((clusterID.size())!=0){
311 std::string
tag =
"clusterIndex multiple=\"" +
DataType(clusterID.size()).toString()+
"\"";
312 double NCellsPerJet = cellID.size()*1./jetCont->size();
313 tag =
"cells multiple=\"" +
DataType(NCellsPerJet).toString()+
"\"";
317 DataMap[
"trackLinkCount"] = trackLinkCount;
◆ initialize()
StatusCode JiveXML::xAODJetRetriever::initialize |
( |
| ) |
|
|
virtual |
Definition at line 28 of file xAODJetRetriever.cxx.
31 ATH_MSG_WARNING(
"Number of btaggers and CDI files do not match. Will not retrieve b-tagging information.");
32 return StatusCode::SUCCESS;
39 btagSelTool.
setTypeAndName(
"BTaggingSelectionTool/btagSelTool_"+taggerName);
49 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.
◆ interfaceID()
const InterfaceID & JiveXML::IDataRetriever::interfaceID |
( |
| ) |
|
|
inlinestaticinherited |
Return the interface identifier.
Definition at line 40 of file IDataRetriever.h.
40 {
return IID_IDataRetriever; }
◆ 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()
◆ retrieve()
StatusCode JiveXML::xAODJetRetriever::retrieve |
( |
ToolHandle< IFormatTool > & |
FormatTool | ) |
|
|
virtual |
For each jet collections retrieve basic parameters.
- Parameters
-
FormatTool | the tool that will create formated output from the DataMap |
Implements JiveXML::IDataRetriever.
Definition at line 57 of file xAODJetRetriever.cxx.
76 return StatusCode::SUCCESS;
◆ 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_bTaggerNames
Gaudi::Property<std::vector<std::string> > JiveXML::xAODJetRetriever::m_bTaggerNames |
|
private |
Initial value:{
this, "BTaggerNames", {"DL1dv01", "GN2v01"}, "Names of the b-taggers to be retrieved"}
Definition at line 62 of file xAODJetRetriever.h.
◆ m_btagSelTools
◆ m_CDIPaths
Gaudi::Property<std::vector<std::string> > JiveXML::xAODJetRetriever::m_CDIPaths |
|
private |
Initial value:{
this, "CDIPaths", {"xAODBTaggingEfficiency/13p6TeV/2023-22-13p6TeV-MC21-CDI_Test_2023-08-1_v1.root", "xAODBTaggingEfficiency/13p6TeV/2023-02_MC23_CDI_GN2v01-noSF.root"}, "Paths to the CDI files storing the b-tagger properties"}
Definition at line 64 of file xAODJetRetriever.h.
◆ m_detStore
◆ m_evtStore
◆ m_keys
SG::ReadHandleKeyArray<xAOD::JetContainer> JiveXML::xAODJetRetriever::m_keys {this,"JetCollections", {"AntiKt4EMPFlowJets","AntiKt4EMTopoJets", "AntiKt4LCTopoJets", "AntiKt10LCTopoJets", "AntiKt10UFOCSSKJets"}, "Name of the Jet containers that will be written out to the JiveXML file"} |
|
private |
◆ m_nTaggers
unsigned int JiveXML::xAODJetRetriever::m_nTaggers =0 |
|
private |
◆ m_tracksName
Gaudi::Property<std::string> JiveXML::xAODJetRetriever::m_tracksName |
|
private |
Initial value:{
this, "TracksName", "InDetTrackParticles_xAOD", "Name of the track container to be retrieved"}
Definition at line 60 of file xAODJetRetriever.h.
◆ m_typeName
const std::string JiveXML::xAODJetRetriever::m_typeName = "Jet" |
|
private |
The data type that is generated by this retriever.
Definition at line 54 of file xAODJetRetriever.h.
◆ m_varHandleArraysDeclared
◆ m_vhka
◆ m_writeJetQuality
Gaudi::Property<bool> JiveXML::xAODJetRetriever::m_writeJetQuality |
|
private |
Initial value:{
this, "WriteJetQuality", false, "Write out extended jet quality information. False by default."}
Definition at line 58 of file xAODJetRetriever.h.
The documentation for this class was generated from the following files:
virtual double phi() const
The azimuthal angle ( ) of the particle.
char data[hepevt_bytes_allocation_ATLAS]
Extra patterns decribing particle interation process.
Scalar phi() const
phi method
std::string find(const std::string &s)
return a remapped string
Scalar eta() const
pseudorapidity method
const TrackParticle * track(unsigned int index) const
Retrieve a const pointer to a Rec::TrackParticle.
std::vector< DataType > DataVect
Defines a map with a key and a vector of DataType objects e.g.
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
std::unordered_map< std::string, asg::AnaToolHandle< IBTaggingSelectionTool > > m_btagSelTools
const DataMap getData(const xAOD::JetContainer *, const std::string &jetkey)
Puts the variables into a DataMap.
const std::string m_typeName
The data type that is generated by this retriever.
#define ATH_MSG_VERBOSE(x)
std::map< std::string, DataVect > DataMap
Class providing the definition of the 4-vector interface.
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
virtual void setOwner(IDataHandleHolder *o)=0
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Description of a calorimeter cluster.
Gaudi::Property< std::string > m_tracksName
virtual double eta() const
The pseudorapidity ( ) of the particle.
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
FlowElement_v1 FlowElement
Definition of the current "pfo version".
signal_t signalType() const
@ ParticleFlow
The object is a particle-flow object.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
size_t index() const
Return the index of this element within its container.
const CaloClusterCellLink * getCellLinks() const
Get a pointer to the CaloClusterCellLink object (const version)
bool isCharged() const
is a charged PFO
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
const xAOD::IParticle * chargedObject(std::size_t i) const
Class describing a particle flow object.
virtual std::string dataTypeName() const
Return the name of the data type that is generated by this retriever.
std::string to_string(const DetectorType &type)
double charge(const T &p)
SG::ReadHandleKeyArray< xAOD::JetContainer > m_keys
@ Ok
The correction was done successfully.
#define ATH_MSG_WARNING(x)
Return value from object correction CP tools.
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Gaudi::Property< std::vector< std::string > > m_CDIPaths
Class describing a TrackParticle.
size_type size() const noexcept
Returns the number of elements in the collection.
Gaudi::Property< std::vector< std::string > > m_bTaggerNames
A detector object made of other lower level object(s)