TwoBitStats.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2003 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/Base/Extensions/Pulsar/TwoBitStats.h
 
#ifndef __TwoBitStats_h
#define __TwoBitStats_h
 
#include "Pulsar/ArchiveExtension.h"
 
namespace Pulsar {


 
  class TwoBitStats : public Pulsar::Archive::Extension {
 
  public:
    
    TwoBitStats ();

    TwoBitStats (const TwoBitStats& extension);

    const TwoBitStats& operator= (const TwoBitStats& extension);

    TwoBitStats* clone () const { return new TwoBitStats( *this ); }

    void resize (unsigned nsample, unsigned ndig);

    unsigned get_ndig () const;

    unsigned get_nsample () const;

    void set_threshold (float threshold);
    
    float get_threshold () const { return threshold; }

    void set_cutoff_sigma (float cutoff_sigma);
    
    float get_cutoff_sigma () const { return cutoff_sigma; }

    const std::vector<float>& get_histogram (unsigned idig) const;

    void set_histogram (const std::vector<float>&, unsigned idig);

    void zero ();

    double get_mean_Phi (unsigned idig) const;

    double get_distortion (unsigned idig) const;
 
  protected:

    unsigned ndig;

    unsigned nsample;

    float threshold;

    float cutoff_sigma;

    std::vector< std::vector<float> > histogram;

    void range_check (unsigned idig, const char* method) const;
 
  };
 
}
 
#endif