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