---

CHIMESpectrum.h

Go to the documentation of this file.

// $Id: CHIMESpectrum.h,v 1.9 2001/10/18 03:30:16 srobinsn Exp $


// File: CHIMESpectrum.h
//
//
//  Patrick Nolan, Stanford University, 1998
//
#ifndef CHIME_SPECTRUM_H
#define CHIME_SPECTRUM_H

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


class CHIMESpectrum : public Spectrum
{
  
 public:
    CHIMESpectrum(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 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(double lat, double 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 "CHIMESpectrum";}
    //   use default destructor, copy constructor, and assignment op.
    
    void setParticleName(std::string name);
 

    typedef std::pair<int,float> Intrp;

 protected:
    void init(std::string params);
 private:
    class InterpVec : public std::vector<float> {

    public:
        InterpVec();  // constructor
        Intrp search(float x) const;
        float interpolate(Intrp) const;
    };

    
    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)
    double /*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< CHIMESpectrum > m_observer; //obsever tag

};

class CHIMEMax : public CHIMESpectrum
{
 public:
     CHIMEMax();
     ~CHIMEMax(){}
};

class CHIMEMin : public CHIMESpectrum
{
 public:
     CHIMEMin();
     ~CHIMEMin(){}
};


class CHIMEAvg : public CHIMESpectrum
{
 public:
     CHIMEAvg();
     ~CHIMEAvg(){}
};


#endif // CHIME_SPECTRUM_H

---