ISF::GenParticleSimQuasiStableFilter Class Reference

#include <GenParticleSimQuasiStableFilter.h>

Inheritance diagram for ISF::GenParticleSimQuasiStableFilter:

Collaboration diagram for ISF::GenParticleSimQuasiStableFilter:

Public Member Functions
	GenParticleSimQuasiStableFilter (const std::string &t, const std::string &n, const IInterface *p)
	Constructor. More...
	~GenParticleSimQuasiStableFilter ()
	Destructor. More...
StatusCode	initialize ()
	Athena algtool's Hooks. More...
StatusCode	finalize ()
bool	pass (const HepMC::GenParticle &particle) const
	passes through to the private version More...

Private Attributes
std::vector< int >	m_statuscode = {1,2}
	returns true if the the particle and all daughters are on the accept list More...
std::vector< double >	m_minProdRadius = {30.19*Gaudi::Units::mm,0.0*Gaudi::Units::mm}
	minimum production radius below which particles should be ignored More...
std::vector< double >	m_minDecayRadius = {30.19*Gaudi::Units::mm,30.19*Gaudi::Units::mm}
	minimum decay radius below which particles should be ignored More...

Detailed Description

Quasi stable particle and their daughter filter for HepMC particles to be used in the stack filling process. Checks this particle and all daughters.

Definition at line 34 of file GenParticleSimQuasiStableFilter.h.

Constructor & Destructor Documentation

GenParticleSimQuasiStableFilter()

ISF::GenParticleSimQuasiStableFilter::GenParticleSimQuasiStableFilter	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Constructor.

Definition at line 19 of file GenParticleSimQuasiStableFilter.cxx.

  : base_class(t,n,p)
{
    // different options
    declareProperty("StatusCode", m_statuscode);
    declareProperty("MinProdRadius", m_minProdRadius);
    declareProperty("MinDecayRadius", m_minDecayRadius);
}

~GenParticleSimQuasiStableFilter()

ISF::GenParticleSimQuasiStableFilter::~GenParticleSimQuasiStableFilter ( )

inline

Destructor.

Definition at line 41 of file GenParticleSimQuasiStableFilter.h.

{}

Member Function Documentation

finalize()

StatusCode ISF::GenParticleSimQuasiStableFilter::finalize

(

)

Definition at line 74 of file GenParticleSimQuasiStableFilter.cxx.

{
    ATH_MSG_VERBOSE("Finalizing ...");
    return StatusCode::SUCCESS;
}

initialize()

StatusCode ISF::GenParticleSimQuasiStableFilter::initialize

(

)

Athena algtool's Hooks.

Definition at line 31 of file GenParticleSimQuasiStableFilter.cxx.

{
    ATH_MSG_VERBOSE("Initializing ...");
    
    m_minProdRadius.resize(m_statuscode.size(),0);
    m_minDecayRadius.resize(m_statuscode.size(),0);
 
    // All done!
    return StatusCode::SUCCESS;
}

pass()

bool ISF::GenParticleSimQuasiStableFilter::pass

(

const HepMC::GenParticle &

particle

)

const

passes through to the private version

Definition at line 45 of file GenParticleSimQuasiStableFilter.cxx.

{
   auto particle=&particleref;
#endif
  ATH_MSG_VERBOSE( "Checking whether " << particle << " passes the filter." );
  
  bool passFilter = false;
  for(unsigned int i=0;i<m_statuscode.size();++i) {
    if(particle->status()==m_statuscode[i]) {
      passFilter=true;
      if(particle->end_vertex()        && (m_minDecayRadius[i] > particle->end_vertex()->position().perp())        ) passFilter = false;
      if(particle->production_vertex() && (m_minProdRadius[i]  > particle->production_vertex()->position().perp()) ) passFilter = false;
      if(passFilter) break;
    }
  }
  
  if(passFilter) {
    ATH_MSG_DEBUG( "Accpeting particle " << particle );
    if(particle->production_vertex()) ATH_MSG_DEBUG("  prod r="<<particle->production_vertex()->position().perp()<<" : "<<particle->production_vertex ());
    if(particle->end_vertex()) ATH_MSG_DEBUG("  decay r="<<particle->end_vertex()->position().perp()<<" : "<<particle->end_vertex ());
  } else {
//      ATH_MSG_DEBUG( "Rejecting particle " << particle );
//      if(particle.production_vertex ()) ATH_MSG_DEBUG("  prod r="<<particle.production_vertex()->position().perp()<<", min="<<m_minProdRadius[i]<<" : "<<*(particle.production_vertex ()));
//      if(particle.end_vertex ()) ATH_MSG_DEBUG("  decay r="<<particle.end_vertex()->position().perp()<<", min="<<m_minDecayRadius[i]<<" : "<<*(particle.end_vertex ()));
  }
 
  return passFilter;
}

Member Data Documentation

m_minDecayRadius

std::vector< double > ISF::GenParticleSimQuasiStableFilter::m_minDecayRadius = {30.19*Gaudi::Units::mm,30.19*Gaudi::Units::mm}

private

minimum decay radius below which particles should be ignored

Definition at line 58 of file GenParticleSimQuasiStableFilter.h.

m_minProdRadius

std::vector< double > ISF::GenParticleSimQuasiStableFilter::m_minProdRadius = {30.19*Gaudi::Units::mm,0.0*Gaudi::Units::mm}

private

minimum production radius below which particles should be ignored

Definition at line 57 of file GenParticleSimQuasiStableFilter.h.

m_statuscode

std::vector< int > ISF::GenParticleSimQuasiStableFilter::m_statuscode = {1,2}

private

returns true if the the particle and all daughters are on the accept list

status codes to accept

Definition at line 56 of file GenParticleSimQuasiStableFilter.h.

---

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

• GenParticleSimQuasiStableFilter.h
• GenParticleSimQuasiStableFilter.cxx