SimplePredictor.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2007 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/Util/resources/Pulsar/SimplePredictor.h
 
#ifndef __SimplePredictor_h
#define __SimplePredictor_h
 
#include "Pulsar/Predictor.h"
#include "sky_coord.h"
 
#include <string>
#include <vector>
 
 
namespace Pulsar {
 
  class SimplePredictor : public Predictor {
 
  public:
 
    static double precision;
 
    SimplePredictor (const char* filename = 0);
 
    ~SimplePredictor ();
 
    Predictor* clone () const;
 
    Generator* generator () const;
 
    void insert (const Pulsar::Predictor*);
 
    void keep (const std::vector<MJD>& epochs);
 
    bool matches (const Pulsar::Predictor*) const;
 
    void set_observing_frequency (long double MHz);
 
    long double get_observing_frequency () const;
 
    Pulsar::Phase phase (const MJD& t) const;
 
    MJD iphase (const Pulsar::Phase& phase, const MJD* guess) const;
 
    long double frequency (const MJD& t) const;
 
    Pulsar::Phase dispersion (const MJD &t, long double MHz) const;
 
    void load (FILE*);
 
    void unload (FILE*) const;
 
    std::string get_name () const { return name; }
 
    sky_coord get_coordinates () const { return coordinates; }
 
    double get_dispersion_measure () const { return dispersion_measure; }
 
    void set_reference_epoch (const MJD& epoch) { reference_epoch = epoch; }
 
    MJD get_reference_epoch () const { return reference_epoch; }
 
    MJD    get_reftime       () const { return reference_epoch; }
    Pulsar::Phase  get_refphase      () const { return 0.0; }
    double get_reffrequency  () const { return coefs.at(0); }
    
  private:
 
    std::vector<long double> coefs;
 
    MJD reference_epoch;
 
    std::string name;
 
    sky_coord coordinates;
 
    double dispersion_measure;
 
    std::string telescope;
 
    void parse(std::string key, std::string value);
 
  };
 
}
 
#endif