Inheritance diagram for reco_tool::PeakFitterGaussElimination:

Public Member Functions
	PeakFitterGaussElimination (const fhicl::ParameterSet &pset)

void	findPeakParameters (const std::vector< float > &, const ICandidateHitFinder::HitCandidateVec &, PeakParamsVec &, double &, int &) const override

Private Types
using	PeakParamsVec = std::vector< PeakFitParams_t >

Private Attributes
float	fStepSize
	Step size used by gaussian elim alg. More...

float	fMax
	Max. More...

std::unique_ptr< util::GaussianEliminationAlg >	fGEAlg

Detailed Description

Definition at line 15 of file PeakFitterGaussElimination_tool.cc.

Member Typedef Documentation

using reco_tool::IPeakFitter::PeakParamsVec = std::vector<PeakFitParams_t>

inherited

Definition at line 34 of file IPeakFitter.h.

Constructor & Destructor Documentation

reco_tool::PeakFitterGaussElimination::PeakFitterGaussElimination ( const fhicl::ParameterSet & pset )

explicit

Definition at line 35 of file PeakFitterGaussElimination_tool.cc.

References fMax, fStepSize, and fhicl::ParameterSet::get().

   {
     // Start by recovering the parameters
     fStepSize = pset.get<float>("StepSize", 0.1);
     fMax = pset.get<float>("Max", 0.5);
 
     //    fGEAlg    = std::make_unique<util::GaussianEliminationAlg>(fStepSize, fMax);
 
     return;
   }

Member Function Documentation

void reco_tool::PeakFitterGaussElimination::findPeakParameters	(	const std::vector< float > &	roiSignalVec,
		const ICandidateHitFinder::HitCandidateVec &	hitCandidateVec,
		PeakParamsVec &	,
		double &	,
		int &
	)		const

overridevirtual

Implements reco_tool::IPeakFitter.

Definition at line 47 of file PeakFitterGaussElimination_tool.cc.

References bin, and DEFINE_ART_CLASS_TOOL.

   {
     // This module tries to use the method for fitting hits found in the RRFHitFinder
     // from Wes Ketchum. It uses the gaussian elimation algorithm he set up.
     //
     // *** NOTE: this algorithm assumes the reference time for input hit candidates is to
     //           the first tick of the input waveform (ie 0)
     //
     if (hitCandidateVec.empty()) return;
 
     std::vector<float> meanVec;
     std::vector<float> sigmaVec;
     std::vector<float> heightVec;
 
     for (const auto& hitCandidate : hitCandidateVec) {
       float candMean = hitCandidate.hitCenter;
       float candSigma = hitCandidate.hitSigma;
       size_t bin = std::floor(candMean);
 
       bin = std::min(bin, roiSignalVec.size() - 1);
 
       float candHeight =
         roiSignalVec[bin] - (candMean - (float)bin) * (roiSignalVec[bin] - roiSignalVec[bin + 1]);
 
       meanVec.push_back(candMean);
       sigmaVec.push_back(candSigma);
       heightVec.push_back(candHeight);
     }
 
     return;
   }