---

GalElSpectrum.h

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/FluxSvc/src/GalElSpectrum.h,v 1.7 2002/07/25 05:18:58 srobinsn Exp $
// Original author: P. L. Nolan, pln@egret1.Stanford.EDU

/*
Notation: probably one of those things is called "flux" and the other
is something else.  Can someone look it up?
*/

#ifndef GAL_EL_SPECTRUM_H
#define GAL_EL_SPECTRUM_H

#include "Spectrum.h"
#include "CHIMESpectrum.h"
#include "facilities/Observer.h"
#include <string>

class GalElSpectrum : public Spectrum
{
    
public:
    GalElSpectrum(const std::string& params);
    virtual ~GalElSpectrum();
    
    virtual double calculate_rate(double old_rate);
    
    virtual double flux(double) const;
    
    virtual double solidAngle()const;
    
    float flux(float cut) const;
    
    virtual float flux(float lat, float lon) const;
    virtual float flux(std::pair<double, double> coords) const;
    
    virtual float operator() (float)const;
    
    float cutoff () const {return m_cutoff;};
    
    virtual void setPosition(float lat, float lon);
    virtual void setPosition(std::pair<double,double> coords);
    
    int askGPS();
    
    float findCutoff(float rflux) const;
    
    float findCutoff(float lat, float lon) const;
    float findCutoff(std::pair<double,double> coords) const;
    
    virtual std::pair<float,float> dir(float energy)const;
    
    
    virtual std::string title() const;
    virtual const char * particleName() const;
    inline  const char * nameOf() const {return "GalElSpectrum";}
    //   use default destructor, copy constructor, and assignment op.
    
    void setParticleName(std::string name);
    
private:
    CHIMESpectrum m_pspec;
    void init(float lat, float lon);
    float m_expo;   // exponent of power-law spectrum (integral)
    float m_norm;
    float m_normfact; // fraction of sky not blocked by earth
    float m_cutoff; // current cutoff energy
    float m_coscutoff;  // zenith angle of horizon
    float m_flux;   // current flux (set when cutoff changes)
    
    static const float m_rearth;    // radius of earth in km
    static const float m_altitude;  // altitude of circular orbit
    
    std::string m_particle_name;
    ObserverAdapter< GalElSpectrum > m_observer; //obsever tag
};

#endif // GAL_EL_SPECTRUM_H

---