#include <VertexDecoratorAlg.h>

Inheritance diagram for InDetGNNHardScatterSelection::VertexDecoratorAlg:

Public Member Functions
	VertexDecoratorAlg (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize () override
StatusCode	execute (const EventContext &ctx) override
	Execute method.
virtual StatusCode	sysInitialize ()
	Initialization method invoked by the framework.
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`.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
	Get filter decision:
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
	Set filter decision:
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
virtual bool	isReEntrant () const override final
	Legacy algorithms are not thread-safe.
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
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.
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.
TH1 *	bookGetPointer (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
TH1 *	bookGetPointer (TH1 *hist, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
TH1 *	bookGetPointer (TH1 &histRef, std::string tDir="", std::string stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
TTree *	bookGetPointer (const TTree &treeRef, std::string tDir="", std::string stream="")
	Simplify the booking and registering (into THistSvc) of TTrees.
TGraph *	bookGetPointer (const TGraph &graphRef, std::string tDir="", std::string stream="")
	Simplify the booking and registering (into THistSvc) of TGraphs.
TEfficiency *	bookGetPointer (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency *	bookGetPointer (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency *	bookGetPointer (TEfficiency &effRef, std::string tDir="", std::string stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode	book (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
StatusCode	book (TH1 *hist, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
StatusCode	book (TH1 &histRef, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
StatusCode	book (const TTree &treeRef, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TTrees.
StatusCode	book (const TGraph &graphRef, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TGraphs.
StatusCode	book (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode	book (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode	book (TEfficiency &effRef, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
TH1 *	hist (const std::string &histName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered histograms of any type.
TH2 *	hist2d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered 2-d histograms.
TH3 *	hist3d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered 3-d histograms.
TTree *	tree (const std::string &treeName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered TTrees.
TGraph *	graph (const std::string &graphName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered TGraphs.
TEfficiency *	efficiency (const std::string &effName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered TEfficiency.

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

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>
void	buildBookingString (std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
	Method to build individual booking string.
void	myReplace (std::string &str, const std::string &oldStr, const std::string &newStr)
	Helper method to replace sub-string.
hash_t	hash (const std::string &histName) const
	Method to calculate a 32-bit hash from a string.

Private Attributes
SG::ReadHandleKey< xAOD::VertexContainer >	m_vertexInKey
SG::ReadHandleKey< xAOD::EventInfo >	m_eventInKey
SG::ReadHandleKey< xAOD::ElectronContainer >	m_electronsInKey
SG::ReadHandleKey< xAOD::MuonContainer >	m_muonsInKey
SG::ReadHandleKey< xAOD::JetContainer >	m_jetsInKey
SG::ReadHandleKey< xAOD::PhotonContainer >	m_photonsInKey
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_ntrk
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_sumPt
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_chi2Over_ndf
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_z_asym
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_weighted_z_asym
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_weighted_z_kurt
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_z_skew
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_photon_deltaz
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_photon_deltaPhi
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_actualInterPerXing
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_photonLinksKey
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_jetLinksKey
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_electronLinksKey
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_muonLinksKey
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_gnnScore
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_nElectrons
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_nMuons
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_nJets
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_mDecor_nPhotons
SG::WriteHandleKey< xAOD::CompositeParticleContainer >	m_multiPhotonsOutKey
SG::WriteDecorHandleKey< xAOD::VertexContainer >	m_multiPhotonLinksKey
SG::ReadDecorHandleKey< xAOD::VertexContainer >	m_deltaZKey
SG::ReadDecorHandleKey< xAOD::VertexContainer >	m_deltaPhiKey
SG::ReadDecorHandleKey< xAOD::PhotonContainer >	m_caloPointingZKey
SG::ReadDecorHandleKey< xAOD::PhotonContainer >	m_zCommonKey
SG::ReadDecorHandleKey< xAOD::PhotonContainer >	m_zCommonErrorKey
ToolHandle< CP::TrackVertexAssociationTool >	m_trkVtxAssociationTool
ToolHandle< GNNTool >	m_gnnTool
ServiceHandle< ITHistSvc >	m_histSvc
	Default constructor: AthHistogramAlgorithm();.
std::string	m_prefix
	Name of the ROOT output stream (file).
std::string	m_rootDir
	Name of the ROOT directory.
std::string	m_histNamePrefix
	The prefix for the histogram THx name.
std::string	m_histNamePostfix
	The postfix for the histogram THx name.
std::string	m_histTitlePrefix
	The prefix for the histogram THx title.
std::string	m_histTitlePostfix
	The postfix for the histogram THx title.
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default).
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared
HistMap_t	m_histMap
	The map of histogram names to their pointers.
EffMap_t	m_effMap
	The map of histogram names to their pointers.
TreeMap_t	m_treeMap
	The map of TTree names to their pointers.
GraphMap_t	m_graphMap
	The map of TGraph names to their pointers.
std::string	m_streamName
	Name of the ROOT output stream (file).
std::string	m_name
	Instance name.
MsgStream	m_msg
	Cached Message Stream.

Detailed Description

Definition at line 40 of file VertexDecoratorAlg.h.

Member Typedef Documentation

◆ EffMap_t

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

privateinherited

Typedef for convenience.

Definition at line 207 of file AthHistogramming.h.

◆ GraphMap_t

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

privateinherited

Typedef for convenience.

Definition at line 221 of file AthHistogramming.h.

◆ hash_t

typedef uint32_t AthHistogramming::hash_t

privateinherited

typedef for the internal hash

Definition at line 171 of file AthHistogramming.h.

◆ HistMap_t

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

privateinherited

Typedef for convenience.

Definition at line 200 of file AthHistogramming.h.

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::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 214 of file AthHistogramming.h.

Constructor & Destructor Documentation

◆ VertexDecoratorAlg()

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

Definition at line 19 of file VertexDecoratorAlg.cxx.

      : AthHistogramAlgorithm (name, pSvcLocator)
  {
  }

Member Function Documentation

◆ 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 337 of file AthHistogramming.h.

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

◆ 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 305 of file AthHistogramming.h.

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

◆ 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 405 of file AthHistogramming.h.

{
  // Call the other Book method and see if it returns a valid pointer
  TTree* treePointer = this->bookGetPointer( treeRef, tDir, stream );
  if ( treePointer )
    {
      return StatusCode::SUCCESS;
    }
  else
    {
      return StatusCode::FAILURE;
    }
}

◆ 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 358 of file AthHistogramming.h.

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

◆ 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 348 of file AthHistogramming.h.

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

◆ 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 326 of file AthHistogramming.h.

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

◆ 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 316 of file AthHistogramming.h.

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

◆ 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 283 of file AthHistogramming.h.

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

◆ 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 388 of file AthHistogramming.cxx.

{
  // Get a pointer
  const TGraph* graphPointer = &graphRef;
 
  // Check that we got a valid pointer
  if ( !graphPointer )
    {
      m_msg << MSG::WARNING
            << "We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
 
  // Modify the name and title according to the prefixes of this classes instance
  std::string graphName  = graphPointer->GetName();
  const std::string graphTitle = graphPointer->GetTitle();
 
  // Check if the hash for this graphName already exists, i.e., if we have a hash collision
  const hash_t graphHash = this->hash(graphName);
  GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
  if ( it != m_graphMap.end() ) // It does exist!
    {
      m_msg << MSG::WARNING
            << "Detected a hash collision. The hash for the TGraph with name=" << graphName
            << " already exists and points to a TGraph with name=" << it->second->GetName()
            << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
      return NULL;
    }
 
  // Create a clone that has the new name
  TGraph* graphClone = dynamic_cast< TGraph* >( graphPointer->Clone((m_histNamePrefix+graphName+m_histNamePostfix).c_str()) );
  if( !graphClone )
    {
      m_msg << MSG::WARNING
            << "We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
  graphClone->SetTitle ((m_histTitlePrefix+graphTitle+m_histTitlePostfix).c_str());
 
  // Massage the final string to book things
  std::string bookingString("");
  this->buildBookingString( bookingString, graphName, tDir, stream );
 
  // Register the TGraph into the THistSvc
  if ( !((histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
    {
      m_msg << MSG::WARNING
            << "Problem registering TGraph with name " << graphName
            << ", title " << graphTitle
            << " in " << m_name << "!" << endmsg;
      return NULL;
    }
 
  // Also register it in the local map of string to pointer
  m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
 
  return graphClone;
}

◆ 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 262 of file AthHistogramming.h.

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

◆ 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 273 of file AthHistogramming.cxx.

{
  // Get a pointer
  const TTree* treePointer = &treeRef;
 
  // Check that we got a valid pointer
  if ( !treePointer )
    {
      m_msg << MSG::WARNING
            << "We got an invalid TTree pointer in the BookGetPointer(TTree*) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
 
  // Modify the name and title according to the prefixes of this classes instance
  std::string treeName  = treePointer->GetName();
  const std::string treeTitle = treePointer->GetTitle();
 
  // Check if the hash for this treeName already exists, i.e., if we have a hash collision
  const hash_t treeHash = this->hash(treeName);
  TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
  if ( it != m_treeMap.end() ) // It does exist!
    {
      m_msg << MSG::WARNING
            << "Detected a hash collision. The hash for the TTree with name=" << treeName
            << " already exists and points to a TTree with name=" << it->second->GetName()
            << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
      return NULL;
    }
 
  // Create a clone that has the new name
  TTree* treeClone = dynamic_cast< TTree* >( treePointer->Clone(treeName.c_str()) );
  if( !treeClone )
    {
      m_msg << MSG::WARNING
            << "We couldn't clone the TTree in the BookGetPointer(TTree&) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
  treeClone->SetTitle (treeTitle.c_str());
 
  // Massage the final string to book things
  std::string bookingString("");
  this->buildBookingString( bookingString, treeName, tDir, stream );
 
  // Register the TTree into the THistSvc
  if ( !((histSvc()->regTree(bookingString, treeClone)).isSuccess()) )
    {
      m_msg << MSG::WARNING
            << "Problem registering TTree with name " << treeName
            << ", title " << treeTitle
            << " in " << m_name << "!" << endmsg;
      return NULL;
    }
 
  // Also register it in the local map of string to pointer
  m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treeClone ) );
 
  return treeClone;
}

◆ 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 123 of file AthHistogramming.cxx.

{
  std::string effName(effRef.GetName());
  const std::string effTitle(effRef.GetTitle());
  std::string bookingString;
  this->buildBookingString(bookingString, effName, tDir, stream);
  const std::string finalEffName = m_histNamePrefix + effName + m_histNamePostfix;
  effRef.SetTitle((m_histTitlePrefix + effTitle + m_histTitlePostfix).c_str());
  effRef.SetName(finalEffName.c_str());
 
  const hash_t effHash = this->hash(effName);
  EffMap_t::const_iterator it = m_effMap.find(effHash);
  if (it != m_effMap.end()) {
    m_msg << MSG::WARNING
          << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
          << " already exists and points to a TEfficiency with name=" << it->second->GetName()
          << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
    return nullptr;
  }
 
  if (!histSvc()->regEfficiency(bookingString, &effRef).isSuccess()) {
    m_msg << MSG::WARNING
          << "Problem registering TEfficiency with name " << effName
          << ", name prefix " << m_histNamePrefix
          << ", title " << effTitle
          << ", title prefix " << m_histTitlePrefix
          << ", and title postfix " << m_histTitlePostfix
          << " in " << m_name << "!" << endmsg;
    return nullptr;
  }
 
  m_effMap.insert(m_effMap.end(), std::pair<const hash_t, TEfficiency*>(effHash, &effRef));
 
  return &effRef;
}

◆ 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 295 of file AthHistogramming.h.

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

◆ 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 83 of file AthHistogramming.cxx.

{
  std::string histName(histRef.GetName());
  const std::string histTitle(histRef.GetTitle());
  std::string bookingString;
 
  this->buildBookingString(bookingString, histName, tDir, stream);
 
  const std::string finalHistName = m_histNamePrefix + histName + m_histNamePostfix;
 
  histRef.SetTitle((m_histTitlePrefix + histTitle + m_histTitlePostfix).c_str());
  histRef.SetName(finalHistName.c_str());
 
  const hash_t histHash = this->hash(histName);
  HistMap_t::const_iterator it = m_histMap.find(histHash);
  if (it != m_histMap.end()) {
    m_msg << MSG::WARNING
          << "Detected a hash collision. The hash for the histogram with name=" << histName
          << " already exists and points to a histogram with name=" << it->second->GetName()
          << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
    return nullptr;
  }
 
  if (!histSvc()->regHist(bookingString, &histRef).isSuccess()) {
    m_msg << MSG::WARNING
          << "Problem registering histogram with name " << histName
          << ", name prefix " << m_histNamePrefix
          << ", title " << histTitle
          << ", title prefix " << m_histTitlePrefix
          << ", and title postfix " << m_histTitlePostfix
          << " in " << m_name << "!" << endmsg;
    return nullptr;
  }
 
  m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, &histRef));
 
  return &histRef;
}

◆ 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 274 of file AthHistogramming.h.

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

◆ 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 521 of file AthHistogramming.cxx.

{
  // Massage the final string to book things
  if(tDir.empty()) tDir = m_rootDir;
  size_t pos = histName.rfind('/');
  if(pos != std::string::npos){
    tDir+='/';
    tDir.append(histName, 0,pos);
    histName.erase(0,pos+1);
  };
  if(stream.empty()) stream = m_streamName;
 
  if(usePrefixPostfix){
    bookingString = "/"+stream+"/"+tDir+"/"+m_histNamePrefix+histName+m_histNamePostfix;
  } else {
    bookingString = "/"+stream+"/"+tDir+"/"+histName;
  }
  while(bookingString.find("//") != std::string::npos){
    this->myReplace(bookingString,"//","/");
  }
 
  return;
}

◆ 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 51 of file AthHistogramming.cxx.

{
  m_histSvc          = histSvc;
  m_streamName       = prefix;
  m_rootDir          = rootDir;
  m_histNamePrefix   = histNamePrefix;
  m_histNamePostfix  = histNamePostfix;
  m_histTitlePrefix  = histTitlePrefix;
  m_histTitlePostfix = histTitlePostfix;
 
  return StatusCode::SUCCESS;
}

◆ declareGaudiProperty()

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

◆ detStore()

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

AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ 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 211 of file AthHistogramming.cxx.

{
  // Build a 32 bit hash out of the name
  const hash_t effHash = this->hash(effName);
 
  // See if this entry exists in the map
  EffMap_t::const_iterator it = m_effMap.find( effHash );
  if ( it == m_effMap.end() ) // It doesn't exist!
    { // Let's see into the THistSvc if somebody else has registered the TEfficiency...
 
      // Need to copy the strings as we will massage them from here on
      std::string effNameCopy = effName;
      std::string tDirCopy     = tDir;
      std::string streamCopy   = stream;
 
      // Massage the final string to book things
      std::string bookingString("");
      this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy ,false);
 
      TEfficiency* effPointer(NULL);
      if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
        {
          // Massage the final string to book things
          std::string bookingString("");
          this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy, true );
 
          if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
            {
              m_msg << MSG::WARNING
                    << "Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
                    << m_histNamePrefix + effNameCopy + m_histNamePostfix
                    << " or with name " << effNameCopy
                    << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
                    << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
              return NULL;
            }
          // If we get to here, we actually found the TEfficiency in the THistSvc.
          // So let's add it to the local cache map and return its pointer
          m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
          return effPointer;
        }
      // If we get to here, we actually found the TEfficiency in the THistSvc.
      // So let's add it to the local cache map and return its pointer
      m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
      return effPointer;
    }
 
  // Return the pointer to the TEfficiency that we got from the local cache map
  return it->second;
}

◆ evtStore()

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

AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ execute()

StatusCode VertexDecoratorAlg::execute ( const EventContext & ctx )

overridevirtual

Execute method.

Implements AthAlgorithm.

Definition at line 104 of file VertexDecoratorAlg.cxx.

  {
    SG::ReadHandle<xAOD::VertexContainer> vertices(m_vertexInKey, ctx);
    ATH_CHECK(vertices.isValid());
 
    SG::ReadHandle<xAOD::PhotonContainer> photonsIn(m_photonsInKey, ctx);
    ATH_CHECK(photonsIn.isValid());
    
    SG::ReadHandle<xAOD::ElectronContainer> electronsIn(m_electronsInKey, ctx);
    ATH_CHECK(electronsIn.isValid());
    SG::ReadHandle<xAOD::MuonContainer> muonsIn(m_muonsInKey, ctx);
    ATH_CHECK(muonsIn.isValid());
    SG::ReadHandle<xAOD::JetContainer> jetsIn(m_jetsInKey, ctx);
    ATH_CHECK(jetsIn.isValid());
    SG::ReadHandle<xAOD::EventInfo> eventInfo(m_eventInKey,ctx);
    ATH_CHECK(eventInfo.isValid());
 
    SG::WriteDecorHandle<xAOD::VertexContainer, std::vector<ElementLink<xAOD::PhotonContainer>>> dec_photonLinks(m_photonLinksKey, ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, std::vector<ElementLink<xAOD::JetContainer>>> dec_jetLinks(m_jetLinksKey, ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, std::vector<ElementLink<xAOD::ElectronContainer>>> dec_electronLinks(m_electronLinksKey, ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, std::vector<ElementLink<xAOD::MuonContainer>>> dec_muonLinks(m_muonLinksKey, ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, std::vector<ElementLink<xAOD::CompositeParticleContainer>>> dec_multiPhotonLinks(m_multiPhotonLinksKey, ctx);
 
    SG::ReadDecorHandle<xAOD::VertexContainer, float> acc_deltaZ(m_deltaZKey, ctx);
    SG::ReadDecorHandle<xAOD::VertexContainer, float> acc_deltaPhi(m_deltaPhiKey, ctx);
    SG::ReadDecorHandle<xAOD::PhotonContainer, float> acc_caloPointingZ(m_caloPointingZKey, ctx);
    SG::ReadDecorHandle<xAOD::PhotonContainer, float> acc_zCommon(m_zCommonKey, ctx);
    SG::ReadDecorHandle<xAOD::PhotonContainer, float> acc_zCommonError(m_zCommonErrorKey, ctx);
 
    SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_gnnScore(m_mDecor_gnnScore,ctx);
    
    // Decorations needed by the GNNTool
    SG::WriteDecorHandle<xAOD::VertexContainer, int> dec_ntrk(m_mDecor_ntrk, ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_sumPt(m_mDecor_sumPt,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_chi2Over_ndf(m_mDecor_chi2Over_ndf,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_z_asym(m_mDecor_z_asym,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_weighted_z_asym(m_mDecor_weighted_z_asym,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_z_kurt(m_mDecor_weighted_z_kurt,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_z_skew(m_mDecor_z_skew,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_photon_deltaz(m_mDecor_photon_deltaz,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_photon_deltaPhi(m_mDecor_photon_deltaPhi,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, float> dec_actualInterPerXing(m_mDecor_actualInterPerXing,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, int> dec_nElectrons(m_mDecor_nElectrons,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, int> dec_nMuons(m_mDecor_nMuons,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, int> dec_nJets(m_mDecor_nJets,ctx);
    SG::WriteDecorHandle<xAOD::VertexContainer, int> dec_nPhotons(m_mDecor_nPhotons,ctx);
    SG::WriteHandle<xAOD::CompositeParticleContainer> mpHandle(m_multiPhotonsOutKey, ctx);
 
    auto mpCont = std::make_unique<xAOD::CompositeParticleContainer>();
    auto mpAux  = std::make_unique<xAOD::CompositeParticleAuxContainer>();
    mpCont->setStore(mpAux.get());
 
    const xAOD::Vertex* bestVtxForMP = nullptr;
    ElementLink<xAOD::CompositeParticleContainer> mpLink; // will be valid only if created
    bool haveMP = false;
 
    // For MT safety, decoration keys must be initialized and available upfront.
    const bool haveZCommonDecor = acc_zCommon.isAvailable() && acc_zCommonError.isAvailable();
    const bool haveCaloPointingDecor = acc_caloPointingZ.isAvailable();
    
    if (photonsIn->size() > 1 && (!haveZCommonDecor || !haveCaloPointingDecor))
    {
      ATH_MSG_ERROR("photonsIn has size > 1 but required decorations are missing: "
                    << m_zCommonKey.key() << " and/or " << m_zCommonErrorKey.key()
                    << " and/or " << m_caloPointingZKey.key());
      return StatusCode::FAILURE;
    }
 
    if (photonsIn->size() > 1)
    {
      // Combined pointing Z (weighted)
      double sumW = 0.0;
      double sumWZ = 0.0;
 
      TLorentzVector p4sum;
      int nUsed = 0;
      for (const xAOD::Photon* ph : *photonsIn)
      {
        const float zCommon = acc_zCommon(*ph);
        const float zCaloPointing = acc_caloPointingZ(*ph);
        // Value-based selection only: prefer zCommon, then fallback to caloPointing value.
        const bool useZCommonValue = std::isfinite(zCommon);
        const float z = useZCommonValue ? zCommon : zCaloPointing;
        if (!std::isfinite(z)) {
          ATH_MSG_ERROR("Non-finite photon pointing values found in both "
                        << m_zCommonKey.key() << " and " << m_caloPointingZKey.key());
          return StatusCode::FAILURE;
        }
 
        // Keep weighting consistent with the z source:
        // use zCommonError only when z comes from zCommon.
        double w = 1.0;
        if (useZCommonValue) {
          const float s = acc_zCommonError(*ph);
          if (std::isfinite(s) && s > 0.0f) w = 1.0 / (double(s) * double(s));
        }
        
        sumW  += w;
        sumWZ += w * z;
        p4sum += ph->p4();
        ++nUsed;
      }
 
      if (nUsed > 1 && sumW > 0.0)
      {
        const float zPoint = sumWZ / sumW;
 
        // Choose vertex closest in z to combined pointing
        float bestAbsDZ = 1e30;
        for (const xAOD::Vertex* vtx : *vertices)
        {
          if (vtx->vertexType() == xAOD::VxType::NoVtx) continue;
 
          const float dz = zPoint - vtx->z();
          const float adz = std::abs(dz);
          if (adz < bestAbsDZ)
          {
            bestAbsDZ = adz;
            bestVtxForMP = vtx;
          }
        }
 
        if (bestVtxForMP)
        {
          // Create exactly one CompositeParticle node
          auto* mp = new xAOD::CompositeParticle();
          mpCont->push_back(mp);
 
          mp->setP4(p4sum);
 
          SG::AuxElement::Decorator<float> dec_mp_deltaZ("deltaZ");
          SG::AuxElement::Decorator<float> dec_mp_deltaPhi("deltaPhi");
          SG::AuxElement::Decorator<int>   dec_mp_nPhotons("nPhotons");
          SG::AuxElement::Decorator<float> dec_mp_zPointing("zPointing");
 
          const float deltaZ = zPoint - bestVtxForMP->z();
          dec_mp_deltaZ(*mp) = deltaZ;
          dec_mp_zPointing(*mp) = zPoint;
          dec_mp_nPhotons(*mp) = static_cast<int>(photonsIn->size());
 
          float dphi = 0.0f;
          if (acc_deltaPhi.isAvailable())
          {
            dphi = acc_deltaPhi(*bestVtxForMP);
            if (!std::isfinite(dphi)) dphi = 0.0f;
          }
          dec_mp_deltaPhi(*mp) = dphi;
          mpLink.setElement(mp);
          mpLink.setStorableObject(*mpCont, true);
 
          haveMP = true;
        }
      }
    }
 
    // Record the MultiPhotons container (even if empty)
    ATH_CHECK(mpHandle.record(std::move(mpCont), std::move(mpAux)));
 
    std::map< const xAOD::Vertex*, std::vector<ElementLink<xAOD::JetContainer>> > jetsInVertex;
    std::map< const xAOD::Jet*, std::map< const xAOD::Vertex*, int> > jetVertexPt;
 
    // initialize jet-vertex map 
    for(const xAOD::Vertex *vertex : *vertices){
      jetsInVertex[vertex] = {};
      for(const xAOD::Jet *jet : *jetsIn){
        jetVertexPt[jet][vertex] = 0;
      }
    }
 
    // pre-fill the jet-vertex map
    for(const xAOD::Jet* jet : *jetsIn){
 
      // for each jet, calculate the track pT associated to each vertex
      std::vector<const xAOD::TrackParticle*> ghostTracks = jet->getAssociatedObjects<xAOD::TrackParticle >(xAOD::JetAttribute::GhostTrack);
      for(const xAOD::TrackParticle* jtrk : ghostTracks){
        if( !jtrk ) continue;
        auto jetTrackVertex = m_trkVtxAssociationTool->getUniqueMatchVertexLink(*jtrk, *vertices);
        if(jetTrackVertex) jetVertexPt[jet][*jetTrackVertex] += jtrk->pt();
      }
 
      // find vertex with the largest fraction of jet track pt
      float maxPtFrac = -1;
      const xAOD::Vertex* uniqueVertexAddress = nullptr;
      for (const xAOD::Vertex *vertex : *vertices) {
        if (vertex->vertexType() == xAOD::VxType::NoVtx) continue;
        if(jetVertexPt[jet][vertex] > maxPtFrac){
          maxPtFrac = jetVertexPt[jet][vertex];
          uniqueVertexAddress = vertex;
        }
      }
 
      // add jet to that vertex's vector of links
      ElementLink<xAOD::JetContainer> jetLink;
      jetLink.setElement(jet);
      jetLink.setStorableObject(*jetsIn.ptr(), true);
      jetsInVertex[uniqueVertexAddress].push_back(jetLink);
    }
 
    for (const xAOD::Vertex *vertex : *vertices)
    {
      if (vertex->vertexType() == xAOD::VxType::NoVtx)
        continue;
 
      if (vertex->nTrackParticles() < 2)
        continue;
 
      dec_actualInterPerXing(*vertex) = eventInfo->actualInteractionsPerCrossing();
 
      // variables for calculation of delta Z asymmetry and delta d asymmetry
      float z_asym = 0;
      float sumDZ = 0;
      float deltaZ = 0;
      float modsumDZ = 0;
      float weighted_sumDZ = 0;
      float weighted_deltaZ = 0;
      float weighted_modsumDZ = 0;
      float weighted_z_asym = 0;
 
      // make vector
      std::vector<float> track_deltaZ;
 
      for (size_t i = 0; i < vertex->nTrackParticles(); i++) {
 
        const xAOD::TrackParticle *trackTmp = vertex->trackParticle(i);
 
        if(!trackTmp) continue;
 
        deltaZ = trackTmp->z0() + trackTmp->vz() - vertex->z();
        track_deltaZ.push_back(deltaZ);
        // get the track weight for each track to get the deltaZ/trk_weight
        float trk_weight = vertex->trackWeight(i);
        weighted_deltaZ = deltaZ * trk_weight;
        // sum of delta z
        sumDZ += deltaZ;
        modsumDZ += std::abs(deltaZ);
        weighted_sumDZ += weighted_deltaZ;
        weighted_modsumDZ += std::abs(weighted_deltaZ);
      } // end loop over tracks
 
      if (modsumDZ > 0) {
        z_asym = sumDZ / modsumDZ;
      }
      if (weighted_modsumDZ > 0) {
        weighted_z_asym = weighted_sumDZ / weighted_modsumDZ;
      }
 
      const float number_tracks = track_deltaZ.size(); // get number of tracks
      const float mean_Dz =
        number_tracks > 0 ? sumDZ / number_tracks : 0.F; // calculate average
 
      float z_skew = 0; // skewness of DeltaZ asymmetry
      float z_kurt = 0; // Kurtosis of DeltaZ asymmetry
      float z_var = 0;  // variance of DeltaZ
 
      for (auto i : track_deltaZ)
      {
        float z_zbar = (i - mean_Dz);
        z_var += std::pow(z_zbar, 2);
        z_skew += std::pow(z_zbar, 3);
        z_kurt += std::pow(z_zbar, 4);
      }
      if (number_tracks > 1 && z_var > 0) {
        z_var /= (number_tracks - 1);
        float z_sd = std::sqrt(z_var);
        const float skew_denom = (number_tracks - 1) * std::pow(z_sd, 3);
        const float kurt_denom = (number_tracks - 1) * std::pow(z_sd, 4);
        if (std::isfinite(skew_denom) && skew_denom != 0.F) {
          z_skew /= skew_denom;
        } else {
          z_skew = 0.F;
        }
        if (std::isfinite(kurt_denom) && kurt_denom != 0.F) {
          z_kurt /= kurt_denom;
        } else {
          z_kurt = 0.F;
        }
      }
      else
      {
        ATH_MSG_WARNING("z momenta are NaN: setting to zero");
        z_skew = 0.;
        z_kurt = 0.;
      }
 
      dec_ntrk(*vertex) = number_tracks;
 
      if(!dec_sumPt.isAvailable()){
        dec_sumPt(*vertex) = xAOD::PVHelpers::getVertexSumPt(vertex, 1, false);
      }
      const float numberDoF = vertex->numberDoF();
      dec_chi2Over_ndf(*vertex) =
        numberDoF > 0.F ? vertex->chiSquared() / numberDoF : 0.F;
      dec_z_asym(*vertex) = z_asym;
      dec_weighted_z_asym(*vertex) = weighted_z_asym;
      dec_z_kurt(*vertex) = z_kurt;
      dec_z_skew(*vertex) = z_skew;
 
      if (acc_deltaZ.isAvailable() && photonsIn->size() > 0 ) {
        //protect against rare NaNs before assigning decorator: setting to 0 (-999 cause NaNs)
        if (!std::isfinite(acc_deltaZ(*vertex))) {
          ATH_MSG_WARNING("photon deltaPhi is NaN: setting to -1!");
          dec_photon_deltaz(*vertex) = -1;
        }
        else{
        dec_photon_deltaz(*vertex) = acc_deltaZ(*vertex);
        }
      }
      else{
       dec_photon_deltaz(*vertex) = -1;
      }
      if (acc_deltaPhi.isAvailable() && photonsIn->size() > 0) {
        
        if (!std::isfinite(acc_deltaPhi(*vertex))) {
          ATH_MSG_WARNING("photon deltaPhi is NaN: setting to 0!");
          dec_photon_deltaPhi(*vertex) = -1;
        }
        else{
        dec_photon_deltaPhi(*vertex) = acc_deltaPhi(*vertex);
        }
      }
      else{
       dec_photon_deltaPhi(*vertex) = -1;
      }
 
      // associate objects to vertices
      auto countValid = [](const auto& links) -> int {
        int n = 0;
        for (const auto& l : links) if (l.isValid()) ++n;
        return n;
      };
 
      std::vector<ElementLink<xAOD::ElectronContainer>> electronLinks;
      for(const xAOD::Electron* electron : *electronsIn){
        const auto *id_trk = xAOD::EgammaHelpers::getOriginalTrackParticle(electron);
        if(!id_trk) continue;
        auto eleVertex = m_trkVtxAssociationTool->getUniqueMatchVertexLink(*id_trk, *vertices);
        if(!eleVertex) continue;
        ElementLink<xAOD::ElectronContainer> elLink;
        if(*eleVertex == vertex){
            elLink.setElement(electron);
            elLink.setStorableObject(*electronsIn.ptr(), true);
            electronLinks.push_back(elLink);
        }
      }
      dec_electronLinks(*vertex) = electronLinks;
      dec_nElectrons(*vertex) = countValid(electronLinks);
 
      std::vector<ElementLink<xAOD::PhotonContainer>> photonLinks;
      for(const xAOD::Photon* photon : *photonsIn){
        ElementLink<xAOD::PhotonContainer> phLink;
        phLink.setElement(photon);
        phLink.setStorableObject(*photonsIn.ptr(), true);
        photonLinks.push_back(phLink);
      }
      dec_photonLinks(*vertex) = photonLinks;
      dec_nPhotons(*vertex) = countValid(photonLinks);
      
      // multi-photon link 
      std::vector<ElementLink<xAOD::CompositeParticleContainer>> mpLinks;
      if (haveMP && vertex == bestVtxForMP)
      {
        mpLinks.push_back(mpLink);
      }
      dec_multiPhotonLinks(*vertex) = mpLinks;
 
      // for jets, use prefilled map
      dec_jetLinks(*vertex) = jetsInVertex[vertex];
      dec_nJets(*vertex) = countValid(jetsInVertex[vertex]);
 
      std::vector<ElementLink<xAOD::MuonContainer>> muonLinks;
      for(const xAOD::Muon* muon : *muonsIn){
        const auto *tp = muon->trackParticle(xAOD::Muon::InnerDetectorTrackParticle);
        if(!tp) continue;
        try{
          auto muonVertex = m_trkVtxAssociationTool->getUniqueMatchVertexLink(*tp, *vertices);
          if(!muonVertex) continue;
          ElementLink<xAOD::MuonContainer> muonLink;
          if(*muonVertex == vertex){
            muonLink.setElement(muon);
            muonLink.setStorableObject(*muonsIn.ptr(), true);
            muonLinks.push_back(muonLink);
          }
        }catch(...) {
          ATH_MSG_DEBUG("Skipping muon as the track is not associated to any PV ");
          ATH_MSG_DEBUG("Muon pT, eta = " << muon->pt() << " " << muon->eta());
        }
 
      }
      dec_muonLinks(*vertex) = muonLinks;
      dec_nMuons(*vertex) = countValid(muonLinks);
 
      // Finally, decorate the vertices with the GNN score
      float score_phsvertex = m_gnnTool->decorate(*vertex);
 
      dec_gnnScore(*vertex) = score_phsvertex;
    }
 
    return StatusCode::SUCCESS;
  }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::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 & AthCommonAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 89 of file AthCommonAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

◆ filterPassed()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

◆ 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 453 of file AthHistogramming.cxx.

{
  // Build a 32 bit hash out of the name
  const hash_t graphHash = this->hash(graphName);
 
  // See if this entry exists in the map
  GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
  if ( it == m_graphMap.end() ) // It doesn't exist!
    { // Let's see into the THistSvc if somebody else has registered the TGraph...
 
      // Need to copy the strings as we will massage them from here on
      std::string graphNameCopy = graphName;
      std::string tDirCopy      = tDir;
      std::string streamCopy    = stream;
 
      // Massage the final string to book things
      std::string bookingString("");
      this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, false);
 
      TGraph* graphPointer(nullptr);
      if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
        {
          // Massage the final string to book things
          std::string bookingString("");
          this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, true );
 
          if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
            {
              m_msg << MSG::WARNING
                    << "Problem retrieving the TGraph with name (including pre- and post-fixes) "
                    << m_histNamePrefix + graphNameCopy + m_histNamePostfix
                    << " or with name " << graphNameCopy
                    << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
                    << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
              return nullptr;
            }
          // If we get to here, we actually found the TGraph in the THistSvc.
          // So let's add it to the local cache map and return its pointer
          m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
          return graphPointer;
        }
      // If we get to here, we actually found the TGraph in the THistSvc.
      // So let's add it to the local cache map and return its pointer
      m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
      return graphPointer;
    }
 
 
  // Return the pointer to the TGraph that we got from the local cache map
  return it->second;
}

◆ 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 430 of file AthHistogramming.h.

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

◆ 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 164 of file AthHistogramming.cxx.

{
  std::string histNameCopy = histName;
  std::string tDirCopy     = tDir;
  std::string streamCopy   = stream;
 
  std::string bookingString;
  this->buildBookingString(bookingString, histNameCopy, tDirCopy, streamCopy, false);
 
  const hash_t histHash = this->hash(histNameCopy);
 
  HistMap_t::const_iterator it = m_histMap.find(histHash);
  if (it == m_histMap.end()) {
    TH1* histPointer(nullptr);
 
    if (!histSvc()->getHist(bookingString, histPointer).isSuccess()) {
      std::string prefixedHistNameCopy = histName;
      std::string prefixedTDirCopy     = tDir;
      std::string prefixedStreamCopy   = stream;
 
      std::string prefixedBookingString;
      this->buildBookingString(prefixedBookingString,
                               prefixedHistNameCopy,
                               prefixedTDirCopy,
                               prefixedStreamCopy,
                               true);
 
      if (!histSvc()->getHist(prefixedBookingString, histPointer).isSuccess()) {
        m_msg << MSG::WARNING
              << "Problem retrieving the histogram with name (including pre- and post-fixes) "
              << m_histNamePrefix + histNameCopy + m_histNamePostfix
              << " or with name " << histNameCopy
              << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
              << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
        return nullptr;
      }
    }
 
    m_histMap.insert(m_histMap.end(), std::pair<const hash_t, TH1*>(histHash, histPointer));
    return histPointer;
  }
 
  return it->second;
}

◆ 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 371 of file AthHistogramming.h.

{
  // Get the TH1 pointer
  TH1* th1Pointer = this->hist(histName, tDir, stream);
  if ( !th1Pointer )
    {
      m_msg << MSG::ERROR
            << "Cannot get a 2-d histogram with name " << histName
            << "... will probably seg-fault!" << endmsg;
      return NULL;
    }
  // If the TH1 pointer is valid, simply return the dynamic_cast
  return dynamic_cast<TH2*>( th1Pointer );
}

◆ 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 388 of file AthHistogramming.h.

{
  // Get the TH1 pointer
  TH1* th1Pointer = this->hist(histName, tDir, stream);
  if ( !th1Pointer )
    {
      m_msg << MSG::ERROR
            << "Cannot get a 3-d histogram with name " << histName
            << "... will probably seg-fault!" << endmsg;
      return NULL;
    }
  // If the TH1 pointer is valid, simply return the dynamic_cast
  return dynamic_cast<TH3*>( th1Pointer );
}

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

{
  return m_histSvc;
}

◆ initialize()

StatusCode VertexDecoratorAlg::initialize ( )

override

Definition at line 25 of file VertexDecoratorAlg.cxx.

  {
    ATH_CHECK(m_vertexInKey.initialize());
    ATH_CHECK(m_eventInKey.initialize());
    ATH_CHECK(m_photonsInKey.initialize());
    ATH_CHECK(m_electronsInKey.initialize());
    ATH_CHECK(m_muonsInKey.initialize());
    ATH_CHECK(m_jetsInKey.initialize());
    ATH_CHECK(m_multiPhotonsOutKey.initialize());
 
    const std::string baseName = m_vertexInKey.key();
 
    // WriteHandleKeys
    m_photonLinksKey = baseName + ".photonLinks";
    m_jetLinksKey = baseName + ".jetLinks";
    m_electronLinksKey = baseName + ".electronLinks";
    m_muonLinksKey = baseName + ".muonLinks";
    m_mDecor_ntrk = baseName + "." + m_mDecor_ntrk.key();
    m_mDecor_sumPt = baseName + "." + m_mDecor_sumPt.key();
    m_mDecor_chi2Over_ndf = baseName + "." + m_mDecor_chi2Over_ndf.key();
    m_mDecor_z_asym = baseName + "." + m_mDecor_z_asym.key();
    m_mDecor_weighted_z_asym = baseName + "." + m_mDecor_weighted_z_asym.key();
    m_mDecor_weighted_z_kurt = baseName + "." + m_mDecor_weighted_z_kurt.key();
    m_mDecor_z_skew = baseName + "." + m_mDecor_z_skew.key();
    m_mDecor_photon_deltaz = baseName + "." + m_mDecor_photon_deltaz.key();
    m_mDecor_photon_deltaPhi = baseName + "." + m_mDecor_photon_deltaPhi.key();
    m_mDecor_actualInterPerXing = baseName + "." + m_mDecor_actualInterPerXing.key();
    m_multiPhotonLinksKey = baseName + ".multiPhotonLinks";
    m_mDecor_gnnScore = baseName + "." + m_mDecor_gnnScore.key();
    m_mDecor_nElectrons = baseName + "." + m_mDecor_nElectrons.key();
    m_mDecor_nMuons = baseName + "." + m_mDecor_nMuons.key();
    m_mDecor_nPhotons = baseName + "." + m_mDecor_nPhotons.key();
    m_mDecor_nJets = baseName + "." + m_mDecor_nJets.key();
 
    ATH_CHECK(m_multiPhotonLinksKey.initialize());
    ATH_CHECK(m_photonLinksKey.initialize());
    ATH_CHECK(m_jetLinksKey.initialize());
    ATH_CHECK(m_electronLinksKey.initialize());
    ATH_CHECK(m_muonLinksKey.initialize());
    ATH_CHECK(m_mDecor_ntrk.initialize());
    ATH_CHECK(m_mDecor_sumPt.initialize());
    ATH_CHECK(m_mDecor_chi2Over_ndf.initialize());
    ATH_CHECK(m_mDecor_z_asym.initialize());
    ATH_CHECK(m_mDecor_weighted_z_asym.initialize());
    ATH_CHECK(m_mDecor_weighted_z_kurt.initialize());
    ATH_CHECK(m_mDecor_z_skew.initialize());
    ATH_CHECK(m_mDecor_photon_deltaz.initialize());
    ATH_CHECK(m_mDecor_photon_deltaPhi.initialize());
    ATH_CHECK(m_mDecor_actualInterPerXing.initialize());
    ATH_CHECK(m_mDecor_gnnScore.initialize());
    ATH_CHECK(m_mDecor_nElectrons.initialize());
    ATH_CHECK(m_mDecor_nMuons.initialize());
    ATH_CHECK(m_mDecor_nJets.initialize());
    ATH_CHECK(m_mDecor_nPhotons.initialize());
 
    // ReadHandleKeys
    m_deltaZKey = baseName + ".deltaZ";
    m_deltaPhiKey = baseName + ".deltaPhi";
    ATH_CHECK(m_deltaZKey.initialize());
    ATH_CHECK(m_deltaPhiKey.initialize());
 
    // additional ReadHandleKeys to declare dependencies to the scheduler
    std::string photonBaseName = m_photonsInKey.key();
 
    m_caloPointingZKey = photonBaseName + ".caloPointingZ";
    m_zCommonKey = photonBaseName + ".zCommon";
    m_zCommonErrorKey = photonBaseName + ".zCommonError";
 
    ATH_CHECK(m_caloPointingZKey.initialize());
    ATH_CHECK(m_zCommonKey.initialize());
    ATH_CHECK(m_zCommonErrorKey.initialize());
 
    // Tools
    ATH_CHECK(m_gnnTool.retrieve());
    ATH_CHECK(m_trkVtxAssociationTool.retrieve());
 
    return StatusCode::SUCCESS;
  }

◆ inputHandles()

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

◆ isClonable()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

inlineoverridevirtualinherited

Specify if the algorithm is clonable.

Only relevant for non-reentrant algorithms. Actual number of clones needs to be set via the "Cardinality" property.

Reimplemented in AFP_DigiTop, AlgB, AlgT, BCM_Digitization, CscDigitBuilder, CscDigitToCscRDO, G4AtlasAlg, G4RunAlg, HGTD_Digitization, HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ISF::SimKernelMT, ITk::StripDigitization, ITkPixelCablingAlg, ITkStripCablingAlg, LArHitEMapMaker, LArTTL1Maker, LUCID_DigiTop, LVL1::L1TopoSimulation, MergeCalibHits, MergeGenericMuonSimHitColl, MergeHijingPars, MergeMcEventCollection, MergeTrackRecordCollection, MergeTruthJets, MergeTruthParticles, MuonDigitizer, PileUpMTAlg, PixelDigitization, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_Digitization, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, SGInputLoader, Simulation::BeamEffectsAlg, TileHitVecToCnt, TileMuonFitter, TilePulseForTileMuonReceiver, TileRawChannelMaker, TRTDigitization, and ZDC_DigiTop.

Definition at line 68 of file AthCommonAlgorithm.h.

                                           {
    return true;
  }

◆ isReEntrant()

virtual bool AthAlgorithm::isReEntrant ( ) const

inlinefinaloverrideprotectedvirtualinherited

Legacy algorithms are not thread-safe.

Definition at line 47 of file AthAlgorithm.h.

47{ return false; }

◆ msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

◆ 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 551 of file AthHistogramming.cxx.

{
  size_t pos = 0;
  while((pos = str.find(oldStr, pos)) != std::string::npos)
    {
      str.replace(pos, oldStr.length(), newStr);
      pos += newStr.length();
    }
}

◆ outputHandles()

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

◆ 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< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

◆ setFilterPassed()

virtual void AthCommonAlgorithm< Gaudi::Algorithm >::setFilterPassed	(	bool	state,
		const EventContext &	ctx ) const

inlinevirtualinherited

Set filter decision:

Reimplemented in AthFilterAlgorithm.

Definition at line 99 of file AthCommonAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

◆ sysExecute()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

◆ 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 AthCommonAlgorithm< Gaudi::Algorithm >.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 75 of file AthHistogramAlgorithm.cxx.

{
  // ---- stolen from GaudiKernel/Algorithm::sysInitialize -------
  // Bypass the initialization if the algorithm
  // has already been initialized.
  if ( Gaudi::StateMachine::INITIALIZED <= FSMState() ) return StatusCode::SUCCESS;
 
  // Set the Algorithm's properties
  bindPropertiesTo( serviceLocator()->getOptsSvc() );
 
  // Bypass the initialization if the algorithm is disabled.
  // Need to do this after bindPropertiesTo.
  if ( !isEnabled( ) ) return StatusCode::SUCCESS;
  
  // ---- stolen from GaudiKernel/Algorithm::sysInitialize ------- END ---
 
 
  // Get the THistSvc
  ATH_CHECK ( histSvc().retrieve() );
 
  // Configure the underlying AthHistogramming helper
  ATH_CHECK ( AthHistogramming::configAthHistogramming( m_histSvc,
                                                        m_prefix,          m_rootDir, 
                                                        m_histNamePrefix,  m_histNamePostfix,
                                                        m_histTitlePrefix, m_histTitlePostfix ) );
  
  // Print some setup information into the log file
  ATH_MSG_DEBUG ("Initializing " << name() << "...");
  ATH_MSG_DEBUG (" using THistService     = " << m_histSvc );
  ATH_MSG_DEBUG (" using RootStreamName   = " << m_prefix );
  ATH_MSG_DEBUG (" using RootDirName      = " << m_rootDir );
  ATH_MSG_DEBUG (" using HistNamePrefix   = " << m_histNamePrefix );
  ATH_MSG_DEBUG (" using HistNamePostfix  = " << m_histNamePostfix );
  ATH_MSG_DEBUG (" using HistTitlePrefix  = " << m_histTitlePrefix );
  ATH_MSG_DEBUG (" using HistTitlePostfix = " << m_histTitlePostfix );
 
 
  // re-direct to base class...
  return AthAlgorithm::sysInitialize();
}

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::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 339 of file AthHistogramming.cxx.

{
  // Build a 32 bit hash out of the name
  const hash_t treeHash = this->hash(treeName);
 
  // See if this entry exists in the map
  TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
  if ( it == m_treeMap.end() ) // It doesn't exist!
    { // Let's see into the THistSvc if somebody else has registered the TTree...
 
      // Need to copy the strings as we will massage them from here on
      std::string treeNameCopy = treeName;
      std::string tDirCopy     = tDir;
      std::string streamCopy   = stream;
 
      // Massage the final string to book things
      std::string bookingString("");
      this->buildBookingString( bookingString, treeNameCopy, tDirCopy, streamCopy );
 
      TTree* treePointer(NULL);
      if ( !((histSvc()->getTree(bookingString, treePointer)).isSuccess()) )
        {
          m_msg << MSG::WARNING
                << "Problem retrieving the TTree with name " << treeNameCopy
                << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
                << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
          return NULL;
        }
      // If we get to here, we actually found the TTree in the THistSvc.
      // So let's add it to the local cache map and return its pointer
      m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
      return treePointer;
    }
 
  // Return the pointer to the TTree that we got from the local cache map
  return it->second;
}

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

◆ m_caloPointingZKey

SG::ReadDecorHandleKey<xAOD::PhotonContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_caloPointingZKey

private

Initial value:

{

this, "caloPointingZKey", "", "" }

Definition at line 117 of file VertexDecoratorAlg.h.

117 {

118 this, "caloPointingZKey", "", "" };

◆ m_deltaPhiKey

SG::ReadDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_deltaPhiKey

private

Initial value:

{

this, "deltaPhiKey", "", "" }

Definition at line 112 of file VertexDecoratorAlg.h.

112 {

113 this, "deltaPhiKey", "", "" };

◆ m_deltaZKey

SG::ReadDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_deltaZKey

private

Initial value:

{

this, "deltaZKey", "", "" }

Definition at line 110 of file VertexDecoratorAlg.h.

110 {

111 this, "deltaZKey", "", "" };

◆ m_detStore

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

privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

◆ m_effMap

EffMap_t AthHistogramming::m_effMap

privateinherited

The map of histogram names to their pointers.

Definition at line 210 of file AthHistogramming.h.

◆ m_electronLinksKey

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_electronLinksKey

private

Initial value:

{

this, "electronLinks", "", "" }

Definition at line 87 of file VertexDecoratorAlg.h.

87 {

88 this, "electronLinks", "", "" };

◆ m_electronsInKey

SG::ReadHandleKey<xAOD::ElectronContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_electronsInKey

private

Initial value:

{

this, "electronsIn", "Electrons", "containerName to read"}

Definition at line 53 of file VertexDecoratorAlg.h.

53 {

54 this, "electronsIn", "Electrons", "containerName to read"};

◆ m_eventInKey

SG::ReadHandleKey<xAOD::EventInfo> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_eventInKey

private

Initial value:

{

this, "eventIn", "EventInfo", "containerName to read"}

Definition at line 51 of file VertexDecoratorAlg.h.

51 {

52 this, "eventIn", "EventInfo", "containerName to read"};

◆ m_evtStore

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

privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

◆ m_gnnTool

ToolHandle<GNNTool> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_gnnTool

private

Initial value:

{
        this, "gnnTool",
        "InDetGNNHardScatterSelection/GNNTool",
        "GNN tool to use"}

Definition at line 130 of file VertexDecoratorAlg.h.

                                 {
        this, "gnnTool",
        "InDetGNNHardScatterSelection/GNNTool",
        "GNN tool to use"};

◆ m_graphMap

GraphMap_t AthHistogramming::m_graphMap

privateinherited

The map of TGraph names to their pointers.

Definition at line 224 of file AthHistogramming.h.

◆ m_histMap

HistMap_t AthHistogramming::m_histMap

privateinherited

The map of histogram names to their pointers.

Definition at line 203 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_jetLinksKey

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_jetLinksKey

private

Initial value:

{

this, "jetLinks", "", "" }

Definition at line 85 of file VertexDecoratorAlg.h.

85 {

86 this, "jetLinks", "", "" };

◆ m_jetsInKey

SG::ReadHandleKey<xAOD::JetContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_jetsInKey

private

Initial value:

{

this, "jetsIn", "AntiKt4EMTopoJets", "containerName to read"}

Definition at line 57 of file VertexDecoratorAlg.h.

57 {

58 this, "jetsIn", "AntiKt4EMTopoJets", "containerName to read"};

◆ m_mDecor_actualInterPerXing

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_actualInterPerXing

private

Initial value:

{

this, "decor_actualInterPerXing", "actualIntPerXing", ""}

Definition at line 80 of file VertexDecoratorAlg.h.

80 {

81 this, "decor_actualInterPerXing", "actualIntPerXing", ""};

◆ m_mDecor_chi2Over_ndf

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_chi2Over_ndf

private

Initial value:

{

this, "decor_chi2Over_ndfKey", "chi2Over_ndf", ""}

Definition at line 66 of file VertexDecoratorAlg.h.

66 {

67 this, "decor_chi2Over_ndfKey", "chi2Over_ndf", ""};

◆ m_mDecor_gnnScore

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_gnnScore

private

Initial value:

{

this, "GNNScoreKey", "HSGN2_phsvertex", "Score from HSGNN"}

Definition at line 92 of file VertexDecoratorAlg.h.

92 {

93 this, "GNNScoreKey", "HSGN2_phsvertex", "Score from HSGNN"};

◆ m_mDecor_nElectrons

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_nElectrons

private

Initial value:

{

this, "decor_nElectrons", "nElectrons", ""}

Definition at line 95 of file VertexDecoratorAlg.h.

95 {

96 this, "decor_nElectrons", "nElectrons", ""};

◆ m_mDecor_nJets

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_nJets

private

Initial value:

{

this, "decor_nJets", "nJets", ""}

Definition at line 99 of file VertexDecoratorAlg.h.

99 {

100 this, "decor_nJets", "nJets", ""};

◆ m_mDecor_nMuons

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_nMuons

private

Initial value:

{

this, "decor_nMuons", "nMuons", ""}

Definition at line 97 of file VertexDecoratorAlg.h.

97 {

98 this, "decor_nMuons", "nMuons", ""};

◆ m_mDecor_nPhotons

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_nPhotons

private

Initial value:

{

this, "decor_nPhotons", "nPhotons", ""}

Definition at line 101 of file VertexDecoratorAlg.h.

101 {

102 this, "decor_nPhotons", "nPhotons", ""};

◆ m_mDecor_ntrk

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_ntrk

private

Initial value:

{

this, "decor_ntrkKey", "ntrk", ""}

Definition at line 62 of file VertexDecoratorAlg.h.

62 {

63 this, "decor_ntrkKey", "ntrk", ""};

◆ m_mDecor_photon_deltaPhi

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_photon_deltaPhi

private

Initial value:

{

this, "decor_photon_deltaPhi", "photon_deltaPhi", ""}

Definition at line 78 of file VertexDecoratorAlg.h.

78 {

79 this, "decor_photon_deltaPhi", "photon_deltaPhi", ""};

◆ m_mDecor_photon_deltaz

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_photon_deltaz

private

Initial value:

{

this, "decor_photon_deltaz", "photon_deltaz", ""}

Definition at line 76 of file VertexDecoratorAlg.h.

76 {

77 this, "decor_photon_deltaz", "photon_deltaz", ""};

◆ m_mDecor_sumPt

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_sumPt

private

Initial value:

{

this, "decor_sumPtKey", "sumPt", ""}

Definition at line 64 of file VertexDecoratorAlg.h.

64 {

65 this, "decor_sumPtKey", "sumPt", ""};

◆ m_mDecor_weighted_z_asym

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_weighted_z_asym

private

Initial value:

{

this, "decor_weighted_z_asym", "weighted_z_asymmetry", ""}

Definition at line 70 of file VertexDecoratorAlg.h.

70 {

71 this, "decor_weighted_z_asym", "weighted_z_asymmetry", ""};

◆ m_mDecor_weighted_z_kurt

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_weighted_z_kurt

private

Initial value:

{

this, "decor_weighted_z_kurt", "z_kurtosis", ""}

Definition at line 72 of file VertexDecoratorAlg.h.

72 {

73 this, "decor_weighted_z_kurt", "z_kurtosis", ""};

◆ m_mDecor_z_asym

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_z_asym

private

Initial value:

{

this, "decor_z_asymmetryKey", "z_asymmetry", ""}

Definition at line 68 of file VertexDecoratorAlg.h.

68 {

69 this, "decor_z_asymmetryKey", "z_asymmetry", ""};

◆ m_mDecor_z_skew

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_mDecor_z_skew

private

Initial value:

{

this, "decor_z_skew", "z_skewness", ""}

Definition at line 74 of file VertexDecoratorAlg.h.

74 {

75 this, "decor_z_skew", "z_skewness", ""};

◆ m_msg

MsgStream AthHistogramming::m_msg

privateinherited

Cached Message Stream.

Definition at line 250 of file AthHistogramming.h.

◆ m_multiPhotonLinksKey

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_multiPhotonLinksKey

private

Initial value:

{

this, "multiPhotonLinks", "multiPhotonLinks", "Per-vertex link(s) to multi-photon node"}

Definition at line 106 of file VertexDecoratorAlg.h.

106 {

107 this, "multiPhotonLinks", "multiPhotonLinks", "Per-vertex link(s) to multi-photon node"};

◆ m_multiPhotonsOutKey

SG::WriteHandleKey<xAOD::CompositeParticleContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_multiPhotonsOutKey

private

Initial value:

{

this, "multiPhotonsOut", "MultiPhotons", "Output container for multi-photon node"}

Definition at line 104 of file VertexDecoratorAlg.h.

104 {

105 this, "multiPhotonsOut", "MultiPhotons", "Output container for multi-photon node"};

◆ m_muonLinksKey

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_muonLinksKey

private

Initial value:

{

this, "muonLinks", "", "" }

Definition at line 89 of file VertexDecoratorAlg.h.

89 {

90 this, "muonLinks", "", "" };

◆ m_muonsInKey

SG::ReadHandleKey<xAOD::MuonContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_muonsInKey

private

Initial value:

{

this, "muonsIn", "Muons", "containerName to read"}

Definition at line 55 of file VertexDecoratorAlg.h.

55 {

56 this, "muonsIn", "Muons", "containerName to read"};

◆ m_name

std::string AthHistogramming::m_name

privateinherited

Instance name.

Definition at line 247 of file AthHistogramming.h.

◆ m_photonLinksKey

SG::WriteDecorHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_photonLinksKey

private

Initial value:

{

this, "photonLinks", "", "" }

Definition at line 83 of file VertexDecoratorAlg.h.

83 {

84 this, "photonLinks", "", "" };

◆ m_photonsInKey

SG::ReadHandleKey<xAOD::PhotonContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_photonsInKey

private

Initial value:

{

this, "photonsIn", "Photons", "containerName to read"}

Definition at line 59 of file VertexDecoratorAlg.h.

59 {

60 this, "photonsIn", "Photons", "containerName to read"};

◆ m_prefix

std::string AthHistogramAlgorithm::m_prefix

privateinherited

Name of the ROOT output stream (file).

Definition at line 88 of file AthHistogramAlgorithm.h.

◆ m_rootDir

std::string AthHistogramAlgorithm::m_rootDir

privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

◆ m_streamName

std::string AthHistogramming::m_streamName

privateinherited

Name of the ROOT output stream (file).

Definition at line 228 of file AthHistogramming.h.

◆ m_treeMap

TreeMap_t AthHistogramming::m_treeMap

privateinherited

The map of TTree names to their pointers.

Definition at line 217 of file AthHistogramming.h.

◆ m_trkVtxAssociationTool

ToolHandle<CP::TrackVertexAssociationTool> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_trkVtxAssociationTool

private

Initial value:

{
        this, "TrackVertexAssociationTool", 
        "CP::TrackVertexAssociationTool/TrackVertexAssociationTool", 
        "Track vertex association tool to use"}

Definition at line 125 of file VertexDecoratorAlg.h.

                                                                    {
        this, "TrackVertexAssociationTool", 
        "CP::TrackVertexAssociationTool/TrackVertexAssociationTool", 
        "Track vertex association tool to use"};

◆ m_varHandleArraysDeclared

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

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vertexInKey

SG::ReadHandleKey<xAOD::VertexContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_vertexInKey

private

Initial value:

{

this, "vertexIn", "PrimaryVertices", "containerName to read"}

Definition at line 49 of file VertexDecoratorAlg.h.

49 {

50 this, "vertexIn", "PrimaryVertices", "containerName to read"};

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_zCommonErrorKey

SG::ReadDecorHandleKey<xAOD::PhotonContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_zCommonErrorKey

private

Initial value:

{

this, "zCommonErrorKey", "", "" }

Definition at line 121 of file VertexDecoratorAlg.h.

121 {

122 this, "zCommonErrorKey", "", "" };

◆ m_zCommonKey

SG::ReadDecorHandleKey<xAOD::PhotonContainer> InDetGNNHardScatterSelection::VertexDecoratorAlg::m_zCommonKey

private

Initial value:

{

this, "zCommonKey", "", "" }

Definition at line 119 of file VertexDecoratorAlg.h.

119 {

120 this, "zCommonKey", "", "" };

The documentation for this class was generated from the following files:

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ EffMap_t

◆ GraphMap_t

◆ hash_t

◆ HistMap_t

◆ StoreGateSvc_t

◆ TreeMap_t

Constructor & Destructor Documentation

◆ VertexDecoratorAlg()

Member Function Documentation

◆ book() [1/8]

◆ book() [2/8]

◆ book() [3/8]

◆ book() [4/8]

◆ book() [5/8]

◆ book() [6/8]

◆ book() [7/8]

◆ book() [8/8]

◆ bookGetPointer() [1/8]

◆ bookGetPointer() [2/8]

◆ bookGetPointer() [3/8]

◆ bookGetPointer() [4/8]

◆ bookGetPointer() [5/8]

◆ bookGetPointer() [6/8]

◆ bookGetPointer() [7/8]

◆ bookGetPointer() [8/8]

◆ buildBookingString()

◆ configAthHistogramming()

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ efficiency()

◆ evtStore()

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ filterPassed()

◆ graph()

◆ hash()

◆ hist()

◆ hist2d()

◆ hist3d()

◆ histSvc()

◆ initialize()

◆ inputHandles()

◆ isClonable()

◆ isReEntrant()

◆ msg()

◆ msgLvl()

◆ myReplace()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ setFilterPassed()

◆ sysExecute()

◆ sysInitialize()

◆ sysStart()

◆ tree()

◆ updateVHKA()

Member Data Documentation

◆ m_caloPointingZKey

◆ m_deltaPhiKey

◆ m_deltaZKey

◆ m_detStore

◆ m_effMap

◆ m_electronLinksKey

◆ m_electronsInKey

◆ m_eventInKey

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_gnnTool

◆ m_graphMap

◆ m_histMap

◆ m_histNamePostfix