---

main.cxx

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/instrument/src/test/main.cxx,v 1.7 2000/11/03 17:30:57 jrb Exp $

// Test program for the instrument classes

#include "instrument/Instrument.h"
#include "instrument/GlastDetector.h"
#include "instrument/DetectorConverter.h"
#include <iostream>

#include "instrument/CsIDetector.h"
#include "instrument/MCTruth.h"
#include "instrument/Scintillator.h"
#include "instrument/SiDetector.h"
#include "instrument/SiTracker.h"


extern void  FATAL ( const char* msg )
{ 
  std::cerr << "FATAL " <<  msg << std::endl;
}

extern void WARNING ( const char* msg )
{
  std::cerr << "WARNING " << msg << std::endl;
}
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//  Example visitor
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
class   ResponseConverter : public DetectorConverter {
public:
    
  void    forward ( const CsIDetector& d )
    {
      if ( !(d.empty()) ) {
        // write this guy out!
        //std::cout << "csi " << std::endl;
      }
    }
    
  void   forward ( const MCTruth& d )
    {
      if ( !(d.empty()) ) {
        // write out
        //std::cout << "MCTRUTH data!" << std::endl;
      }
    }
    
  void   forward ( const Scintillator& d )
    {
      if ( !(d.empty()) ) {
        // write
        //std::cout<< "ACD data!" << std::endl;
      }
    }
    
  void  forward ( const SiDetector& d )
    {
      if ( !(d.empty()) ) {
        // write out
        //std::cout << "we have some Si Data " << std::endl;
      }
    }
    
};

//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

int main(int argc, char* argv[]) {
    
  Instrument g;
    
  // initialize, by reading default ini and xml files
  int rc = g.initialize();
    
  if( rc) {
    return rc;
  }
  else std::cout << "Loaded " << g.detector_count()
                 << " elements from detector definition file " 
                 << g.xmlFileName() << std::endl;
    
  if( argc <2 ) { 
    std::cerr << "Need to specify an IRF file name!" << std::endl; return 1;
  }
    
  rc = g.openIRF(argv[1]);
    
  if (rc==0 ) { 
    std::cout << "Opened " << g.irfFileName() << " for reading" << std::endl;
  } else return rc;
    
  int irec=0;
  while( rc==0 ) {
    rc = g.readIRF();
    if(rc==0) std::cout << "Record " << (++irec) 
                        <<": found " << g.detectors_with_data() 
                        << " detectors with hits" << std::endl;
    
    // Look at the resulting detectors, brute force
    int count=0;
    for (GlastDetector::MemberList::const_iterator it = 
           GlastDetector::allbegin(); it !=GlastDetector::allend(); ++it) {
      const GlastDetector * d =(*it).second;
      if (!(d->empty())) {
        ++ count;
      }
    }
    // now use a visitor, defined above
    g.accept(ResponseConverter());
  }
    
    
  return 0;
}

---