McParticle.cxx

Go to the documentation of this file.

#include "McTruth/McParticle.h"
#include <iostream>
 
namespace Event {
 
  /// Retrieve particle property
  McParticle::StdHepId McParticle::particleProperty() const { return m_particleID; }
 
  /// Retrieve whether this is a primary particle
  bool McParticle::primaryParticle() const { return ( m_statusFlags & PRIMARY ) == PRIMARY; }
 
  /// Retrieve whether this is a leave particle
  bool McParticle::leafParticle() const { return ( m_statusFlags & LEAF ) == LEAF; }
 
  /// Retrieve whether decayed from generator
  bool McParticle::decayFromGenerator() const {
    return ( m_statusFlags & DECAYED ) == DECAYED;
  }
 
  /// Retrieve wheter decayed in flight
  bool McParticle::decayInFlight() const { return ( m_statusFlags & DECAYFLT ) == DECAYFLT; }
 
  void McParticle::init( McParticle* mother, StdHepId id, unsigned int flags,
                         const HepLorentzVector& initialMomentum,
                         const HepLorentzVector& finalMomentum,
                         const HepLorentzVector& initialPosition,
                         const HepLorentzVector& finalPosition, const std::string process ) {
    initialize( mother, id, flags, initialMomentum, initialPosition, process );
    // finalize(finalMomentum, finalPosition);
    finalize( finalPosition );
  }
 
  void McParticle::initialize( StdHepId id, unsigned int flags,
                               const HepLorentzVector& initialMomentum,
                               const HepLorentzVector& initialPosition,
                               const std::string       process ) {
    m_particleID          = id;
    m_statusFlags         = flags;
    m_initialFourMomentum = initialMomentum;
    m_initialPosition     = initialPosition;
    m_process             = process;
  }
 
  void McParticle::initialize( McParticle* mother, StdHepId id, unsigned int flags,
                               const HepLorentzVector& initialMomentum,
                               const HepLorentzVector& initialPosition,
                               const std::string       process ) {
    m_mother              = mother;
    m_particleID          = id;
    m_statusFlags         = flags;
    m_initialFourMomentum = initialMomentum;
    m_initialPosition     = initialPosition;
    m_process             = process;
    if ( mother != this ) mother->m_daughters.push_back( this );
  }
 
  /*void McParticle::finalize(const HepLorentzVector& finalMomentum,
          const HepLorentzVector& finalPosition)
  {
      m_finalFourMomentum = finalMomentum;
      m_finalPosition = finalPosition;
  }*/
 
  void McParticle::finalize( const HepLorentzVector& finalPosition ) {
    m_finalPosition = finalPosition;
  }
 
  void McParticle::setMother( const SmartRef<McParticle> m ) { m_mother = m; }
 
  const HepLorentzVector& McParticle::initialPosition() const { return m_initialPosition; }
  const HepLorentzVector& McParticle::finalPosition() const { return m_finalPosition; }
  const HepLorentzVector& McParticle::initialFourMomentum() const {
    return m_initialFourMomentum;
  }
  /*const HepLorentzVector&  McParticle::finalFourMomentum()const
  {
      return m_finalFourMomentum;
  }*/
 
  unsigned int McParticle::statusFlags() const { return m_statusFlags; }
 
  /// access to the mother particle
  const McParticle& McParticle::mother() const { return *m_mother; }
 
  /// Remove daughters when in prune mode
  void McParticle::removeDaughter( const SmartRef<McParticle> mcPart ) {
    SmartRefVector<Event::McParticle>::iterator daughtIter;
    for ( daughtIter = m_daughters.begin(); daughtIter != m_daughters.end(); daughtIter++ )
    {
      if ( mcPart == *daughtIter )
      {
        m_daughters.erase( daughtIter );
        break;
      }
    }
    return;
  }
 
} // namespace Event