PMTPulseRecoBase.h

Go to the documentation of this file.

#ifndef PMTPULSERECOBASE_H
#define PMTPULSERECOBASE_H

// STL
#include <vector>
#include <cmath>
#include <functional>
#include <numeric>
#include <iostream>

#include "OpticalRecoTypes.h"
namespace pmtana
{
  
  struct pulse_param{
  public:
    double peak, area, ped_mean, ped_sigma;
    double t_start, t_max, t_end;

    //for vic
    double t_cfdcross;
    
    pulse_param(){
      reset_param();
    }

    ~pulse_param(){}
  
    void reset_param(){
      area = 0;
      peak = -1;
      ped_mean = ped_sigma = -1;
      t_start = t_max = t_end = -1;
      t_cfdcross = -1;
    }
    
  };

  typedef std::vector<pmtana::pulse_param> pulse_param_array;

  class PMTPulseRecoBase {
    
  public:

    PMTPulseRecoBase(const std::string name="noname");

    virtual ~PMTPulseRecoBase();

    const std::string& Name() const;

    const bool Status() const;

    virtual void Reset();

    bool Reconstruct( const pmtana::Waveform_t&,
                      const pmtana::PedestalMean_t&,
                      const pmtana::PedestalSigma_t& );

    const pulse_param& GetPulse(size_t index=0) const;

    const pulse_param_array& GetPulses() const;

    size_t GetNPulse() const {return _pulse_v.size();};

  private:

    std::string _name;

    bool _status;
    
  protected:

    virtual bool RecoPulse( const pmtana::Waveform_t&,
                            const pmtana::PedestalMean_t&,
                            const pmtana::PedestalSigma_t&     ) = 0;
    
    pulse_param_array _pulse_v;

    pulse_param _pulse;

  protected:

    bool Integral   (const std::vector<short> &wf, double &result, size_t begin=0, size_t end=0) const;

    bool Derivative (const std::vector<short> &wf, std::vector<int32_t> &diff, size_t begin=0, size_t end=0) const;

    size_t Max(const std::vector<short> &wf, double &result, size_t begin=0, size_t end=0) const;

    size_t Min(const std::vector<short> &wf, double &result, size_t begin=0, size_t end=0) const;

  };

}
#endif
 // end of doxygen group