DispersionDelay.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2004 - 2024 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/General/Pulsar/DispersionDelay.h
 
#ifndef __Pulsar_DispersionDelay_H
#define __Pulsar_DispersionDelay_H
 
#include "Estimate.h"
 
namespace Pulsar {




  const double dispersion_constant = 2.41e-4;
 
  class DispersionDelay {
 
  public:

    DispersionDelay ();

    void set_dispersion_measure (const double dispersion_measure);
    double get_dispersion_measure () const;

    void set_reference_frequency (double MHz);
    double get_reference_frequency () const;

    void set_frequency (double MHz);
    double get_frequency () const;

    void set_reference_wavelength (double metres);
    double get_reference_wavelength () const;

    void set_wavelength (double metres);
    double get_wavelength () const;

    double evaluate () const;

    typedef double Return;
    void set_measure (double dispersion_measure);
    double get_measure () const;
 
    // given an estimate of the slope \f$ dt / d lambda^2 \f$ in \f$ s m^{-2} \f$, return the DM
    Estimate<double> get_dispersion_measure(const Estimate<double>& slope) const;
 
  protected:

    double dispersion_measure;

    double reference_wavelength;

    double wavelength;
 
  };
 
}
 
#endif