sim::EmEveIdCalculator Class Reference

#include <EmEveIdCalculator.h>

Inheritance diagram for sim::EmEveIdCalculator:

Public Member Functions
	EmEveIdCalculator ()
virtual	~EmEveIdCalculator ()
Public Member Functions inherited from sim::EveIdCalculator
	EveIdCalculator ()
	Constructor and destructor. More...
virtual	~EveIdCalculator ()
void	Init (const sim::ParticleList *list)
	Initialize this calculator for a particular ParticleList. More...
const sim::ParticleList *	ParticleList () const
int	CalculateEveId (const int trackID)

Private Member Functions
virtual int	DoCalculateEveId (const int trackID)

Additional Inherited Members
Protected Attributes inherited from sim::EveIdCalculator
const sim::ParticleList *	m_particleList

Detailed Description

Definition at line 63 of file EmEveIdCalculator.h.

Constructor & Destructor Documentation

sim::EmEveIdCalculator::EmEveIdCalculator ( )

inline

Definition at line 67 of file EmEveIdCalculator.h.

      : EveIdCalculator()   // Make sure the parent class constructor is called
    {}

virtual sim::EmEveIdCalculator::~EmEveIdCalculator ( )

inlinevirtual

Definition at line 70 of file EmEveIdCalculator.h.

{}

Member Function Documentation

int sim::EmEveIdCalculator::DoCalculateEveId ( const int  trackID )

privatevirtual

This is the core method to calculate the eve ID. If another class is going to override the default calculation, this the method that must be implemented.

Reimplemented from sim::EveIdCalculator.

Definition at line 22 of file EmEveIdCalculator.cxx.

  {
    // Almost any eve ID calculation will use this: Get the entire
    // history of the particle and its ancestors in the simulated
    // event. (m_particleList is defined in EveIdCalculator.h)
    const sim::ParticleHistory particleHistory( m_particleList, trackID );

    // You can treat particleHistory like an array, and do something
    // like:

    // for ( int i = particleHistory.size(); i >= 0; --i ) 
    // { const simb::MCParticle* particle = particleHistory[i]; ... }

    // But we know how to use the Standard Template Library (Yes, you
    // do! Don't doubt yourself!) so let's treat particleHistory in
    // the most efficient manner, as an STL container. We want to scan
    // the container from its last element to its first, from the base
    // of the particle production chain towards the primary particle.

    for(auto i = particleHistory.rbegin(); i != particleHistory.rend(); ++i){
      // Get the text string that describes the process that created
      // the particle.
      std::string process = (*i)->Process();

      // Skip it if it was created by pair production, compton
      // scattering, photoelectric effect, bremstrahlung,
      // annihilation, or any ionization. (The ultimate source of
      // the process names are the physics lists used in Geant4.)
      
      if ( process.find("conv")              != std::string::npos ||
           process.find("LowEnConversion")   != std::string::npos ||
           process.find("Pair")              != std::string::npos ||
           process.find("compt")             != std::string::npos ||
           process.find("Compt")             != std::string::npos ||
           process.find("Brem")              != std::string::npos ||
           process.find("phot")              != std::string::npos ||
           process.find("Photo")             != std::string::npos ||
           process.find("Ion")               != std::string::npos ||
           process.find("annihil")           != std::string::npos) continue;         

      // If we get here, the particle wasn't created by any of the
      // above processes. Return its ID.
      return (*i)->TrackId();
    }

    // If we get here, we've skipped every particle in the
    // chain. Perhaps it was empty.
    return 0;
  }

---

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

• nutools/old/ParticleNavigation/EmEveIdCalculator.h
• nutools/old/ParticleNavigation/EmEveIdCalculator.cxx