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 Attributes
float	fStepSize
	Step size used by gaussian elim alg. More...

float	fMax
	Max. More...

std::unique_ptr< util::GaussianEliminationAlg >	fGEAlg

Additional Inherited Members
Private Types inherited from reco_tool::IPeakFitter
using	PeakParamsVec = std::vector< PeakFitParams_t >

Private Member Functions inherited from reco_tool::IPeakFitter
virtual	~IPeakFitter ()=default

Detailed Description

Definition at line 16 of file PeakFitterGaussElimination_tool.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 37 of file PeakFitterGaussElimination_tool.cc.

 {
     // 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 &	peakParamsVec,
		double &	chi2PerNDF,
		int &	NDF
	)		const

overridevirtual

Implements reco_tool::IPeakFitter.

Definition at line 49 of file PeakFitterGaussElimination_tool.cc.

 {
     // 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;
 }