ATLAS Offline Software
Public Types | Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
top::TopxAODAlg Class Referencefinal

#include <TopxAODAlg.h>

Inheritance diagram for top::TopxAODAlg:
Collaboration diagram for top::TopxAODAlg:

Public Types

typedef ServiceHandle< StoreGateSvc > & MetaStorePtr_t
 Type of the metadata store pointer in standalone mode. More...
 
typedef const ServiceHandle< StoreGateSvc > & ConstMetaStorePtr_t
 
typedef ServiceHandle< StoreGateSvc > MetaStore_t
 Type of the metadata store variable in Athena. More...
 

Public Member Functions

 TopxAODAlg (const std::string &name, ISvcLocator *pSvcLocator)
 
StatusCode initialize () override
 
StatusCode beginInputFile () override
 
StatusCode execute () override
 
StatusCode endInputFile () override
 
StatusCode finalize () override
 
::StatusCode requestFileExecute ()
 register this algorithm to have an implementation of fileexecute More...
 
::StatusCode requestBeginInputFile ()
 register this algorithm to have an implementation of beginInputFile More...
 
::StatusCode requestEndInputFile ()
 register this algorithm to have an implementation of endInputFile More...
 
void handle (const Incident &inc)
 receive the given incident More...
 
virtual StatusCode sysInitialize ()
 Initialization method invoked by the framework. More...
 
const ServiceHandle< ITHistSvc > & histSvc () const
 The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc. More...
 
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies. 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 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
 

Protected Member Functions

virtual void print () const
 print the state of the algorithm More...
 
virtual ::StatusCode fileExecute ()
 perform the action exactly once for each file in the dataset More...
 
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...
 
StatusCode configAthHistogramming (const ServiceHandle< ITHistSvc > &histSvc, const std::string &prefix, const std::string &rootDir, const std::string &histNamePrefix, const std::string &histNamePostfix, const std::string &histTitlePrefix, const std::string &histTitlePostfix)
 To be called by the derived classes to fill the internal configuration. More...
 
TH1 * bookGetPointer (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms. More...
 
TH1 * bookGetPointer (TH1 *hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms. More...
 
TH1 * bookGetPointer (TH1 &histRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of histograms. More...
 
TTree * bookGetPointer (const TTree &treeRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TTrees. More...
 
TGraph * bookGetPointer (const TGraph &graphRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs. More...
 
TEfficiency * bookGetPointer (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency. More...
 
TEfficiency * bookGetPointer (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency. More...
 
TEfficiency * bookGetPointer (TEfficiency &effRef, std::string tDir="", std::string stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency. More...
 
StatusCode book (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms. More...
 
StatusCode book (TH1 *hist, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms. More...
 
StatusCode book (TH1 &histRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of histograms. More...
 
StatusCode book (const TTree &treeRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TTrees. More...
 
StatusCode book (const TGraph &graphRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TGraphs. More...
 
StatusCode book (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency. More...
 
StatusCode book (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency. More...
 
StatusCode book (TEfficiency &effRef, const std::string &tDir="", const std::string &stream="")
 Simplify the booking and registering (into THistSvc) of TEfficiency. More...
 
TH1 * hist (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered histograms of any type. More...
 
TH2 * hist2d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered 2-d histograms. More...
 
TH3 * hist3d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered 3-d histograms. More...
 
TTree * tree (const std::string &treeName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TTrees. More...
 
TGraph * graph (const std::string &graphName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TGraphs. More...
 
TEfficiency * efficiency (const std::string &effName, const std::string &tDir="", const std::string &stream="")
 Simplify the retrieval of registered TEfficiency. More...
 

Private Types

typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
 
typedef uint32_t hash_t
 typedef for the internal hash More...
 
typedef std::map< const hash_t, TH1 * > HistMap_t
 Typedef for convenience. More...
 
typedef std::map< const hash_t, TEfficiency * > EffMap_t
 Typedef for convenience. More...
 
typedef std::map< const hash_t, TTree * > TreeMap_t
 Typedef for convenience. More...
 
typedef std::map< const hash_t, TGraph * > GraphMap_t
 Typedef for convenience. More...
 

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 
void buildBookingString (std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
 Method to build individual booking string. More...
 
void myReplace (std::string &str, const std::string &oldStr, const std::string &newStr)
 Helper method to replace sub-string. More...
 
hash_t hash (const std::string &histName) const
 Method to calculate a 32-bit hash from a string. More...
 

Private Attributes

std::shared_ptr< ConfigurationSettings > m_settings {nullptr}
 
std::shared_ptr< TopConfig > m_config {nullptr}
 
ToolHandle< top::TopToolStore > m_topTools
 
ToolHandle< top::EventCleaningSelection > m_eventCleaning
 
ToolHandle< top::ObjectCollectionMaker > m_systObjMaker
 
ToolHandle< top::TopObjectSelection > m_objectSelection
 
std::unique_ptr< top::EventSelectionManager > m_eventSelectionManager
 
std::unique_ptr< top::AnalysisTrackingHelper > m_tracker
 
ToolHandle< top::ParticleLevelLoader > m_particleLevelLoader
 
ToolHandle< top::CalcTtbarPartonHistory > m_topPartonHistory
 
ToolHandle< top::PDFScaleFactorCalculator > m_PDF_SF
 
ToolHandle< top::TopEventMaker > m_topEventMaker
 
ToolHandle< top::ScaleFactorCalculator > m_topScaleFactors
 
std::unique_ptr< top::EventSaverBase > m_eventSaver
 
ToolHandle< PMGTools::IPMGTruthWeightTool > m_pmg_weightTool {"PMGTruthWeightTool"}
 
unsigned int m_doPerfStats = 0u
 
TFile * m_outputFile = nullptr
 
TTree * m_sumPdfWeights = nullptr
 
std::unordered_map< std::string, std::unique_ptr< std::vector< float > > > m_totalEventsPdfWeighted
 
bool m_pdfMetadataExists = false
 
double m_totalEventsWeighted = 0
 
double m_totalEventsWeighted_temp = 0
 
ULong64_t m_totalEvents = 0
 
ULong64_t m_totalEventsInFiles = 0
 
TTree * m_sumWeights = nullptr
 
bool m_recalc_LHE3 = false
 
bool m_recalculateNominalWeightSum = false
 
std::vector< std::string > m_names_LHE3
 
std::vector< double > m_totalEventsWeighted_LHE3
 
std::vector< double > m_totalEventsWeighted_LHE3_temp
 
unsigned int m_entry = 0u
 
unsigned int m_firstEvent = 0u
 
unsigned int m_totalYieldSoFar = 0u
 
unsigned int m_eventSavedReco = 0u
 
unsigned int m_eventSavedRecoLoose = 0u
 
unsigned int m_eventSavedTruth = 0u
 
unsigned int m_eventSavedParticle = 0u
 
bool m_isFirst = true
 
unsigned int m_DSID
 
int m_isAFII
 
std::unordered_map< std::string, std::vector< std::string > > m_boostedTaggersSFSysNames
 
std::string m_configPath {""}
 properties More...
 
std::string m_libraryNames {""}
 
std::string m_derivationStream {""}
 
bool m_crashOnMetaDataFail {true}
 
std::string m_streamName {"output"}
 
unsigned int m_FTAGMapIndex
 
SampleXsection::showering m_showeringAlgo
 
std::string m_year
 
MetaStore_t m_inputMetaStore
 Object accessing the input metadata store. More...
 
MetaStore_t m_outputMetaStore
 Object accessing the output metadata store. More...
 
bool m_hasFileExecute {false}
 the value of hasFileExecute More...
 
bool m_hasBeginInputFile {false}
 the value of hasBeginInputFile More...
 
bool m_hasEndInputFile {false}
 the value of hasEndInputFile More...
 
ServiceHandle< ITHistSvc > m_histSvc
 Default constructor: AthHistogramAlgorithm();. More...
 
std::string m_prefix
 Name of the ROOT output stream (file) More...
 
std::string m_rootDir
 Name of the ROOT directory. More...
 
std::string m_histNamePrefix
 The prefix for the histogram THx name. More...
 
std::string m_histNamePostfix
 The postfix for the histogram THx name. More...
 
std::string m_histTitlePrefix
 The prefix for the histogram THx title. More...
 
std::string m_histTitlePostfix
 The postfix for the histogram THx title. More...
 
DataObjIDColl m_extendedExtraObjects
 
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
 
HistMap_t m_histMap
 The map of histogram names to their pointers. More...
 
EffMap_t m_effMap
 The map of histogram names to their pointers. More...
 
TreeMap_t m_treeMap
 The map of TTree names to their pointers. More...
 
GraphMap_t m_graphMap
 The map of TGraph names to their pointers. More...
 
std::string m_name
 Instance name. More...
 
MsgStream m_msg
 Cached Message Stream. More...
 

Detailed Description

Definition at line 39 of file TopxAODAlg.h.

Member Typedef Documentation

◆ ConstMetaStorePtr_t

typedef const ServiceHandle< StoreGateSvc >& EL::AnaAlgorithm::ConstMetaStorePtr_t
inherited

Definition at line 111 of file AnaAlgorithm.h.

◆ EffMap_t

typedef std::map< const hash_t, TEfficiency* > AthHistogramming::EffMap_t
privateinherited

Typedef for convenience.

Definition at line 205 of file AthHistogramming.h.

◆ GraphMap_t

typedef std::map< const hash_t, TGraph* > AthHistogramming::GraphMap_t
privateinherited

Typedef for convenience.

Definition at line 219 of file AthHistogramming.h.

◆ hash_t

typedef uint32_t AthHistogramming::hash_t
privateinherited

typedef for the internal hash

Definition at line 169 of file AthHistogramming.h.

◆ HistMap_t

typedef std::map< const hash_t, TH1* > AthHistogramming::HistMap_t
privateinherited

Typedef for convenience.

Definition at line 198 of file AthHistogramming.h.

◆ MetaStore_t

typedef ServiceHandle< StoreGateSvc > EL::AnaAlgorithm::MetaStore_t
inherited

Type of the metadata store variable in Athena.

Definition at line 546 of file AnaAlgorithm.h.

◆ MetaStorePtr_t

typedef ServiceHandle< StoreGateSvc >& EL::AnaAlgorithm::MetaStorePtr_t
inherited

Type of the metadata store pointer in standalone mode.

Definition at line 110 of file AnaAlgorithm.h.

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

◆ TreeMap_t

typedef std::map< const hash_t, TTree* > AthHistogramming::TreeMap_t
privateinherited

Typedef for convenience.

Definition at line 212 of file AthHistogramming.h.

Constructor & Destructor Documentation

◆ TopxAODAlg()

top::TopxAODAlg::TopxAODAlg ( const std::string &  name,
ISvcLocator *  pSvcLocator 
)

Definition at line 22 of file TopxAODAlg.cxx.

23  : AnaAlgorithm(name, pSvcLocator) {
24  declareProperty("ConfigPath", m_configPath, "Path to config file to parse by top-el");
25  declareProperty("LibraryNames", m_libraryNames, "Names of loaded .so libraries");
26  declareProperty("DerivationStream", m_derivationStream, "DerivationStream of this sample");
27  declareProperty("CrashOnMetaDataFaile", m_crashOnMetaDataFail,
28  "Crash if metadata reading failed");
29  declareProperty("StreamName", m_streamName, "The name of the output stream");
30  declareProperty("FTAGMapIndex", m_FTAGMapIndex,
31  "Index of the showering algorithm for FTAG MC/MC SFs");
32  declareProperty("ShoweringAlgorithm", m_showeringAlgo,
33  "Name of the showering algorithm for FTAG MC/MC SFs");
34  declareProperty("DataTakingYear", m_year,
35  "Year of data taking (primarily for trigger configuration)");
36  declareProperty("SkipEvents", m_firstEvent, "From which event entry to start (skip before)");
37 }

Member Function Documentation

◆ beginInputFile()

StatusCode top::TopxAODAlg::beginInputFile ( )
override

Definition at line 232 of file TopxAODAlg.cxx.

232  {
233  // Sum of weights and raw number of entries before skimming in current file
234  double sumW_file = 0;
235  ULong64_t initialEvents = 0;
236 
237  // vector of MC generator weights and their names in current file
238  std::vector<float> LHE3_sumW_file;
239  std::vector<std::string> LHE3_names_file;
240  const xAOD::CutBookkeeperContainer *cutBookKeepers = nullptr;
241 
242  if (m_config->isTruthDxAOD()) {
243  ATH_MSG_INFO("Bookkeepers are not read for TRUTH derivations");
244  } else {
245  top::check(inputMetaStore()->retrieve(cutBookKeepers, "CutBookkeepers"),
246  "Cannot retrieve CutBookkeepers");
247  if (m_config->isMC()) {
248 
249  // here we attempt to name the CutBookkeepers based on the MC weight names
250  const std::vector<std::string> &weight_names = m_pmg_weightTool->getWeightNames();
251  // try to retrieve CutBookKeepers for LHE3Weights first
252  top::parseCutBookkeepers(*evtStore()->event(), weight_names.size(), LHE3_names_file,
253  LHE3_sumW_file);
254  // if we have MC generator weights, we rename the bookkeepers in sumWeights TTree to match the
255  // weight names from MetaData
256  top::renameCutBookkeepers(LHE3_names_file, weight_names);
257 
258  // raw number of events taken from "AllExecutedEvents" bookkeeper, which corresponds to 0th MC
259  // weight but these are raw entries, so doesn't matter if 0th MC weight is nominal or not
260  initialEvents = top::getRawEventsBookkeeper(cutBookKeepers);
261 
262  // determine the nominal sum of weight -- we already found the nominal weight in
263  // ScaleFactorCalculator
264  const size_t nominalWeightIndex = m_config->nominalWeightIndex();
265  sumW_file = LHE3_sumW_file.at(nominalWeightIndex);
266  } else {
267  initialEvents = top::getRawEventsBookkeeper(cutBookKeepers);
268  sumW_file = initialEvents; // this is data, it's the same number...
269  }
270  }
271 
272  m_totalEventsWeighted += sumW_file;
273  m_totalEvents += initialEvents;
274 
275  // now we must fill two vectors in sync for MCGeneratorWeights sum of weights
276  if (m_config->doMCGeneratorWeights() && m_config->isMC()) {
277  if (m_totalEventsWeighted_LHE3.size() != 0) {
278  if (m_totalEventsWeighted_LHE3.size() != LHE3_sumW_file.size() ||
279  m_names_LHE3.size() != LHE3_names_file.size() ||
280  m_names_LHE3.size() != m_totalEventsWeighted_LHE3.size()) {
281  ATH_MSG_ERROR("Strange inconsistency of vector sizes in sum of LHE3 weights calculation.");
282  return StatusCode::FAILURE;
283  }
284  for (unsigned int i_genweights = 0; i_genweights < LHE3_names_file.size(); i_genweights++) {
285  if (m_names_LHE3.at(i_genweights) != LHE3_names_file.at(i_genweights)) {
286  ATH_MSG_ERROR("Strange inconsistency in the vector of weight names in sum of LHE3 "
287  "weights calculation.");
288  return StatusCode::FAILURE;
289  } else {
290  m_totalEventsWeighted_LHE3.at(i_genweights) =
291  m_totalEventsWeighted_LHE3.at(i_genweights) + LHE3_sumW_file.at(i_genweights);
292  }
293  }
294  } else {
295  for (unsigned int i_genweights = 0; i_genweights < LHE3_names_file.size(); i_genweights++) {
296  m_names_LHE3.push_back(LHE3_names_file.at(i_genweights));
297  m_totalEventsWeighted_LHE3.push_back(LHE3_sumW_file.at(i_genweights));
298  }
299  }
300  if (!m_names_LHE3.empty()) {
301  ATH_MSG_INFO(
302  "The sum of weights for the following LHE3 weights were retrieved from the input file:");
303  MsgStream &msgInfo = msg(MSG::Level::INFO);
304  for (const std::string& s : m_names_LHE3)
305  msgInfo << s << " ";
306  msgInfo << std::endl;
307  msgInfo.doOutput();
308  } else {
309  ATH_MSG_INFO(
310  "No sum of LHE3 weights could be found in meta-data. Will try to recompute these sums.\n"
311  "This only works on un-skimmed derivations, and the names of these weights may be "
312  "unknown (but we'll try to read them from the PMG tool");
313  m_recalc_LHE3 = true;
314  }
315  }
316 
317  if (m_config->isTruthDxAOD())
318  m_recalculateNominalWeightSum = true;
319 
320  if (m_config->printCDIpathWarning()) {
321  ATH_MSG_WARNING(
322  "\n*************************************************************************\n"
323  << "YOU ARE USING A CUSTOM PATH TO THE CDI FILE WHICH IS NOT THE DEFAULT PATH\n"
324  << " YOU MANY NOT BE USING THE LATEST BTAGGING RECOMMENDATIONS \n"
325  << "*************************************************************************\n\n");
326  }
327  if (m_config->printEgammaCalibModelWarning()) {
328  ATH_MSG_WARNING(
329  "\n*************************************************************************\n"
330  << " YOU HAVE CHANGED DEFAULT EGAMMA CALIBRATION MODEL \n"
331  << " TO USE DEFAULT MODEL, REMOVE 'EGammaCalibrationModel' FROM CONFIG FILE \n"
332  << "*************************************************************************\n\n");
333  }
334 
335  if (m_config->printEIDFileWarning()) {
336  ATH_MSG_WARNING(
337  "\n*************************************************************************\n"
338  << " YOU ARE USING THIS CUSTOM PATH TO THE ELECTRON ID SF FILE: \n\n"
339  << m_config->electronIDSFFilePath() << "(Tight) \n\n"
340  << m_config->electronIDSFFileLoosePath() << "(Loose) \n\n"
341  << " INSTEAD OF THE MOST RECENT RECOMMENDED MAP \n"
342  << " YOU MANY NOT BE USING THE LATEST ELECTRON ID RECOMMENDATIONS \n"
343  << "*************************************************************************\n\n");
344  }
345 
346  return StatusCode::SUCCESS;
347 }

◆ book() [1/8]

StatusCode AthHistogramming::book ( const TEfficiency &  eff,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 335 of file AthHistogramming.h.

336 {
337  // We need to create a non-const clone
338  TEfficiency* effClone = dynamic_cast< TEfficiency* >( eff.Clone() );
339  if ( !effClone ) {
340  m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone" << endmsg;
341  return StatusCode::FAILURE;
342  }
343  return this->book( *effClone, tDir, stream );
344 }

◆ book() [2/8]

StatusCode AthHistogramming::book ( const TGraph &  graphRef,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TGraphs.

◆ book() [3/8]

StatusCode AthHistogramming::book ( const TH1 &  hist,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 303 of file AthHistogramming.h.

304 {
305  // We need to create a non-const clone
306  TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
307  if ( !histClone ) {
308  m_msg << MSG::ERROR << "Couldn't create a TH1 clone" << endmsg;
309  return StatusCode::FAILURE;
310  }
311  return this->book( *histClone, tDir, stream );
312 }

◆ book() [4/8]

StatusCode AthHistogramming::book ( const TTree &  treeRef,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TTrees.

Definition at line 403 of file AthHistogramming.h.

404 {
405  // Call the other Book method and see if it returns a valid pointer
406  TTree* treePointer = this->bookGetPointer( treeRef, tDir, stream );
407  if ( treePointer )
408  {
409  return StatusCode::SUCCESS;
410  }
411  else
412  {
413  return StatusCode::FAILURE;
414  }
415 }

◆ book() [5/8]

StatusCode AthHistogramming::book ( TEfficiency &  effRef,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 356 of file AthHistogramming.h.

357 {
358  // Call the other Book method and see if it returns a valid pointer
359  TEfficiency* effPointer = this->bookGetPointer( effRef, tDir, stream );
360  if ( !effPointer ) {
361  m_msg << MSG::ERROR << "Couldn't book a TEfficiency" << endmsg;
362  return StatusCode::FAILURE;
363  }
364  return StatusCode::SUCCESS;
365 }

◆ book() [6/8]

StatusCode AthHistogramming::book ( TEfficiency *  eff,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 346 of file AthHistogramming.h.

347 {
348  if ( !eff ) {
349  m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency" << endmsg;
350  return StatusCode::FAILURE;
351  }
352  return this->book( *eff, tDir, stream );
353 }

◆ book() [7/8]

StatusCode AthHistogramming::book ( TH1 &  histRef,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 324 of file AthHistogramming.h.

325 {
326  // Call the other Book method and see if it returns a valid pointer
327  TH1* histPointer = this->bookGetPointer( histRef, tDir, stream );
328  if ( !histPointer ) {
329  m_msg << MSG::ERROR << "Couldn't book a TH1" << endmsg;
330  return StatusCode::FAILURE;
331  }
332  return StatusCode::SUCCESS;
333 }

◆ book() [8/8]

StatusCode AthHistogramming::book ( TH1 *  hist,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 314 of file AthHistogramming.h.

315 {
316  if ( !hist ) {
317  m_msg << MSG::ERROR << "Got a zero pointer to a TH1" << endmsg;
318  return StatusCode::FAILURE;
319  }
320  return this->book( *hist, tDir, stream );
321 }

◆ bookGetPointer() [1/8]

TEfficiency * AthHistogramming::bookGetPointer ( const TEfficiency &  eff,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 281 of file AthHistogramming.h.

282 {
283  // We need to create a non-const clone
284  TEfficiency* histClone = dynamic_cast< TEfficiency* >( hist.Clone() );
285  if ( !histClone ) {
286  m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone in bookGetPointer" << endmsg;
287  return 0;
288  }
289  return this->bookGetPointer( *histClone, tDir, stream );
290 
291 }

◆ bookGetPointer() [2/8]

TGraph * AthHistogramming::bookGetPointer ( const TGraph &  graphRef,
std::string  tDir = "",
std::string  stream = "" 
)
protectedinherited

Simplify the booking and registering (into THistSvc) of TGraphs.

Definition at line 427 of file AthHistogramming.cxx.

428 {
429  // Get a pointer
430  const TGraph* graphPointer = &graphRef;
431 
432  // Check that we got a valid pointer
433  if ( !graphPointer )
434  {
435  m_msg << MSG::WARNING
436  << "We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" << m_name
437  << "!" << endmsg;
438  return NULL;
439  }
440 
441  // Modify the name and title according to the prefixes of this classes instance
442  std::string graphName = graphPointer->GetName();
443  const std::string graphTitle = graphPointer->GetTitle();
444 
445  // Check if the hash for this graphName already exists, i.e., if we have a hash collision
446  const hash_t graphHash = this->hash(graphName);
447  GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
448  if ( it != m_graphMap.end() ) // It does exist!
449  {
450  m_msg << MSG::WARNING
451  << "Detected a hash collision. The hash for the TGraph with name=" << graphName
452  << " already exists and points to a TGraph with name=" << it->second->GetName()
453  << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
454  return NULL;
455  }
456 
457  // Create a clone that has the new name
458  TGraph* graphClone = dynamic_cast< TGraph* >( graphPointer->Clone((m_histNamePrefix+graphName+m_histNamePostfix).c_str()) );
459  if( !graphClone )
460  {
461  m_msg << MSG::WARNING
462  << "We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" << m_name
463  << "!" << endmsg;
464  return NULL;
465  }
466  graphClone->SetTitle ((m_histTitlePrefix+graphTitle+m_histTitlePostfix).c_str());
467 
468  // Massage the final string to book things
469  std::string bookingString("");
470  this->buildBookingString( bookingString, graphName, tDir, stream );
471 
472  // Register the TGraph into the THistSvc
473  if ( !((histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
474  {
475  m_msg << MSG::WARNING
476  << "Problem registering TGraph with name " << graphName
477  << ", title " << graphTitle
478  << " in " << m_name << "!" << endmsg;
479  return NULL;
480  }
481 
482  // Also register it in the local map of string to pointer
483  m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
484 
485  return graphClone;
486 }

◆ bookGetPointer() [3/8]

TH1 * AthHistogramming::bookGetPointer ( const TH1 &  hist,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 260 of file AthHistogramming.h.

261 {
262  // We need to create a non-const clone
263  TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
264  if ( !histClone ) {
265  m_msg << MSG::ERROR << "Couldn't create a TH1 clone in bookGetPointer" << endmsg;
266  return 0;
267  }
268  return this->bookGetPointer( *histClone, tDir, stream );
269 
270 }

◆ bookGetPointer() [4/8]

TTree * AthHistogramming::bookGetPointer ( const TTree &  treeRef,
std::string  tDir = "",
std::string  stream = "" 
)
protectedinherited

Simplify the booking and registering (into THistSvc) of TTrees.

Definition at line 312 of file AthHistogramming.cxx.

313 {
314  // Get a pointer
315  const TTree* treePointer = &treeRef;
316 
317  // Check that we got a valid pointer
318  if ( !treePointer )
319  {
320  m_msg << MSG::WARNING
321  << "We got an invalid TTree pointer in the BookGetPointer(TTree*) method of the class" << m_name
322  << "!" << endmsg;
323  return NULL;
324  }
325 
326  // Modify the name and title according to the prefixes of this classes instance
327  std::string treeName = treePointer->GetName();
328  const std::string treeTitle = treePointer->GetTitle();
329 
330  // Check if the hash for this treeName already exists, i.e., if we have a hash collision
331  const hash_t treeHash = this->hash(treeName);
332  TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
333  if ( it != m_treeMap.end() ) // It does exist!
334  {
335  m_msg << MSG::WARNING
336  << "Detected a hash collision. The hash for the TTree with name=" << treeName
337  << " already exists and points to a TTree with name=" << it->second->GetName()
338  << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
339  return NULL;
340  }
341 
342  // Create a clone that has the new name
343  TTree* treeClone = dynamic_cast< TTree* >( treePointer->Clone(treeName.c_str()) );
344  if( !treeClone )
345  {
346  m_msg << MSG::WARNING
347  << "We couldn't clone the TTree in the BookGetPointer(TTree&) method of the class" << m_name
348  << "!" << endmsg;
349  return NULL;
350  }
351  treeClone->SetTitle (treeTitle.c_str());
352 
353  // Massage the final string to book things
354  std::string bookingString("");
355  this->buildBookingString( bookingString, treeName, tDir, stream );
356 
357  // Register the TTree into the THistSvc
358  if ( !((histSvc()->regTree(bookingString, treeClone)).isSuccess()) )
359  {
360  m_msg << MSG::WARNING
361  << "Problem registering TTree with name " << treeName
362  << ", title " << treeTitle
363  << " in " << m_name << "!" << endmsg;
364  return NULL;
365  }
366 
367  // Also register it in the local map of string to pointer
368  m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treeClone ) );
369 
370  return treeClone;
371 }

◆ bookGetPointer() [5/8]

TEfficiency * AthHistogramming::bookGetPointer ( TEfficiency &  effRef,
std::string  tDir = "",
std::string  stream = "" 
)
protectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 146 of file AthHistogramming.cxx.

147 {
148  // Modify the name and title according to the prefixes of this classes instance
149  std::string effName(effRef.GetName());
150  const std::string effTitle(effRef.GetTitle());
151  std::string bookingString("");
152 
153  this->buildBookingString( bookingString, effName, tDir, stream );
154  effRef.SetTitle((m_histTitlePrefix+effTitle+m_histTitlePostfix).c_str() );
155  effRef.SetName(effName.c_str());
156 
157  // Check if the hash for this effName already exists, i.e., if we have a hash collision
158  const hash_t effHash = this->hash(effName);
159  EffMap_t::const_iterator it = m_effMap.find( effHash );
160  if ( it != m_effMap.end() ) // It does exist!
161  {
162  m_msg << MSG::WARNING
163  << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
164  << " already exists and points to a TEfficiency with name=" << it->second->GetName()
165  << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
166  return NULL;
167  }
168 
169  // Set the new name and title for the TEfficiency, based on the prefixes that the user set for this class instance
170  // Create a clone that has the new name
171 
172  // Massage the final string to book things
173 
174  // Register the TEfficiency into the THistSvc
175  if ( !((histSvc()->regEfficiency(bookingString, &effRef)).isSuccess()) )
176  {
177  m_msg << MSG::WARNING
178  << "Problem registering TEfficiency with name " << effName
179  << ", name prefix " << m_histNamePrefix
180  << ", title " << effTitle
181  << ", tile prefix " << m_histTitlePrefix
182  << ", and tile postfix " << m_histTitlePostfix
183  << " in " << m_name << "!" << endmsg;
184  return NULL;
185  }
186 
187  // Also register it in the local map of string to pointer
188  m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, &effRef ) );
189 
190  return &effRef;
191 }

◆ bookGetPointer() [6/8]

TEfficiency * AthHistogramming::bookGetPointer ( TEfficiency *  eff,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 293 of file AthHistogramming.h.

294 {
295  if ( !hist ) {
296  m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency in bookGetPointer" << endmsg;
297  return 0;
298  }
299  return this->bookGetPointer( *hist, tDir, stream );
300 }

◆ bookGetPointer() [7/8]

TH1 * AthHistogramming::bookGetPointer ( TH1 &  histRef,
std::string  tDir = "",
std::string  stream = "" 
)
protectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 98 of file AthHistogramming.cxx.

99 {
100  // Modify the name and title according to the prefixes of this classes instance
101  std::string histName(histRef.GetName());
102  const std::string histTitle(histRef.GetTitle());
103  std::string bookingString("");
104 
105  this->buildBookingString( bookingString, histName, tDir, stream );
106  histRef.SetTitle((m_histTitlePrefix+histTitle+m_histTitlePostfix).c_str() );
107  histRef.SetName(histName.c_str());
108 
109  // Check if the hash for this histName already exists, i.e., if we have a hash collision
110  const hash_t histHash = this->hash(histName);
111  HistMap_t::const_iterator it = m_histMap.find( histHash );
112  if ( it != m_histMap.end() ) // It does exist!
113  {
114  m_msg << MSG::WARNING
115  << "Detected a hash collision. The hash for the histogram with name=" << histName
116  << " already exists and points to a histogram with name=" << it->second->GetName()
117  << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
118  return NULL;
119  }
120 
121  // Set the new name and title for the histogram, based on the prefixes that the user set for this class instance
122  // Create a clone that has the new name
123 
124  // Massage the final string to book things
125 
126  // Register the histogram into the THistSvc
127  if ( !((histSvc()->regHist(bookingString, &histRef)).isSuccess()) )
128  {
129  m_msg << MSG::WARNING
130  << "Problem registering histogram with name " << histName
131  << ", name prefix " << m_histNamePrefix
132  << ", title " << histTitle
133  << ", tile prefix " << m_histTitlePrefix
134  << ", and tile postfix " << m_histTitlePostfix
135  << " in " << m_name << "!" << endmsg;
136  return NULL;
137  }
138 
139  // Also register it in the local map of string to pointer
140  m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, &histRef ) );
141 
142  return &histRef;
143 }

◆ bookGetPointer() [8/8]

TH1 * AthHistogramming::bookGetPointer ( TH1 *  hist,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 272 of file AthHistogramming.h.

273 {
274  if ( !hist ) {
275  m_msg << MSG::ERROR << "Got a zero pointer to a TH1 in bookGetPointer" << endmsg;
276  return 0;
277  }
278  return this->bookGetPointer( *hist, tDir, stream );
279 }

◆ buildBookingString()

void AthHistogramming::buildBookingString ( std::string &  bookingString,
std::string &  histName,
std::string &  tDir,
std::string &  stream,
bool  usePrefixPostfix = false 
)
privateinherited

Method to build individual booking string.

Definition at line 560 of file AthHistogramming.cxx.

565 {
566  // Massage the final string to book things
567  if(tDir.empty()) tDir = m_rootDir;
568  size_t pos = histName.rfind('/');
569  if(pos != std::string::npos){
570  tDir+='/';
571  tDir.append(histName, 0,pos);
572  histName.erase(0,pos+1);
573  };
574  if(stream.empty()) stream = m_streamName;
575 
576  if(usePrefixPostfix){
577  bookingString = "/"+stream+"/"+tDir+"/"+m_histNamePrefix+histName+m_histNamePostfix;
578  } else {
579  bookingString = "/"+stream+"/"+tDir+"/"+histName;
580  }
581  while(bookingString.find("//") != std::string::npos){
582  this->myReplace(bookingString,"//","/");
583  }
584 
585  return;
586 }

◆ configAthHistogramming()

StatusCode AthHistogramming::configAthHistogramming ( const ServiceHandle< ITHistSvc > &  histSvc,
const std::string &  prefix,
const std::string &  rootDir,
const std::string &  histNamePrefix,
const std::string &  histNamePostfix,
const std::string &  histTitlePrefix,
const std::string &  histTitlePostfix 
)
protectedinherited

To be called by the derived classes to fill the internal configuration.

Definition at line 66 of file AthHistogramming.cxx.

70 {
71  m_histSvc = histSvc;
72  m_streamName = prefix;
73  m_rootDir = rootDir;
74  m_histNamePrefix = histNamePrefix;
75  m_histNamePostfix = histNamePostfix;
76  m_histTitlePrefix = histTitlePrefix;
77  m_histTitlePostfix = histTitlePostfix;
78 
79  return StatusCode::SUCCESS;
80 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  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< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  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< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  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< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::declareProperty ( Gaudi::Property< T > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
147  return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
148  }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::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; }

◆ efficiency()

TEfficiency * AthHistogramming::efficiency ( const std::string &  effName,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
protectedinherited

Simplify the retrieval of registered TEfficiency.

Definition at line 250 of file AthHistogramming.cxx.

251 {
252  // Build a 32 bit hash out of the name
253  const hash_t effHash = this->hash(effName);
254 
255  // See if this entry exists in the map
256  EffMap_t::const_iterator it = m_effMap.find( effHash );
257  if ( it == m_effMap.end() ) // It doesn't exist!
258  { // Let's see into the THistSvc if somebody else has registered the TEfficiency...
259 
260  // Need to copy the strings as we will massage them from here on
261  std::string effNameCopy = effName;
262  std::string tDirCopy = tDir;
263  std::string streamCopy = stream;
264 
265  // Massage the final string to book things
266  std::string bookingString("");
267  this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy ,false);
268 
269  TEfficiency* effPointer(NULL);
270  if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
271  {
272  // Massage the final string to book things
273  std::string bookingString("");
274  this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy, true );
275 
276  if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
277  {
278  m_msg << MSG::WARNING
279  << "Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
280  << m_histNamePrefix + effNameCopy + m_histNamePostfix
281  << " or with name " << effNameCopy
282  << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
283  << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
284  return NULL;
285  }
286  // If we get to here, we actually found the TEfficiency in the THistSvc.
287  // So let's add it to the local cache map and return its pointer
288  m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
289  return effPointer;
290  }
291  // If we get to here, we actually found the TEfficiency in the THistSvc.
292  // So let's add it to the local cache map and return its pointer
293  m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
294  return effPointer;
295  }
296 
297  // Return the pointer to the TEfficiency that we got from the local cache map
298  return it->second;
299 }

◆ endInputFile()

StatusCode top::TopxAODAlg::endInputFile ( )
override

Definition at line 656 of file TopxAODAlg.cxx.

656  {
657 
658  // TODO this calculation may be incorrect if EL is instructed to skip some events
659  if (m_tracker)
660  m_tracker->addInputFile(wk()->inputFile()->GetName(), m_entry);
661 
662  // do it at the end, so we can get the DS id from the first event
663  // notice this is different from the normal sum of weights: all entries matter, not only the
664  // highest per file
665  if (m_config->doLHAPDF()) {
666  bool pdfInCBK = false;
667  const xAOD::CutBookkeeperContainer *cutBookKeepers = nullptr;
668  if (inputMetaStore()->retrieve(cutBookKeepers, "CutBookkeepers")) {
669  // try to get the first entry
670  std::string pdf_set = m_totalEventsPdfWeighted.begin()->first;
671  std::string p = pdf_set + "_0";
672  for (const xAOD::CutBookkeeper *cbk : *cutBookKeepers) {
673  std::string pdfName = cbk->name();
674  if (p != pdfName)
675  continue;
676  m_pdfMetadataExists = true;
677  pdfInCBK = true;
678  // reload to set iterator back to original state
679  if (!inputMetaStore()->retrieve(cutBookKeepers, "CutBookkeepers")) {
680  }
681  break;
682  }
683  }
684  // if the pdf info was in the CBK the PDFSumOfWeights retrieve never happens,
685  // and m_cutBookKeepers, will keep its reference to the main CBK
686  if (!pdfInCBK) { // only try this if the info is not in the CBK -- this is legacy stuff
687  if (!inputMetaStore()->retrieve(cutBookKeepers, "PDFSumOfWeights")) {
688  m_pdfMetadataExists = false;
689  }
690  }
691  if (m_pdfMetadataExists && m_config->saveLHAPDFEvent() && m_config->baseLHAPDF().empty()) {
692  // just make a copy of it, but not all derivations have it, so just ignore it if it is not
693  // there
694  for (auto &pdfentry : m_totalEventsPdfWeighted) {
695  std::string pdf_set = pdfentry.first;
696  for (size_t n = 0; n < m_totalEventsPdfWeighted[pdf_set]->size(); ++n) {
697  std::string p = pdf_set + "_" + std::to_string(n);
698  bool foundPdf = false;
699  for (const xAOD::CutBookkeeper *cbk : *cutBookKeepers) {
700  std::string pdfName = cbk->name();
701  if (p != pdfName)
702  continue;
703  m_totalEventsPdfWeighted[pdf_set]->at(n) = cbk->sumOfEventWeights();
704  foundPdf = true;
705  break;
706  }
707  if (!foundPdf) { // asked for PDF weighting, the PDF metadata exists, but this particular
708  // PDF is missing
709  // crash hard
710  ATH_MSG_ERROR("The force is not strong with us, young Padawan ...\n"
711  "You want to save weights to do PDF reweighting using '"
712  << pdf_set << "', which I figure has "
713  << m_totalEventsPdfWeighted[pdf_set]->size()
714  << " sets on it.\n"
715  "There is metadata information for the sum of MC*PDF weights in "
716  "PDFSumOfWeights, "
717  "but none seem to refer to '"
718  << p
719  << "' therefore I do not know how to estimate "
720  "the sum of weights before acceptance for this configuration.\n"
721  "May the force be with you in your next attempt.");
722  return StatusCode::FAILURE;
723  }
724  }
725  }
726  m_sumPdfWeights->Fill();
727  }
728  } // doLHAPDF
729 
730  return StatusCode::SUCCESS;
731 }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::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< Algorithm > >::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; }

◆ execute()

StatusCode top::TopxAODAlg::execute ( )
override

– Truth events –///

– End of Truth events – start of reconstruction level events –///

– We veto events that do not pass the GRL, trigger or have bad calorimeters –/// – No point calibrating and cutting on events that fail these cuts –/// – We only apply a veto if ALL selectors request a cut –/// – And if the user selects "OutputEvents SelectedEvents" –///

– Count initial events –///

– Does event pass the GRL? (always true for MC) –/// – Only veto events when ALL user selectors request GRL –///

– Are the Tile and LAr calorimeters in a good state? (always true for MC) –/// – Only veto events when ALL user selectors request GOODCALO –///

– Do we have a Primary Vertex? –/// – Only veto events when ALL user selectors request PRIVTX – ///

– Wondering which triggers are available ??? –/// – Uncomment this line and get ready for a LOT of output –///

– Does event pass any trigger ? –/// – Only veto events when ALL user selectors request TRIGDEC –/// – And all trigger fail –///

– Nominal objects – /// – Calibrate objects and make all required systematic copies –///

– Object selection (e.g. good electrons, muons, jets etc.). Event selection cuts comes later –///

– Recalculate MissingET based on object selection –///

– Systematic objects – /// – Calibrate objects and make all required systematic copies –///

– Object selection (e.g. good electrons, muons, jets etc.). Event selection cuts comes later –///

– Recalculate MissingET based on object selection –///

– Scale Factor calculation –///

– Loop over all systematics and make a "top::Event" for each –///

– Make a top::Event –///

– Apply event selection –///

– Save event - we defer to eventSaver the decision to write or not –///

– For printout at the end of job, only count nominal events –///

– End loop over all calibration systematics –///

– End of tight events –///

– Loop over all Loose systematics and make a "top::Event" for each –///

– Make a top::Event –///

– Apply event selection –///

– weights for matrix-method fakes estimate from IFF tools, only for nominal –///

– Save event - we defer to eventSaver the decision to write or not –///

– For printout at the end of job, only count nominal events –///

– End loop over all loose calibration systematics –///

– End of loose events –///

– Needed for xAOD output, all systematics go into the same TTree –///

Definition at line 349 of file TopxAODAlg.cxx.

349  {
350  ++m_entry;
351  ++m_totalYieldSoFar;
352 
353  // Pile up and MC event weight - used to normalize the cut flows
354  float mcEventWeight(1.), pileupWeight(1.);
355  if (m_config->isMC())
356  mcEventWeight = m_topScaleFactors->mcEventWeight();
357 
358  if (m_config->doPileupReweighting() && !m_config->isTruthDxAOD()) {
359  top::check(m_topScaleFactors->executePileup(), "Failed to execute pileup reweighting");
360  pileupWeight = m_topScaleFactors->pileupWeight();
361  }
362 
363  if (m_config->isMC()) {
364  // if requested, pre-calculate TopPartons and PDF info
365  // Run topPartonHistory
366  if (m_config->doTopPartonHistory())
367  top::check(m_topPartonHistory->execute(), "Failed to execute m_topPartonHistory");
368 
369  // calculate PDF weights
370  if (m_config->doLHAPDF())
371  top::check(m_PDF_SF->execute(), "Failed to execute PDF SF");
372  }
373 
374  // perform any operation common to both reco and truth level
375  // currently we load the MC generator weights and PDFs if requested
376  m_eventSaver->execute();
377 
379  if (m_config->isMC()) {
380 
381  // Save, if requested, MC truth block, PDFInfo, TopPartons
382  m_eventSaver->saveTruthEvent();
383  if (m_config->doTopPartonLevel())
384  ++m_eventSavedTruth;
385 
386  // Particle level analysis, saved only for truth events passing fiducial selection
387 
388  // --------------------------------------------------
389  // If the truth loader is active, perform truth loading.
390  if (m_particleLevelLoader->active()) {
391  // --------------------------------------------------
392  // Get the top::TruthEvent for the current event
393  top::ParticleLevelEvent particleLevelEvent = m_particleLevelLoader->load();
394 
395  // event selection
396  const bool saveEventInOutputFile =
397  m_eventSelectionManager->applyParticleLevel(particleLevelEvent);
398 
399  if (saveEventInOutputFile) {
400  m_eventSaver->saveParticleLevelEvent(particleLevelEvent);
401  ++m_eventSavedParticle;
402  }
403  }
404  }
405  if (m_totalYieldSoFar == 0 && m_config->isMC() && m_config->doLHAPDF()) {
406  const xAOD::TruthEventContainer *truthEvent(nullptr);
407  top::check(evtStore()->retrieve(truthEvent, m_config->sgKeyTruthEvent()),
408  "Failed to retrieve TruthEvent container for LHAPDF");
409  top::check(truthEvent->size() == 1,
410  "TruthEvent container size != 1, not sure what to do with PDF reweighting");
411  for (const auto *tePtr : *truthEvent) {
412  for (auto &pdf : m_totalEventsPdfWeighted) {
413  if (tePtr->isAvailable<std::vector<float>>("AnalysisTop_" + pdf.first + "_Weights")) {
414  pdf.second->resize(
415  tePtr->auxdata<std::vector<float>>("AnalysisTop_" + pdf.first + "_Weights").size());
416  }
417  }
418  }
419  }
420  // on the first event, set the size of the vector of sum of LHE3 weights in case it needs to be
421  // calculated on the fly
422 
423  if (m_config->isMC()) {
424  const xAOD::EventInfo *ei(nullptr);
425  top::check(evtStore()->retrieve(ei, m_config->sgKeyEventInfo()),
426  "Failed to retrieve LHE3 weights from EventInfo");
427 
428  if (m_recalculateNominalWeightSum) {
429  if (m_totalYieldSoFar == 0)
430  ATH_MSG_INFO("Trying to recalculate nominal weights sum for TRUTH derivation");
431  const size_t nominalWeightIndex = m_config->nominalWeightIndex();
432  m_totalEventsWeighted_temp += ei->mcEventWeights().at(nominalWeightIndex);
433  m_totalEvents++;
434  }
435 
436  if (m_config->doMCGeneratorWeights()) {
437  unsigned int weightsSize = ei->mcEventWeights().size();
438  if (m_recalc_LHE3) {
439  if (m_totalYieldSoFar == 0) {
440  m_totalEventsWeighted_LHE3_temp.resize(weightsSize);
441  for (unsigned int i_LHE3 = 0; i_LHE3 < weightsSize; i_LHE3++) {
442  m_totalEventsWeighted_LHE3_temp.at(i_LHE3) = ei->mcEventWeights().at(i_LHE3);
443  }
444  m_names_LHE3.resize(weightsSize);
445 
446  const std::vector<std::string> &weight_names = m_pmg_weightTool->getWeightNames();
447  if (weight_names.size() != weightsSize) {
448  ATH_MSG_INFO(
449  "In top-xaod, while calculating mc weights sums on the fly, names from PMG tools "
450  "have different size wrt weight vector, we'll not retrieve weight names");
451  std::fill(m_names_LHE3.begin(), m_names_LHE3.end(), "?");
452  } else {
453  for (unsigned int i_wgt = 0; i_wgt < weight_names.size(); i_wgt++)
454  m_names_LHE3[i_wgt] = weight_names[i_wgt];
455  }
456  } else {
457  for (unsigned int i_LHE3 = 0; i_LHE3 < weightsSize; i_LHE3++) {
458  m_totalEventsWeighted_LHE3_temp.at(i_LHE3) =
459  m_totalEventsWeighted_LHE3_temp.at(i_LHE3) + ei->mcEventWeights().at(i_LHE3);
460  }
461  }
462  // if recalc sum of weights not requested, we still cross-check the size of the vectors
463  } else if (weightsSize != m_names_LHE3.size()) {
464  ATH_MSG_ERROR("Strange size inconsistency in the AllExecutedEvents* "
465  "sum of weights bookkeepers from the meta-data and the vector of "
466  "LHE3 weights in the EventInfo container.");
467  return StatusCode::FAILURE;
468  }
469  }
470  }
472  if (m_config->isTruthDxAOD())
473  return StatusCode::SUCCESS;
474 
479 
481  m_eventSelectionManager->countInitial(mcEventWeight, pileupWeight);
482 
485  bool passGRLVeto = m_eventCleaning->applyGRL();
486  if (!passGRLVeto) {
487  setFilterPassed(false);
488  return StatusCode::SUCCESS;
489  }
490  m_eventSelectionManager->countGRL(mcEventWeight, pileupWeight);
491 
494  bool passGoodCalo = m_eventCleaning->applyGoodCalo();
495  if (!passGoodCalo) {
496  setFilterPassed(false);
497  return StatusCode::SUCCESS;
498  }
499  m_eventSelectionManager->countGoodCalo(mcEventWeight, pileupWeight);
500 
503  bool passPriVtx = m_eventCleaning->applyPrimaryVertex();
504  if (!passPriVtx) {
505  setFilterPassed(false);
506  return StatusCode::SUCCESS;
507  }
508  m_eventSelectionManager->countPrimaryVertex(mcEventWeight, pileupWeight);
509 
512  // m_eventCleaning->printTriggers();
513 
517  bool passAnyTriggerVeto = m_eventCleaning->applyTrigger();
518  if (!passAnyTriggerVeto) {
519  setFilterPassed(false);
520  return StatusCode::SUCCESS;
521  }
522 
525  top::check(m_systObjMaker->execute(true), "Failed to execute systObjMaker");
526 
529  top::check(m_objectSelection->execute(true), "Failed to execute objectSelection");
530 
532  top::check(m_systObjMaker->recalculateMET(true), "Failed to recalculateMET with systObjMaker");
533 
536  top::check(m_systObjMaker->execute(false), "Failed to execute systObjMaker");
537 
540  top::check(m_objectSelection->execute(false), "Failed to execute objectSelection");
541 
543  top::check(m_systObjMaker->recalculateMET(false), "Failed to recalculateMET with systObjMaker");
544 
546  if (m_config->isMC())
547  top::check(m_topScaleFactors->execute(), "Failed to calculate scale factors");
548 
549  if (m_config->doTightEvents()) {
551  const xAOD::SystematicEventContainer *allSystematics =
552  m_topEventMaker->systematicEvents(m_config->sgKeyTopSystematicEvents());
553  for (const auto *currentSystematic : *allSystematics) {
554  if (!(currentSystematic->hashValue() == m_config->nominalHashValue() ||
555  m_config->doTightSysts()))
556  continue;
557 
559  top::Event topEvent = m_topEventMaker->makeTopEvent(currentSystematic);
561  const bool passAnyEventSelection =
562  m_eventSelectionManager->apply(topEvent, *currentSystematic);
563  // check if we are using actual mu for mc16d or mc16e
564  if (m_isFirst && m_config->isMC()) {
565  const int runNumber = topEvent.m_info->runNumber();
566  if (runNumber >= 300000) {
567  if ((!m_config->isAFII() && m_config->PileupActualMu_FS().size() == 0) ||
568  (m_config->isAFII() && m_config->PileupActualMu_AF().size() == 0)) {
569  ATH_MSG_WARNING(
570  "\n********************************************************************************"
571  "*******\nYou are running over MC campaigns which support actual mu rewighting, but you are not using it!"
572  "\nYou are strongly adviced to use it.\nCheck "
573  "https://twiki.cern.ch/twiki/bin/view/AtlasProtected/"
574  "TopxAODStartGuideR21#PRW_and_Lumicalc_files\n*************************************"
575  "**************************************************\n");
576  }
577  }
578  }
579  currentSystematic->auxdecor<char>(m_config->passEventSelectionDecoration()) =
580  passAnyEventSelection ? 1 : 0;
581  topEvent.m_saveEvent = passAnyEventSelection;
583  m_eventSaver->saveEvent(topEvent);
584 
586  if (passAnyEventSelection && currentSystematic->hashValue() == m_config->nominalHashValue())
587  ++m_eventSavedReco;
588  }
590  }
592 
593  if (m_config->doLooseEvents()) {
595  const xAOD::SystematicEventContainer *allSystematicsLoose =
596  m_topEventMaker->systematicEvents(m_config->sgKeyTopSystematicEventsLoose());
597  for (const auto *currentSystematic : *allSystematicsLoose) {
598  if (!(currentSystematic->hashValue() == m_config->nominalHashValue() ||
599  m_config->doLooseSysts()))
600  continue;
601 
603  top::Event topEvent = m_topEventMaker->makeTopEvent(currentSystematic);
605  const bool passAnyEventSelection =
606  m_eventSelectionManager->apply(topEvent, *currentSystematic);
607 
608  // check if we are using actual mu for mc16d or mc16e
609  if (m_isFirst && m_config->isMC()) {
610  const int runNumber = topEvent.m_info->runNumber();
611  if (runNumber >= 300000) {
612  if ((!m_config->isAFII() && m_config->PileupActualMu_FS().size() == 0) ||
613  (m_config->isAFII() && m_config->PileupActualMu_AF().size() == 0)) {
614  ATH_MSG_WARNING(
615  "\n********************************************************************************"
616  "*******\nYou are running over mc16d or mc16e sample but you are not using actual "
617  "mu reweighting!\nYou are strongly adviced to use it.\nCheck "
618  "https://twiki.cern.ch/twiki/bin/view/AtlasProtected/"
619  "TopxAODStartGuideR21#PRW_and_Lumicalc_files\n*************************************"
620  "**************************************************\n");
621  }
622  }
623  }
624  currentSystematic->auxdecor<char>(m_config->passEventSelectionDecoration()) =
625  passAnyEventSelection ? 1 : 0;
626  topEvent.m_saveEvent = passAnyEventSelection;
628  if (!m_config->isMC() && m_config->doFakesMMWeightsIFF() &&
629  currentSystematic->hashValue() == m_config->nominalHashValue()) {
630  xAOD::IParticleContainer lepton(SG::VIEW_ELEMENTS);
631  for (xAOD::Electron *t : topEvent.m_electrons)
632  lepton.push_back(static_cast<xAOD::Electron *>(t));
633  for (xAOD::Muon *t : topEvent.m_muons)
634  lepton.push_back(static_cast<xAOD::Muon *>(t));
635 
636  // TODO implement asymptotic matrix method
637  }
639  m_eventSaver->saveEvent(topEvent);
640 
642  if (passAnyEventSelection && currentSystematic->hashValue() == m_config->nominalHashValue())
643  ++m_eventSavedRecoLoose;
644  }
646  }
648 
650  m_eventSaver->saveEventToxAOD();
651  m_isFirst = false;
652 
653  return StatusCode::SUCCESS;
654 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::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

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

51 {
52  // If we didn't find any symlinks to add, just return the collection
53  // from the base class. Otherwise, return the extended collection.
54  if (!m_extendedExtraObjects.empty()) {
55  return m_extendedExtraObjects;
56  }
57  return Algorithm::extraOutputDeps();
58 }

◆ fileExecute()

StatusCode EL::AnaAlgorithm::fileExecute ( )
protectedinherited

perform the action exactly once for each file in the dataset

Ideally you don't use this, but instead rely on meta-data tools instead. However, there are enough people asking for it that I decided to implement it anyways.

\warn To use this you have to call requestFileExecute to use this.

\warn The user should not expect this to be called at any particular point in execution. If a file is split between multiple jobs this will be called in only one of these jobs, and not the others. It usually gets called before the first event in a file, but that is not guaranteed and relying on this is a bug.

\warn The execution order of beginInputFile and fileExecute is currently unspecified.

\warn fileExecute does not work with sub-file splitting in Athena, i.e. processing half the events of a file in one job the other half in another job. this should not normally happen, unless you do crazy things like run AthenaMP or explicitly select sub-file splitting in panda. in that case you are on your own.

Definition at line 342 of file AnaAlgorithm.cxx.

344  {
345  return StatusCode::SUCCESS;
346  }

◆ finalize()

StatusCode top::TopxAODAlg::finalize ( )
override

Definition at line 733 of file TopxAODAlg.cxx.

733  {
734 
735  if (m_doPerfStats > 0)
736  xAOD::PerfStats::instance().stop(); // Stop the PerfStats timer
737 
738  if (m_config->doLHAPDF())
739  top::check(m_PDF_SF->finalize(), "Failed to finalize PDF SF");
740 
741  // this fills the sum of weights
742  // if the list of input files has more than
743  // one type of dataset (ie: the first file is ttbar and the second one
744  // is W+jets), then the Fill statement below needs to come in the end of
745  // each file's processing instead ... but why would you do that?!
746  if (m_recalc_LHE3) { // in case the sum of LHE3 weight has been re-calculated with double (but we
747  // need floats in the end)
748  for (double d : m_totalEventsWeighted_LHE3_temp)
749  m_totalEventsWeighted_LHE3.push_back(d);
750  }
751  if (m_recalculateNominalWeightSum) {
752  m_totalEventsWeighted = m_totalEventsWeighted_temp;
753  }
754  m_sumWeights->Fill();
755  m_outputFile->cd();
756 
757  if (m_config->doLHAPDF()) {
758  // Save sum of weights of PDF variations
759  // Only do this if the metadata is not available
760  if (m_totalEventsInFiles != m_totalEvents && !m_pdfMetadataExists) {
761  ATH_MSG_ERROR("These are not the droids you are looking for.\n"
762  "You are running over skimmed derivations. We cannot determine "
763  "the sum of MC*PDF weights before skimming "
764  "because no PDF metadata is available in the file!");
765  return StatusCode::FAILURE;
766  }
767  // save recomputed sum weights
768  if ((!m_config->baseLHAPDF().empty() || !m_pdfMetadataExists)) {
769  for (auto &pdf_result : m_config->getLHAPDFResults())
770  *(m_totalEventsPdfWeighted[pdf_result.first]) = pdf_result.second;
771  m_sumPdfWeights->Fill();
772  }
773  } // doLHAPDF
774 
775  m_outputFile->cd();
776  m_eventSelectionManager->finalise();
777  m_eventSaver->finalize();
778  m_outputFile->cd();
779  if (m_tracker)
780  m_tracker->writeTree("AnalysisTracking");
781 
782  if (!m_config->isTruthDxAOD()) {
783  if (m_config->doTightEvents())
784  ATH_MSG_INFO("Events saved to output file nominal reconstruction tree: " << m_eventSavedReco);
785  if (m_config->doLooseEvents())
786  ATH_MSG_INFO("Events saved to output file nominal loose reconstruction tree: "
787  << m_eventSavedRecoLoose);
788  }
789  if (m_config->isMC()) {
790  ATH_MSG_INFO("Events saved to output file truth tree : " << m_eventSavedTruth);
791  if (m_particleLevelLoader->active()) {
792  ATH_MSG_INFO("Events saved to output file particle level tree : " << m_eventSavedParticle);
793  }
794  }
795  ATH_MSG_INFO("Total sum-of-weights (for normalization) : " << m_totalEventsWeighted);
796 
797  // print some xAOD performance info
798  if (m_doPerfStats == 1)
799  xAOD::IOStats::instance().stats().Print("Summary");
800  if (m_doPerfStats == 2)
801  xAOD::IOStats::instance().stats().Print();
802 
803  return StatusCode::SUCCESS;
804 }

◆ graph()

TGraph * AthHistogramming::graph ( const std::string &  graphName,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
protectedinherited

Simplify the retrieval of registered TGraphs.

Definition at line 492 of file AthHistogramming.cxx.

493 {
494  // Build a 32 bit hash out of the name
495  const hash_t graphHash = this->hash(graphName);
496 
497  // See if this entry exists in the map
498  GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
499  if ( it == m_graphMap.end() ) // It doesn't exist!
500  { // Let's see into the THistSvc if somebody else has registered the TGraph...
501 
502  // Need to copy the strings as we will massage them from here on
503  std::string graphNameCopy = graphName;
504  std::string tDirCopy = tDir;
505  std::string streamCopy = stream;
506 
507  // Massage the final string to book things
508  std::string bookingString("");
509  this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, true);
510 
511  TGraph* graphPointer(NULL);
512  if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
513  {
514  // Massage the final string to book things
515  std::string bookingString("");
516  this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, false );
517 
518  if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
519  {
520  m_msg << MSG::WARNING
521  << "Problem retrieving the TGraph with name (including pre- and post-fixes) "
522  << m_histNamePrefix + graphNameCopy + m_histNamePostfix
523  << " or with name " << graphNameCopy
524  << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
525  << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
526  return NULL;
527  }
528  // If we get to here, we actually found the TGraph in the THistSvc.
529  // So let's add it to the local cache map and return its pointer
530  m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
531  return graphPointer;
532  }
533  // If we get to here, we actually found the TGraph in the THistSvc.
534  // So let's add it to the local cache map and return its pointer
535  m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
536  return graphPointer;
537  }
538 
539 
540  // Return the pointer to the TGraph that we got from the local cache map
541  return it->second;
542 }

◆ handle()

void EL::AnaAlgorithm::handle ( const Incident &  inc)
inherited

receive the given incident

Guarantee
basic
Failures
incident handling errors

Definition at line 520 of file AnaAlgorithm.cxx.

522  {
523  if (inc.type() == IncidentType::BeginInputFile)
524  {
525  if (m_hasBeginInputFile)
526  ANA_CHECK_THROW (beginInputFile ());
527  if (m_hasFileExecute)
528  ANA_CHECK_THROW (fileExecute ());
529  } else if (inc.type() == IncidentType::EndInputFile)
530  {
531  if (m_hasEndInputFile)
532  ANA_CHECK_THROW (endInputFile ());
533  } else
534  {
535  ATH_MSG_WARNING( "Unknown incident type received: " << inc.type() );
536  }
537  }

◆ hash()

AthHistogramming::hash_t AthHistogramming::hash ( const std::string &  histName) const
inlineprivateinherited

Method to calculate a 32-bit hash from a string.

Definition at line 428 of file AthHistogramming.h.

429 {
430  const uint64_t hash64 = CxxUtils::crc64( histName );
431  return (hash_t)(hash64 & 0xFFFFFFFF);
432 }

◆ hist()

TH1 * AthHistogramming::hist ( const std::string &  histName,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
protectedinherited

Simplify the retrieval of registered histograms of any type.

Definition at line 198 of file AthHistogramming.cxx.

199 {
200  // Build a 32 bit hash out of the name
201  const hash_t histHash = this->hash(histName);
202 
203  // See if this entry exists in the map
204  HistMap_t::const_iterator it = m_histMap.find( histHash );
205  if ( it == m_histMap.end() ) // It doesn't exist!
206  { // Let's see into the THistSvc if somebody else has registered the histogram...
207 
208  // Need to copy the strings as we will massage them from here on
209  std::string histNameCopy = histName;
210  std::string tDirCopy = tDir;
211  std::string streamCopy = stream;
212 
213  // Massage the final string to book things
214  std::string bookingString("");
215  this->buildBookingString( bookingString, histNameCopy, tDirCopy, streamCopy ,false);
216 
217  TH1* histPointer(NULL);
218  if ( !((histSvc()->getHist(bookingString, histPointer)).isSuccess()) )
219  {
220  // Massage the final string to book things
221  std::string bookingString("");
222  this->buildBookingString( bookingString, histNameCopy, tDirCopy, streamCopy, true );
223 
224  if ( !((histSvc()->getHist(bookingString, histPointer)).isSuccess()) )
225  {
226  m_msg << MSG::WARNING
227  << "Problem retrieving the histogram with name (including pre- and post-fixes) "
228  << m_histNamePrefix + histNameCopy + m_histNamePostfix
229  << " or with name " << histNameCopy
230  << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
231  << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
232  return NULL;
233  }
234  // If we get to here, we actually found the histogram in the THistSvc.
235  // So let's add it to the local cache map and return its pointer
236  m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
237  return histPointer;
238  }
239  // If we get to here, we actually found the histogram in the THistSvc.
240  // So let's add it to the local cache map and return its pointer
241  m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
242  return histPointer;
243  }
244 
245 
246  // Return the pointer to the histogram that we got from the local cache map
247  return it->second;
248 }

◆ hist2d()

TH2 * AthHistogramming::hist2d ( const std::string &  histName,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the retrieval of registered 2-d histograms.

Definition at line 369 of file AthHistogramming.h.

370 {
371  // Get the TH1 pointer
372  TH1* th1Pointer = this->hist(histName, tDir, stream);
373  if ( !th1Pointer )
374  {
375  m_msg << MSG::ERROR
376  << "Cannot get a 2-d histogram with name " << histName
377  << "... will probably seg-fault!" << endmsg;
378  return NULL;
379  }
380  // If the TH1 pointer is valid, simply return the dynamic_cast
381  return dynamic_cast<TH2*>( th1Pointer );
382 }

◆ hist3d()

TH3 * AthHistogramming::hist3d ( const std::string &  histName,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
inlineprotectedinherited

Simplify the retrieval of registered 3-d histograms.

Definition at line 386 of file AthHistogramming.h.

387 {
388  // Get the TH1 pointer
389  TH1* th1Pointer = this->hist(histName, tDir, stream);
390  if ( !th1Pointer )
391  {
392  m_msg << MSG::ERROR
393  << "Cannot get a 3-d histogram with name " << histName
394  << "... will probably seg-fault!" << endmsg;
395  return NULL;
396  }
397  // If the TH1 pointer is valid, simply return the dynamic_cast
398  return dynamic_cast<TH3*>( th1Pointer );
399 }

◆ histSvc()

const ServiceHandle< ITHistSvc > & AthHistogramAlgorithm::histSvc ( ) const
inlineinherited

The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc.

Definition at line 113 of file AthHistogramAlgorithm.h.

114 {
115  return m_histSvc;
116 }

◆ initialize()

StatusCode top::TopxAODAlg::initialize ( )
override

Definition at line 39 of file TopxAODAlg.cxx.

39  {
40  ANA_CHECK(requestBeginInputFile());
41  ANA_CHECK(requestEndInputFile());
42 
43  m_settings.reset(top::ConfigurationSettings::get(true));
44  m_settings->loadFromFile(m_configPath);
45 
46  m_config.reset(new top::TopConfig());
47  // read FileMetadata for our TopConfig class
48  const xAOD::FileMetaData *fmd = nullptr;
49  if (inputMetaStore()->contains<xAOD::FileMetaData>("FileMetaData")) {
50  ANA_CHECK(inputMetaStore()->retrieve(fmd, "FileMetaData"));
51  top::readMetaData(fmd, m_config);
52  } else {
53  ATH_MSG_ERROR("Could not retrieve FileMetadata. We thus cannot reliably pre-determine DSID, "
54  "whether this is MC, and simulation type.");
55  if (m_crashOnMetaDataFail) {
56  return StatusCode::FAILURE;
57  }
58  }
59 
60  m_config->setConfigSettings(m_settings.get());
61 
62  if (m_config->isMC() && !m_config->isTruthDxAOD()) {
63  // at this point we should also configure showering algorithm for FTAG MC/MC SFs
64  m_config->setMapIndex(m_FTAGMapIndex);
65  m_config->setShoweringAlgorithm(m_showeringAlgo);
66  }
67 
68  // trigger configuration
69  m_config->SetYear(m_year);
70  m_config->SetTriggersToYear(m_config->isMC());
71 
72  m_firstEvent = m_config->numberOfEventsToSkip();
73 
74  // analysis tracking
75  {
76  bool useTracking = true;
77  m_settings->retrieve("WriteTrackingData", useTracking);
78  if (useTracking) {
79  m_tracker.reset(new top::AnalysisTrackingHelper());
80  }
81  }
82 
83  // I/O Performance stats?
84  // Summary gives a summary
85  // Full gives detailed info on each collection in the file
86  if (m_settings->value("PerfStats") == "Summary")
87  m_doPerfStats = 1;
88 
89  if (m_settings->value("PerfStats") == "Full")
90  m_doPerfStats = 2;
91 
92  // Setup all asg::AsgTools
93  m_topTools = new top::TopToolStore("top::TopToolStore");
94  top::check(m_topTools->setProperty("config", m_config), "Failed to setProperty of topTools");
95  top::check(m_topTools->initialize(), "Failed to initialize topTools");
96 
97  if (m_config->isMC() && !m_pmg_weightTool.retrieve()) {
98  ATH_MSG_ERROR("Cannot retrieve PMGTruthWeightTool");
99  return StatusCode::FAILURE;
100  }
101 
102  // EventCleaningSelection
103  // Decorates EventInfo with GRL decision (always true for MC)
104  // Evaluates all triggers requested from all selections
105  // Trigger matching
106  // Check for good Tile and LAr calorimeters
107  // Let's check this passes before running expensive calibrations
108  // If ANY selection does not request a trigger, we won't veto any events
109  m_eventCleaning = new top::EventCleaningSelection("top::EventCleaningSelection");
110  top::check(m_eventCleaning->setProperty("config", m_config),
111  "Failed to setProperty of triggerGRLSelection");
112  top::check(m_eventCleaning->initialize(), "Failed to initialize triggerGRLSelection");
113  m_eventCleaning->setEventSelections(m_settings->selections());
114 
115  // Systematic object collection making
116  m_systObjMaker = new top::ObjectCollectionMaker("top::ObjectCollectionMaker");
117  top::check(m_systObjMaker->setProperty("config", m_config),
118  "Failed to setProperty of systObjMaker");
119  if (!m_config->isTruthDxAOD())
120  top::check(m_systObjMaker->initialize(), "Failed to initialize systObjMaker");
121 
122  // object selection
123  m_objectSelection = top::loadObjectSelection(m_config);
124  m_objectSelection->print(msg(MSG::Level::INFO));
125 
126  // setup event-level cuts
127  m_outputFile = wk()->getOutputFile(m_streamName);
128  m_eventSelectionManager.reset(new top::EventSelectionManager(
129  m_settings->selections(), m_outputFile, m_libraryNames, m_config));
130 
131  // The loader tool for top::ParticleLevelEvent objects
132  m_particleLevelLoader = new top::ParticleLevelLoader(m_config);
133 
134  // The loader tool for Upgrade objects
135  // top::UpgradeObjectLoader upgradeLoader(m_config);
136 
137  // Fix the configuration - it now knows about:
138  // * all objects collections to work with
139  // * all systematic variations
140  m_config->fixConfiguration();
141 
142  // OK let's printout the TopConfig
143  ATH_MSG_INFO(*m_config << "\n");
144 
145  if (m_tracker)
146  m_tracker->setTopConfig(m_config);
147 
148  // Top parton history for MC events
149  // This is quite ugly and simple, it will be harmonized with in the future
150  // along with all other factory methods (it's not a factory method right now)
151  if (m_settings->value("TopPartonHistory") == "ttbar") {
152  m_topPartonHistory = new top::CalcTtbarPartonHistory("top::CalcTtbarPartonHistory");
153  top::check(m_topPartonHistory->setProperty("config", m_config),
154  "Failed to setProperty of top::CalcTtbarPartonHistory");
155  } else {
156  ATH_MSG_ERROR("Only ttbar TopPartonHistory is supported for the moment!");
157  return StatusCode::FAILURE;
158  }
159 
160  // LHAPDF SF calculation
161  m_PDF_SF = new top::PDFScaleFactorCalculator("top::PDFScaleFactorCalculator");
162  if (m_config->doLHAPDF()) {
163  top::check(m_PDF_SF->setProperty("config", m_config),
164  "Failed to set config for PDF SF Calculator");
165  top::check(m_PDF_SF->initialize(), "Failed to initialize PDF SF calculator");
166  }
167 
168  // make top::Event objects
169  m_topEventMaker = new top::TopEventMaker("top::TopEventMaker");
170  top::check(m_topEventMaker->setProperty("config", m_config),
171  "Failed to setProperty of top::TopEventMaker");
172  top::check(m_topEventMaker->initialize(), "Failed to initialize top::TopEventMaker");
173  // Debug messages?
174  // m_topEventMaker.msg().setLevel(MSG::DEBUG);
175 
176  m_topScaleFactors = new top::ScaleFactorCalculator("top::ScaleFactorCalculator");
177  top::check(m_topScaleFactors->setProperty("config", m_config),
178  "Failed to setProperty of top::ScaleFactorCalculator");
179  top::check(m_topScaleFactors->initialize(), "Failed to initialize top::ScaleFactorCalculator");
180 
181  // TODO asymptotic matrix method
182 
183  // A list of extra branches that we create in this code and want to save
184  std::vector<std::string> extraBranches;
185  m_eventSelectionManager->addExtraBranches(extraBranches);
186 
187  // output file format (and any user-modified code)
188  m_eventSaver.reset(top::loadEventSaver(m_config));
189  m_eventSaver->initialize(m_config, m_outputFile, extraBranches);
190 
191  // save sum of weights before derivations
192  m_outputFile->cd();
193  m_sumWeights = new TTree("sumWeights", "");
194  m_boostedTaggersSFSysNames = m_config->boostedTaggersSFSysNames();
195  m_DSID = m_config->getDSID();
196  m_isAFII = m_config->isAFII();
197 
198  m_sumWeights->Branch("dsid", &m_DSID);
199  m_sumWeights->Branch("isAFII", &m_isAFII);
200  m_sumWeights->Branch("totalEventsWeighted", &m_totalEventsWeighted);
201  if (m_config->doMCGeneratorWeights()) { // the main problem is that we don't have the list of
202  // names a priori
203  m_sumWeights->Branch("totalEventsWeighted_mc_generator_weights", &m_totalEventsWeighted_LHE3);
204  m_sumWeights->Branch("names_mc_generator_weights", &m_names_LHE3);
205  }
206  m_sumWeights->Branch("totalEvents", &m_totalEvents, "totalEvents/l");
207 
208  for (auto &it : m_boostedTaggersSFSysNames) {
209  m_sumWeights->Branch(("sysNames_" + it.first).c_str(), &it.second);
210  }
211 
212  unsigned int dsidPdf = m_config->getDSID();
213  if (m_config->doLHAPDF()) {
214  m_sumPdfWeights = new TTree("PDFsumWeights", "");
215  m_sumPdfWeights->Branch("dsid", &dsidPdf);
216  for (const auto &pdf_set : m_config->LHAPDFSets()) {
217  m_totalEventsPdfWeighted[pdf_set] = std::make_unique<std::vector<float>>();
218  m_sumPdfWeights->Branch(pdf_set.c_str(), m_totalEventsPdfWeighted[pdf_set].get());
219  }
220  }
221 
222  // the analysis loop
223  ATH_MSG_INFO("Starting event loop\n");
224 
225  // Performance stats
226  if (m_doPerfStats > 0)
227  xAOD::PerfStats::instance().start(); // start Perfstats timer
228 
229  return StatusCode::SUCCESS;
230 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::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.

◆ inputMetaStore() [1/2]

AnaAlgorithm::MetaStorePtr_t EL::AnaAlgorithm::inputMetaStore ( )
inherited

Definition at line 81 of file AnaAlgorithm.cxx.

82  {
83 #ifdef XAOD_STANDALONE
84  return &m_inputMetaStore;
85 #else
86  return m_inputMetaStore;
87 #endif // XAOD_STANDALONE
88  }

◆ inputMetaStore() [2/2]

AnaAlgorithm::ConstMetaStorePtr_t EL::AnaAlgorithm::inputMetaStore ( ) const
inherited

Accessor for the input metadata store

Definition at line 72 of file AnaAlgorithm.cxx.

73  {
74 #ifdef XAOD_STANDALONE
75  return &m_inputMetaStore;
76 #else
77  return m_inputMetaStore;
78 #endif // XAOD_STANDALONE
79  }

◆ msg() [1/2]

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

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::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< Algorithm >::msgLvl ( const MSG::Level  lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30  {
31  return this->msgLevel(lvl);
32  }

◆ myReplace()

void AthHistogramming::myReplace ( std::string &  str,
const std::string &  oldStr,
const std::string &  newStr 
)
privateinherited

Helper method to replace sub-string.

Definition at line 590 of file AthHistogramming.cxx.

593 {
594  size_t pos = 0;
595  while((pos = str.find(oldStr, pos)) != std::string::npos)
596  {
597  str.replace(pos, oldStr.length(), newStr);
598  pos += newStr.length();
599  }
600 }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::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.

◆ outputMetaStore() [1/2]

AnaAlgorithm::MetaStorePtr_t EL::AnaAlgorithm::outputMetaStore ( )
inherited

Definition at line 101 of file AnaAlgorithm.cxx.

102  {
103 #ifdef XAOD_STANDALONE
104  return &m_outputMetaStore;
105 #else
106  return m_outputMetaStore;
107 #endif // XAOD_STANDALONE
108  }

◆ outputMetaStore() [2/2]

AnaAlgorithm::ConstMetaStorePtr_t EL::AnaAlgorithm::outputMetaStore ( ) const
inherited

Accessor for the output metadata store

Definition at line 92 of file AnaAlgorithm.cxx.

93  {
94 #ifdef XAOD_STANDALONE
95  return &m_outputMetaStore;
96 #else
97  return m_outputMetaStore;
98 #endif // XAOD_STANDALONE
99  }

◆ print()

void EL::AnaAlgorithm::print ( ) const
protectedvirtualinherited

print the state of the algorithm

This is mostly to allow algorithms to add a little debugging information if they feel like it.

Definition at line 336 of file AnaAlgorithm.cxx.

338  {}

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

◆ requestBeginInputFile()

StatusCode EL::AnaAlgorithm::requestBeginInputFile ( )
inherited

register this algorithm to have an implementation of beginInputFile

Guarantee
strong
Failures
beginInputFile not supported

Definition at line 270 of file AnaAlgorithm.cxx.

272  {
273  m_hasBeginInputFile = true;
274 
275 #ifndef XAOD_STANDALONE
276  // Connect to the IncidentSvc:
277  ServiceHandle< IIncidentSvc > incSvc( "IncidentSvc", name() );
278  ATH_CHECK( incSvc.retrieve() );
279 
280  // Set up the right callback, but ensure we don't double-register
281  // if we are called twice
282  incSvc->removeListener( this, IncidentType::BeginInputFile );
283  incSvc->addListener( this, IncidentType::BeginInputFile, 0, true );
284 #endif
285 
286  return StatusCode::SUCCESS;
287  }

◆ requestEndInputFile()

StatusCode EL::AnaAlgorithm::requestEndInputFile ( )
inherited

register this algorithm to have an implementation of endInputFile

Guarantee
strong
Failures
endInputFile not supported

Definition at line 291 of file AnaAlgorithm.cxx.

293  {
294  m_hasEndInputFile = true;
295 
296 #ifndef XAOD_STANDALONE
297  // Connect to the IncidentSvc:
298  ServiceHandle< IIncidentSvc > incSvc( "IncidentSvc", name() );
299  ATH_CHECK( incSvc.retrieve() );
300 
301  // Set up the right callback, but ensure we don't double-register
302  // if we are called twice
303  incSvc->removeListener( this, IncidentType::EndInputFile );
304  incSvc->addListener( this, IncidentType::EndInputFile, 0, true );
305 #endif
306 
307  return StatusCode::SUCCESS;
308  }

◆ requestFileExecute()

StatusCode EL::AnaAlgorithm::requestFileExecute ( )
inherited

register this algorithm to have an implementation of fileexecute

Guarantee
strong
Failures
fileExecute not supported

Definition at line 249 of file AnaAlgorithm.cxx.

251  {
252  m_hasFileExecute = true;
253 
254 #ifndef XAOD_STANDALONE
255  // Connect to the IncidentSvc:
256  ServiceHandle< IIncidentSvc > incSvc( "IncidentSvc", name() );
257  ATH_CHECK( incSvc.retrieve() );
258 
259  // Set up the right callback, but ensure we don't double-register
260  // if we are called twice
261  incSvc->removeListener( this, IncidentType::BeginInputFile );
262  incSvc->addListener( this, IncidentType::BeginInputFile, 0, true );
263 #endif
264 
265  return StatusCode::SUCCESS;
266  }

◆ sysInitialize()

StatusCode AthHistogramAlgorithm::sysInitialize ( )
virtualinherited

Initialization method invoked by the framework.

This method is responsible for any bookkeeping of initialization required by the framework itself. It will in turn invoke the initialize() method of the derived algorithm, and of any sub-algorithms which it creates.

Reimplemented from AthAlgorithm.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 75 of file AthHistogramAlgorithm.cxx.

76 {
77  // ---- stolen from GaudiKernel/Algorithm::sysInitialize -------
78  // Bypass the initialization if the algorithm
79  // has already been initialized.
80  if ( Gaudi::StateMachine::INITIALIZED <= FSMState() ) return StatusCode::SUCCESS;
81 
82  // Set the Algorithm's properties
83  bindPropertiesTo( serviceLocator()->getOptsSvc() );
84 
85  // Bypass the initialization if the algorithm is disabled.
86  // Need to do this after bindPropertiesTo.
87  if ( !isEnabled( ) ) return StatusCode::SUCCESS;
88 
89  // ---- stolen from GaudiKernel/Algorithm::sysInitialize ------- END ---
90 
91 
92  // Get the THistSvc
93  ATH_CHECK ( histSvc().retrieve() );
94 
95  // Configure the underlying AthHistogramming helper
96  ATH_CHECK ( AthHistogramming::configAthHistogramming( m_histSvc,
97  m_prefix, m_rootDir,
98  m_histNamePrefix, m_histNamePostfix,
99  m_histTitlePrefix, m_histTitlePostfix ) );
100 
101  // Print some setup information into the log file
102  ATH_MSG_DEBUG ("Initializing " << name() << "...");
103  ATH_MSG_DEBUG (" using THistService = " << m_histSvc );
104  ATH_MSG_DEBUG (" using RootStreamName = " << m_prefix );
105  ATH_MSG_DEBUG (" using RootDirName = " << m_rootDir );
106  ATH_MSG_DEBUG (" using HistNamePrefix = " << m_histNamePrefix );
107  ATH_MSG_DEBUG (" using HistNamePostfix = " << m_histNamePostfix );
108  ATH_MSG_DEBUG (" using HistTitlePrefix = " << m_histTitlePrefix );
109  ATH_MSG_DEBUG (" using HistTitlePostfix = " << m_histTitlePostfix );
110 
111 
112  // re-direct to base class...
113  return AthAlgorithm::sysInitialize();
114 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::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.

◆ tree()

TTree * AthHistogramming::tree ( const std::string &  treeName,
const std::string &  tDir = "",
const std::string &  stream = "" 
)
protectedinherited

Simplify the retrieval of registered TTrees.

Definition at line 378 of file AthHistogramming.cxx.

379 {
380  // Build a 32 bit hash out of the name
381  const hash_t treeHash = this->hash(treeName);
382 
383  // See if this entry exists in the map
384  TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
385  if ( it == m_treeMap.end() ) // It doesn't exist!
386  { // Let's see into the THistSvc if somebody else has registered the TTree...
387 
388  // Need to copy the strings as we will massage them from here on
389  std::string treeNameCopy = treeName;
390  std::string tDirCopy = tDir;
391  std::string streamCopy = stream;
392 
393  // Massage the final string to book things
394  std::string bookingString("");
395  this->buildBookingString( bookingString, treeNameCopy, tDirCopy, streamCopy );
396 
397  TTree* treePointer(NULL);
398  if ( !((histSvc()->getTree(bookingString, treePointer)).isSuccess()) )
399  {
400  m_msg << MSG::WARNING
401  << "Problem retrieving the TTree with name " << treeNameCopy
402  << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
403  << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
404  return NULL;
405  }
406  // If we get to here, we actually found the TTree in the THistSvc.
407  // So let's add it to the local cache map and return its pointer
408  m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
409  return treePointer;
410  }
411 
412  // Return the pointer to the TTree that we got from the local cache map
413  return it->second;
414 }

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::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_boostedTaggersSFSysNames

std::unordered_map<std::string, std::vector<std::string> > top::TopxAODAlg::m_boostedTaggersSFSysNames
private

Definition at line 93 of file TopxAODAlg.h.

◆ m_config

std::shared_ptr<TopConfig> top::TopxAODAlg::m_config {nullptr}
private

Definition at line 50 of file TopxAODAlg.h.

◆ m_configPath

std::string top::TopxAODAlg::m_configPath {""}
private

properties

Definition at line 96 of file TopxAODAlg.h.

◆ m_crashOnMetaDataFail

bool top::TopxAODAlg::m_crashOnMetaDataFail {true}
private

Definition at line 99 of file TopxAODAlg.h.

◆ m_derivationStream

std::string top::TopxAODAlg::m_derivationStream {""}
private

Definition at line 98 of file TopxAODAlg.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_doPerfStats

unsigned int top::TopxAODAlg::m_doPerfStats = 0u
private

Definition at line 65 of file TopxAODAlg.h.

◆ m_DSID

unsigned int top::TopxAODAlg::m_DSID
private

Definition at line 91 of file TopxAODAlg.h.

◆ m_effMap

EffMap_t AthHistogramming::m_effMap
privateinherited

The map of histogram names to their pointers.

Definition at line 208 of file AthHistogramming.h.

◆ m_entry

unsigned int top::TopxAODAlg::m_entry = 0u
private

Definition at line 82 of file TopxAODAlg.h.

◆ m_eventCleaning

ToolHandle<top::EventCleaningSelection> top::TopxAODAlg::m_eventCleaning
private

Definition at line 52 of file TopxAODAlg.h.

◆ m_eventSavedParticle

unsigned int top::TopxAODAlg::m_eventSavedParticle = 0u
private

Definition at line 88 of file TopxAODAlg.h.

◆ m_eventSavedReco

unsigned int top::TopxAODAlg::m_eventSavedReco = 0u
private

Definition at line 85 of file TopxAODAlg.h.

◆ m_eventSavedRecoLoose

unsigned int top::TopxAODAlg::m_eventSavedRecoLoose = 0u
private

Definition at line 86 of file TopxAODAlg.h.

◆ m_eventSavedTruth

unsigned int top::TopxAODAlg::m_eventSavedTruth = 0u
private

Definition at line 87 of file TopxAODAlg.h.

◆ m_eventSaver

std::unique_ptr<top::EventSaverBase> top::TopxAODAlg::m_eventSaver
private

Definition at line 62 of file TopxAODAlg.h.

◆ m_eventSelectionManager

std::unique_ptr<top::EventSelectionManager> top::TopxAODAlg::m_eventSelectionManager
private

Definition at line 55 of file TopxAODAlg.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_firstEvent

unsigned int top::TopxAODAlg::m_firstEvent = 0u
private

Definition at line 83 of file TopxAODAlg.h.

◆ m_FTAGMapIndex

unsigned int top::TopxAODAlg::m_FTAGMapIndex
private

Definition at line 101 of file TopxAODAlg.h.

◆ m_graphMap

GraphMap_t AthHistogramming::m_graphMap
privateinherited

The map of TGraph names to their pointers.

Definition at line 222 of file AthHistogramming.h.

◆ m_hasBeginInputFile

bool EL::AnaAlgorithm::m_hasBeginInputFile {false}
privateinherited

the value of hasBeginInputFile

Definition at line 590 of file AnaAlgorithm.h.

◆ m_hasEndInputFile

bool EL::AnaAlgorithm::m_hasEndInputFile {false}
privateinherited

the value of hasEndInputFile

Definition at line 594 of file AnaAlgorithm.h.

◆ m_hasFileExecute

bool EL::AnaAlgorithm::m_hasFileExecute {false}
privateinherited

the value of hasFileExecute

Definition at line 586 of file AnaAlgorithm.h.

◆ m_histMap

HistMap_t AthHistogramming::m_histMap
privateinherited

The map of histogram names to their pointers.

Definition at line 201 of file AthHistogramming.h.

◆ m_histNamePostfix

std::string AthHistogramAlgorithm::m_histNamePostfix
privateinherited

The postfix for the histogram THx name.

Definition at line 97 of file AthHistogramAlgorithm.h.

◆ m_histNamePrefix

std::string AthHistogramAlgorithm::m_histNamePrefix
privateinherited

The prefix for the histogram THx name.

Definition at line 94 of file AthHistogramAlgorithm.h.

◆ m_histSvc

ServiceHandle<ITHistSvc> AthHistogramAlgorithm::m_histSvc
privateinherited

Default constructor: AthHistogramAlgorithm();.

a handle on the Hist/TTree registration service

Definition at line 83 of file AthHistogramAlgorithm.h.

◆ m_histTitlePostfix

std::string AthHistogramAlgorithm::m_histTitlePostfix
privateinherited

The postfix for the histogram THx title.

Definition at line 103 of file AthHistogramAlgorithm.h.

◆ m_histTitlePrefix

std::string AthHistogramAlgorithm::m_histTitlePrefix
privateinherited

The prefix for the histogram THx title.

Definition at line 100 of file AthHistogramAlgorithm.h.

◆ m_inputMetaStore

MetaStore_t EL::AnaAlgorithm::m_inputMetaStore
privateinherited

Object accessing the input metadata store.

Definition at line 551 of file AnaAlgorithm.h.

◆ m_isAFII

int top::TopxAODAlg::m_isAFII
private

Definition at line 92 of file TopxAODAlg.h.

◆ m_isFirst

bool top::TopxAODAlg::m_isFirst = true
private

Definition at line 89 of file TopxAODAlg.h.

◆ m_libraryNames

std::string top::TopxAODAlg::m_libraryNames {""}
private

Definition at line 97 of file TopxAODAlg.h.

◆ m_msg

MsgStream AthHistogramming::m_msg
privateinherited

Cached Message Stream.

Definition at line 248 of file AthHistogramming.h.

◆ m_name

std::string AthHistogramming::m_name
privateinherited

Instance name.

Definition at line 245 of file AthHistogramming.h.

◆ m_names_LHE3

std::vector<std::string> top::TopxAODAlg::m_names_LHE3
private

Definition at line 77 of file TopxAODAlg.h.

◆ m_objectSelection

ToolHandle<top::TopObjectSelection> top::TopxAODAlg::m_objectSelection
private

Definition at line 54 of file TopxAODAlg.h.

◆ m_outputFile

TFile* top::TopxAODAlg::m_outputFile = nullptr
private

Definition at line 66 of file TopxAODAlg.h.

◆ m_outputMetaStore

MetaStore_t EL::AnaAlgorithm::m_outputMetaStore
privateinherited

Object accessing the output metadata store.

Definition at line 555 of file AnaAlgorithm.h.

◆ m_particleLevelLoader

ToolHandle<top::ParticleLevelLoader> top::TopxAODAlg::m_particleLevelLoader
private

Definition at line 57 of file TopxAODAlg.h.

◆ m_PDF_SF

ToolHandle<top::PDFScaleFactorCalculator> top::TopxAODAlg::m_PDF_SF
private

Definition at line 59 of file TopxAODAlg.h.

◆ m_pdfMetadataExists

bool top::TopxAODAlg::m_pdfMetadataExists = false
private

Definition at line 69 of file TopxAODAlg.h.

◆ m_pmg_weightTool

ToolHandle<PMGTools::IPMGTruthWeightTool> top::TopxAODAlg::m_pmg_weightTool {"PMGTruthWeightTool"}
private

Definition at line 63 of file TopxAODAlg.h.

◆ m_prefix

std::string AthHistogramAlgorithm::m_prefix
privateinherited

Name of the ROOT output stream (file)

Definition at line 88 of file AthHistogramAlgorithm.h.

◆ m_recalc_LHE3

bool top::TopxAODAlg::m_recalc_LHE3 = false
private

Definition at line 75 of file TopxAODAlg.h.

◆ m_recalculateNominalWeightSum

bool top::TopxAODAlg::m_recalculateNominalWeightSum = false
private

Definition at line 76 of file TopxAODAlg.h.

◆ m_rootDir

std::string AthHistogramAlgorithm::m_rootDir
privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

◆ m_settings

std::shared_ptr<ConfigurationSettings> top::TopxAODAlg::m_settings {nullptr}
private

Definition at line 49 of file TopxAODAlg.h.

◆ m_showeringAlgo

SampleXsection::showering top::TopxAODAlg::m_showeringAlgo
private

Definition at line 102 of file TopxAODAlg.h.

◆ m_streamName

std::string top::TopxAODAlg::m_streamName {"output"}
private

Definition at line 100 of file TopxAODAlg.h.

◆ m_sumPdfWeights

TTree* top::TopxAODAlg::m_sumPdfWeights = nullptr
private

Definition at line 67 of file TopxAODAlg.h.

◆ m_sumWeights

TTree* top::TopxAODAlg::m_sumWeights = nullptr
private

Definition at line 74 of file TopxAODAlg.h.

◆ m_systObjMaker

ToolHandle<top::ObjectCollectionMaker> top::TopxAODAlg::m_systObjMaker
private

Definition at line 53 of file TopxAODAlg.h.

◆ m_topEventMaker

ToolHandle<top::TopEventMaker> top::TopxAODAlg::m_topEventMaker
private

Definition at line 60 of file TopxAODAlg.h.

◆ m_topPartonHistory

ToolHandle<top::CalcTtbarPartonHistory> top::TopxAODAlg::m_topPartonHistory
private

Definition at line 58 of file TopxAODAlg.h.

◆ m_topScaleFactors

ToolHandle<top::ScaleFactorCalculator> top::TopxAODAlg::m_topScaleFactors
private

Definition at line 61 of file TopxAODAlg.h.

◆ m_topTools

ToolHandle<top::TopToolStore> top::TopxAODAlg::m_topTools
private

Definition at line 51 of file TopxAODAlg.h.

◆ m_totalEvents

ULong64_t top::TopxAODAlg::m_totalEvents = 0
private

Definition at line 72 of file TopxAODAlg.h.

◆ m_totalEventsInFiles

ULong64_t top::TopxAODAlg::m_totalEventsInFiles = 0
private

Definition at line 73 of file TopxAODAlg.h.

◆ m_totalEventsPdfWeighted

std::unordered_map<std::string, std::unique_ptr<std::vector<float> > > top::TopxAODAlg::m_totalEventsPdfWeighted
private

Definition at line 68 of file TopxAODAlg.h.

◆ m_totalEventsWeighted

double top::TopxAODAlg::m_totalEventsWeighted = 0
private

Definition at line 70 of file TopxAODAlg.h.

◆ m_totalEventsWeighted_LHE3

std::vector<double> top::TopxAODAlg::m_totalEventsWeighted_LHE3
private

Definition at line 78 of file TopxAODAlg.h.

◆ m_totalEventsWeighted_LHE3_temp

std::vector<double> top::TopxAODAlg::m_totalEventsWeighted_LHE3_temp
private

Definition at line 80 of file TopxAODAlg.h.

◆ m_totalEventsWeighted_temp

double top::TopxAODAlg::m_totalEventsWeighted_temp = 0
private

Definition at line 71 of file TopxAODAlg.h.

◆ m_totalYieldSoFar

unsigned int top::TopxAODAlg::m_totalYieldSoFar = 0u
private

Definition at line 84 of file TopxAODAlg.h.

◆ m_tracker

std::unique_ptr<top::AnalysisTrackingHelper> top::TopxAODAlg::m_tracker
private

Definition at line 56 of file TopxAODAlg.h.

◆ m_treeMap

TreeMap_t AthHistogramming::m_treeMap
privateinherited

The map of TTree names to their pointers.

Definition at line 215 of file AthHistogramming.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_year

std::string top::TopxAODAlg::m_year
private

Definition at line 103 of file TopxAODAlg.h.

The documentation for this class was generated from the following files:
AthHistogramming::m_rootDir
std::string m_rootDir
Name of the ROOT directory.
Definition: AthHistogramming.h:229
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
AthHistogramming::histSvc
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
Definition: AthHistogramming.h:421
AthHistogramAlgorithm::histSvc
const ServiceHandle< ITHistSvc > & histSvc() const
The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a...
Definition: AthHistogramAlgorithm.h:113
top::TopxAODAlg::m_eventSavedRecoLoose
unsigned int m_eventSavedRecoLoose
Definition: TopxAODAlg.h:86
EL::AnaAlgorithm::m_hasBeginInputFile
bool m_hasBeginInputFile
the value of hasBeginInputFile
Definition: AnaAlgorithm.h:590
EL::AnaAlgorithm::endInputFile
virtual ::StatusCode endInputFile()
perform the action for the end of an input file
Definition: AnaAlgorithm.cxx:359
top::Event::m_info
const xAOD::EventInfo * m_info
Pointer to the event info - nullptr if not loaded, but that's probably a bad idea.
Definition: Event.h:90
AthHistogramming::myReplace
void myReplace(std::string &str, const std::string &oldStr, const std::string &newStr)
Helper method to replace sub-string.
Definition: AthHistogramming.cxx:590
AthHistogramming::m_msg
MsgStream m_msg
Cached Message Stream.
Definition: AthHistogramming.h:248
python.SystemOfUnits.s
int s
Definition: SystemOfUnits.py:131
AthHistogramming::m_histMap
HistMap_t m_histMap
The map of histogram names to their pointers.
Definition: AthHistogramming.h:201
python.PerfMonSerializer.p
def p
Definition: PerfMonSerializer.py:743
AddEmptyComponent.histName
string histName
Definition: AddEmptyComponent.py:64
top::TopxAODAlg::m_particleLevelLoader
ToolHandle< top::ParticleLevelLoader > m_particleLevelLoader
Definition: TopxAODAlg.h:57
AthHistogramAlgorithm::m_histSvc
ServiceHandle< ITHistSvc > m_histSvc
Default constructor: AthHistogramAlgorithm();.
Definition: AthHistogramAlgorithm.h:83
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
top::TopxAODAlg::m_recalculateNominalWeightSum
bool m_recalculateNominalWeightSum
Definition: TopxAODAlg.h:76
xAOD::IOStats::stats
ReadStats & stats()
Access the object belonging to the current thread.
Definition: IOStats.cxx:17
xAOD::PerfStats::start
void start(bool clear=true)
Start the statistics collection.
top::TopxAODAlg::m_eventSavedParticle
unsigned int m_eventSavedParticle
Definition: TopxAODAlg.h:88
AthHistogramming::book
StatusCode book(const TH1 &hist, const std::string &tDir="", const std::string &stream="")
Simplify the booking and registering (into THistSvc) of histograms.
Definition: AthHistogramming.h:303
top::EventSelectionManager
Maybe you want to run multiple selections (e+jets, mu+jets) on the same input files at the same time.
Definition: EventSelectionManager.h:40
SG::VIEW_ELEMENTS
@ VIEW_ELEMENTS
this data object is a view, it does not own its elmts
Definition: OwnershipPolicy.h:18
top::TopxAODAlg::m_topScaleFactors
ToolHandle< top::ScaleFactorCalculator > m_topScaleFactors
Definition: TopxAODAlg.h:61
AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
top::TopxAODAlg::m_systObjMaker
ToolHandle< top::ObjectCollectionMaker > m_systObjMaker
Definition: TopxAODAlg.h:53
hist_file_dump.d
d
Definition: hist_file_dump.py:137
top::TopxAODAlg::m_boostedTaggersSFSysNames
std::unordered_map< std::string, std::vector< std::string > > m_boostedTaggersSFSysNames
Definition: TopxAODAlg.h:93
EL::AnaAlgorithm::m_inputMetaStore
MetaStore_t m_inputMetaStore
Object accessing the input metadata store.
Definition: AnaAlgorithm.h:551
xAOD::CutBookkeeper_v1
Description of the class that is used to keep track of event counts.
Definition: CutBookkeeper_v1.h:29
top::TopxAODAlg::m_totalYieldSoFar
unsigned int m_totalYieldSoFar
Definition: TopxAODAlg.h:84
EL::AnaAlgorithm::m_outputMetaStore
MetaStore_t m_outputMetaStore
Object accessing the output metadata store.
Definition: AnaAlgorithm.h:555
EL::AnaAlgorithm::inputMetaStore
ConstMetaStorePtr_t inputMetaStore() const
Definition: AnaAlgorithm.cxx:72
xAOD::ReadStats::Print
void Print(::Option_t *option="") const
Print information about the collected statistics.
top::TopxAODAlg::m_totalEvents
ULong64_t m_totalEvents
Definition: TopxAODAlg.h:72
top::TopxAODAlg::m_tracker
std::unique_ptr< top::AnalysisTrackingHelper > m_tracker
Definition: TopxAODAlg.h:56
skel.it
it
Definition: skel.GENtoEVGEN.py:423
top::TopxAODAlg::m_showeringAlgo
SampleXsection::showering m_showeringAlgo
Definition: TopxAODAlg.h:102
top::TopxAODAlg::m_totalEventsPdfWeighted
std::unordered_map< std::string, std::unique_ptr< std::vector< float > > > m_totalEventsPdfWeighted
Definition: TopxAODAlg.h:68
top::TopxAODAlg::m_eventSavedReco
unsigned int m_eventSavedReco
Definition: TopxAODAlg.h:85
AthHistogramAlgorithm::m_prefix
std::string m_prefix
Name of the ROOT output stream (file)
Definition: AthHistogramAlgorithm.h:88
EL::AnaAlgorithm::beginInputFile
virtual ::StatusCode beginInputFile()
perform the action for the beginning of an input file
Definition: AnaAlgorithm.cxx:351
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
top::TopxAODAlg::m_totalEventsWeighted_LHE3
std::vector< double > m_totalEventsWeighted_LHE3
Definition: TopxAODAlg.h:78
AthHistogramAlgorithm::m_histNamePrefix
std::string m_histNamePrefix
The prefix for the histogram THx name.
Definition: AthHistogramAlgorithm.h:94
top::TopxAODAlg::m_settings
std::shared_ptr< ConfigurationSettings > m_settings
Definition: TopxAODAlg.h:49
ANA_CHECK
#define ANA_CHECK(EXP)
check whether the given expression was successful
Definition: Control/AthToolSupport/AsgMessaging/AsgMessaging/MessageCheck.h:324
top::loadObjectSelection
top::TopObjectSelection * loadObjectSelection(const std::shared_ptr< top::TopConfig > &config)
Users may want to define (and load!) their own object selection.
Definition: Tools.cxx:291
AthHistogramAlgorithm::m_histTitlePrefix
std::string m_histTitlePrefix
The prefix for the histogram THx title.
Definition: AthHistogramAlgorithm.h:100
top::TopxAODAlg::m_eventCleaning
ToolHandle< top::EventCleaningSelection > m_eventCleaning
Definition: TopxAODAlg.h:52
top::readMetaData
bool readMetaData(TFile *inputFile, const std::shared_ptr< top::TopConfig > &config)
This function will be used to load the metadata object and pull information from it.
Definition: Tools.cxx:331
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
dumpFileToPlots.pdfName
string pdfName
Definition: dumpFileToPlots.py:21
AthHistogramming::m_histNamePostfix
std::string m_histNamePostfix
The postfix for the histogram THx name.
Definition: AthHistogramming.h:235
EL::AnaAlgorithm::AnaAlgorithm
AnaAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
constructor with parameters
Definition: AnaAlgorithm.cxx:40
EL::AnaAlgorithm::fileExecute
virtual ::StatusCode fileExecute()
perform the action exactly once for each file in the dataset
Definition: AnaAlgorithm.cxx:343
top::TopxAODAlg::m_recalc_LHE3
bool m_recalc_LHE3
Definition: TopxAODAlg.h:75
AthenaPoolTestWrite.stream
string stream
Definition: AthenaPoolTestWrite.py:12
top::TopxAODAlg::m_isAFII
int m_isAFII
Definition: TopxAODAlg.h:92
H5Utils::getTree
std::string getTree(const std::string &file_name)
Definition: getTree.cxx:36
top::TopConfig
Definition: TopConfig.h:41
xAOD::Muon_v1
Class describing a Muon.
Definition: Muon_v1.h:38
xAOD::EventInfo_v1::runNumber
uint32_t runNumber() const
The current event's run number.
top::TopxAODAlg::m_topTools
ToolHandle< top::TopToolStore > m_topTools
Definition: TopxAODAlg.h:51
top::ParticleLevelLoader
Loading tool which creates a particle level event object.
Definition: ParticleLevelLoader.h:53
top::TopxAODAlg::m_year
std::string m_year
Definition: TopxAODAlg.h:103
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
top::TopEventMaker
TopEventMaker produce top::Event objects from xAOD::SystematicEvent objects.
Definition: TopEventMaker.h:32
ANA_CHECK_THROW
#define ANA_CHECK_THROW(EXP)
check whether the given expression was successful, throwing an exception on failure
Definition: Control/AthToolSupport/AsgMessaging/AsgMessaging/MessageCheck.h:339
python.setupRTTAlg.size
int size
Definition: setupRTTAlg.py:39
top::TopxAODAlg::m_sumWeights
TTree * m_sumWeights
Definition: TopxAODAlg.h:74
AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
top::TopxAODAlg::m_totalEventsWeighted_LHE3_temp
std::vector< double > m_totalEventsWeighted_LHE3_temp
Definition: TopxAODAlg.h:80
top::Event::m_muons
xAOD::MuonContainer m_muons
Container of muons (can be sorted)
Definition: Event.h:105
AthCommonDataStore
Definition: AthCommonDataStore.h:52
AthAlgorithm::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthAlgorithm.cxx:66
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
top::TopxAODAlg::m_entry
unsigned int m_entry
Definition: TopxAODAlg.h:82
event
POOL::TEvent event(POOL::TEvent::kClassAccess)
CaloCondBlobAlgs_fillNoiseFromASCII.inputFile
string inputFile
Definition: CaloCondBlobAlgs_fillNoiseFromASCII.py:17
top::Event::m_electrons
xAOD::ElectronContainer m_electrons
Container of electrons (can be sorted)
Definition: Event.h:99
beamspotman.n
n
Definition: beamspotman.py:731
top::TopxAODAlg::m_FTAGMapIndex
unsigned int m_FTAGMapIndex
Definition: TopxAODAlg.h:101
top::TopxAODAlg::m_DSID
unsigned int m_DSID
Definition: TopxAODAlg.h:91
endmsg
#define endmsg
Definition: AnalysisConfig_Ntuple.cxx:63
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
top::ConfigurationSettings::get
static ConfigurationSettings * get(bool reset=false)
Design patterns 101.
Definition: ConfigurationSettings.cxx:714
EL::AnaAlgorithm::requestBeginInputFile
::StatusCode requestBeginInputFile()
register this algorithm to have an implementation of beginInputFile
Definition: AnaAlgorithm.cxx:271
top::TopxAODAlg::m_streamName
std::string m_streamName
Definition: TopxAODAlg.h:100
TH3
Definition: rootspy.cxx:440
top::AnalysisTrackingHelper
Helper for collecting data for analysis tracking.
Definition: AnalysisTrackingHelper.h:20
top::TopxAODAlg::m_firstEvent
unsigned int m_firstEvent
Definition: TopxAODAlg.h:83
checkCorrelInHIST.prefix
dictionary prefix
Definition: checkCorrelInHIST.py:391
dumpFileToPlots.treeName
string treeName
Definition: dumpFileToPlots.py:20
top::check
void check(bool thingToCheck, const std::string &usefulFailureMessage)
Print an error message and terminate if thingToCheck is false.
Definition: EventTools.cxx:15
xAOD::uint64_t
uint64_t
Definition: EventInfo_v1.cxx:123
top::parseCutBookkeepers
void parseCutBookkeepers(xAOD::TEvent &xaodEvent, const std::size_t size, std::vector< std::string > &names, std::vector< float > &sumW)
Search bookkeepers for ones matching AllExecutedEvents, and which originate from AOD before skimming.
Definition: Tools.cxx:89
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
top::ParticleLevelEvent
Definition: ParticleLevelEvent.h:24
xAOD::CutBookkeeperContainer_v1
Container that holds the Container of all CutBookkeepers.
Definition: CutBookkeeperContainer_v1.h:27
EL::AnaAlgorithm::m_hasFileExecute
bool m_hasFileExecute
the value of hasFileExecute
Definition: AnaAlgorithm.h:586
AthHistogramming::m_graphMap
GraphMap_t m_graphMap
The map of TGraph names to their pointers.
Definition: AthHistogramming.h:222
top::TopxAODAlg::m_PDF_SF
ToolHandle< top::PDFScaleFactorCalculator > m_PDF_SF
Definition: TopxAODAlg.h:59
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
AthHistogramming::m_streamName
std::string m_streamName
Name of the ROOT output stream (file)
Definition: AthHistogramming.h:226
AthHistogramming::configAthHistogramming
StatusCode configAthHistogramming(const ServiceHandle< ITHistSvc > &histSvc, const std::string &prefix, const std::string &rootDir, const std::string &histNamePrefix, const std::string &histNamePostfix, const std::string &histTitlePrefix, const std::string &histTitlePostfix)
To be called by the derived classes to fill the internal configuration.
Definition: AthHistogramming.cxx:66
xAOD::IOStats::instance
static IOStats & instance()
Singleton object accessor.
Definition: IOStats.cxx:11
top::CalcTtbarPartonHistory
Definition: CalcTtbarPartonHistory.h:33
DataVector
Derived DataVector<T>.
Definition: DataVector.h:581
AthHistogramming::hash_t
uint32_t hash_t
typedef for the internal hash
Definition: AthHistogramming.h:169
TH2
Definition: rootspy.cxx:373
CxxUtils::crc64
uint64_t crc64(const CRCTable &table, const char *data, size_t data_len)
Find the CRC-64 of a string,.
Definition: crc64.cxx:696
xAOD::FileMetaData_v1
Class holding file-level metadata about an xAOD file.
Definition: FileMetaData_v1.h:34
AthHistogramming::m_histTitlePostfix
std::string m_histTitlePostfix
The postfix for the histogram THx title.
Definition: AthHistogramming.h:241
top::TopxAODAlg::m_doPerfStats
unsigned int m_doPerfStats
Definition: TopxAODAlg.h:65
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
top::TopToolStore
Definition: TopToolStore.h:49
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
top::TopxAODAlg::m_names_LHE3
std::vector< std::string > m_names_LHE3
Definition: TopxAODAlg.h:77
top::getRawEventsBookkeeper
ULong64_t getRawEventsBookkeeper(const xAOD::CutBookkeeperContainer *cutBookKeepers)
Get raw number of events before skimming from "AllExecutedEvents" bookkeeper.
Definition: Tools.cxx:125
top::TopxAODAlg::m_objectSelection
ToolHandle< top::TopObjectSelection > m_objectSelection
Definition: TopxAODAlg.h:54
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
top::TopxAODAlg::m_outputFile
TFile * m_outputFile
Definition: TopxAODAlg.h:66
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:192
ActsTrk::to_string
std::string to_string(const DetectorType &type)
Definition: GeometryDefs.h:34
top::TopxAODAlg::m_crashOnMetaDataFail
bool m_crashOnMetaDataFail
Definition: TopxAODAlg.h:99
EL::AnaAlgorithm::requestEndInputFile
::StatusCode requestEndInputFile()
register this algorithm to have an implementation of endInputFile
Definition: AnaAlgorithm.cxx:292
AthHistogramming::hash
hash_t hash(const std::string &histName) const
Method to calculate a 32-bit hash from a string.
Definition: AthHistogramming.h:428
top::TopxAODAlg::m_totalEventsInFiles
ULong64_t m_totalEventsInFiles
Definition: TopxAODAlg.h:73
module_driven_slicing.getHist
def getHist(name, makehist)
Definition: module_driven_slicing.py:121
xAOD::PerfStats::instance
static PerfStats & instance()
Function accessing the singleton instance.
xAOD::Electron_v1
Definition: Electron_v1.h:34
top::TopxAODAlg::m_totalEventsWeighted_temp
double m_totalEventsWeighted_temp
Definition: TopxAODAlg.h:71
top::TopxAODAlg::m_eventSavedTruth
unsigned int m_eventSavedTruth
Definition: TopxAODAlg.h:87
python.LumiBlobConversion.pos
pos
Definition: LumiBlobConversion.py:18
xAOD::EventInfo_v1
Class describing the basic event information.
Definition: EventInfo_v1.h:43
top::TopxAODAlg::m_eventSelectionManager
std::unique_ptr< top::EventSelectionManager > m_eventSelectionManager
Definition: TopxAODAlg.h:55
xxh3::hash64
std::uint64_t hash64(const void *data, std::size_t size)
Passthrough to XXH3_64bits.
Definition: XXH.cxx:9
DeMoAtlasDataLoss.runNumber
string runNumber
Definition: DeMoAtlasDataLoss.py:64
AthAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Definition: AthAlgorithm.h:79
top::TopxAODAlg::m_eventSaver
std::unique_ptr< top::EventSaverBase > m_eventSaver
Definition: TopxAODAlg.h:62
a
TList * a
Definition: liststreamerinfos.cxx:10
AthHistogramming::m_effMap
EffMap_t m_effMap
The map of histogram names to their pointers.
Definition: AthHistogramming.h:208
h
AthHistogramming::m_histTitlePrefix
std::string m_histTitlePrefix
The prefix for the histogram THx title.
Definition: AthHistogramming.h:238
lumiFormat.fill
fill
Definition: lumiFormat.py:111
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
top::TopxAODAlg::m_totalEventsWeighted
double m_totalEventsWeighted
Definition: TopxAODAlg.h:70
python.Constants.INFO
int INFO
Definition: Control/AthenaCommon/python/Constants.py:16
AthHistogramAlgorithm::m_histNamePostfix
std::string m_histNamePostfix
The postfix for the histogram THx name.
Definition: AthHistogramAlgorithm.h:97
TH1
Definition: rootspy.cxx:268
AthHistogramming::m_treeMap
TreeMap_t m_treeMap
The map of TTree names to their pointers.
Definition: AthHistogramming.h:215
top::TopxAODAlg::m_derivationStream
std::string m_derivationStream
Definition: TopxAODAlg.h:98
AthCommonMsg< Algorithm >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
AthHistogramming::m_histNamePrefix
std::string m_histNamePrefix
The prefix for the histogram THx name.
Definition: AthHistogramming.h:232
EL::AnaAlgorithm::m_hasEndInputFile
bool m_hasEndInputFile
the value of hasEndInputFile
Definition: AnaAlgorithm.h:594
top::TopxAODAlg::m_isFirst
bool m_isFirst
Definition: TopxAODAlg.h:89
top::TopxAODAlg::m_pmg_weightTool
ToolHandle< PMGTools::IPMGTruthWeightTool > m_pmg_weightTool
Definition: TopxAODAlg.h:63
PowhegPythia8EvtGen_jetjet.pdf
pdf
Definition: PowhegPythia8EvtGen_jetjet.py:4
top::ObjectCollectionMaker
Definition: ObjectCollectionMaker.h:45
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:616
top::TopxAODAlg::m_pdfMetadataExists
bool m_pdfMetadataExists
Definition: TopxAODAlg.h:69
AthHistogramAlgorithm::m_histTitlePostfix
std::string m_histTitlePostfix
The postfix for the histogram THx title.
Definition: AthHistogramAlgorithm.h:103
declareProperty
#define declareProperty(n, p, h)
Definition: BaseFakeBkgTool.cxx:15
top::renameCutBookkeepers
void renameCutBookkeepers(std::vector< std::string > &bookkeeper_names, const std::vector< std::string > &pmg_weight_names)
Rename CutBookkeeper names according to MC generator weight names reported by PMGTruthWeightTool.
Definition: Tools.cxx:143
AthHistogramming::hist
TH1 * hist(const std::string &histName, const std::string &tDir="", const std::string &stream="")
Simplify the retrieval of registered histograms of any type.
Definition: AthHistogramming.cxx:198
xAOD::PerfStats::stop
void stop()
Stop the statistics collection.
top::TopxAODAlg::m_topEventMaker
ToolHandle< top::TopEventMaker > m_topEventMaker
Definition: TopxAODAlg.h:60
AthHistogramming::m_histSvc
ServiceHandle< ITHistSvc > m_histSvc
Pointer to the THistSvc (event store by default)
Definition: AthHistogramming.h:194
dqt_zlumi_alleff_HIST.eff
int eff
Definition: dqt_zlumi_alleff_HIST.py:113
str
Definition: BTagTrackIpAccessor.cxx:11
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:790
top::Event
Very simple class to hold event data after reading from a file.
Definition: Event.h:49
top::TopxAODAlg::m_topPartonHistory
ToolHandle< top::CalcTtbarPartonHistory > m_topPartonHistory
Definition: TopxAODAlg.h:58
top::TopxAODAlg::m_sumPdfWeights
TTree * m_sumPdfWeights
Definition: TopxAODAlg.h:67
top::TopxAODAlg::m_libraryNames
std::string m_libraryNames
Definition: TopxAODAlg.h:97
AthHistogramming::m_name
std::string m_name
Instance name.
Definition: AthHistogramming.h:245
top::ScaleFactorCalculator
Definition: ScaleFactorCalculator.h:44
AthHistogramAlgorithm::m_rootDir
std::string m_rootDir
Name of the ROOT directory.
Definition: AthHistogramAlgorithm.h:91
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
top::loadEventSaver
top::EventSaverBase * loadEventSaver(const std::shared_ptr< top::TopConfig > &config)
Users may also want to write out custom ntuples / xAODs.
Definition: Tools.cxx:314
top::PDFScaleFactorCalculator
For testing PDF reweighting with LHAPDF6.
Definition: PDFScaleFactorCalculator.h:37
top::TopxAODAlg::m_configPath
std::string m_configPath
properties
Definition: TopxAODAlg.h:96
top::EventCleaningSelection
Definition: EventCleaningSelection.h:32
fitman.k
k
Definition: fitman.py:528
top::TopxAODAlg::m_config
std::shared_ptr< TopConfig > m_config
Definition: TopxAODAlg.h:50
AthHistogramming::buildBookingString
void buildBookingString(std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
Method to build individual booking string.
Definition: AthHistogramming.cxx:560
ServiceHandle< IIncidentSvc >
AthHistogramming::bookGetPointer
TH1 * bookGetPointer(const TH1 &hist, const std::string &tDir="", const std::string &stream="")
Simplify the booking and registering (into THistSvc) of histograms.
Definition: AthHistogramming.h:260
top::Event::m_saveEvent
bool m_saveEvent
save the event?
Definition: Event.h:87
Generated on Sun Jun 30 2024 22:28:39 for ATLAS Offline Software by   doxygen 1.8.18