HitFilterAlg.cxx

Go to the documentation of this file.

#include "HitFilterAlg.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

namespace hit{

  HitFilterAlg::HitFilterAlg(fhicl::ParameterSet const & p) {   
    this->reconfigure(p);
  }
  
  void HitFilterAlg::reconfigure(fhicl::ParameterSet const & p) {
    // Implementation of optional member function here.
    fMinPulseHeight   = p.get< std::vector<float> >("MinPulseHeight");
    fMinPulseSigma    = p.get< std::vector<float> >("MinPulseSigma");
  }

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

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

    const geo::WireID& wireID   = hit.WireID();
    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