---

rootplot.h

Go to the documentation of this file.

// Flux test program that generates a ROOT macro to plot the flux
//

// $Header: /nfs/slac/g/glast/ground/cvs/FluxSvc/src/test/flux/rootplot.h,v 1.2 2002/05/01 22:33:06 srobinsn Exp $

// Original author: Theodore Hierath

#include "../../FluxMgr.h"
#include "FluxSvc/EventSource.h"
//#include "FluxSvc/ISpectrumFactory.h"
#include "../../SpectrumFactoryTable.h"
#include "rootEnergyHist.h"
#include "rootAngleHist.h"

#include <fstream>

class rootplot
{
public:
    
    const int NUM_BINS;
    const int LOOP;
    
    const double TIME;
    
    const double ENERGY_MIN;
    const double ENERGY_MAX;
    
    void help() {
        std::cout << 
            "   Simple test program for a particle source.\n"
            "   Command line args are \n"
            "      '-sum' sums up the histograms'\n"
            "      '-bins <number of bins for histogram>'\n"
            "      '-events <number of events to create>'\n"
            "      '-min' or '-energy_min' sets minimum energy\n"
            "      '-max' or '-energy_max' sets maximum energy\n"
            "      '-list' lists the available spectra\n"
            "      '-file <file name>' writes energy vs. flux to filename.txt\n"
            "      '-trueflux' graphs using flux using per steradian\n"
            "      '-flux' graphs the flux vs. E instead of E*flux vs E\n"
            "      '-flux_min' manually set lower limit for graph\n"
            "      '-flux_max' manually set upper limit for graph\n"
            "      '-graph <log | semilogx | semilogy | linear>'\n"
            "      '-longsrc <sourcename>' for long-term energy averaging\n"
            "      '-time <time in seconds>' for the flux at time\n"
            "      '-help' for this help"
            << std::endl;
    }
    
    void listSources(const std::list<std::string>& source_list ) {
        std::cout << "List of available sources:" << std::endl;
        for( std::list<std::string>::const_iterator it = source_list.begin(); 
        it != source_list.end(); ++it) { 
            std::cout << '\t'<< *it << std::endl;
        }
    }
    
    void listSpectra() {
        std::cout << "List of loaded Spectrum objects: " << std::endl;
        std::list<std::string> spectra(SpectrumFactoryTable::instance()->spectrumList());
        for( std::list<std::string>::const_iterator it = spectra.begin(); 
        it != spectra.end(); ++it) { 
            std::cout << '\t'<< *it << std::endl;
        }
    }
    
#define DLL_DECL_SPECTRUM(x)   extern const ISpectrumFactory& x##Factory; x##Factory.addRef();
    
    void flux_load() {
        
        // these are the spectra that we want to make available
        DLL_DECL_SPECTRUM( CHIMESpectrum);
        DLL_DECL_SPECTRUM( AlbedoPSpectrum);
        DLL_DECL_SPECTRUM( HeSpectrum);
        DLL_DECL_SPECTRUM( GalElSpectrum);   
//        DLL_DECL_SPECTRUM( CrElectron);
//        DLL_DECL_SPECTRUM( CrProton);
        DLL_DECL_SPECTRUM( FILESpectrum);
        //   DLL_DECL_SPECTRUM( GRBSpectrum);
        
    }
      
    
    rootplot(std::vector<char*> argv);

};

---