HitFilterAlg.cxx

Go to the documentation of this file.

#include "HitFilterAlg.h"

#include "larcoreobj/SimpleTypesAndConstants/geo_types.h"
#include "lardataobj/RecoBase/Hit.h"

#include "fhiclcpp/ParameterSet.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

namespace hit{

  HitFilterAlg::HitFilterAlg(fhicl::ParameterSet const & p) :  
    fMinPulseHeight(p.get< std::vector<float> >("MinPulseHeight")),
    fMinPulseSigma(p.get< std::vector<float> >("MinPulseSigma"))
{
}

  bool HitFilterAlg::IsGoodHit(const recob::Hit& hit) const {

    const float              hitPH    = hit.PeakAmplitude();
    const float              hitSigma = hit.RMS();

    const geo::WireID& wireID   = hit.WireID();
    const size_t             view     = wireID.Plane;

    if (view >= fMinPulseSigma.size() || view >= fMinPulseHeight.size()) {
      mf::LogError("HitFilterAlg") << "Filtering settings not configured for all views! Will not filter hits in unconfigured views!";
      return true;
    }

    if ( hitPH    > fMinPulseHeight[view] &&
         hitSigma > fMinPulseSigma[view] ) {
      return true;
    }
    else return false;
  }
}//end namespace hit