AlgoFixedWindow.cxx

Go to the documentation of this file.

//
//  AlgoFixedWindow source
//

#include "fhiclcpp/ParameterSet.h"

#include "AlgoFixedWindow.h"

namespace pmtana {

  //*******************************************************************************
  AlgoFixedWindow::AlgoFixedWindow(const std::string name) : PMTPulseRecoBase(name)
  //*******************************************************************************
  {
    Reset();

    _index_start = 0;

    _index_end = 0;
  }

  //****************************************************************************************
  AlgoFixedWindow::AlgoFixedWindow(
    const fhicl::ParameterSet& pset,
    std::unique_ptr<pmtana::RiseTimeCalculatorBase> risetimecalculator,
    //AlgoFixedWindow::AlgoFixedWindow(const ::fcllite::PSet& pset,
    const std::string name)
    : PMTPulseRecoBase(name)
  //****************************************************************************************
  {
    _risetime_calc_ptr = std::move(risetimecalculator);

    Reset();

    _index_start = pset.get<size_t>("StartIndex");

    _index_end = pset.get<size_t>("EndIndex");
  }

  //***************************************************************
  void AlgoFixedWindow::Reset()
  //***************************************************************
  {
    if (!(_pulse_v.size())) _pulse_v.push_back(_pulse);

    _pulse_v[0].reset_param();
  }

  //***************************************************************
  bool AlgoFixedWindow::RecoPulse(const Waveform_t& wf,
                                  const PedestalMean_t& mean_v,
                                  const PedestalSigma_t& sigma_v)
  //***************************************************************
  {
    this->Reset();

    if (_index_start >= wf.size()) return true;

    _pulse_v[0].t_start = (double)(_index_start);

    _pulse_v[0].ped_mean = mean_v.front();

    _pulse_v[0].ped_sigma = sigma_v.front();

    if (!_index_end)

      _pulse_v[0].t_end = (double)(wf.size() - 1);

    else if (_index_end < wf.size())

      _pulse_v[0].t_end = (double)_index_end;

    else

      _pulse_v[0].t_end = wf.size() - 1;

    _pulse_v[0].t_max =
      PMTPulseRecoBase::Max(wf, _pulse_v[0].peak, _index_start, _pulse_v[0].t_end);

    _pulse_v[0].peak -= mean_v.front();

    PMTPulseRecoBase::Integral(wf, _pulse_v[0].area, _index_start, _pulse_v[0].t_end);

    _pulse_v[0].area =
      _pulse_v[0].area - (_pulse_v[0].t_end - _pulse_v[0].t_start + 1) * mean_v.front();

    if (_risetime_calc_ptr)
      _pulse_v[0].t_rise = _risetime_calc_ptr->RiseTime(
        {wf.begin() + _pulse.t_start, wf.begin() + _pulse.t_end},
        {mean_v.begin() + _pulse.t_start, mean_v.begin() + _pulse.t_end},
        true);

    return true;
  }

}