ExcisionUnpacker.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2008 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// dspsr/Kernel/Classes/dsp/ExcisionUnpacker.h
 
#ifndef __ExcisionUnpacker_h
#define __ExcisionUnpacker_h
 
#include "dsp/HistUnpacker.h"
#include "JenetAnderson98.h"
 
namespace dsp
{
  class WeightedTimeSeries;
  class Input;
 
  class ExcisionUnpacker: public HistUnpacker
  {
 
  public:
 
    ExcisionUnpacker (const char* name = "ExcisionUnpacker");
 
    void reserve ();
 
    void set_output (TimeSeries*);
 
    void match_resolution (Input*);
 
    unsigned get_resolution () const;
 
    virtual unsigned get_input_offset (unsigned idig) const;
 
    virtual unsigned get_input_incr () const;
 
    virtual unsigned get_output_incr () const;
 
    virtual void set_ndat_per_weight (unsigned ndat_per_weight);
 
    unsigned get_ndat_per_weight () const { return ndat_per_weight; }
 
    virtual void set_nstate (unsigned nstate) { set_ndat_per_weight (nstate); }
 
    virtual void set_cutoff_sigma (float cutoff_sigma);
    
    virtual float get_cutoff_sigma() const { return cutoff_sigma; }
 
    virtual void set_threshold (float threshold) { }
 
    virtual float get_threshold () const { return 0.0; }
 
    //
    //
    //
 
    float get_fraction_low () const { return ja98.get_mean_Phi(); }
 
    unsigned get_nlow_min() const { return nlow_min; }
 
    unsigned get_nlow_max() const { return nlow_max; }
 
  protected:
 
    Reference::To<WeightedTimeSeries> weighted_output;
 
    void set_limits ();
 
    virtual void build ();
 
    void unpack ();
 
    virtual void dig_unpack (const unsigned char* input_data, 
                             float* output_data,
                             uint64_t ndat,
                             unsigned long* hist,
                             unsigned* weights = 0,
                             unsigned nweights = 0) = 0;
 
    template<class U, class Iterator>
    void excision_unpack (U& unpack, Iterator& input,
                          float* output_data, uint64_t ndat,
                          unsigned long* hist,
                          unsigned* weights, unsigned nweights);
 
    float cutoff_sigma;
 
    unsigned nlow_min;
 
    unsigned nlow_max;
 
    JenetAnderson98 ja98;
 
    unsigned current_digitizer;
 
    void not_built ();
 
  private:
 
    bool built;
 
    unsigned ndat_per_weight;
 
  };
}
 
#endif