---

HeSpectrum.h

Go to the documentation of this file.

// $Id: HeSpectrum.h,v 1.8 2001/10/20 07:20:34 srobinsn Exp $

// File: HeSpectrum.h
//
//
#ifndef HE_SPECTRUM_H
#define HE_SPECTRUM_H

#include "flux/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() 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

---