---

MCTruth.cxx

Go to the documentation of this file.

// $Id: MCTruth.cxx,v 1.6 2001/08/23 06:05:27 burnett Exp $
//
//
//    Implementation of class MCTruth  functions
//

#include "instrument/MCTruth.h"

#include "instrument/DetectorVisitor.h"
#include "instrument/DetectorConverter.h"
#include "instrument/FluxGenerator.h"

// static variable implementations
MCTruth* MCTruth::s_instance;
const GlastDetector::_PersistKey MCTruth::classPersistKey = 
GlastDetector::_PersistKey("MCTruth");

namespace {
  bool MCTruthPersistFlag = 
    GlastDetector::addPersistence(
                                  MCTruth::classPersistKey, 
                                  new MCTruth::_Factory);

};


MCTruth::MCTruth (const DOM_Element& elem)
  : GlastDetector(elem) 
{ 
  s_instance = this; 
}


MCTruth::~MCTruth()
{}

void
MCTruth::accept(DetectorVisitor& v) {
  v.visit(this);
  GlastDetector::accept(v);
}

void
MCTruth::accept(DetectorConverter& c) {
  c.backward(*this);
}

void
MCTruth::accept(DetectorConverter& c) const {
  c.forward(*this);
}
// Some local routines to read in MCParticle objects
//
template< class T>
static inline T read_3_vector(std::istream& in, double f)
{
    double x, y, z;
    in >> x >> y >> z;
    return T(x/f,y/f,z/f);
}

    

void    MCTruth::ParticleData::read (std::istream& in)
{
        char name[10];
        double x,y,z;
        in >> name >> m_energy; 
        m_name=std::string(name);
#if 1
        in >> x >> y >> z; 
        m_pos=Point(x,y,z);
        in >> x >> y >> z;
        m_p = 1.e-4*Vector(x,y,z);
#else
    double one = 1.0; // keep Sun CC happy
    m_pos = read_3_vector<Point>( in, one );
    m_dir = read_3_vector<Vector>(in, 1.0e6);
    in >> m_code;
    in >> m_matZ;

    // now read in the particle data
    delete m_particles;
    m_particles = new ParticleData();
    m_particles->read (in);
#endif
}

void MCTruth::ParticleData::write ( std::ostream& out) const
{
    out << m_name << " " << static_cast<int>(m_energy) << " ";
    out << static_cast<int>(m_pos.x()+0.5) << " ";
    out << static_cast<int>(m_pos.y()+0.5) << " ";
    out << static_cast<int>(m_pos.z()+0.5) << " ";
    out << static_cast<int>(m_p.x()*1e4) << " ";
    out << static_cast<int>(m_p.y()*1e4) << " ";
    out << static_cast<int>(m_p.z()*1e4) << " ";

}

//default constructor: no particle?
MCTruth::MCTruth(_Particle* p)
  : m_particle(p? p : new FluxGenerator), m_id(0) 
{
  s_instance = this; 
}

---