ReceptionModel.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2004 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/Polarimetry/Pulsar/ReceptionModel.h
 
#ifndef __ReceptionModel_H
#define __ReceptionModel_H
 
#include "Pulsar/MeasurementEquation.h"
#include "MEAL/SumRule.h"
#include "Estimate.h"
 
namespace Calibration {
 
  class CoherencyMeasurementSet;


 
  class ReceptionModel : public MeasurementEquation
  {
 
  public:

    ReceptionModel ();

    ~ReceptionModel ();
 
    // ///////////////////////////////////////////////////////////////////
    //
    // Add data to the model
    //
    // ///////////////////////////////////////////////////////////////////

    virtual void add_data (CoherencyMeasurementSet& data);

    unsigned get_ndata () const;

    const CoherencyMeasurementSet& get_data (unsigned idata) const;
 
    // ///////////////////////////////////////////////////////////////////
    //
    // Remove data from the model
    //
    // ///////////////////////////////////////////////////////////////////

    void delete_data ();

    void erase_input (unsigned index) override;

    void erase_transformation (unsigned index) override;
 
    // ///////////////////////////////////////////////////////////////////
    //
    // Report on the model
    //
    // ///////////////////////////////////////////////////////////////////
 
    class Report;

    void add_prefit_report (Report*);

    void add_postfit_report (Report*);
 
    // ///////////////////////////////////////////////////////////////////
    //
    // Fit the model
    //
    // ///////////////////////////////////////////////////////////////////
 
    class Solver;

    static Solver* new_default_Solver ();
    
    void set_solver (Solver*);

    Solver* get_solver ();
    const Solver* get_solver () const;

    void solve ();

    bool get_solved () const;

    void copy_fit (const ReceptionModel*);
 
    // ///////////////////////////////////////////////////////////////////
    //
    // Model implementation
    //
    // ///////////////////////////////////////////////////////////////////

    std::string get_name () const;
 
  protected:

    void range_check (unsigned idata, const char* method) const;

    Reference::To<Solver> solver;

    std::vector< Reference::To<Report> > prefit_reports;

    std::vector< Reference::To<Report> > postfit_reports;
 
 private:

    friend class Solver;

    std::vector<CoherencyMeasurementSet> data;
 
  };
 
 
}
 
#endif