HistUnpacker.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2005 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// dspsr/Kernel/Classes/dsp/HistUnpacker.h
 
#ifndef __HistUnpacker_h
#define __HistUnpacker_h
 
#include "dsp/Unpacker.h"
#include <vector>
 
namespace dsp {
 
  class HistUnpacker: public Unpacker
  {
 
  public:
 
    static bool keep_histogram;
 
    HistUnpacker (const char* name = "HistUnpacker");
 
    virtual ~HistUnpacker ();
 
    void combine (const Operation*);
 
    void reset ();
 
    virtual double get_optimal_variance ();
 
    virtual void set_ndig (unsigned ndig);
    virtual unsigned get_ndig () const;
 
    virtual unsigned get_ndim_per_digitizer () const;
 
    virtual void set_nstate (unsigned nstate);
    unsigned get_nstate () const { return nstate; }
 
    virtual unsigned get_output_offset (unsigned idig) const;
 
    virtual unsigned get_output_ipol (unsigned idig) const;
 
    virtual unsigned get_output_ichan (unsigned idig) const;
 
    virtual void get_histogram (std::vector<unsigned long>&, unsigned) const;
 
    template <typename T>
    void get_histogram (std::vector<T>& data, unsigned idig) const;
 
    unsigned long* get_histogram (unsigned idig, unsigned expect = 0);
    const unsigned long* get_histogram (unsigned idig) const;
 
    double get_histogram_mean (unsigned idig) const;
 
    unsigned long get_histogram_total (unsigned idig) const;
 
    void zero_histogram ();
 
  protected:
 
    virtual void set_default_ndig ();
 
    void set_nstate_internal (unsigned _nstate);
    unsigned get_nstate_internal () const;
 
  private:
 
    unsigned nstate;
 
    unsigned nstate_internal;
 
    unsigned ndig;
 
    std::vector< std::vector< unsigned long > > histograms;
 
    void resize ();
 
    bool resize_needed;
  };
  
}
 
template <typename T> void
dsp::HistUnpacker::get_histogram (std::vector<T>& data, unsigned idig) const
{
  std::vector<unsigned long> hist;
  get_histogram (hist, idig);
 
  data.resize( hist.size() );
 
  for (unsigned i=0; i<hist.size(); i++)
    data[i] = T(hist[i]);
}
 
#endif