VanVleck.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2008 by Paul Demorest
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
#ifndef __VanVleck
#define __VanVleck
 
#include <vector>
#include "NormalDistribution.h"


class VanVleck {
 
 public:

  VanVleck();

  void set_nlevel(int n);

  int get_nlevel() const { return nlevel; }

  void set_uniform_threshold();

  void set_threshold(float *thr);

  void set_uniform_levels();

  void set_canonical_levels();

  void set_levels(float *lev);

  double out_pow(double in_pow, double in_dc=0.0);

  double d_out_pow(double in_pow, double in_dc=0.0);

  double in_pow(double out_pow);


  void set_in_power(double px, double py);

  void correct_acf(std::vector<float> data);

  void correct_ccf(std::vector<float> data);
  
 protected:

  int nlevel;

  double pow_x; 
  double pow_y;

  std::vector<float> thresh;

  std::vector<float> levels;

  void compute_coeffs();

  double c1x, c1y, c3x, c3y;
 
};
 
#endif