JenetAnderson98A5.h

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

class JenetAnderson98::EquationA5 : public Reference::Able {
 
 public:

  void set_nsamp (unsigned nsamp);

  unsigned get_nsamp () const { return losq.size(); }
  
  template<typename T>
  void set_A6 (const std::vector<T>& prob)
  {
    double sum = 0.0;
    for (unsigned i=0; i<prob.size(); i++)
      sum += prob[i];
 
    prob_Phi.resize( prob.size() );
    for (unsigned i=0; i<prob.size(); i++)
      prob_Phi[i] = prob[i] / sum;
 
    set_nsamp( prob.size() );
  }

  double evaluate (double mean_Phi);

  double invert (double sigma_hat);
 
 private:
 
  std::vector<double> losq;
  std::vector<double> hisq;
  std::vector<double> fact;
  double dA5_dmean_Phi;
 
  std::vector<double> prob_Phi;
 
};
 
#endif