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