#include <OldSpclMcFilterTool.h>

Inheritance diagram for OldSpclMcFilterTool:

Collaboration diagram for OldSpclMcFilterTool:

Public Member Functions
	OldSpclMcFilterTool (const std::string &type, const std::string &name, const IInterface *parent)
	Constructor with parameters: More...

virtual	~OldSpclMcFilterTool ()
	Destructor: More...

StatusCode	buildMcAod (const McEventCollection in, McEventCollection out)
	This method will check the validity of the input McEventCollection and build a filtered one from the strategy implemented by this concrete tool. More...

StatusCode	initialize ()
	Athena Algorithm's Hooks. More...

StatusCode	execute ()

StatusCode	finalize ()

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

Static Public Member Functions
static const InterfaceID &	interfaceID ()
	Inline methods: More...

Protected Types
typedef ToolHandle< IMcVtxFilterTool >	McVtxFilterTool_t

typedef ToolHandle< ITruthIsolationTool >	IsolTool_t

Protected Member Functions
	OldSpclMcFilterTool ()
	Default constructor: More...

StatusCode	selectSpclMcBarcodes ()
	Retrieve the GenParticles from the GenEvent object (located into the McEventCollection container), filter them and store their barcode into the list of barcodes of particles to keep. More...

StatusCode	shapeGenEvent (McEventCollection *outEvt)
	Build the out McEventCollection by removing all the GenParticles whose barcode is not in the list of barcodes to keep. More...

StatusCode	reconnectParticles (const McEventCollection inEvt, McEventCollection out)
	Reconnect the particles of the filtered McEventCollection. More...

StatusCode	rebuildLinks (const HepMC::GenEvent mcEvent, HepMC::GenEvent outEvt, const HepMC::GenParticlePtr &mcPart)
	Rebuild the links between particles which were connected, eg by a long FSR chain. More...

StatusCode	initializeTool ()
	to get tesIO More...

virtual StatusCode	finalizeTool ()

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

Protected Attributes
std::set< int >	m_barcodes
	List of particles (in fact their unique barcode) which have been labelled as "special". More...

DoubleProperty	m_rIsol
	deltaR isolation energy cut for electrons, muons, taus and photons More...

DoubleProperty	m_ptGamMin
	Minimum threshold for transverse momentum of photons. More...

DoubleProperty	m_ptMin
	Minimum threshold for transverse momentum for all particles. More...

DoubleProperty	m_etaRange
	Eta acceptance cut applied on all stable particles. More...

BooleanProperty	m_includeSimul
	Switch to include or not particles from detector simulation (Geant4) More...

BooleanProperty	m_includePartonShowers
	Switch to include or not parton showers. More...

TruthHelper::GenAccessIO *	m_tesIO
	Pointer to a StoreGateSvc helper (Mc particle retrieval with predicates) More...

BooleanProperty	m_removeDecayToSelf
	Switch to remove particles which decay into themselves (t->tg) but only for generated particles, not the ones from the Geant4 interactions. More...

McVtxFilterTool_t	m_mcVtxFilterTool
	Pointer to the McVtxFilterTool to be able to select additional vertices on some decay pattern criterion. More...

IsolTool_t	m_isolationTool
	Pointer to the TruthIsolationTool to be able to compute transverse energy isolations for various isolation cones cuts. More...

BooleanProperty	m_doEtIsolation
	Switch to compute or not the Et-isolations for `TruthParticle` (and their underlying `HepMC::GenParticle`). More...

SG::ReadHandleKey< McEventCollection >	m_mcEventsReadHandleKey {this,"McEvents","TruthEvent","ReadHandleKey for input McEventCollection one wants to filter"}
	ReadHandleKey for the input McEventCollection one wants to filter. More...

SG::WriteHandleKey< McEventCollection >	m_mcEventsOutputWriteHandleKey {this,"McEventsOutput","GEN_AOD","WriteHandleKey for the output McEventCollection which has been filtered"}
	Location of the output McEventCollection which has been filtered. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
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

Definition at line 29 of file OldSpclMcFilterTool.h.

Member Typedef Documentation

◆ IsolTool_t

typedef ToolHandle<ITruthIsolationTool> TruthParticleFilterBaseTool::IsolTool_t

protectedinherited

Definition at line 80 of file TruthParticleFilterBaseTool.h.

◆ McVtxFilterTool_t

typedef ToolHandle<IMcVtxFilterTool> TruthParticleFilterBaseTool::McVtxFilterTool_t

protectedinherited

Definition at line 74 of file TruthParticleFilterBaseTool.h.

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ OldSpclMcFilterTool() [1/2]

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

Constructor with parameters:

Definition at line 41 of file OldSpclMcFilterTool.cxx.

                                                  : 
   TruthParticleFilterBaseTool( type, name, parent ),
   m_tesIO( nullptr )
 {
   //
   // Property declaration
   // 
   //declareProperty( "Property", m_nProperty );
  
   declareProperty( "rIsolation",
            m_rIsol = 0.45,
            "DeltaR isolation energy cut for electrons, muons, "
            "taus and photons" );
  
   declareProperty( "ptGammaMin",
            m_ptGamMin = 0.5*GeV,
            "Minimum threshold for transverse momentum of photons" );
  
   declareProperty( "ptMin",             
            m_ptMin = 0.5*GeV,
            "Minimum threshold for transverse momentum for all particles.\n"
            "Warning: this cut is applied *before* Pt photon cut !" );
  
   declareProperty( "etaRange", 
            m_etaRange = 5.0,
            "Eta acceptance cut applied on all stable particles" );
  
   // switches
  
   declareProperty( "IncludeSimul",         
            m_includeSimul = true,
            "Switch to include or not particles from detector simulation "
            "(Geant 4)" );
  
   declareProperty( "IncludePartonShowers", 
            m_includePartonShowers = true,
            "Switch to include or not parton showers" );
  
   declareProperty( "RemoveDecayToSelf",    
            m_removeDecayToSelf = true,
            "Switch to remove particles which decay into themselves (t->tg) "
            "*but* only for generated particles, not the ones from the " 
            "Geant4 interactions" );
  
 }

◆ ~OldSpclMcFilterTool()

OldSpclMcFilterTool::~OldSpclMcFilterTool ( )

virtual

Destructor:

Definition at line 91 of file OldSpclMcFilterTool.cxx.

 { 
   delete m_tesIO; m_tesIO = nullptr;
   ATH_MSG_DEBUG("Calling destructor");
 }

◆ OldSpclMcFilterTool() [2/2]

OldSpclMcFilterTool::OldSpclMcFilterTool ( )

protected

Default constructor:

Member Function Documentation

◆ buildMcAod()

StatusCode OldSpclMcFilterTool::buildMcAod	(	const McEventCollection *	in,
		McEventCollection *	out
	)

virtual

This method will check the validity of the input McEventCollection and build a filtered one from the strategy implemented by this concrete tool.

Implements ITruthParticleFilterTool.

Definition at line 101 of file OldSpclMcFilterTool.cxx.

 {
   if ( nullptr == in || nullptr == out ) {
     ATH_MSG_ERROR("Invalid pointer to McEventCollection !" << endmsg
           << "  in: " << in << endmsg
           << " out: " << out);
     return StatusCode::FAILURE;
   }
  
   // we just copy the input McEventCollection and put it into the output one
   out->operator=( *in );
  
   // select the barcodes of the "special" particles from the input GenEvent
   if ( selectSpclMcBarcodes().isFailure() ) {
     ATH_MSG_ERROR("Could not select the \"special\" barcodes !!");
     return StatusCode::FAILURE;
   }
  
   // remove the not "special" particles from the filtered McEventCollection
   if ( shapeGenEvent(out).isFailure() ) {
     ATH_MSG_ERROR("Could not remove the not \"special\" particles from the "\
           "filtered McEventCollection !!");
     return StatusCode::FAILURE;
   }
  
   // reconnect the particles
   if ( reconnectParticles(in, out).isFailure() ) {
     ATH_MSG_ERROR("Could not reconnect the particles in the filtered "\
           "McEventCollection !!");
     return StatusCode::FAILURE;
   }
  
   return StatusCode::SUCCESS;
 }

◆ 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()

StatusCode TruthParticleFilterBaseTool::execute ( )

virtualinherited

Implements ITruthParticleFilterTool.

Definition at line 120 of file TruthParticleFilterBaseTool.cxx.

 {  
   ATH_MSG_DEBUG("Executing " << name() << "...");
  
   //Setup Handle to read input container
   SG::ReadHandle<McEventCollection> mcEventsReadHandle(m_mcEventsReadHandleKey);
  
   if (!mcEventsReadHandle.isValid()){
     ATH_MSG_ERROR("Invalid ReadHandle to McEventColleciton with key: " << m_mcEventsReadHandleKey.key());
     return StatusCode::FAILURE;
   }
  
   //Setup WriteHandle, and then record new McEventCollection.
   SG::WriteHandle<McEventCollection> mcEventsOutputWriteHandle(m_mcEventsOutputWriteHandleKey);
   ATH_CHECK(mcEventsOutputWriteHandle.record(std::make_unique<McEventCollection>()));  
   
   if (!mcEventsOutputWriteHandle.isValid()){
     ATH_MSG_ERROR("Invalid WriteHamdle for McEventCollection with key ["
                   <<m_mcEventsOutputWriteHandleKey.key()  << "] !!");
     return StatusCode::FAILURE;
   }
   
   // Compute isolation for gamma/lepton.
   if ( m_doEtIsolation.value() ) {
     ATH_MSG_VERBOSE("Computing Et isolations...");
     if ( m_isolationTool->buildEtIsolations(m_mcEventsReadHandleKey.key()).isFailure() ) {
       ATH_MSG_ERROR("Could not compute Et isolations !!");
       return StatusCode::FAILURE;
     }
   } //> end do Et-isolation
  
     if ( this->buildMcAod( mcEventsReadHandle.ptr(), mcEventsOutputWriteHandle.ptr() ).isFailure() ) {
     ATH_MSG_ERROR("Could not buildMcAod(in,out) !!");
     return StatusCode::FAILURE;
   }
  
   // We have slimmed the mcEventsOutputWriteHandle.
   // To not bias the map of isolation energies for this GenEvent, we alias
   // it to its original one
   if ( m_doEtIsolation.value() &&
        !m_isolationTool->registerAlias( m_mcEventsOutputWriteHandleKey.key(), 
                                         m_mcEventsReadHandleKey.key() 
                                         ).isSuccess() ) {
     ATH_MSG_WARNING("Could not create an alias in the map of "\
             "isolation energies !");
   }
  
   return StatusCode::SUCCESS;
 }

◆ 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()

StatusCode TruthParticleFilterBaseTool::finalize ( )

virtualinherited

Implements ITruthParticleFilterTool.

Definition at line 105 of file TruthParticleFilterBaseTool.cxx.

 {
   ATH_MSG_INFO("Finalizing " << name() << "...");
  
   m_mcVtxFilterTool->stats();
  
   // Give the concrete (derived) tool a chance to finalize itself
   if ( finalizeTool().isFailure() ) {
     ATH_MSG_ERROR("Could not finalize concrete tool !");
     return StatusCode::FAILURE;
   }
  
   return StatusCode::SUCCESS;
 }

◆ finalizeTool()

StatusCode TruthParticleFilterBaseTool::finalizeTool ( )

inlineprotectedvirtualinherited

Definition at line 124 of file TruthParticleFilterBaseTool.h.

 {
   return StatusCode::SUCCESS;
 }

◆ initialize()

StatusCode TruthParticleFilterBaseTool::initialize ( )

virtualinherited

Athena Algorithm's Hooks.

Retrieves a private AlgTool to filter a McEventCollection

Retrieves the isolation tool (public, to be used also in TruthParticleCnvTool)

Implements ITruthParticleFilterTool.

Definition at line 65 of file TruthParticleFilterBaseTool.cxx.

 {
   ATH_MSG_INFO("Initializing " << name() << "...");
  
   if ( !m_mcVtxFilterTool.retrieve().isSuccess() ) {
     ATH_MSG_ERROR("Creation of algTool IMcVtxFilterTool FAILED !");
     return StatusCode::FAILURE;
   }
  
   if( m_doEtIsolation.value()) {
     ATH_CHECK(m_isolationTool.retrieve());
   }
   else {
     m_isolationTool.disable();
   }
  
   //initialize DataHandleKeys
   ATH_CHECK(m_mcEventsReadHandleKey.initialize());
   ATH_CHECK(m_mcEventsOutputWriteHandleKey.initialize());
   
   ATH_MSG_INFO
     (" DoEtIsolations: [" << std::boolalpha << m_doEtIsolation.value()
      << "]" << endmsg
      << " McEvents:       [" << m_mcEventsReadHandleKey.key() << "]" << endmsg
      << " McEventsOutput: [" << m_mcEventsOutputWriteHandleKey.key() << "]");
   
   // Give the concrete (derived) tool a chance to initialize itself
   if ( initializeTool().isFailure() ) {
     ATH_MSG_ERROR("Could not initialize concrete tool !");
     return StatusCode::FAILURE;
   }
  
   ATH_MSG_INFO("Options of the McVtxFilterTool:");
   m_mcVtxFilterTool->displayOptions();
  
   return StatusCode::SUCCESS;
 }

◆ initializeTool()

StatusCode OldSpclMcFilterTool::initializeTool ( )

protectedvirtual

to get tesIO

Reimplemented from TruthParticleFilterBaseTool.

Definition at line 626 of file OldSpclMcFilterTool.cxx.

 {
   ATH_MSG_DEBUG("Calling initializeTool");
   delete m_tesIO; m_tesIO = nullptr;
   // accessor for particles  
   m_tesIO = new GenAccessIO();
   if( m_tesIO == nullptr ) {
     ATH_MSG_ERROR("Unable to retrieve GenAccessIO pointer");
     return StatusCode::FAILURE;
   }
   return StatusCode::SUCCESS;
 }

◆ inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

◆ interfaceID()

const InterfaceID & ITruthParticleFilterTool::interfaceID ( )

inlinestaticinherited

Inline methods:

Definition at line 70 of file ITruthParticleFilterTool.h.

 { 
    return IID_ITruthParticleFilterTool; 
 }

◆ 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);
   }

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

◆ rebuildLinks()

StatusCode OldSpclMcFilterTool::rebuildLinks	(	const HepMC::GenEvent *	mcEvent,
		HepMC::GenEvent *	outEvt,
		const HepMC::GenParticlePtr &	mcPart
	)

protected

Rebuild the links between particles which were connected, eg by a long FSR chain.

Definition at line 485 of file OldSpclMcFilterTool.cxx.

 {
  
   if ( !mcPart ) {
     ATH_MSG_WARNING("Null GenParticle: can not rebuildLinks");
     return StatusCode::FAILURE;
   }
  
   if ( mcPart->end_vertex() ) {
     ATH_MSG_VERBOSE("GenParticle has already a decay vertex : nothing to do");
     return StatusCode::SUCCESS;
   }
  
   if ( !mcEvt ) {
     ATH_MSG_WARNING("Null input HepMC::GenEvent : can not rebuildLinks");
     return StatusCode::FAILURE;
   }
  
   if ( !outEvt ) {
     ATH_MSG_WARNING("Null output HepMC::GenEvent: can not rebuildLinks");
     return StatusCode::FAILURE;
   }
 #ifdef HEPMC3
   return StatusCode::SUCCESS;
 #else
   // Cache some useful infos
   const int pdgId = mcPart->pdg_id();
   const int bc    = HepMC::barcode(mcPart);
  
 //AV: Const correctness is broken for HepMC2.
   HepMC::GenParticlePtr inPart = HepMC::barcode_to_particle(mcEvt,bc);
   HepMC::GenVertexPtr   dcyVtx = inPart->end_vertex();
  
   if ( !dcyVtx ) {
     ATH_MSG_VERBOSE("No decay vertex for the particle #" << bc << " : " << "No link to rebuild...");
     return StatusCode::SUCCESS;
   }
  
   std::list<int> bcChildPart;
   std::list<int> bcChildVert;
  
   //
   // Loop over all descendants of the GenParticle
   // Store the barcode of the GenParticles entering into each GenVertex
   //
   const HepMC::GenVertex::vertex_iterator endVtx = dcyVtx->vertices_end(HepMC::descendants);
   for ( HepMC::GenVertex::vertex_iterator itrVtx = dcyVtx->vertices_begin( HepMC::descendants );
     itrVtx != endVtx;
     ++itrVtx ) {
     bool foundPdgId = false;
     HepMC::GenVertex::particles_in_const_iterator endPart = (*itrVtx)->particles_in_const_end();
     for ( HepMC::GenVertex::particles_in_const_iterator itrPart = (*itrVtx)->particles_in_const_begin();
       itrPart != endPart;
       ++itrPart ) {
  
       // because the vertices are traversed in POST ORDER !!
       // (quoting GenVertex::vertex_iterator)
       bcChildPart.push_front( (*itrPart)->barcode() );
  
       if ( (*itrPart)->pdg_id() == pdgId ) {
     foundPdgId = true;
       }
     }//> loop over in-going particles of this vertex
  
     if ( foundPdgId ) {
       bcChildVert.push_front( (*itrVtx)->barcode() );
     }
  
   }//> loop over descendants of decay vertex
  
   //
   // Now we loop over the previously stored barcodes and
   // we connect our GenParticle to the first found barcode
   // 
   std::list<int>::const_iterator bcVtxEnd = bcChildVert.end();
   for ( std::list<int>::const_iterator itrBcVtx = bcChildVert.begin();
     itrBcVtx != bcVtxEnd;
     ++itrBcVtx ) {
     HepMC::GenVertex * childVtx = outEvt->barcode_to_vertex(*itrBcVtx);
     if ( childVtx ) {
       if ( childVtx->particles_in_size() > 0 ) {
     HepMC::GenVertex::particles_in_const_iterator endPart = childVtx->particles_in_const_end();
     for ( HepMC::GenVertex::particles_in_const_iterator itrPart = childVtx->particles_in_const_begin();
           itrPart != endPart;
           ++itrPart ) {
       if ( (*itrPart)->pdg_id() == pdgId ) {
         HepMC::GenVertex * prodVtx = (*itrPart)->production_vertex();
         if ( prodVtx ) {
           if ( prodVtx->particles_in_size() > 0 ) {
         // Humm... This is not what we'd have expected
         // so we skip it
         if ( msgLvl(MSG::VERBOSE) ) {
           msg(MSG::VERBOSE)
             << "found a particle = "
             << (*itrPart) << ", "
             << "but its production vertex has incoming particles !"
             << endmsg;
           continue;
         }
         // create a GenVertex which will be the decay vertex of our
         // GenParticle and the production vertex of the GenParticle
         // we just found
         HepMC::GenVertexPtr linkVtx = HepMC::newGenVertexPtr();
         outEvt->add_vertex( linkVtx );
         linkVtx->add_particle_in( mcPart );
         linkVtx->add_particle_out( *itrPart );
         
         msg(MSG::ERROR)
           << "====================================================="
           << endmsg
           << "Created a GenVertex - link !"
           << std::endl;
         std::stringstream vtxLink("");
         linkVtx->print(vtxLink);
         msg(MSG::ERROR)
           << vtxLink.str()
           << endmsg
           << "====================================================="
           << endmsg;
           }
         }
       }
     }//> loop over incoming particles
       } else { 
     // no incoming particle : so we just add this particle
     // a bit odd though : FIXME ?
     childVtx->add_particle_in(mcPart);
     msg(MSG::WARNING) << "Odd situation:" << std::endl;
     std::stringstream vtxDump( "" );
     childVtx->print(vtxDump);
     msg(MSG::WARNING) << vtxDump.str() << endmsg;
     return StatusCode::SUCCESS;
       }//> end if incoming particles
     }//> found a child-vertex
   }//> loop over child-vertex-barcodes
   return StatusCode::FAILURE;
 #endif
 }

◆ reconnectParticles()

StatusCode OldSpclMcFilterTool::reconnectParticles	(	const McEventCollection *	inEvt,
		McEventCollection *	out
	)

protected

Reconnect the particles of the filtered McEventCollection.

We need the original McEventCollection to know if 2 particles were connected

Definition at line 429 of file OldSpclMcFilterTool.cxx.

 {
   if ( nullptr == in || nullptr == out ) {
     ATH_MSG_ERROR("Invalid pointer to McEventCollection !!" << endmsg
           << "  in: " << in << endmsg
           << " out: " << out);
     return StatusCode::FAILURE;
   }
 #ifdef HEPMC3
   return StatusCode::SUCCESS;
 #else
   for ( unsigned int iEvt = 0; iEvt != in->size(); ++iEvt) {
     const HepMC::GenEvent * evt    = (*in)[iEvt];
     HepMC::GenEvent       * outEvt = (*out)[iEvt];
     
     // Reconnect the particles
     ATH_MSG_VERBOSE("Reconnecting particles...");
     for (const auto& itrPart: *outEvt) {
       if ( itrPart->end_vertex() ) {
     continue;
       }
       if ( rebuildLinks( evt, outEvt, itrPart ).isFailure() ) {
     ATH_MSG_WARNING("Could not rebuild links for this particle = "<< itrPart);
       } else if ( msgLvl(MSG::VERBOSE) ) {
     msg(MSG::VERBOSE)
       << "==========================================================="
       << endmsg << "Production vertex for particle "  << itrPart << " : ";
     if ( itrPart->production_vertex() ) {
       std::stringstream prodVtx("");
       HepMC::Print::line(prodVtx,itrPart->production_vertex());
       msg(MSG::VERBOSE) << std::endl << prodVtx.str() << endmsg;
     } else {
       msg(MSG::VERBOSE) << "[No production vertex]" << endmsg;
     }
     
     msg(MSG::VERBOSE) << "Decay vertex for particle " 
               << itrPart << " : ";
     if ( itrPart->end_vertex() ) {
       std::stringstream dcyVtx("");
       HepMC::Print::line(dcyVtx, itrPart->end_vertex());
       msg(MSG::VERBOSE) << std::endl
                 << dcyVtx.str()
                 << endmsg;
     } else {
       msg(MSG::VERBOSE) << endmsg << "[No decay vertex]" << endmsg;
     }
       }//> end VERBOSE messages
       
     }//> loop over particles
   }//> loop over GenEvents
   
   return StatusCode::SUCCESS;
 #endif
 }

◆ 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();
   }

◆ selectSpclMcBarcodes()

StatusCode OldSpclMcFilterTool::selectSpclMcBarcodes ( )

protected

Retrieve the GenParticles from the GenEvent object (located into the McEventCollection container), filter them and store their barcode into the list of barcodes of particles to keep.

skip stuff with no end-vertex

stable particles

Save photons

Save photon only if it does not decay (ie no end_vertex)

ask McVtxFilterTool for its decision ie: final word is given to the user

Definition at line 140 of file OldSpclMcFilterTool.cxx.

 {
   StatusCode sc = StatusCode::SUCCESS;
  
   // Get all of the generated particles (does not have Geant secondaries)
  
   std::vector<HepMC::ConstGenParticlePtr> particles;
   if ( m_includeSimul ) {
     sc = m_tesIO->getMC(particles, false,  m_mcEventsReadHandleKey.key());
   } else {
     sc = m_tesIO->getMC(particles, true, m_mcEventsReadHandleKey.key());
   }
   if ( sc.isFailure() ) {
     ATH_MSG_ERROR("Could not get Monte Carlo particles from TDS at : "
           << m_mcEventsReadHandleKey.key());
     return StatusCode::FAILURE;
   }
  
   m_barcodes.clear();
  
   // Loop over all particles, selecting special ones
   // keep track of them using their barcodes
  
   for (const auto& part : particles) {
  
     const int id      = part->pdg_id();
     const int ida     = std::abs(id);
     const HepMC::FourVector hlv = part->momentum();
     const double pt   = hlv.perp();
     const double eta  = hlv.pseudoRapidity();
     const double mass = hlv.m();
     const int barcode = HepMC::barcode(part);
  
     HepMC::ConstGenVertexPtr  decayVtx = part->end_vertex();
     HepMC::ConstGenVertexPtr  prodVtx  = part->production_vertex();
  
     bool isSpcl = false;
  
     if( !MC::isStable(part) && !decayVtx ) continue;
  
     const bool accept = pt > m_ptMin;
  
     // Select special particles
  
     if (m_includeSimul ) {
         if ( MC::isSimStable(part) && accept && std::abs(eta) < m_etaRange ) isSpcl = true;
     } else {
         if ( MC::isGenStable(part) && accept && std::abs(eta) < m_etaRange ) isSpcl = true;
     }
     // e, mu, tau, neutrino 
     // Hard coded cut pt>2GeV to remove beam-jet Dalitz decays, etc.
     if( MC::isSMLepton(part) && pt>2.*GeV && std::abs(eta) < m_etaRange ) isSpcl = true;
  
    if ( MC::isPhoton(part) && pt>m_ptGamMin && std::abs(eta) < m_etaRange ) {
      if ( nullptr == part->end_vertex() ) {
         isSpcl=true;
      } 
    } 
  
     // Any long-lived B meson or baryon 5x1 or 5xx2
     int ifl=0;
     if( ida>100 && ida<1000) ifl=ida/100;
     if( ida>1000 && ida<10000) ifl=ida/1000;
     int jj1 = ida%10;
     if( (ifl==5 && jj1<3 && mass<9.0*GeV) || ( ifl==4 && jj1<3 && mass<2.4*GeV) ) {
       if (accept) isSpcl = true;
     } 
  
     // Any heavy particle other than quarks or gluons.
     // M > 9.0GeV excludes B hadrons but includes Upsilon, W, Z, t, ...
     if( ida>5 && ida!=21 && ida!=9 && mass>9.*GeV && accept ) isSpcl=true;
  
     // save the quarks
     if (ida >= 1 && ida <= 5 && accept ) isSpcl = true;
  
     // save the gluons
     if ( ida == 21 && accept ) isSpcl = true;
  
     // save the double charged Higgs
     if ( (ida == 9900041 || ida == 9900042) && accept ) isSpcl = true; 
  
     // SUSY particles -- relies on PDG codes
     if( ida>1000000 && ida<3000000 && accept ) isSpcl = true;
  
     // Exclude invalid special particles
     
     // Monte Carlo internal particles per PDG
     if( ida>80 && ida<101 ) isSpcl=false;
  
     // Particles decaying into itself
     // This prevents documentaries from being reported twice and
     // eliminates (if asked to) parton showers like q -> qg.
     // But should not delete particle if child is from GEANT.
     if( isSpcl && decayVtx ) {
         auto dcyVtx = part->end_vertex();
         for(const auto& child: *dcyVtx) {
     if( child->pdg_id()==id && //> looking for parton showers or documentaries
         HepMC::barcode(child) !=barcode  && //> just to be sure that merging GEN_EVENT/G4Truth is OK
         !HepMC::is_simulation_particle(child)  //> child is not from GEANT
         ) {
       if ( m_includePartonShowers ) {
         // we keep the parent particle
         // this could be vertices of the form :
         // q -> q + g + b + bbar(hence the outgoing quark isn't decaying:documentary)
         // t -> t + W + b (with outgoing top being a documentary)
         // q -> q + g (real parton shower)
         isSpcl = true;
       } else {
         isSpcl = false;
       }
     }
       }
     }
  
     if ( ( decayVtx && m_mcVtxFilterTool->isAccepted(decayVtx) ) || 
      ( prodVtx  && m_mcVtxFilterTool->isAccepted(prodVtx) ) ) {
       isSpcl = true;
     }
  
     // Save special particles and children
     if( !isSpcl ) continue;
     m_barcodes.insert(barcode); // add it to list
  
     // Children
     if( isSpcl && decayVtx ) {
       for(const auto& child: *(part->end_vertex())) {
         if( HepMC::is_simulation_particle(child) && !m_removeDecayToSelf) { 
       m_barcodes.insert(HepMC::barcode(child));// its not there already
         }
       }
     }
  
   }
   ATH_MSG_DEBUG("Read " << particles.size() 
         << " and selected  " <<  m_barcodes.size() << " particles");
  
   return StatusCode::SUCCESS;
 }

◆ shapeGenEvent()

StatusCode OldSpclMcFilterTool::shapeGenEvent ( McEventCollection * outEvt )

protected

Build the out McEventCollection by removing all the GenParticles whose barcode is not in the list of barcodes to keep.

Remove also the dead vertices.

Definition at line 294 of file OldSpclMcFilterTool.cxx.

 {
   //now remove all the particles except those whose barcodes are marked
   for ( McEventCollection::iterator evt = genAod->begin(); evt != genAod->end();++evt) {
     std::vector<HepMC::GenParticlePtr> going_out;
  
 #ifdef HEPMC3
     for ( auto & p :  ((HepMC::GenEvent*)(*evt))->particles()) {
 //AV: We modify event      
       int pBC = HepMC::barcode(p);
       ATH_MSG_DEBUG("[pdg,bc]= " << p->pdg_id() << ", " << pBC);
       if ( m_barcodes.count(pBC) != 0 ) continue; 
       going_out.push_back(p); // list of useless particles
       auto pvtx = p->production_vertex();
       auto evtx = p->end_vertex();
       if (pvtx) pvtx->remove_particle_out(p); //remove from production vertex from useless partilcle
       if (evtx) { // if it has end vertex, may need to move the out partilces
         if(pvtx){ // move the partilces back
           for (auto& pp: evtx->particles_out()) {
             pvtx->add_particle_out(pp);
           }
         }
         evtx->remove_particle_out(p); // disconnect from end vertex
       }
     }//> loop over particles
 #else
     std::list<int> evtBarcodes;
     for ( const auto& p: **evt ) {
       evtBarcodes.push_back( HepMC::barcode(p) );
     }
     for ( std::list<int>::const_iterator itrBc = evtBarcodes.begin();itrBc != evtBarcodes.end(); ++itrBc ) {
 //AV: We modify event      
       HepMC::GenParticlePtr p = HepMC::barcode_to_particle((HepMC::GenEvent*)(*evt),*itrBc);
       int pBC = HepMC::barcode(p);
       ATH_MSG_DEBUG("[pdg,bc]= " << p->pdg_id() << ", " << pBC);
       if ( m_barcodes.count(pBC) == 0 ) { 
         going_out.push_back(p); // list of useless particles
         auto pvtx = p->production_vertex();
         auto evtx = p->end_vertex();
     if (pvtx) pvtx->remove_particle(p); //remove from production vertex from useless partilcle
     if (evtx) { // if it has end vertex, may need to move the out partilces
       if(pvtx){ // move the partilces back
         if ( msgLvl(MSG::DEBUG) ) {
           msg(MSG::DEBUG) << "\tin endVtx   "<< endmsg;
         }
         while ( evtx->particles_out_const_begin() !=  evtx->particles_out_const_end()) {
           HepMC::GenVertex::particles_out_const_iterator np = evtx->particles_out_const_begin();
           pvtx->add_particle_out(*np); // note that this really is a MOVE!!! it get taken off evtx by magic
         }
       }//> end if [prod vertex]
       evtx->remove_particle(p); // disconnect from end vertex     
     }//> end if [decay vertex]
       }//> particle has to be removed
     }//> loop over particles (via their barcode)
 #endif
  
  
 #ifdef HEPMC3
     // there may be a bunch of vertices with no particles connected to them:
     // ==> Get rid of them //AV: Not sure if this is needed
     std::vector<HepMC::ConstGenVertexPtr> going_out_again;
     for ( auto& v: (*evt)->vertices()) {
       if ( v->particles_in().empty() && v->particles_out().empty() ){
         going_out_again.push_back(v);
       }
     }//> loop over vertices
 //HepMC3 uses smart pointers
 #else
     // now get rid of all dead particles
     for ( std::vector<HepMC::GenParticle*>::iterator d = going_out.begin(); 
       d != going_out.end(); 
       ++d ){
       delete *d;
     }
  
     // there may be a bunch of vertices with no particles connected to them:
     // ==> Get rid of them
     std::vector<HepMC::GenVertex*> going_out_again;
     for ( HepMC::GenEvent::vertex_const_iterator v = (*evt)->vertices_begin();
       v != (*evt)->vertices_end(); ++v ) {
       if ( (*v)->particles_in_size() == 0 && (*v)->particles_out_size() == 0 ){
     going_out_again.push_back(*v);
       }
     }//> loop over vertices
  
    // now get rid of all dead vertices
     for ( std::vector<HepMC::GenVertex*>::iterator d = going_out_again.begin();
       d != going_out_again.end(); 
       ++d ){
       delete *d;
     }
 #endif
     
   }//> loop over GenEvents in McEventCollection
   
   // Set the signal_process_vertex to NULL if not to be recorded 
   for ( McEventCollection::iterator evt = genAod->begin();  evt != genAod->end(); ++evt) {
 #ifdef HEPMC3
     auto sigProcVtx = HepMC::signal_process_vertex(*evt); 
     if (!sigProcVtx) continue;
       const int sigProcBC = HepMC::barcode(sigProcVtx); 
       bool isInColl = false; 
       for ( const auto& itrVtx: (*evt)->vertices() ) { 
         if ( sigProcBC == HepMC::barcode(itrVtx) ) { 
           isInColl = true; 
           break; 
         } 
       }  //> loop over vertices 
 //AV: We don't set nullptr as signal vertex in HepMC3
     if ( !isInColl ) { 
        (*evt)->remove_attribute("signal_process_vertex");
     }
 #else
     const HepMC::GenVertex * sigProcVtx = (*evt)->signal_process_vertex(); 
     if ( 0 != sigProcVtx ) { 
       const int sigProcBC = sigProcVtx->barcode(); 
       bool isInColl = false; 
       for ( HepMC::GenEvent::vertex_const_iterator itrVtx = (*evt)->vertices_begin(); 
         itrVtx != (*evt)->vertices_end(); 
         ++itrVtx ) { 
     if ( sigProcBC == (*itrVtx)->barcode() ) { 
       isInColl = true; 
       break; 
     } 
       }  //> loop over vertices 
       if ( !isInColl ) { 
     (*evt)->set_signal_process_vertex(0); 
       } 
     }//> a signal_process_vertex has been setup 
 #endif
   }//> loop over GenEvent's 
  
   return StatusCode::SUCCESS;
 }

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

std::set<int> OldSpclMcFilterTool::m_barcodes

protected

List of particles (in fact their unique barcode) which have been labelled as "special".

Definition at line 103 of file OldSpclMcFilterTool.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_doEtIsolation

BooleanProperty TruthParticleFilterBaseTool::m_doEtIsolation

protectedinherited

Switch to compute or not the Et-isolations for TruthParticle (and their underlying HepMC::GenParticle).

Default is to not compute these Et-isolations (and save CPU)

Definition at line 94 of file TruthParticleFilterBaseTool.h.

◆ m_etaRange

DoubleProperty OldSpclMcFilterTool::m_etaRange

protected

Eta acceptance cut applied on all stable particles.

Definition at line 120 of file OldSpclMcFilterTool.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_includePartonShowers

BooleanProperty OldSpclMcFilterTool::m_includePartonShowers

protected

Switch to include or not parton showers.

Definition at line 128 of file OldSpclMcFilterTool.h.

◆ m_includeSimul

BooleanProperty OldSpclMcFilterTool::m_includeSimul

protected

Switch to include or not particles from detector simulation (Geant4)

Definition at line 124 of file OldSpclMcFilterTool.h.

◆ m_isolationTool

IsolTool_t TruthParticleFilterBaseTool::m_isolationTool

protectedinherited

Pointer to the TruthIsolationTool to be able to compute transverse energy isolations for various isolation cones cuts.

See McParticleEvent/TruthParticleParameters.h for cone cuts.

Definition at line 85 of file TruthParticleFilterBaseTool.h.

◆ m_mcEventsOutputWriteHandleKey

SG::WriteHandleKey<McEventCollection> TruthParticleFilterBaseTool::m_mcEventsOutputWriteHandleKey {this,"McEventsOutput","GEN_AOD","WriteHandleKey for the output McEventCollection which has been filtered"}

protectedinherited

Location of the output McEventCollection which has been filtered.

Definition at line 104 of file TruthParticleFilterBaseTool.h.

◆ m_mcEventsReadHandleKey

SG::ReadHandleKey<McEventCollection> TruthParticleFilterBaseTool::m_mcEventsReadHandleKey {this,"McEvents","TruthEvent","ReadHandleKey for input McEventCollection one wants to filter"}

protectedinherited

ReadHandleKey for the input McEventCollection one wants to filter.

Definition at line 100 of file TruthParticleFilterBaseTool.h.

◆ m_mcVtxFilterTool

McVtxFilterTool_t TruthParticleFilterBaseTool::m_mcVtxFilterTool

protectedinherited

Pointer to the McVtxFilterTool to be able to select additional vertices on some decay pattern criterion.

Definition at line 78 of file TruthParticleFilterBaseTool.h.

◆ m_ptGamMin

DoubleProperty OldSpclMcFilterTool::m_ptGamMin

protected

Minimum threshold for transverse momentum of photons.

Definition at line 111 of file OldSpclMcFilterTool.h.

◆ m_ptMin

DoubleProperty OldSpclMcFilterTool::m_ptMin

protected

Minimum threshold for transverse momentum for all particles.

Warning: this cut is applied before Pt photon cut !

Definition at line 116 of file OldSpclMcFilterTool.h.

◆ m_removeDecayToSelf

BooleanProperty OldSpclMcFilterTool::m_removeDecayToSelf

protected

Switch to remove particles which decay into themselves (t->tg) but only for generated particles, not the ones from the Geant4 interactions.

Definition at line 137 of file OldSpclMcFilterTool.h.

◆ m_rIsol

DoubleProperty OldSpclMcFilterTool::m_rIsol

protected

deltaR isolation energy cut for electrons, muons, taus and photons

Definition at line 107 of file OldSpclMcFilterTool.h.

◆ m_tesIO

TruthHelper::GenAccessIO* OldSpclMcFilterTool::m_tesIO

protected

Pointer to a StoreGateSvc helper (Mc particle retrieval with predicates)

Definition at line 132 of file OldSpclMcFilterTool.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

Static Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ IsolTool_t

◆ McVtxFilterTool_t

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ OldSpclMcFilterTool() [1/2]

◆ ~OldSpclMcFilterTool()

◆ OldSpclMcFilterTool() [2/2]

Member Function Documentation

◆ buildMcAod()

◆ 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()

◆ finalizeTool()

◆ initialize()

◆ initializeTool()

◆ inputHandles()

◆ interfaceID()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ rebuildLinks()

◆ reconnectParticles()

◆ renounce()

◆ renounceArray()

◆ selectSpclMcBarcodes()

◆ shapeGenEvent()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_barcodes

◆ m_detStore

◆ m_doEtIsolation

◆ m_etaRange

◆ m_evtStore

◆ m_includePartonShowers

◆ m_includeSimul

◆ m_isolationTool

◆ m_mcEventsOutputWriteHandleKey

◆ m_mcEventsReadHandleKey

◆ m_mcVtxFilterTool

◆ m_ptGamMin

◆ m_ptMin

◆ m_removeDecayToSelf

◆ m_rIsol

◆ m_tesIO

◆ m_varHandleArraysDeclared

◆ m_vhka