OldSpclMcFilterTool Class Reference

#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:
virtual	~OldSpclMcFilterTool ()
	Destructor:
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.
StatusCode	initialize ()
	Athena Algorithm's Hooks.
StatusCode	execute ()
StatusCode	finalize ()
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	sysInitialize () override
	Perform system initialization for an algorithm.
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

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

Protected Types
typedef ToolHandle< IMcVtxFilterTool >	McVtxFilterTool_t
typedef ToolHandle< ITruthIsolationTool >	IsolTool_t

Protected Member Functions
	OldSpclMcFilterTool ()
	Default constructor:
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.
StatusCode	shapeGenEvent (McEventCollection *outEvt)
	Build the out McEventCollection by removing all the GenParticles whose barcode is not in the list of barcodes to keep.
StatusCode	reconnectParticles (const McEventCollection *inEvt, McEventCollection *out)
	Reconnect the particles of the filtered McEventCollection.
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.
StatusCode	initializeTool ()
	to get tesIO
virtual StatusCode	finalizeTool ()
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.

Protected Attributes
std::set< int >	m_barcodes
	List of particles (in fact their unique barcode) which have been labelled as "special".
DoubleProperty	m_rIsol
	deltaR isolation energy cut for electrons, muons, taus and photons
DoubleProperty	m_ptGamMin
	Minimum threshold for transverse momentum of photons.
DoubleProperty	m_ptMin
	Minimum threshold for transverse momentum for all particles.
DoubleProperty	m_etaRange
	Eta acceptance cut applied on all stable particles.
BooleanProperty	m_includeSimul
	Switch to include or not particles from detector simulation (Geant4)
BooleanProperty	m_includePartonShowers
	Switch to include or not parton showers.
TruthHelper::GenAccessIO *	m_tesIO
	Pointer to a StoreGateSvc helper (Mc particle retrieval with predicates)
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.
McVtxFilterTool_t	m_mcVtxFilterTool
	Pointer to the McVtxFilterTool to be able to select additional vertices on some decay pattern criterion.
IsolTool_t	m_isolationTool
	Pointer to the TruthIsolationTool to be able to compute transverse energy isolations for various isolation cones cuts.
BooleanProperty	m_doEtIsolation
	Switch to compute or not the Et-isolations for TruthParticle (and their underlying HepMC::GenParticle).
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.
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.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

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

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

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

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

declareProperty()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

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.

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

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

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

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.

{this,"McEventsOutput","GEN_AOD","WriteHandleKey for the output McEventCollection which has been filtered"};

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.

{this,"McEvents","TruthEvent","ReadHandleKey for input  McEventCollection one wants to filter"};

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:

• OldSpclMcFilterTool.h
• OldSpclMcFilterTool.cxx