---

McVertex.cxx

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/GlastEvent/src/MonteCarlo/McVertex.cxx,v 1.3 2001/01/29 23:20:47 igable Exp $

#include <iostream>
#include "GlastEvent/MonteCarlo/McVertex.h"
#include "GlastEvent/Utilities/CLHEPStreams.h"


// FIXME!!:
// The next `using' directive is an ad-hoc declaration for the transition
// phase to the namespace `GlastEvent'.
// After the transition is completed, it should be removed and each
// function definition should have the namespace identifier.
using namespace GlastEvent;

StreamBuffer& McVertex::serialize( StreamBuffer& s ) const
{
  int tmp = m_vertexType;
  ContainedObject::serialize(s);
  return s
    << m_subEvtID
    << m_initialPosition
    << m_finalPosition
    << m_timeOfFlight
    << tmp // m_vertexType
    << m_initialFourMomentum
    << m_finalFourMomentum
    << m_mcParticle(this)
    << m_motherMcParticle(this)
    << m_daughterMcParticles(this);
}

StreamBuffer& McVertex::serialize( StreamBuffer& s )
{
  int tmp;
  ContainedObject::serialize(s);
  s
    >> m_subEvtID
    >> m_initialPosition
    >> m_finalPosition
    >> m_timeOfFlight
    >> tmp // m_vertexType
    >> m_initialFourMomentum
    >> m_finalFourMomentum
    >> m_mcParticle(this)
    >> m_motherMcParticle(this)
    >> m_daughterMcParticles(this);

  m_vertexType = originType(tmp);
  return s;
}

std::ostream& McVertex::fillStream( std::ostream& s ) const
{
  s << "class McVertex"
    << "  (SubEvent:" << m_subEvtID << ")"
    << " :"
    << "\n    initialPosition (x, y, z) = ( "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_initialPosition.x() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_initialPosition.y() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_initialPosition.z() << " )"
    << "\n    finalPosition (x, y, z) = ( "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_finalPosition.x() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_finalPosition.y() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_finalPosition.z() << " )"
    << "\n    Time of flight     = "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_timeOfFlight
    << "\n    Initial 4-momentum (px, py, pz, E) = ( "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_initialFourMomentum.px() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_initialFourMomentum.py() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_initialFourMomentum.pz() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_initialFourMomentum.e()  << " )"
    << "\n    Final 4-momentum (px, py, pz, E) = ( "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_finalFourMomentum.px() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_finalFourMomentum.py() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_finalFourMomentum.pz() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_finalFourMomentum.e()  << " )"
    << "\n    Vertex Type           = " << int(m_vertexType)
    << "\n    Pair McParticle       = " << m_mcParticle(this)
    << "\n    Mother McParticle     = " << m_motherMcParticle(this)
    << "\n    Daughter McParticles  = ";
    SmartRefVector<McParticle>::const_iterator it;
    for (it = m_daughterMcParticles.begin(); it != m_daughterMcParticles.end(); it++){
      s << *it << "\n                            ";
    }
  return s;
}


const HepPoint3D& McVertex::initialPosition () const
{
    return m_initialPosition;
}
HepPoint3D& McVertex::initialPosition ()
{
    return m_initialPosition;
}
void McVertex::setInitialPosition (const HepPoint3D& value)
{
    m_initialPosition = value;
}

const HepPoint3D& McVertex::finalPosition () const
{
    return m_finalPosition;
}
HepPoint3D& McVertex::finalPosition ()
{
    return m_finalPosition;
}
void McVertex::setFinalPosition (const HepPoint3D& value)
{
    m_finalPosition = value;
}

double McVertex::timeOfFlight () const
{
    return m_timeOfFlight;
}
void McVertex::setTimeOfFlight (double value)
{
    m_timeOfFlight = value;
}
McVertex::originType McVertex::vertexType () const
{
    return m_vertexType;
}
void McVertex::setVertexType (McVertex::originType value)
{
    m_vertexType = value;
}


const HepLorentzVector& McVertex::initialFourMomentum() const
{
    return m_initialFourMomentum;
}
HepLorentzVector& McVertex::initialFourMomentum()
{
    return m_initialFourMomentum;
}
void McVertex::setInitialFourMomentum( const HepLorentzVector& value )
{
    m_initialFourMomentum = value;
}


const HepLorentzVector& McVertex::finalFourMomentum() const
{
    return m_finalFourMomentum;
}
HepLorentzVector& McVertex::finalFourMomentum()
{
    return m_finalFourMomentum;
}
void McVertex::setFinalFourMomentum( const HepLorentzVector& value )
{
    m_finalFourMomentum = value;
}


const McParticle* McVertex::mcParticle() const
{
    return m_mcParticle;
}
McParticle* McVertex::mcParticle()
{
    return m_mcParticle;
}
void McVertex::setMcParticle( McParticle* value )
{
    m_mcParticle = value;
}
void McVertex::setMcParticle( SmartRef<McParticle> value )
{
    m_mcParticle = value;
}


const McParticle* McVertex::motherMcParticle() const
{
    return m_motherMcParticle;
}
McParticle* McVertex::motherMcParticle()
{
    return m_motherMcParticle;
}

void McVertex::setMotherMcParticle( McParticle* value )
{
    m_motherMcParticle = value;
}
void McVertex::setMotherMcParticle( SmartRef<McParticle> value )
{
    m_motherMcParticle = value;
}

const SmartRefVector<McParticle>& McVertex::daughterMcParticles() const
{
    return m_daughterMcParticles;
}
SmartRefVector<McParticle>& McVertex::daughterMcParticles()
{
    return m_daughterMcParticles;
}

void McVertex::setDaughterMcParticles( const SmartRefVector<McParticle>& value )
{
    m_daughterMcParticles = value;
}

void McVertex::removeDaughterMcParticles()
{
    m_daughterMcParticles.clear();
}

void McVertex::addDaughterMcParticle( McParticle* value )
{
    m_daughterMcParticles.push_back(value);
}
void McVertex::addDaughterMcParticle( SmartRef<McParticle> value )
{
    m_daughterMcParticles.push_back(value);
}


short McVertex::subEvtID() const
{
    return m_subEvtID;
}
void McVertex::setSubEvtID( short value )
{
    m_subEvtID = value;
}

---