JenetAnderson98.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2006 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/Util/genutil/JenetAnderson98.h
 
#ifndef __Jenet_Anderson_98
#define __Jenet_Anderson_98
 
#include "Reference.h"
#include <vector>


class JenetAnderson98 : public Reference::Able {
 
 public:


  static double get_optimal_spacing (unsigned bits);

  JenetAnderson98 ();

  void set_Phi (double Phi);

  void set_sigma_n (double sigma_n);

  double get_hi () const { return hi; }

  double get_lo () const { return lo; }

  double get_A () const { return A; }

  void set_threshold ( double t = get_optimal_spacing(2) );
  double get_threshold () const { return threshold; }

  double get_mean_Phi () const { return mean_Phi; }
  void set_mean_Phi (double mean_Phi);

  double get_var_Phi () const { return var_Phi; }

  void get_prob_Phi (unsigned L, std::vector<float>& prob_Phi);

  double A4 (double Phi);

  double invert_A4 (double sigma_hat);

  double A14 (double mean_Phi);

  double invert_A14 (double sigma_hat);

  class Probability;

  class Plot;

  class EquationA5;

  void set_measured_prob_Phi (const Probability* data);

  const Probability* get_measured_prob_Phi () const;
 
 
 protected:

  void set_alpha (double alpha, double Phi);

  double threshold;

  Reference::To<const Probability> measured;
 
  double lo;
  double hi;
  double A;
  double mean_Phi;
  double var_Phi;
 
};

class JenetAnderson98::Probability : public Reference::Able  {
  
public:
  
  virtual unsigned get_nsample() const = 0;
  
  virtual float get_threshold () const = 0;
  
  virtual float get_cutoff_sigma () const = 0;
  
  virtual unsigned get_ndig () const = 0;

  virtual void get_histogram (std::vector<float>& h, unsigned dig) const = 0;

  double evaluate (unsigned idig, double Phi, double* dprob_dPhi = 0);
 
private:
  unsigned last_idig;
  std::vector<float> last_hist;
 
};
 
#endif