DynamicSpectrum.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2010 by Paul Demorest
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
#ifndef __Pulsar_DynamicSpectrum_h
#define __Pulsar_DynamicSpectrum_h
 
#include "Pulsar/Algorithm.h"
#include "Estimate.h"
#include "MJD.h"
 
namespace Pulsar
{
  class Archive;
  class Flux;

  class DynamicSpectrum : public Algorithm
  {   
  public:

    DynamicSpectrum();

    ~DynamicSpectrum();

    void set_Archive(const Archive* _arch);

    void set_flux_method(Flux *flux);

    Reference::To<Flux> get_flux_method() { return flux; }

    void compute();

    void reset();

    Estimate<double> get_ds(unsigned ichan, unsigned isub);

    const double *get_ds();

    const double *get_ds_err();

    const unsigned get_nchan() const { return nchan; }

    const unsigned get_nsub() const { return nsub; }

    const MJD get_rel_epoch(unsigned isub);

    const double get_freq(unsigned ichan, unsigned isub=0);

    void unload(const std::string& filename, const std::string& command="");

    static DynamicSpectrum *load(const std::string& filename);
 
  protected:

    Reference::To<const Archive> arch;

    Reference::To<Flux> flux;

    bool computed;

    unsigned nchan;

    unsigned nsub;

    double *ds_data;

    double *ds_data_err;

    double ds_max;

    double ds_min;
 
  };
 
}
 
 
#endif