RegionAboveThresholdFinder.cxx

Go to the documentation of this file.

#include "RegionAboveThresholdFinder.h"
#include <stdexcept>

void hit::RegionAboveThresholdFinder::FillStartAndEndTicks(const std::vector<float>& signal,
                                                           std::vector<unsigned int>& start_ticks,
                                                           std::vector<unsigned int>& end_ticks)
{

  start_ticks.clear();
  end_ticks.clear();

  bool in_RAT = false;
  for (unsigned int i_tick = 0; i_tick < signal.size(); i_tick++) {

    if (!in_RAT && signal[i_tick] >= fThreshold) {
      start_ticks.push_back(i_tick);
      in_RAT = true;
    }
    else if (in_RAT && signal[i_tick] < fThreshold) {
      end_ticks.push_back(i_tick);
      in_RAT = false;
    }
  }

  if (in_RAT) end_ticks.push_back(signal.size());

  if (end_ticks.size() != start_ticks.size())
    throw std::runtime_error(
      "ERROR in RegionAboveThresholdFinder: start and end tick vectors not equal.");
}