---

Glast.cxx

Go to the documentation of this file.

// $Id: Glast.cxx,v 1.3 2000/11/22 16:19:16 jrb Exp $

#include "instrument/Glast.h"

#include "instrument/Tower.h"
#include "instrument/MCTruth.h"
#include "instrument/Scintillator.h"
#include "xml/IFile.h"
#include "xml/IFileManager.h"

#include <strstream>
#include <fstream>
#include <cassert>

int Glast::xNum = 1;
int Glast::yNum = 1;
double Glast::wall_gap = 0;
double Glast::ke_cutoff = .0005;
double Glast::zModDepth;
double Glast::zModCntr;
double Glast::Width;

Glast *Glast::s_instance = 0;

const GlastDetector::_PersistKey
Glast::classPersistKey = GlastDetector::_PersistKey("Glast");

namespace {
    
  // register with the IFileManager
  bool glastLoadParams = xml::IFileManager::loadParams(Glast::loadParameters);
  bool glastpersistValid = 
    GlastDetector::addPersistence(Glast::classPersistKey, new Glast::_Factory);
    
};


// constructor
Glast::Glast () : GlastDetector() {
  s_instance = this;
}
// destructor
Glast::~Glast () {
  s_instance = 0;
}

// persistence constructor - XML (DOM api)
Glast::Glast (const DOM_Element& elem) : GlastDetector(elem) {
  s_instance = this;
}

// "persist" using XML (DOM api)
// DOM_Element&    Glast::persist (DOM_Element& parent) const {
//  DOM_Element& selem = GlastDetector::persist(parent);
DOM_Element    Glast::persist (DOM_Element& parent) const {
  DOM_Element selem = GlastDetector::persist(parent);
    
  return selem;
}

//  setup parameters, done once
void Glast::loadParameters(xml::IFile& ini)
{
  // get basic parameters from the ini file
  xNum = ini.getInt("glast","xNum"),
    yNum = ini.getInt("glast","yNum");
  wall_gap=   ini.getDouble("walls", "gap");
  ke_cutoff=  ini.getDouble("glast", "ke_cutoff", .0005);
    
}

void
Glast::printOn (std::ostream& out) const
{
  GlastDetector::printOn (out);
}

void Glast::accept (DetectorVisitor& v)
{
  m_visit_flag = true;
  v.visit(this); // allows visitor to initialize, since this would be first
  v.visit(&m_truth);  // allow visitor to get the MCTruth

  // now go thru children unless aborted
  if(m_visit_flag) GlastDetector::accept(v);
}

void Glast::write(std::ostream & os) const
{
  GlastDetector::write (os);
}

void Glast::read(std::istream & is)
{
  GlastDetector::read (is);
}

void
Glast::visit(Tower* t) 
{
  addChild(t); // keep track of towers
}

void
Glast::visit(Scintillator* s)
{
  addChild (s); // keep track of scintillators
}

Glast::iterator Glast::begin () {
  return iterator(begin());
}
Glast::iterator Glast::end () {
  return iterator(end());
}
Glast::const_iterator Glast::begin () const {
  return const_iterator(begin());
}
Glast::const_iterator Glast::end () const {
  return const_iterator(end());
}
const Tower*    Glast::tower (const unsigned int& i) const {
  if (i >= nChildren())   return 0;
  return static_cast<const Tower*>(*(begin() + i));
}

---