Inheritance diagram for reco_tool::CandHitStandard:

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

	~CandHitStandard ()

void	configure (const fhicl::ParameterSet &pset) override

void	findHitCandidates (const std::vector< float > &, size_t, size_t, size_t, HitCandidateVec &) const override

void	MergeHitCandidates (const std::vector< float > &, const HitCandidateVec &, MergeHitCandidateVec &) const override

Private Types
using	HitCandidate_t = struct HitCandidate{size_t startTick

using	HitCandidateVec = std::vector< HitCandidate_t >

using	MergeHitCandidateVec = std::vector< HitCandidateVec >

using	Waveform = std::vector< float >

Private Member Functions
void	findHitCandidates (std::vector< float >::const_iterator, std::vector< float >::const_iterator, size_t, size_t, HitCandidateVec &) const

Private Attributes
int	fPlane
	Plane we are meant to work with. More...

float	fRoiThreshold
	minimum maximum to minimum peak distance More...

const geo::GeometryCore *	fGeometry = lar::providerFrom<geo::Geometry>()

size_t	stopTick

size_t	maxTick

size_t	minTick

float	maxDerivative

float	minDerivative

float	hitCenter

float	hitSigma

float	hitHeight

Detailed Description

Definition at line 21 of file CandHitStandard_tool.cc.

Member Typedef Documentation

using reco_tool::ICandidateHitFinder::HitCandidate_t = struct HitCandidate { size_t startTick

inherited

Definition at line 30 of file ICandidateHitFinder.h.

using reco_tool::ICandidateHitFinder::HitCandidateVec = std::vector<HitCandidate_t>

inherited

Definition at line 41 of file ICandidateHitFinder.h.

using reco_tool::ICandidateHitFinder::MergeHitCandidateVec = std::vector<HitCandidateVec>

inherited

Definition at line 42 of file ICandidateHitFinder.h.

using reco_tool::ICandidateHitFinder::Waveform = std::vector<float>

inherited

Definition at line 44 of file ICandidateHitFinder.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 57 of file CandHitStandard_tool.cc.

References configure().

 {
     configure(pset);
 }

reco_tool::CandHitStandard::~CandHitStandard ( )

Definition at line 62 of file CandHitStandard_tool.cc.

63 {

64 }

Member Function Documentation

void reco_tool::CandHitStandard::configure ( const fhicl::ParameterSet & pset )

overridevirtual

Implements reco_tool::ICandidateHitFinder.

Definition at line 66 of file CandHitStandard_tool.cc.

References fPlane, fRoiThreshold, and fhicl::ParameterSet::get().

Referenced by CandHitStandard().

 {
     // Start by recovering the parameters
     fPlane        = pset.get< int   >("Plane",        0);
     fRoiThreshold = pset.get< float >("RoiThreshold", 5.);
 
     return;
 }

void reco_tool::CandHitStandard::findHitCandidates	(	const std::vector< float > &	waveform,
		size_t	roiStartTick,
		size_t	channel,
		size_t	eventCount,
		HitCandidateVec &	hitCandidateVec
	)		const

overridevirtual

Implements reco_tool::ICandidateHitFinder.

Definition at line 75 of file CandHitStandard_tool.cc.

References geo::GeometryCore::ChannelToWire(), and fGeometry.

Referenced by findHitCandidates().

 {
     // Recover the plane index for this method
     std::vector<geo::WireID> wids  = fGeometry->ChannelToWire(channel);
     size_t                   plane = wids[0].Plane;
     
     // Use the recursive version to find the candidate hits
     findHitCandidates(waveform.begin(),waveform.end(),roiStartTick,plane,hitCandidateVec);
     
     return;
 }

void reco_tool::CandHitStandard::findHitCandidates	(	std::vector< float >::const_iterator	startItr,
		std::vector< float >::const_iterator	stopItr,
		size_t	roiStartTick,
		size_t	planeIdx,
		HitCandidateVec &	hitCandidateVec
	)		const

private

Definition at line 91 of file CandHitStandard_tool.cc.

References findHitCandidates(), fRoiThreshold, and max.

 {
     // Need a minimum number of ticks to do any work here
     if (std::distance(startItr,stopItr) > 4)
     {
         // Find the highest peak in the range given
         auto maxItr = std::max_element(startItr, stopItr);
         
         float maxValue = *maxItr;
         int   maxTime  = std::distance(startItr,maxItr);
         
         if (maxValue > fRoiThreshold)
         {
             // backwards to find first bin for this candidate hit
             auto firstItr = std::distance(startItr,maxItr) > 2 ? maxItr - 1 : startItr;
             
             while(firstItr != startItr)
             {
                 // Check for pathology where waveform goes too negative
                 if (*firstItr < -fRoiThreshold) break;
                 
                 // Check both sides of firstItr and look for min/inflection point
                 if (*firstItr < *(firstItr+1) && *firstItr <= *(firstItr-1)) break;
                 
                 firstItr--;
             }
             
             int firstTime = std::distance(startItr,firstItr);
             
             // Recursive call to find all candidate hits earlier than this peak
             findHitCandidates(startItr, firstItr + 1, roiStartTick, planeIdx, hitCandidateVec);
             
             // forwards to find last bin for this candidate hit
             auto lastItr = std::distance(maxItr,stopItr) > 2 ? maxItr + 1 : stopItr - 1;
             
             while(lastItr != stopItr - 1)
             {
                 // Check for pathology where waveform goes too negative
                 if (*lastItr < -fRoiThreshold) break;
                 
                 // Check both sides of firstItr and look for min/inflection point
                 if (*lastItr <= *(lastItr+1) && *lastItr < *(lastItr-1)) break;
                 
                 lastItr++;
             }
             
             int lastTime = std::distance(startItr,lastItr);
             
             // Now save this candidate's start and max time info
             HitCandidate hitCandidate;
             hitCandidate.startTick     = roiStartTick + firstTime;
             hitCandidate.stopTick      = roiStartTick + lastTime;
             hitCandidate.maxTick       = roiStartTick + firstTime;
             hitCandidate.minTick       = roiStartTick + lastTime;
             hitCandidate.maxDerivative = *(startItr + firstTime);
             hitCandidate.minDerivative = *(startItr + lastTime);
             hitCandidate.hitCenter     = roiStartTick + maxTime;
             hitCandidate.hitSigma      = std::max(2.,float(lastTime - firstTime) / 6.);
             hitCandidate.hitHeight     = maxValue;
             
             hitCandidateVec.push_back(hitCandidate);
             
             // Recursive call to find all candidate hits later than this peak
             findHitCandidates(lastItr + 1, stopItr, roiStartTick + std::distance(startItr,lastItr + 1), planeIdx, hitCandidateVec);
         }
     }
     
     return;
 }

void reco_tool::CandHitStandard::MergeHitCandidates	(	const std::vector< float > &	signalVec,
		const HitCandidateVec &	hitCandidateVec,
		MergeHitCandidateVec &	mergedHitsVec
	)		const

overridevirtual

Implements reco_tool::ICandidateHitFinder.

Definition at line 165 of file CandHitStandard_tool.cc.

References DEFINE_ART_CLASS_TOOL.

 {
     // If no hits then nothing to do here
     if (hitCandidateVec.empty()) return;
     
     // The idea is to group hits that "touch" so they can be part of common fit, those that
     // don't "touch" are fit independently. So here we build the output vector to achieve that
     HitCandidateVec groupedHitVec;
     int             lastTick = hitCandidateVec.front().stopTick;
     
     // Step through the input hit candidates and group them by proximity
     for(const auto& hitCandidate : hitCandidateVec)
     {
         // Check condition that we have a new grouping
         if (int(hitCandidate.startTick) - lastTick > 1)
         {
             mergedHitsVec.emplace_back(groupedHitVec);
             
             groupedHitVec.clear();
         }
         
         // Add the current hit to the current group
         groupedHitVec.emplace_back(hitCandidate);
         
         lastTick = hitCandidate.stopTick;
     }
     
     // Check end condition
     if (!groupedHitVec.empty()) mergedHitsVec.emplace_back(groupedHitVec);
     
     return;
 }

Member Data Documentation

const geo::GeometryCore* reco_tool::CandHitStandard::fGeometry = lar::providerFrom<geo::Geometry>()

private

Definition at line 52 of file CandHitStandard_tool.cc.

Referenced by findHitCandidates().

int reco_tool::CandHitStandard::fPlane

private

Plane we are meant to work with.

Definition at line 49 of file CandHitStandard_tool.cc.

Referenced by configure().

float reco_tool::CandHitStandard::fRoiThreshold

private

minimum maximum to minimum peak distance

Definition at line 50 of file CandHitStandard_tool.cc.

Referenced by configure(), and findHitCandidates().

float reco_tool::ICandidateHitFinder::hitCenter

inherited

Definition at line 36 of file ICandidateHitFinder.h.

float reco_tool::ICandidateHitFinder::hitHeight

inherited

Definition at line 38 of file ICandidateHitFinder.h.

float reco_tool::ICandidateHitFinder::hitSigma

inherited

Definition at line 37 of file ICandidateHitFinder.h.

float reco_tool::ICandidateHitFinder::maxDerivative

inherited

Definition at line 34 of file ICandidateHitFinder.h.

size_t reco_tool::ICandidateHitFinder::maxTick

inherited

Definition at line 32 of file ICandidateHitFinder.h.

float reco_tool::ICandidateHitFinder::minDerivative

inherited

Definition at line 35 of file ICandidateHitFinder.h.

size_t reco_tool::ICandidateHitFinder::minTick

inherited

Definition at line 33 of file ICandidateHitFinder.h.

size_t reco_tool::ICandidateHitFinder::stopTick

inherited

Definition at line 31 of file ICandidateHitFinder.h.

Referenced by reco_tool::CandHitDerivative::findHitCandidates(), and reco_tool::CandHitMorphological::findHitCandidates().

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

larreco/v06_64_02/source/larreco/HitFinder/HitFinderTools/CandHitStandard_tool.cc

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation