hit::RFFHitFinderAlg Class Reference

#include "RFFHitFinderAlg.h"

Public Member Functions
	RFFHitFinderAlg (fhicl::ParameterSet const &)
void	SetFitterParamsVectors (unsigned int num_planes)
void	Run (std::vector< recob::Wire > const &, std::vector< recob::Hit > &, geo::WireReadoutGeom const &)

Private Member Functions
void	SetFitterParams (unsigned int)
void	EmplaceHit (std::vector< recob::Hit > &, recob::Wire const &, float const &, raw::TDCtick_t const &, raw::TDCtick_t const &, geo::SigType_t const &, geo::WireID const &)

Private Attributes
const float	SQRT_TWO_PI = 2.506628
std::vector< float >	fMatchThresholdVec
std::vector< unsigned int >	fMergeMultiplicityVec
std::vector< float >	fAmpThresholdVec
RFFHitFitter	fFitter

Detailed Description

Definition at line 31 of file RFFHitFinderAlg.h.

Constructor & Destructor Documentation

hit::RFFHitFinderAlg::RFFHitFinderAlg

(

fhicl::ParameterSet const &

p

)

Title: RFFHitFinderAlg Class Author: Wes Ketchum (wketc.nosp@m.hum@.nosp@m.lanl..nosp@m.gov)

Description: Class that runs the RFF HitFinder. Implements an RFFHitFitter, and takes the result and stores it in recob::Hit objects.

Input: recob::Wire Output: recob::Hit

Definition at line 21 of file RFFHitFinderAlg.cxx.

References fAmpThresholdVec, fMatchThresholdVec, fMergeMultiplicityVec, and fhicl::ParameterSet::get().

{
  fMatchThresholdVec = p.get<std::vector<float>>("MeanMatchThreshold");
  fMergeMultiplicityVec = p.get<std::vector<unsigned int>>("MinMergeMultiplicity");
  fAmpThresholdVec = p.get<std::vector<float>>("AmplitudeThreshold", std::vector<float>(1, 0.0));
}

Member Function Documentation

void hit::RFFHitFinderAlg::EmplaceHit ( std::vector< recob::Hit > &  hitVector, recob::Wire const &  wire, float const &  summedADCTotal, raw::TDCtick_t const &  startTick, raw::TDCtick_t const &  endTick, geo::SigType_t const &  sigtype, geo::WireID const &  wireID  )

private

Definition at line 79 of file RFFHitFinderAlg.cxx.

References hit::RFFHitFitter::AmplitudeErrorVector(), hit::RFFHitFitter::AmplitudeVector(), recob::Wire::Channel(), fFitter, hit::RFFHitFitter::MeanErrorVector(), hit::RFFHitFitter::MeanVector(), hit::RFFHitFitter::NHits(), hit::RFFHitFitter::SigmaErrorVector(), hit::RFFHitFitter::SigmaVector(), SQRT_TWO_PI, and recob::Wire::View().

Referenced by Run().

{

  float totalArea = 0.0;
  std::vector<float> areaVector(fFitter.NHits());
  std::vector<float> areaErrorVector(fFitter.NHits());
  std::vector<float> areaFracVector(fFitter.NHits());

  for (size_t ihit = 0; ihit < fFitter.NHits(); ihit++) {
    areaVector[ihit] = fFitter.AmplitudeVector()[ihit] * fFitter.SigmaVector()[ihit] * SQRT_TWO_PI;
    areaErrorVector[ihit] =
      SQRT_TWO_PI * std::sqrt(fFitter.AmplitudeVector()[ihit] * fFitter.SigmaErrorVector()[ihit] *
                                fFitter.AmplitudeVector()[ihit] * fFitter.SigmaErrorVector()[ihit] +
                              fFitter.AmplitudeErrorVector()[ihit] * fFitter.SigmaVector()[ihit] *
                                fFitter.AmplitudeErrorVector()[ihit] * fFitter.SigmaVector()[ihit]);
    totalArea += areaVector[ihit];
  }

  for (size_t ihit = 0; ihit < fFitter.NHits(); ihit++) {
    areaFracVector[ihit] = areaVector[ihit] / totalArea;

    hitVector.emplace_back(wire.Channel(),
                           startTick,
                           endTick,
                           fFitter.MeanVector()[ihit] + (float)startTick,
                           fFitter.MeanErrorVector()[ihit],
                           fFitter.SigmaVector()[ihit],
                           fFitter.AmplitudeVector()[ihit],
                           fFitter.AmplitudeErrorVector()[ihit],
                           summedADCTotal * areaFracVector[ihit],
                           summedADCTotal * areaFracVector[ihit], //NOT CORRECTLY FILLED
                           areaVector[ihit],
                           areaErrorVector[ihit],
                           fFitter.NHits(),
                           ihit,
                           -999.,
                           -999,
                           wire.View(),
                           sigtype,
                           wireID);
  }
}

void hit::RFFHitFinderAlg::Run	(	std::vector< recob::Wire > const &	wireVector,
		std::vector< recob::Hit > &	hitVector,
		geo::WireReadoutGeom const &	wireReadoutGeom
	)

Definition at line 55 of file RFFHitFinderAlg.cxx.

References geo::WireReadoutGeom::ChannelToWire(), EmplaceHit(), fFitter, hit::RFFHitFitter::RunFitter(), SetFitterParams(), and geo::WireReadoutGeom::SignalType().

Referenced by hit::RFFHitFinder::produce().

{
  hitVector.reserve(wireVector.size());
  for (auto const& wire : wireVector) {
    geo::SigType_t const& sigtype = wireReadoutGeom.SignalType(wire.Channel());
    geo::WireID const wireID = wireReadoutGeom.ChannelToWire(wire.Channel()).at(0);

    SetFitterParams(wire.View());

    for (auto const& roi : wire.SignalROI().get_ranges()) {
      fFitter.RunFitter(roi.data());

      const float summedADCTotal = std::accumulate(roi.data().begin(), roi.data().end(), 0.0);
      const raw::TDCtick_t startTick = roi.begin_index();
      const raw::TDCtick_t endTick = roi.begin_index() + roi.size();

      EmplaceHit(hitVector, wire, summedADCTotal, startTick, endTick, sigtype, wireID);
    } //end loop over ROIs on wire

  } //end loop over wires
}

void hit::RFFHitFinderAlg::SetFitterParams ( unsigned int  p )

private

Definition at line 50 of file RFFHitFinderAlg.cxx.

References fAmpThresholdVec, fFitter, fMatchThresholdVec, fMergeMultiplicityVec, and hit::RFFHitFitter::SetFitterParams().

Referenced by Run().

{
  fFitter.SetFitterParams(fMatchThresholdVec[p], fMergeMultiplicityVec[p], fAmpThresholdVec[p]);
}

void hit::RFFHitFinderAlg::SetFitterParamsVectors

(

unsigned int

num_planes

)

Definition at line 28 of file RFFHitFinderAlg.cxx.

References fAmpThresholdVec, fMatchThresholdVec, and fMergeMultiplicityVec.

Referenced by hit::RFFHitFinder::beginJob().

{
  //If size zero, throw. If size one, assume same for all planes.
  //If size > 1 but < n_planes, throw. If size = n_plane, good.

  if (fMatchThresholdVec.size() == 0 || fMergeMultiplicityVec.size() == 0 ||
      fAmpThresholdVec.size() == 0)
    throw std::runtime_error("Error in RFFHitFinderAlg: Configured with zero planes.");

  if ((fMatchThresholdVec.size() > 1 && fMatchThresholdVec.size() < n_planes) ||
      (fMergeMultiplicityVec.size() > 1 && fMergeMultiplicityVec.size() < n_planes) ||
      (fAmpThresholdVec.size() > 1 && fAmpThresholdVec.size() < n_planes))
    throw std::runtime_error("Error in RFFHitFinderAlg: Configured with incorrect n_planes.");

  if (fMatchThresholdVec.size() == 1) fMatchThresholdVec.resize(n_planes, fMatchThresholdVec[0]);

  if (fMergeMultiplicityVec.size() == 1)
    fMergeMultiplicityVec.resize(n_planes, fMergeMultiplicityVec[0]);

  if (fAmpThresholdVec.size() == 1) fAmpThresholdVec.resize(n_planes, fAmpThresholdVec[0]);
}

Member Data Documentation

std::vector<float> hit::RFFHitFinderAlg::fAmpThresholdVec

private

Definition at line 46 of file RFFHitFinderAlg.h.

Referenced by RFFHitFinderAlg(), SetFitterParams(), and SetFitterParamsVectors().

RFFHitFitter hit::RFFHitFinderAlg::fFitter

private

Definition at line 58 of file RFFHitFinderAlg.h.

Referenced by EmplaceHit(), Run(), and SetFitterParams().

std::vector<float> hit::RFFHitFinderAlg::fMatchThresholdVec

private

Definition at line 44 of file RFFHitFinderAlg.h.

Referenced by RFFHitFinderAlg(), SetFitterParams(), and SetFitterParamsVectors().

std::vector<unsigned int> hit::RFFHitFinderAlg::fMergeMultiplicityVec

private

Definition at line 45 of file RFFHitFinderAlg.h.

Referenced by RFFHitFinderAlg(), SetFitterParams(), and SetFitterParamsVectors().

const float hit::RFFHitFinderAlg::SQRT_TWO_PI = 2.506628

private

Definition at line 33 of file RFFHitFinderAlg.h.

Referenced by EmplaceHit().

---

The documentation for this class was generated from the following files:

• larreco/v10_01_07/source/larreco/HitFinder/RFFHitFinderAlg.h
• larreco/v10_01_07/source/larreco/HitFinder/RFFHitFinderAlg.cxx