ATLAS Offline Software
Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
JiveXML::xAODJetRetriever Class Reference

Retrieves the Jet objects (AntiKt4TopoEMJets etc.) More...

#include <xAODJetRetriever.h>

Inheritance diagram for JiveXML::xAODJetRetriever:
Collaboration diagram for JiveXML::xAODJetRetriever:

Public Member Functions

 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
 

Static Public Member Functions

static const InterfaceID & interfaceID ()
 Return the interface identifier. More...
 

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
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. More...
 

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

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

Private Attributes

const std::string m_typeName = "Jet"
 The data type that is generated by this retriever. More...
 
SG::ReadHandleKeyArray< xAOD::JetContainerm_keys {this,"JetCollections", {"AntiKt4EMPFlowJets","AntiKt4EMTopoJets", "AntiKt4LCTopoJets", "AntiKt10LCTopoJets", "AntiKt10UFOCSSKJets"}, "Name of the Jet containers that will be written out to the JiveXML file"}
 
Gaudi::Property< bool > m_writeJetQuality
 
Gaudi::Property< std::string > m_tracksName
 
Gaudi::Property< std::vector< std::string > > m_bTaggerNames
 
Gaudi::Property< std::vector< std::string > > m_CDIPaths
 
std::unordered_map< std::string, asg::AnaToolHandle< IBTaggingSelectionTool > > m_btagSelTools
 
unsigned int m_nTaggers =0
 
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
 

Detailed Description

Retrieves the Jet objects (AntiKt4TopoEMJets etc.)

Definition at line 36 of file xAODJetRetriever.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

◆ xAODJetRetriever()

JiveXML::xAODJetRetriever::xAODJetRetriever ( const std::string &  type,
const std::string &  name,
const IInterface *  parent 
)

Standard Constructor.

This is the standard AthAlgTool constructor.

Parameters
typeAlgTool type name
nameAlgTool instance name
parentAlgTools parent owning this tool

Definition at line 25 of file xAODJetRetriever.cxx.

25  :

Member Function Documentation

◆ dataTypeName()

virtual std::string JiveXML::xAODJetRetriever::dataTypeName ( ) const
inlinevirtual

Return the name of the data type that is generated by this retriever.

Implements JiveXML::IDataRetriever.

Definition at line 49 of file xAODJetRetriever.h.

49 { return m_typeName; };

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  hndl,
const SG::VarHandleKeyArrayType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

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  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  hndl,
const SG::VarHandleType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  t,
const SG::NotHandleType  
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleBase hndl,
const std::string &  doc,
const SG::VarHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation 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.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKey hndl,
const std::string &  doc,
const SG::VarHandleKeyType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation 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.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType  
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation 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.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation 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.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

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  }

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

95 { return m_detStore; }

◆ evtStore() [1/2]

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.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ 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

◆ getData()

const DataMap JiveXML::xAODJetRetriever::getData ( const xAOD::JetContainer jetCont,
const std::string &  jetkey 
)

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.

84  {
85 
86  ATH_MSG_DEBUG( "in getData()" );
87 
89 
90  DataVect et; et.reserve(jetCont->size());
91  DataVect phi; phi.reserve(jetCont->size());
92  DataVect eta; eta.reserve(jetCont->size());
93  DataVect mass; mass.reserve(jetCont->size());
94  DataVect energy; energy.reserve(jetCont->size());
95  DataVect bTagName; bTagName.reserve(jetCont->size());
96  DataVect bTagValue; bTagValue.reserve(jetCont->size());
97  DataVect charge; energy.reserve(jetCont->size());
98  DataVect idVec; idVec.reserve(jetCont->size());
99  DataVect px; px.reserve(jetCont->size());
100  DataVect py; py.reserve(jetCont->size());
101  DataVect pz; pz.reserve(jetCont->size());
102  DataVect jvf; jvf.reserve(jetCont->size());
103  DataVect jvt; jvt.reserve(jetCont->size());
104  DataVect emfrac; emfrac.reserve(jetCont->size());
105 
106  DataVect trackKey; trackKey.reserve(jetCont->size());
107  DataVect trackContKey; trackContKey.reserve(jetCont->size());
108  DataVect trackLinkCount; trackLinkCount.reserve(jetCont->size());
109  DataVect clusterID; clusterID.reserve(jetCont->size());
110 
111  DataVect cellID; cellID.reserve(jetCont->size());
112  DataVect numCells; numCells.reserve(jetCont->size());
113 
114  int id = 0;
115 
116  int counter = 0;
117  for (const auto jet : *jetCont) {
118  if(!jet)continue;
119  ATH_MSG_DEBUG( " Jet #" << counter++ << " : eta = " << jet->eta() << ", phi = " << jet->phi() << ", pt = " << jet->pt() );
120 
121  /* retrieve associated tracks and calo clusters */
122  size_t numConstit = jet->numConstituents();
123  std::vector<std::string> tempCellID;
124  size_t trackcounter = 0;
125  //TODO FIXME jet->rawConstituent(0) fails for some keys, including (but not exclusively) HLT_AntiKt10EMPFlowCSSKJets_nojcalib_ftf, HLT_AntiKt10EMPFlowCSSKSoftDropBeta100Zcut10Jets_jes_ftf, HLT_AntiKt4EMPFlowJets_nojcalib_ftf, HLT_AntiKt4EMPFlowJets_subjesgscIS_ftf, HLT_AntiKt4EMPFlowJets_subresjesgscIS_ftf
126  if (numConstit > 0){
127  const xAOD::IParticle* constituent = jet->rawConstituent(0);
128  if (constituent){
129  xAOD::Type::ObjectType ctype=jet->rawConstituent(0)->type();
130  // PFlow and Flow jets from athena/Reconstruction/Jet/JetMomentTools/Root/JetTrackMomentsTool.cxx
131  if (ctype == xAOD::Type::ParticleFlow) {
132  // This jet is either a PFlow jet (constituent type: xAOD::FlowElement::PFlow) or UFO jets
133  for (size_t i = 0; i < numConstit; i++) {
134  const xAOD::PFO *constit =
135  dynamic_cast<const xAOD::PFO *>(jet->rawConstituent(i));
136  if (constit->isCharged()) {
137  const xAOD::TrackParticle *thisTrack = constit->track(0); // by construction xAOD::PFO can only have one track, in eflowRec usage
138  trackKey.emplace_back(DataType(thisTrack->index()));
139  trackContKey.emplace_back(m_tracksName.value());
140  trackcounter++;
141  } // We have a charged PFO
142  } // Loop on jet constituents
143  } else if (ctype == xAOD::Type::FlowElement) {
144  // This jet is made from xAOD::FlowElement, so we calculate the pflow moments if they're PFOs
145  size_t numConstit = jet->numConstituents();
146  for (size_t i = 0; i < numConstit; i++) {
147  const xAOD::FlowElement *constit = dynamic_cast<const xAOD::FlowElement *>(jet->rawConstituent(i));
148  // UFO jet constituents have signalType xAOD::FlowElement::Charged or xAOD::FlowElement::Neutral
149  // PFlow jet constituents have signalType xAOD::FlowElement::ChargedPFlow or xAOD::FlowElement::NeutralPFlow
150  if (constit != nullptr && ((constit->signalType() & xAOD::FlowElement::PFlow) || constit->signalType() == xAOD::FlowElement::Charged)) {
151  if (constit->isCharged()) {
152  const xAOD::TrackParticle *thisTrack = dynamic_cast<const xAOD::TrackParticle *>( constit->chargedObject( 0)); // PFO should have only 1 track
153  if (thisTrack != nullptr) {
154  trackKey.emplace_back(DataType(thisTrack->index()));
155  trackContKey.emplace_back(m_tracksName.value());
156  trackcounter++;
157  }
158  else
159  ATH_MSG_WARNING( "Charged PFO had no associated TrackParticle");
160  } // We have a charged PFO
161  } // The FlowElement is a PFO
162  } // Loop on jet constituents
163  } else if (ctype == xAOD::Type::CaloCluster) {
164  // get associated cluster
165  for (size_t j = 0; j < numConstit; ++j) {
166  const xAOD::CaloCluster *cluster = dynamic_cast<const xAOD::CaloCluster *>(jet->rawConstituent(j));
167  clusterID.emplace_back(DataType(cluster->index()));
168  if(!cluster->getCellLinks())continue;
169  for (const auto cc : *(cluster->getCellLinks())) {
170  if (std::find(tempCellID.begin(), tempCellID.end(), std::to_string( cc->caloDDE()->identify().get_compact())) != tempCellID.end()) {
171  continue;
172  } else {
173  cellID.emplace_back( DataType(cc->caloDDE()->identify().get_compact()));
174  tempCellID.emplace_back( std::to_string(cc->caloDDE()->identify().get_compact()));
175  }
176  }
177  ATH_MSG_VERBOSE(" Associated cluster: eta = " << cluster->eta() << ", phi = " << cluster->phi());
178  }
179  // get ghost associated tracks
180  std::vector<const xAOD::TrackParticle*> ghosttracks;
181  jet->getAssociatedObjects<xAOD::TrackParticle>(xAOD::JetAttribute::GhostTrack, ghosttracks);
182  if (ghosttracks.empty()) {
183  ATH_MSG_VERBOSE(" Associated track: ERROR");
184  } else {
185  for (size_t i = 0; i < ghosttracks.size(); i++) {
186 
187  // can access the base track class, should be able to get tracker hits ?
188  // const Trk::Track* baseTrack = dynamic_cast< const Trk::Track* >( ghosttracks[i]->track());
189 
190  trackKey.emplace_back(DataType(ghosttracks[i]->index()));
191  trackContKey.emplace_back(m_tracksName.value());
192 
193  ATH_MSG_VERBOSE(" Associated track: d0 = " << ghosttracks[i]->d0() << ", pt = " << ghosttracks[i]->pt());
194  }
195  trackcounter = ghosttracks.size();
196  }
197  } else if (ctype == xAOD::Type::TrackParticle) {
198  for (size_t j = 0; j < numConstit; ++j) {
199  const xAOD::TrackParticle *track =
200  dynamic_cast<const xAOD::TrackParticle *>(jet->rawConstituent(j));
201  if (!track) {
202  ATH_MSG_VERBOSE(" Associated track: ERROR");
203  } else {
204  trackKey.emplace_back(DataType(track->index()));
205  trackContKey.emplace_back(m_tracksName.value());
206  trackcounter++;
207  ATH_MSG_VERBOSE(" Associated track: d0 = " << track->d0() << ", pt = " << track->pt());
208  }
209  }
210  }
211  }
212  }
213  trackLinkCount.emplace_back(DataType(trackcounter));
214  numCells.emplace_back(DataType(tempCellID.size()));
215 
216 
217  phi.emplace_back(DataType(jet->phi()));
218  eta.emplace_back(DataType(jet->eta()));
219  et.emplace_back(DataType(jet->pt()/GeV)); // hack ! no et in xAOD_Jet_v1 currently
220  idVec.emplace_back( DataType( ++id ));
221 
222  mass.emplace_back(DataType(jet->m()/GeV));
223  energy.emplace_back( DataType(jet->e()/GeV ) );
224 
225  px.emplace_back(DataType(jet->px()/GeV));
226  py.emplace_back(DataType(jet->py()/GeV));
227  pz.emplace_back(DataType(jet->pz()/GeV));
228 
229  // bjet tagger values
230  if (jetkey!="AntiKt4EMPFlowJets" || (m_nTaggers==0)){
231  for (auto taggerName : m_bTaggerNames) {
232  bTagName.emplace_back(DataType("None"));
233  bTagValue.emplace_back(DataType(0.));
234  }
235  }else{
236  double btagValue;
237  for (auto taggerName : m_bTaggerNames) {
238  CP::CorrectionCode code = m_btagSelTools[taggerName]->getTaggerWeight(*jet, btagValue);
239  if (code != CP::CorrectionCode::Ok) {
240  ATH_MSG_DEBUG("Failed to get btagging weight for tagger " << taggerName);
241  btagValue = 0;
242  }
243  bTagName.emplace_back(DataType(taggerName));
244  bTagValue.emplace_back(DataType(btagValue));
245  }
246  }
247 
248  float chargeread;
249  if (!jet->getAttribute<float>(xAOD::JetAttribute::Charge, chargeread)) {
250  ATH_MSG_DEBUG("Jet charge unavailable!");
251  charge.emplace_back( DataType( 0. ));
252  }else{
253  charge.emplace_back( DataType( chargeread ));
254  }
255 
256  // updated for data15
257  // from: Reconstruction/MET/METReconstruction/Root/METJetFilterTool.cxx
258  std::vector<float> jvfread;
259  if(!jet->getAttribute<std::vector<float> >(xAOD::JetAttribute::JVF,jvfread)) {
260  ATH_MSG_DEBUG("Jet JVF unavailable!");
261  jvf.emplace_back( DataType( 1. ));
262  }else{
263  jvf.emplace_back( DataType( jvfread[0] ));
264  }
265 
266  float jvtread;
267  if(!jet->getAttribute<float>(xAOD::JetAttribute::Jvt,jvtread)) {
268  ATH_MSG_DEBUG("Jet JVT unavailable!");
269  jvt.emplace_back(DataType(0.));
270  } else {
271  jvt.emplace_back(DataType(jvtread));
272  }
273 
274  float emfracread = 0;
275  if(!jet->getAttribute(xAOD::JetAttribute::EMFrac,emfracread)) {
276  ATH_MSG_DEBUG("Jet EMFrac unavailable!");
277  emfrac.emplace_back( DataType( 0. ));
278  }else{
279  emfrac.emplace_back( DataType( emfracread ));
280  }
281 
282  } // end loop
283 
284  // four-vectors
285  DataMap["phi"] = phi;
286  DataMap["eta"] = eta;
287  DataMap["et"] = et;
288  DataMap["energy"] = energy;
289  DataMap["mass"] = mass;
290  std::string str_nTaggers = m_nTaggers>0 ? std::to_string(m_nTaggers) : "1"; // default to 1 if no btaggers so that atlantis can process the jets properly
291  DataMap["bTagName multiple=\""+str_nTaggers+"\""] = bTagName; // assigned by hand !
292  DataMap["bTagValue multiple=\""+str_nTaggers+"\""] = bTagValue;
293  DataMap["charge"] = charge;
294  DataMap["id"] = idVec;
295  DataMap["px"] = px;
296  DataMap["py"] = py;
297  DataMap["pz"] = pz;
298  DataMap["jvf"] = jvf;
299  DataMap["jvt"] = jvt;
300  DataMap["emfrac"] = emfrac;
301 
302  if ((trackKey.size()) != 0){
303  double NTracksPerVertex = trackKey.size()*1./jetCont->size();
304  std::string tag = "trackIndex multiple=\"" +DataType(NTracksPerVertex).toString()+"\"";
305  DataMap[tag] = trackKey;
306  tag = "trackKey multiple=\"" +DataType(NTracksPerVertex).toString()+"\"";
307  DataMap[tag] = trackContKey;
308  }
309 
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()+"\"";
314  DataMap[tag]=cellID;
315  }
316 
317  DataMap["trackLinkCount"] = trackLinkCount;
318  DataMap["numCells"] = numCells;
319 
320  ATH_MSG_DEBUG( dataTypeName() << " retrieved with " << phi.size() << " entries" );
321 
322  return DataMap;
323  }

◆ initialize()

StatusCode JiveXML::xAODJetRetriever::initialize ( )
virtual

Definition at line 28 of file xAODJetRetriever.cxx.

28  {
29  ATH_CHECK(m_keys.initialize());
30  if (m_bTaggerNames.size()!=m_CDIPaths.size()){
31  ATH_MSG_WARNING("Number of btaggers and CDI files do not match. Will not retrieve b-tagging information.");
32  return StatusCode::SUCCESS;
33  } else {
34  m_nTaggers = m_bTaggerNames.size();
35  }
36  for (unsigned int i=0; i<m_nTaggers; i++){
37  std::string taggerName = m_bTaggerNames[i];
39  btagSelTool.setTypeAndName("BTaggingSelectionTool/btagSelTool_"+taggerName);
40  ATH_CHECK(btagSelTool.setProperty("TaggerName", taggerName));
41  ATH_CHECK(btagSelTool.setProperty("JetAuthor", "AntiKt4EMPFlowJets")); // only AntiKt4EMPFlowJets is supported
42  ATH_CHECK(btagSelTool.setProperty("OperatingPoint", "FixedCutBEff_70")); // the working point doesn't matter because we don't cut on the tagger discriminant, but it must exist in CDI.
43  ATH_CHECK(btagSelTool.setProperty("FlvTagCutDefinitionsFileName", m_CDIPaths[i]));
44  ATH_CHECK(btagSelTool.setProperty( "MinPt", 0.0));
45  ATH_CHECK(btagSelTool.initialize());
46  m_btagSelTools.emplace(taggerName, btagSelTool);
47  }
48 
49  return StatusCode::SUCCESS;
50  }

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

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

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

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

◆ 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();
383  PBASE::renounce (h);
384  }

◆ 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  {
365  handlesArray.renounce();
366  }

◆ retrieve()

StatusCode JiveXML::xAODJetRetriever::retrieve ( ToolHandle< IFormatTool > &  FormatTool)
virtual

For each jet collections retrieve basic parameters.

Parameters
FormatToolthe tool that will create formated output from the DataMap

Implements JiveXML::IDataRetriever.

Definition at line 57 of file xAODJetRetriever.cxx.

57  {
58 
59  ATH_MSG_DEBUG( "in retrieve()" );
60 
61  // Loop through the keys and retrieve the corresponding data
62  for (const auto& key : m_keys) {
64  if (cont.isValid()) {
65  DataMap data = getData(&(*cont), key.key());
66  if (FormatTool->AddToEvent(dataTypeName(), key.key() + "_xAOD", &data).isFailure()) {
67  ATH_MSG_WARNING("Failed to add collection " << key.key());
68  } else {
69  ATH_MSG_DEBUG(" (" << key.key() << ") retrieved");
70  }
71  } else {
72  ATH_MSG_WARNING("Collection " << key.key() << " not found in SG");
73  }
74  }
75 
76  return StatusCode::SUCCESS;
77  }

◆ 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 DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ 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) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ 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

std::unordered_map<std::string, asg::AnaToolHandle<IBTaggingSelectionTool> > JiveXML::xAODJetRetriever::m_btagSelTools
private

Definition at line 67 of file xAODJetRetriever.h.

◆ 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

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_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

Definition at line 56 of file xAODJetRetriever.h.

◆ m_nTaggers

unsigned int JiveXML::xAODJetRetriever::m_nTaggers =0
private

Definition at line 68 of file xAODJetRetriever.h.

◆ 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

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.

◆ 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:
asg::AnaToolHandle::initialize
StatusCode initialize()
initialize the tool
xAOD::CaloCluster_v1::phi
virtual double phi() const
The azimuthal angle ( ) of the particle.
Definition: CaloCluster_v1.cxx:256
data
char data[hepevt_bytes_allocation_ATLAS]
Definition: HepEvt.cxx:11
et
Extra patterns decribing particle interation process.
GeV
#define GeV
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:18
asg::AnaToolHandle< IBTaggingSelectionTool >
JiveXML::xAODJetRetriever::m_nTaggers
unsigned int m_nTaggers
Definition: xAODJetRetriever.h:68
test_pyathena.px
px
Definition: test_pyathena.py:18
phi
Scalar phi() const
phi method
Definition: AmgMatrixBasePlugin.h:67
xAOD::FlowElement_v1::Charged
@ Charged
Definition: FlowElement_v1.h:38
find
std::string find(const std::string &s)
return a remapped string
Definition: hcg.cxx:135
Base_Fragment.mass
mass
Definition: Sherpa_i/share/common/Base_Fragment.py:59
xAOD::et
et
Definition: TrigEMCluster_v1.cxx:25
ObjectType
ObjectType
Definition: BaseObject.h:11
eta
Scalar eta() const
pseudorapidity method
Definition: AmgMatrixBasePlugin.h:83
SG::ReadHandle
Definition: StoreGate/StoreGate/ReadHandle.h:67
xAOD::PFO_v1::track
const TrackParticle * track(unsigned int index) const
Retrieve a const pointer to a Rec::TrackParticle.
Definition: PFO_v1.cxx:691
JiveXML::DataVect
std::vector< DataType > DataVect
Defines a map with a key and a vector of DataType objects e.g.
Definition: DataType.h:58
DataType
OFFLINE_FRAGMENTS_NAMESPACE::PointerType DataType
Definition: RoIBResultByteStreamTool.cxx:25
asg::AnaToolHandle::setTypeAndName
void setTypeAndName(const std::string &val_typeAndName)
set the value of type and name
test_pyathena.pt
pt
Definition: test_pyathena.py:11
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
JiveXML::xAODJetRetriever::m_btagSelTools
std::unordered_map< std::string, asg::AnaToolHandle< IBTaggingSelectionTool > > m_btagSelTools
Definition: xAODJetRetriever.h:67
JiveXML::xAODJetRetriever::getData
const DataMap getData(const xAOD::JetContainer *, const std::string &jetkey)
Puts the variables into a DataMap.
Definition: xAODJetRetriever.cxx:84
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
JiveXML::xAODJetRetriever::m_typeName
const std::string m_typeName
The data type that is generated by this retriever.
Definition: xAODJetRetriever.h:54
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
JiveXML::DataMap
std::map< std::string, DataVect > DataMap
Definition: DataType.h:59
xAOD::IParticle
Class providing the definition of the 4-vector interface.
Definition: Event/xAOD/xAODBase/xAODBase/IParticle.h:41
asg::AnaToolHandle::setProperty
StatusCode setProperty(const std::string &property, const T2 &value)
set the given property of the tool.
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
xAOD::CaloCluster
CaloCluster_v1 CaloCluster
Define the latest version of the calorimeter cluster class.
Definition: Event/xAOD/xAODCaloEvent/xAODCaloEvent/CaloCluster.h:19
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition: Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
xAOD::FlowElement_v1::PFlow
@ PFlow
Definition: FlowElement_v1.h:45
xAOD::FlowElement_v1::isCharged
bool isCharged() const
Definition: FlowElement_v1.cxx:56
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
xAOD::CaloCluster_v1
Description of a calorimeter cluster.
Definition: CaloCluster_v1.h:62
histSizes.code
code
Definition: histSizes.py:129
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
jet
Definition: JetCalibTools_PlotJESFactors.cxx:23
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
ParticleGun_FastCalo_ChargeFlip_Config.energy
energy
Definition: ParticleGun_FastCalo_ChargeFlip_Config.py:78
JiveXML::xAODJetRetriever::m_tracksName
Gaudi::Property< std::string > m_tracksName
Definition: xAODJetRetriever.h:60
xAOD::CaloCluster_v1::eta
virtual double eta() const
The pseudorapidity ( ) of the particle.
Definition: CaloCluster_v1.cxx:251
lumiFormat.i
int i
Definition: lumiFormat.py:85
xAOD::JetAttribute::EMFrac
@ EMFrac
Definition: JetAttributes.h:112
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition: AthCommonDataStore.h:145
xAOD::FlowElement
FlowElement_v1 FlowElement
Definition of the current "pfo version".
Definition: FlowElement.h:16
xAOD::FlowElement_v1::signalType
signal_t signalType() const
TRT::Track::d0
@ d0
Definition: InnerDetector/InDetCalibEvent/TRT_CalibData/TRT_CalibData/TrackInfo.h:62
Amg::pz
@ pz
Definition: GeoPrimitives.h:40
test_pyathena.parent
parent
Definition: test_pyathena.py:15
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
xAODType::ParticleFlow
@ ParticleFlow
The object is a particle-flow object.
Definition: ObjectType.h:41
xAOD::JetAttribute::JVF
@ JVF
Definition: JetAttributes.h:99
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
SG::AuxElement::index
size_t index() const
Return the index of this element within its container.
xAOD::CaloCluster_v1::getCellLinks
const CaloClusterCellLink * getCellLinks() const
Get a pointer to the CaloClusterCellLink object (const version)
Definition: CaloCluster_v1.cxx:859
AthAlgTool::AthAlgTool
AthAlgTool()
Default constructor:
xAOD::PFO_v1::isCharged
bool isCharged() const
is a charged PFO
Definition: PFO_v1.cxx:251
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
xAOD::FlowElement_v1::chargedObject
const xAOD::IParticle * chargedObject(std::size_t i) const
Definition: FlowElement_v1.cxx:127
Amg::py
@ py
Definition: GeoPrimitives.h:39
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
xAOD::PFO_v1
Class describing a particle flow object.
Definition: PFO_v1.h:35
JiveXML::xAODJetRetriever::dataTypeName
virtual std::string dataTypeName() const
Return the name of the data type that is generated by this retriever.
Definition: xAODJetRetriever.h:49
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
ActsTrk::to_string
std::string to_string(const DetectorType &type)
Definition: GeometryDefs.h:34
charge
double charge(const T &p)
Definition: AtlasPID.h:986
xAOD::JetAttribute::Charge
@ Charge
Definition: JetAttributes.h:182
JiveXML::xAODJetRetriever::m_keys
SG::ReadHandleKeyArray< xAOD::JetContainer > m_keys
Definition: xAODJetRetriever.h:56
CP::CorrectionCode::Ok
@ Ok
The correction was done successfully.
Definition: CorrectionCode.h:38
DeMoScan.index
string index
Definition: DeMoScan.py:362
a
TList * a
Definition: liststreamerinfos.cxx:10
xAOD::TauJetParameters::numCells
@ numCells
Definition: TauDefs.h:171
h
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
CP::CorrectionCode
Return value from object correction CP tools.
Definition: CorrectionCode.h:31
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:629
xAOD::JetAttribute::Jvt
@ Jvt
Definition: JetAttributes.h:101
JiveXML::xAODJetRetriever::m_CDIPaths
Gaudi::Property< std::vector< std::string > > m_CDIPaths
Definition: xAODJetRetriever.h:64
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:801
xAOD::track
@ track
Definition: TrackingPrimitives.h:513
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
CaloCondBlobAlgs_fillNoiseFromASCII.tag
string tag
Definition: CaloCondBlobAlgs_fillNoiseFromASCII.py:23
test_pyathena.counter
counter
Definition: test_pyathena.py:15
xAOD::JetAttribute::GhostTrack
@ GhostTrack
Definition: JetAttributes.h:252
DataVector::size
size_type size() const noexcept
Returns the number of elements in the collection.
JiveXML::xAODJetRetriever::m_bTaggerNames
Gaudi::Property< std::vector< std::string > > m_bTaggerNames
Definition: xAODJetRetriever.h:62
fitman.k
k
Definition: fitman.py:528
python.handimod.cc
int cc
Definition: handimod.py:522
xAOD::FlowElement_v1
A detector object made of other lower level object(s)
Definition: FlowElement_v1.h:25
mapkey::key
key
Definition: TElectronEfficiencyCorrectionTool.cxx:37