---

Spectrum.h

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/flux/flux/Spectrum.h,v 1.10 2001/10/20 07:20:35 srobinsn Exp $
//
//
// Spectrum: base class for energy spectrum objects
// SimpleSpectrum: define a particle and spectral index
//

#ifndef GLAST_SPECTRUM_H
#define GLAST_SPECTRUM_H

#include <string>
#include <utility>
#include <vector>

//forward declaration
class HepRandomEngine;


class Spectrum {
public:
    
    //    class  Direction : public std::pair<float,float> {
    //    public:
    //        double costh()const{return this.first;}
    //        double phi()const {return this.second;}
    //    };
    
    virtual float operator()(float /*r*/)const{return 0;};
    // this is all that these objects do. Must be virtual for
    // polymorphism
    // returns kinetic energy for random number r in [0,1). 
    // requried that it be monatonic
    // NB. Default is to return zero, an indicator that the actual Spectrum object
    //     implements a method that makes direct use of the random generator
    
    //virtual double calculate_rate (double old_rate) = 0;
    
    virtual const char * particleName()const=0;
    
    virtual double    flux ( ) const;
    
    
    virtual double solidAngle()const;
    
    virtual void        setPosition ( double lat, double lon , double time=0);
    
    virtual std::string title()const=0;
    
    float fraction(float energy);
    
    
    //   virtual Direction dir(float energy)const;
    // return solid angle pair (costh, phi) for the given energy
    
    virtual ~Spectrum();
    
    
    virtual std::pair<float,float> dir(float energy)const;
    
    virtual double energySrc(HepRandomEngine* engine);
    virtual std::pair<double,double> dir(double energy, HepRandomEngine* engine);
    
    
protected:
    Spectrum(const std::vector<float>& /*params*/){};
    Spectrum(/* double lat = 0, double lon = 0, double time=0*/) 
        : m_lat(0), m_lon(0), m_time(0){}
    // all constructors protected to ensure an abstract class
    
    virtual void parseParamList(std::string input, std::vector<float>& output) const;
    
    double    m_lat, m_lon;   // latitude and longitudinal coordinates
    double    m_time;
};

#endif    

---