---

main.cxx

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/reconstruction/src/test/main.cxx,v 1.2 2000/12/08 23:01:09 igable Exp $

// Test program for the reconstruction classes

#include "instrument/Instrument.h"

#include "reconstruction/data/RecoGlastData.h"
#include "reconstruction/GlastRecon.h"
#include "reconstruction/ReconData.h"
#include "reconstruction/GlastTuple.h"
#include "reconstruction/PrintReconData.h"

#include <iostream>

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

extern void WARNING ( const char* msg )
{
    std::cerr << "WARNING " << msg << 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 ) { 
        // get the path to the instrument from its root environment variable: if not there,
        // assume that we are in the root
        const char * local_path = ::getenv("INSTRUMENTROOT");
        
        std::string input_file = std::string(local_path? local_path: "");
        
        input_file +=  std::string("/src/test/muontest.irf");
        
        rc = g.openIRF(input_file);
        std::cout << "Opening test irf file "<< input_file << std::endl; 
    }else   {
        // open file in args
        rc = g.openIRF(argv[1]);
    }
    if (rc==0 ) { 
        std::cout << "Opened " << g.irfFileName() << " for reading" << std::endl;
    } else {
        std::cout << "Failed to open " << g.irfFileName() << " for reading" << std::endl;
        return rc;
    }
        // create the GlastData and GlastRecon objects
    RecoGlastData data;
    GlastRecon recon;
    
    
    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;
        
        data.clear();
        recon.clear();

        // fill  data from the IRF
        g.accept(data);
        
        // see what is there
        data.printOn(std::cout);
        
        //  pass data to the GlastRecon object.
        recon.accept(ReconData(&data));
        
        // print out the full tuple
        recon.accept(PrintReconData(std::cout));
        
        
    }
    
    
    return 0;
}

---