PMTPulseRecoBase.h

Go to the documentation of this file.

#ifndef PMTPULSERECOBASE_H
#define PMTPULSERECOBASE_H

// STL
#include <string>
#include <vector>

#include "OpticalRecoTypes.h"
#include "larana/OpticalDetector/OpHitFinder/RiseTimeTools/RiseTimeCalculatorBase.h"

#include <memory>

namespace pmtana {

  struct pulse_param {
  public:
    double peak, area, ped_mean, ped_sigma;
    double t_start, t_max, t_end, t_rise;

    //for vic
    double t_cfdcross;

    pulse_param() { reset_param(); }

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

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

  class PMTPulseRecoBase {

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

    virtual ~PMTPulseRecoBase() = default;

    const std::string& Name() 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;

    std::unique_ptr<pmtana::RiseTimeCalculatorBase> _risetime_calc_ptr = nullptr;

  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