PedAlgoUB.cxx

Go to the documentation of this file.

//
//  PedAlgoUB source
//

#include "PedAlgoUB.h"

#include "fhiclcpp/ParameterSet.h"

namespace pmtana {

  //************************************************
  PedAlgoUB::PedAlgoUB(const std::string name) : PMTPedestalBase(name), _beamgatealgo(name)
  //************************************************
  {}

  //*************************************************************
  PedAlgoUB::PedAlgoUB(const fhicl::ParameterSet& pset,
                       //PedAlgoUB::PedAlgoUB(const ::fcllite::PSet &pset,
                       const std::string name)
    : PMTPedestalBase(name)
    //, _beamgatealgo(pset.get_pset("BeamGateAlgo"),"BeamGateAlgo")
    , _beamgatealgo(pset, "BeamGateAlgo")
  //*************************************************************
  {
    _beam_gate_samples = pset.get<unsigned int>("BeamGateSamples");
  }

  //*********************************************************************
  bool PedAlgoUB::ComputePedestal(const pmtana::Waveform_t& wf,
                                  pmtana::PedestalMean_t& mean_v,
                                  pmtana::PedestalSigma_t& sigma_v)
  //*********************************************************************
  {

    if (wf.size() < _beam_gate_samples) {

      double ped_mean = wf.front(); //first sample
      double ped_sigma = 0;

      for (auto& v : mean_v)
        v = ped_mean;
      for (auto& v : sigma_v)
        v = ped_sigma;

      return true;
    }

    else {

      _beamgatealgo.Evaluate(wf);
      mean_v = _beamgatealgo.Mean();
      sigma_v = _beamgatealgo.Sigma();

      return true;
    }
  }

}