---

HeSpectrum.h

Go to the documentation of this file.

// $Id: HeSpectrum.h,v 1.6 2002/07/23 19:00:57 srobinsn Exp $
#ifndef HE_SPECTRUM_H
#define HE_SPECTRUM_H

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


class HeSpectrum : public Spectrum
{
    
public:
    HeSpectrum(const std::string& params);
    
    virtual double calculate_rate(double old_rate);
    
    virtual double flux(double) const;
    
    virtual double solidAngle()const;
    
    float flux(float cut) const;
    
    
    virtual float operator() (float)const;
    
    float cutoff () const {return m_cutoff;};
    
    virtual void setPosition(double lat, double lon);
    
    int askGPS();
    
    float findCutoff(float rflux) const;
    
    float findCutoff(float lat, float lon) 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 "HeSpectrum";}
    //   use default destructor, copy constructor, and assignment op.
    
    void setParticleName(std::string name);
    
    typedef std::pair<int,float> Intrp;
private:
    double HeSpectrum::flux(double lat, double lon) const;
    
    class InterpVec : public std::vector<float> {
        
    public:
        InterpVec();  // constructor
        Intrp search(float x) const;
        float interpolate(Intrp) const;
    };
    
    void init(float lat, float lon);
    InterpVec m_en; // values of energy in the table
    InterpVec m_fl; // values of integral flux in the table
    InterpVec m_fluxes; // integral flux, unmodulated by earth's field
    float m_expo;   // exponent of power-law spectrum above table (integral)
    float m_normfact; // fraction of sky not blocked by earth
    float m_tot;    // total flux above the current cutoff energy
    float m_upper;  // total flux in the power-law part of spectrum
    float m_etop;   // energy at which table transitions to power law
    float m_cutoff; // current cutoff energy
    float m_coscutoff;  // zenith angle of horizon
    float m_fluxTbl[73][13];  // table of total flux vs. latitude and longitude
    float m_flux;   // current flux (set when cutoff changes)
    
    float rad2() const; // square of distance (in km) from dipole center
    float exposure(float E) const;  // geomagnetic cutoff factor
    float cosomega(float E) const;  // Opening angle of Størmer cone
    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< HeSpectrum > m_observer; //obsever tag
    
};

#endif // HE_SPECTRUM_H

---