/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Event/McTruth/McTruth-00-02-19/src/McParticle.cxx

Go to the documentation of this file.

#include <iostream>
#include "McTruth/McParticle.h"

namespace Event {

McParticle::StdHepId McParticle::particleProperty() const
{
  return m_particleID;
}

bool McParticle::primaryParticle() const
{
  return (m_statusFlags & PRIMARY)==PRIMARY;
}

bool McParticle::leafParticle() const
{
  return (m_statusFlags & LEAF)==LEAF;
}

bool McParticle::decayFromGenerator() const
{
  return (m_statusFlags & DECAYED)==DECAYED;
}

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

const McParticle& McParticle::mother()const
{
    return *m_mother;
}

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

---