Write out the selected events in xAOD format. More...

#include <EventSaverxAOD.h>

Inheritance diagram for top::EventSaverxAOD:

Collaboration diagram for top::EventSaverxAOD:

Public Member Functions
	EventSaverxAOD ()
	Kept simple, set the output container prefix (so they have a different name to the input containers). More...

virtual	~EventSaverxAOD ()

virtual void	initialize (std::shared_ptr< top::TopConfig > config, TFile *file, const std::vector< std::string > &extraBranches)
	This does some configuration of the output file. More...

virtual StatusCode	initialize ()
	Dummy implementation of the initialisation function. More...

virtual void	saveEvent (const top::Event &event)
	Save an event. More...

virtual void	finalize ()
	xAOD needs to write some more stuff to the file at the end of a job More...

virtual void	execute ()
	Method to execute any calculations common for reco/particle/parton level. More...

virtual void	saveEventToxAOD ()

virtual void	saveTruthEvent ()
	@breif Run for every event More...

virtual void	saveParticleLevelEvent (const top::ParticleLevelEvent &)
	Function to save the particle level event. More...

virtual void	saveUpgradeEvent (const top::ParticleLevelEvent &)
	Function to save the upgrade event. More...

	ClassDef (top::EventSaverBase, 0)

virtual void	print () const
	Print the state of the tool. More...

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)

void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
	ClassDef (top::EventSaverxAOD, 0)

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

Private Attributes
const std::string	m_prefix
	Name to prepend to the output containers. More...

std::shared_ptr< top::TopConfig >	m_config
	We need access to the configuration file to get the container names. More...

TFile *	m_outputFile
	We need to hold on to the output file. More...

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

Write out the selected events in xAOD format.

At some point this will be the future. But probably not right now. Still too many missing features.

Definition at line 20 of file EventSaverxAOD.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ EventSaverxAOD()

top::EventSaverxAOD::EventSaverxAOD ( )

Kept simple, set the output container prefix (so they have a different name to the input containers).

Definition at line 20 of file EventSaverxAOD.cxx.

                                  :
     asg::AsgTool("top::EventSaverxAOD"),
     m_prefix("CorrectedSelected"),
     m_config(nullptr),
     m_outputFile(nullptr) {
   }

◆ ~EventSaverxAOD()

top::EventSaverxAOD::~EventSaverxAOD ( )

virtual

Definition at line 27 of file EventSaverxAOD.cxx.

27 {

28 }

Member Function Documentation

◆ ClassDef() [1/2]

top::EventSaverBase::ClassDef	(	top::EventSaverBase	,
		0
	)

inherited

◆ ClassDef() [2/2]

top::EventSaverxAOD::ClassDef	(	top::EventSaverxAOD	,
		0
	)

private

◆ declareGaudiProperty() [1/4]

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

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	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());
  
   }

◆ declareGaudiProperty() [3/4]

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

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::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.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

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

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

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

◆ declareProperty() [5/6]

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

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

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

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > & 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< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

virtual void top::EventSaverBase::execute ( )

inlinevirtualinherited

Method to execute any calculations common for reco/particle/parton level.

Reimplemented in top::EventSaverFlatNtuple.

Definition at line 48 of file EventSaverBase.h.

48 {return;}

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ finalize()

void top::EventSaverxAOD::finalize ( )

virtual

xAOD needs to write some more stuff to the file at the end of a job

Reimplemented from top::EventSaverBase.

Definition at line 208 of file EventSaverxAOD.cxx.

                                 {
     m_outputFile->Write();
     top::check(evtStore()->event()->finishWritingTo(
                  m_outputFile), "EventSaverxAOD::finalise trying to call finishWritingTo");
   }

◆ getKey()

SG::sgkey_t asg::AsgTool::getKey ( const void * ptr ) const

inherited

Get the (hashed) key of an object that is in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::TEvent both provide ways for getting the SG::sgkey_t key for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also: asg::AsgTool::getName

Parameters

ptr	The bare pointer to the object that the event store should know about

Returns: The hashed key of the object in the store. If not found, an invalid (zero) key.

Definition at line 119 of file AsgTool.cxx.

                                                     {
  
 #ifdef XAOD_STANDALONE
       // In case we use @c xAOD::TEvent, we have a direct function call
       // for this.
       return evtStore()->event()->getKey( ptr );
 #else
       const SG::DataProxy* proxy = evtStore()->proxy( ptr );
       return ( proxy == nullptr ? 0 : proxy->sgkey() );
 #endif // XAOD_STANDALONE
    }

◆ getName()

const std::string & asg::AsgTool::getName ( const void * ptr ) const

inherited

Get the name of an object that is / should be in the event store.

This is a bit of a special one. StoreGateSvc and xAOD::TEvent both provide ways for getting the std::string name for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.

In order to allow tools to efficiently perform this operation, they can use this helper function.

See also: asg::AsgTool::getKey

Parameters

ptr	The bare pointer to the object that the event store should know about

Returns: The string name of the object in the store. If not found, an empty string.

Definition at line 106 of file AsgTool.cxx.

                                                             {
  
 #ifdef XAOD_STANDALONE
       // In case we use @c xAOD::TEvent, we have a direct function call
       // for this.
       return evtStore()->event()->getName( ptr );
 #else
       const SG::DataProxy* proxy = evtStore()->proxy( ptr );
       static const std::string dummy = "";
       return ( proxy == nullptr ? dummy : proxy->name() );
 #endif // XAOD_STANDALONE
    }

◆ getProperty()

template<class T >

const T* asg::AsgTool::getProperty ( const std::string & name ) const

inherited

Get one of the tool's properties.

◆ initialize() [1/2]

virtual StatusCode top::EventSaverxAOD::initialize ( )

inlinevirtual

Dummy implementation of the initialisation function.

It's here to allow the dual-use tools to skip defining an initialisation function. Since many are doing so...

Reimplemented from asg::AsgTool.

Definition at line 51 of file EventSaverxAOD.h.

51 {return StatusCode::SUCCESS;}

◆ initialize() [2/2]

void top::EventSaverxAOD::initialize	(	std::shared_ptr< top::TopConfig >	config,
		TFile *	file,
		const std::vector< std::string > &	extraBranches
	)

virtual

This does some configuration of the output file.

It also keeps track of which branches should be saved.

You surely don't need all the branches, but if you do remove the corresponding line.

We can't do much in the constructor, like I would like, because we're letting root make this. So instead we have a function here.

Parameters

config	The top::TopConfig object that contains lots of settings and stuff.
file	The output TFile that you want to save things to.
extraBranches	List of extra branches that (currently) are attached to EventInfo. e.g. the branch that tells you if an event passed a certain selection. If you add it here, it'll be written to the output file.

Implements top::EventSaverBase.

Definition at line 30 of file EventSaverxAOD.cxx.

                                                                              {
     m_config = config;
     m_outputFile = file;
     top::check(evtStore()->event()->writeTo(m_outputFile), "EventSaverxAOD trying to call writeTo on output file");
  
     //EventInfo
     std::string eventInfoList = "runNumber.eventNumber.eventTypeBitmask.averageInteractionsPerCrossing";
     if (config->isMC()) eventInfoList += ".mcChannelNumber.mcEventWeights.PileupWeight";
  
     //Add extra branches - e.g. if it passed the ee / mumu / emu selection per event
     for (const auto& branchName : extraBranches)
       eventInfoList += "." + branchName;
  
     evtStore()->event()->setAuxItemList("EventInfoAux.", eventInfoList);
  
     //Electrons
     evtStore()->event()->setAuxItemList(m_prefix + m_config->sgKeyElectrons() + "Aux.", "pt.eta.phi.m.charge");
  
     //Muons
     evtStore()->event()->setAuxItemList(
       m_prefix + m_config->sgKeyMuons() + "Aux.", "pt.eta.phi.m.charge.muonType.combinedTrackParticleLink");
  
     //Once mu has a built-in charge we can get rid of this
     //Muon charge is a property of the associated track - e.g. d0.z0.phi.theta.qOverP
     evtStore()->event()->setAuxItemList("InDetTrackParticlesAux.", "qOverP");
     evtStore()->event()->setAuxItemList("ExtrapolatedMuonTrackParticlesAux.", "qOverP");
     evtStore()->event()->setAuxItemList("CombinedMuonTrackParticlesAux.", "qOverP");
  
     //Jets
     evtStore()->event()->setAuxItemList(m_prefix + m_config->sgKeyJets() + "Aux.", "pt.eta.phi.m.btaggingLink");
  
   }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msg_level_name()

const std::string & asg::AsgTool::msg_level_name ( ) const

inherited

A deprecated function for getting the message level's name.

Instead of using this, weirdly named function, user code should get the string name of the current minimum message level (in case they really need it...), with:

MSG::name( msg().level() )

This function's name doesn't follow the ATLAS coding rules, and as such will be removed in the not too distant future.

Returns: The string name of the current minimum message level that's printed

Definition at line 101 of file AsgTool.cxx.

                                                   {
  
       return MSG::name( msg().level() );
    }

◆ msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ print()

void asg::AsgTool::print ( ) const

virtualinherited

Print the state of the tool.

Implements asg::IAsgTool.

Reimplemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, top::TopObjectSelection, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.

Definition at line 131 of file AsgTool.cxx.

                              {
  
       ATH_MSG_INFO( "AsgTool " << name() << " @ " << this );
       return;
    }

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ saveEvent()

void top::EventSaverxAOD::saveEvent ( const top::Event & event )

virtual

Save an event.

Why am I only writing out the nominal tree.

I thought the fill for the xaod would allow me to specify the tree name, but it seems I can only do that in the writeTo. Maybe their is a clever way around it, but I'm too tired of this.

Why not save shallow copies? This would be nice. But aren't shallow copies unsorted? We'd have to remember to sort them.

It looks like branch removal from things in containers (jets, leptons) is still not working. So the xAOD files are still huge.

Parameters

event	The top event for this particular systematic variation.
xaodEvent	The xAOD event object, so we can actually save.

Reimplemented from top::EventSaverBase.

Definition at line 64 of file EventSaverxAOD.cxx.

                                                       {
     //only save the nominal variation for now - speak to Attila about TEvent.fill not taking a tree name
     if (event.m_hashValue != m_config->nominalHashValue()) return;
  
     // record the event?
     if (m_config->saveOnlySelectedEvents() && !event.m_saveEvent) return;
  
     // Systematic Event
     const xAOD::SystematicEventContainer* allSystematics(nullptr);
     top::check(evtStore()->retrieve(allSystematics,
                                     m_config->sgKeyTopSystematicEvents()),
                "Failed to retrieve xAOD::SystematicEventContainer");
     //cppcheck-suppress uninitvar
     xAOD::SystematicEventContainer* allSystematics_output = new xAOD::SystematicEventContainer {};
     SG::IAuxStore* allSystematics_aux = evtStore()->event()->recordAux(m_config->sgKeyTopSystematicEvents() + "Aux.");
     allSystematics_output->setStore(allSystematics_aux);
     for (const auto *currentSystematic : *allSystematics) {
       //cppcheck-suppress uninitvar
       xAOD::SystematicEvent* out = new xAOD::SystematicEvent {};
       out->makePrivateStore(*currentSystematic);
       allSystematics_output->push_back(out);
     }
     top::check(evtStore()->event()->record(allSystematics_output,
                                            m_config->sgKeyTopSystematicEvents()),
                "Failed to record xAOD::SystematicEventContainer");
  
  
     // Top Parton History
     if (m_config->doTopPartonHistory()) {
       if (evtStore()->contains<xAOD::PartonHistoryContainer>(m_config->sgKeyTopPartonHistory())) {
         const xAOD::PartonHistoryContainer* partonHistory(nullptr);
         top::check(evtStore()->retrieve(partonHistory,
                                         m_config->sgKeyTopPartonHistory()), "Failed to retrieve Top Parton History");
         //cppcheck-suppress uninitvar
         xAOD::PartonHistoryContainer* partonHistory_output = new xAOD::PartonHistoryContainer {};
         SG::IAuxStore* partonHistory_aux = evtStore()->event()->recordAux(m_config->sgKeyTopPartonHistory() + "Aux.");
         partonHistory_output->setStore(partonHistory_aux);
         for (const auto *x : *partonHistory) {
           //cppcheck-suppress uninitvar
           xAOD::PartonHistory* out = new xAOD::PartonHistory {};
           out->makePrivateStore(*x);
           partonHistory_output->push_back(out);
         }
         top::check(evtStore()->event()->record(partonHistory_output,
                                                m_config->sgKeyTopPartonHistory()), "Failed to record Parton History");
       }
     }
  
  
  
     //Event Info
     top::check(evtStore()->event()->copy("EventInfo"), "copying event info failed");
  
     //Electrons
     if (m_config->useElectrons()) {
       xAOD::ElectronContainer* electrons_output = new xAOD::ElectronContainer();
       SG::IAuxStore* electrons_aux = evtStore()->event()->recordAux(m_prefix + m_config->sgKeyElectrons() + "Aux.");
       electrons_output->setStore(electrons_aux);
  
       for (const auto* const elPtr : event.m_electrons) {
         xAOD::Electron* outputElPtr = new xAOD::Electron();
         outputElPtr->makePrivateStore(*elPtr);
         electrons_output->push_back(outputElPtr);
       }
  
       top::check(evtStore()->event()->record(electrons_output,
                                              m_prefix + m_config->sgKeyElectrons()),
                  "xAOD::TEvent record " + m_prefix + m_config->sgKeyElectrons() + " failed");
     }
  
     //Muons
     if (m_config->useMuons()) {
       xAOD::MuonContainer* muons_output = new xAOD::MuonContainer();
       SG::IAuxStore* muons_aux = evtStore()->event()->recordAux(m_prefix + m_config->sgKeyMuons() + "Aux.");
       muons_output->setStore(muons_aux);
  
       for (const auto* const muPtr : event.m_muons) {
         xAOD::Muon* outputMuPtr = new xAOD::Muon();
         outputMuPtr->makePrivateStore(*muPtr);
         muons_output->push_back(outputMuPtr);
       }
  
       top::check(evtStore()->event()->record(muons_output,
                                              m_prefix + m_config->sgKeyMuons()),
                  "xAOD::TEvent record " + m_prefix + m_config->sgKeyMuons() + " failed");
  
       //for now...
       //Muon charge is a property of the associated track
       //Copy tracks (for muon charge)
       check(evtStore()->event()->copy("InDetTrackParticles"), "xAOD::TEvent copy InDetTrackParticles failed");
       check(evtStore()->event()->copy(
               "ExtrapolatedMuonTrackParticles"),
             "xAOD::TEvent copy InDetTrackParticles ExtrapolatedMuonTrackParticles");
       check(evtStore()->event()->copy(
               "CombinedMuonTrackParticles"), "xAOD::TEvent copy InDetTrackParticles CombinedMuonTrackParticles");
     }
  
     //Jets
     if (m_config->useJets()) {
       xAOD::JetContainer* jets_output = new xAOD::JetContainer();
       SG::IAuxStore* jets_aux = evtStore()->event()->recordAux(m_prefix + m_config->sgKeyJets() + "Aux.");
       jets_output->setStore(jets_aux);
  
       for (const auto* const jetPtr : event.m_jets) {
         xAOD::Jet* outputJetPtr = new xAOD::Jet();
         outputJetPtr->makePrivateStore(*jetPtr);
         jets_output->push_back(outputJetPtr);
       }
  
       top::check(evtStore()->event()->record(jets_output,
                                              m_prefix + m_config->sgKeyJets()),
                  "xAOD::TEvent record " + m_prefix + m_config->sgKeyJets() + " failed");
  
       //b-tagging variables need the BTagging_JetNameWithoutJetsOnTheEnd container
       //remove Jets from the end of the container name and add BTagging_
       const std::string btagName = "BTagging_" + m_config->sgKeyJets().substr(0, m_config->sgKeyJets().size() - 4);
       check(evtStore()->event()->copy(btagName), "xAOD::TEvent copy " + btagName + " failed");
     }
  
     //Large Jets
     if (m_config->useLargeRJets()) {
       xAOD::JetContainer* largeJets_output = new xAOD::JetContainer();
       SG::IAuxStore* largeJets_aux = evtStore()->event()->recordAux(m_prefix + m_config->sgKeyLargeRJets() + "Aux.");
       largeJets_output->setStore(largeJets_aux);
  
       for (const auto* const jetPtr : event.m_largeJets) {
         xAOD::Jet* outputJetPtr = new xAOD::Jet();
         outputJetPtr->makePrivateStore(*jetPtr);
         largeJets_output->push_back(outputJetPtr);
       }
  
       check(evtStore()->event()->record(largeJets_output,
                                         m_prefix + m_config->sgKeyLargeRJets()),
             "xAOD::TEvent record " + m_prefix + m_config->sgKeyLargeRJets() + " failed");
     }
  
     //MET
     top::check(evtStore()->event()->copy(m_config->sgKeyMissingEt()),
                "xAOD::TEvent copy " + m_config->sgKeyMissingEt() + " failed");
  
     //do it!
     top::check(evtStore()->event()->fill(), "xAOD::TEvent fill failed");
   }

◆ saveEventToxAOD()

virtual void top::EventSaverBase::saveEventToxAOD ( )

inlinevirtualinherited

Reimplemented in top::EventSaverxAODNext.

Definition at line 60 of file EventSaverBase.h.

60 {return;}

◆ saveParticleLevelEvent()

virtual void top::EventSaverBase::saveParticleLevelEvent ( const top::ParticleLevelEvent & )

inlinevirtualinherited

Function to save the particle level event.

This function has a default implementation that does ... nothing! This is deliberate since doing so does not force developers of custom EventSaver tools to implement particle level saving.

Parameters

plEvent The particle level event whose contents will be saved.

Reimplemented in top::EventSaverFlatNtuple.

Definition at line 79 of file EventSaverBase.h.

                                                                         {
       /* Deliberately No-Op */
       return;
     }

◆ saveTruthEvent()

virtual void top::EventSaverBase::saveTruthEvent ( )

inlinevirtualinherited

@breif Run for every event

This fills, if requested, the: MC truth record PDF info TopPartons

Reimplemented in top::EventSaverFlatNtuple.

Definition at line 70 of file EventSaverBase.h.

70 {return;}

◆ saveUpgradeEvent()

virtual void top::EventSaverBase::saveUpgradeEvent ( const top::ParticleLevelEvent & )

inlinevirtualinherited

Function to save the upgrade event.

This function has a default implementation that does ... nothing! This is deliberate since doing so does not force developers of custom EventSaver tools to implement upgrade saving.

Parameters

upgradeEvent The upgrade event whose contents will be saved.

Definition at line 91 of file EventSaverBase.h.

                                                                   {
       /* Deliberately No-Op */
       return;
     }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                                {
     // 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_config

std::shared_ptr<top::TopConfig> top::EventSaverxAOD::m_config

private

We need access to the configuration file to get the container names.

Definition at line 80 of file EventSaverxAOD.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_outputFile

TFile* top::EventSaverxAOD::m_outputFile

private

We need to hold on to the output file.

Definition at line 83 of file EventSaverxAOD.h.

◆ m_prefix

const std::string top::EventSaverxAOD::m_prefix

private

Name to prepend to the output containers.

Definition at line 77 of file EventSaverxAOD.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ EventSaverxAOD()

◆ ~EventSaverxAOD()

Member Function Documentation

◆ ClassDef() [1/2]

◆ ClassDef() [2/2]

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ execute()

◆ extraDeps_update_handler()

◆ finalize()

◆ getKey()

◆ getName()

◆ getProperty()

◆ initialize() [1/2]

◆ initialize() [2/2]

◆ inputHandles()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msg_level_name()

◆ msgLvl()

◆ outputHandles()

◆ print()

◆ renounce()

◆ renounceArray()

◆ saveEvent()

◆ saveEventToxAOD()

◆ saveParticleLevelEvent()

◆ saveTruthEvent()

◆ saveUpgradeEvent()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_config

◆ m_detStore

◆ m_evtStore

◆ m_outputFile

◆ m_prefix

◆ m_varHandleArraysDeclared

◆ m_vhka