StandardCandles.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2004 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/General/Pulsar/StandardCandles.h
 
#ifndef __Pulsar_StandardCandles_h
#define __Pulsar_StandardCandles_h
 
#include "Pulsar/Config.h"
 
#include "sky_coord.h"
#include "Types.h"
 
#include <string>
#include <vector>
 
namespace Pulsar {
 
  class StandardCandles : public Reference::Able {
 
  public:

    static bool verbose;

    static Option<std::string> default_filename;

    static Option<float> on_radius;

    static Option<float> off_radius;

    StandardCandles ();
    
    StandardCandles (const std::string& filename);
    
    ~StandardCandles ();
    
    void unload (const std::string& filename);

    void load (const std::string& filename);
    
    unsigned size () const { return entries.size(); }
 
    std::string get_filename () const { return loaded_filename; }


    class Entry {
      
    public:
      
      std::vector<std::string> source_name;

      sky_coord position;

      double reference_frequency;

      double reference_flux;

      double spectral_index;

      std::vector<double> spectral_coeffs;

      Entry () { init(); }

      explicit Entry (std::string& str) { load(str); }

      ~Entry();
      
      void load (const std::string& str);

      void unload (std::string& str);
      
      bool matches (const std::string& name) const;

      double get_flux_mJy (double MHz);
 
    protected:
      
      void init ();
      
    };

    Entry match (const std::string& source, double MHz) const;

    Signal::Source guess (std::string& name, sky_coord& pos) const;

    Entry get_entry (const unsigned idx) const { return entries[idx]; }
 
  protected:
    
    std::vector<Entry> entries;

    std::string loaded_filename;
    
  };
 
 
}
 
#endif