Public Member Functions | Private Member Functions | Private Attributes | List of all members
reco_tool::CandHitStandard Class Reference
Inheritance diagram for reco_tool::CandHitStandard:
reco_tool::ICandidateHitFinder

Public Member Functions

 CandHitStandard (const fhicl::ParameterSet &pset)
 
void findHitCandidates (const recob::Wire::RegionsOfInterest_t::datarange_t &, const size_t, const size_t, const size_t, HitCandidateVec &) const override
 
void MergeHitCandidates (const recob::Wire::RegionsOfInterest_t::datarange_t &, const HitCandidateVec &, MergeHitCandidateVec &) const override
 

Private Member Functions

void findHitCandidates (std::vector< float >::const_iterator, std::vector< float >::const_iterator, const size_t, const size_t, HitCandidateVec &) const
 
- Private Member Functions inherited from reco_tool::ICandidateHitFinder
virtual ~ICandidateHitFinder () noexcept=default
 

Private Attributes

const float fRoiThreshold
 minimum maximum to minimum peak distance More...
 
const geo::GeometryCore * fGeometry = lar::providerFrom<geo::Geometry>()
 

Additional Inherited Members

- Private Types inherited from reco_tool::ICandidateHitFinder
using HitCandidateVec = std::vector< HitCandidate >
 
using MergeHitCandidateVec = std::vector< HitCandidateVec >
 
using Waveform = std::vector< float >
 

Detailed Description

Definition at line 15 of file CandHitStandard_tool.cc.

Constructor & Destructor Documentation

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

Definition at line 44 of file CandHitStandard_tool.cc.

45  : fRoiThreshold(pset.get<float>("RoiThreshold", 5.))
46  {}
reco_tool::CandHitStandard::fRoiThreshold
const float fRoiThreshold
minimum maximum to minimum peak distance
Definition: CandHitStandard_tool.cc:37
fhicl::ParameterSet::get
T get(std::string const &key) const
Definition: ParameterSet.h:271

Member Function Documentation

void reco_tool::CandHitStandard::findHitCandidates ( const recob::Wire::RegionsOfInterest_t::datarange_t &  dataRange,
const size_t  roiStartTick,
const size_t  channel,
const size_t  eventCount,
HitCandidateVec &  hitCandidateVec 
) const
overridevirtual

Implements reco_tool::ICandidateHitFinder.

Definition at line 49 of file CandHitStandard_tool.cc.

54  {
55  // Recover the actual waveform
56  const Waveform& waveform = dataRange.data();
57 
58  // Recover the plane index for this method
59  std::vector<geo::WireID> wids = fGeometry->ChannelToWire(channel);
60  const size_t plane = wids[0].Plane;
61 
62  // Use the recursive version to find the candidate hits
63  findHitCandidates(waveform.begin(), waveform.end(), roiStartTick, plane, hitCandidateVec);
64 
65  return;
66  }
reco_tool::CandHitStandard::findHitCandidates
void findHitCandidates(const recob::Wire::RegionsOfInterest_t::datarange_t &, const size_t, const size_t, const size_t, HitCandidateVec &) const override
Definition: CandHitStandard_tool.cc:49
reco_tool::ICandidateHitFinder::Waveform
std::vector< float > Waveform
Definition: ICandidateHitFinder.h:39
geo::GeometryCore::ChannelToWire
std::vector< geo::WireID > ChannelToWire(raw::ChannelID_t const channel) const
Returns a list of wires connected to the specified TPC channel.
Definition: GeometryCore.cxx:1107
channel
uint8_t channel
Definition: CRTFragment.hh:201
reco_tool::CandHitStandard::fGeometry
const geo::GeometryCore * fGeometry
Definition: CandHitStandard_tool.cc:39
void reco_tool::CandHitStandard::findHitCandidates ( std::vector< float >::const_iterator  startItr,
std::vector< float >::const_iterator  stopItr,
const size_t  roiStartTick,
const size_t  planeIdx,
HitCandidateVec &  hitCandidateVec 
) const
private

Definition at line 69 of file CandHitStandard_tool.cc.

74  {
75  // Need a minimum number of ticks to do any work here
76  if (std::distance(startItr, stopItr) > 4) {
77  // Find the highest peak in the range given
78  auto maxItr = std::max_element(startItr, stopItr);
79 
80  float maxValue = *maxItr;
81  int maxTime = std::distance(startItr, maxItr);
82 
83  if (maxValue > fRoiThreshold) {
84  // backwards to find first bin for this candidate hit
85  auto firstItr = std::distance(startItr, maxItr) > 2 ? maxItr - 1 : startItr;
86 
87  while (firstItr != startItr) {
88  // Check for pathology where waveform goes too negative
89  if (*firstItr < -fRoiThreshold) break;
90 
91  // Check both sides of firstItr and look for min/inflection point
92  if (*firstItr < *(firstItr + 1) && *firstItr <= *(firstItr - 1)) break;
93 
94  firstItr--;
95  }
96 
97  int firstTime = std::distance(startItr, firstItr);
98 
99  // Recursive call to find all candidate hits earlier than this peak
100  findHitCandidates(startItr, firstItr + 1, roiStartTick, planeIdx, hitCandidateVec);
101 
102  // forwards to find last bin for this candidate hit
103  auto lastItr = std::distance(maxItr, stopItr) > 2 ? maxItr + 1 : stopItr - 1;
104 
105  while (lastItr != stopItr - 1) {
106  // Check for pathology where waveform goes too negative
107  if (*lastItr < -fRoiThreshold) break;
108 
109  // Check both sides of firstItr and look for min/inflection point
110  if (*lastItr <= *(lastItr + 1) && *lastItr < *(lastItr - 1)) break;
111 
112  lastItr++;
113  }
114 
115  int lastTime = std::distance(startItr, lastItr);
116 
117  // Now save this candidate's start and max time info
118  HitCandidate hitCandidate;
119  hitCandidate.startTick = roiStartTick + firstTime;
120  hitCandidate.stopTick = roiStartTick + lastTime;
121  hitCandidate.maxTick = roiStartTick + firstTime;
122  hitCandidate.minTick = roiStartTick + lastTime;
123  hitCandidate.maxDerivative = *(startItr + firstTime);
124  hitCandidate.minDerivative = *(startItr + lastTime);
125  hitCandidate.hitCenter = roiStartTick + maxTime;
126  hitCandidate.hitSigma = std::max(2., float(lastTime - firstTime) / 6.);
127  hitCandidate.hitHeight = maxValue;
128 
129  hitCandidateVec.push_back(hitCandidate);
130 
131  // Recursive call to find all candidate hits later than this peak
132  findHitCandidates(lastItr + 1,
133  stopItr,
134  roiStartTick + std::distance(startItr, lastItr + 1),
135  planeIdx,
136  hitCandidateVec);
137  }
138  }
139 
140  return;
141  }
reco_tool::CandHitStandard::findHitCandidates
void findHitCandidates(const recob::Wire::RegionsOfInterest_t::datarange_t &, const size_t, const size_t, const size_t, HitCandidateVec &) const override
Definition: CandHitStandard_tool.cc:49
reco_tool::CandHitStandard::fRoiThreshold
const float fRoiThreshold
minimum maximum to minimum peak distance
Definition: CandHitStandard_tool.cc:37
protoana::distance
double distance(double x1, double y1, double z1, double x2, double y2, double z2)
Definition: truepionXsection_module.cc:217
max
static int max(int a, int b)
Definition: fortranscanner.cpp:60366
void reco_tool::CandHitStandard::MergeHitCandidates ( const recob::Wire::RegionsOfInterest_t::datarange_t &  rangeData,
const HitCandidateVec &  hitCandidateVec,
MergeHitCandidateVec &  mergedHitsVec 
) const
overridevirtual

Implements reco_tool::ICandidateHitFinder.

Definition at line 144 of file CandHitStandard_tool.cc.

148  {
149  // If no hits then nothing to do here
150  if (hitCandidateVec.empty()) return;
151 
152  // The idea is to group hits that "touch" so they can be part of common fit, those that
153  // don't "touch" are fit independently. So here we build the output vector to achieve that
154  HitCandidateVec groupedHitVec;
155  int lastTick = hitCandidateVec.front().stopTick;
156 
157  // Step through the input hit candidates and group them by proximity
158  for (const auto& hitCandidate : hitCandidateVec) {
159  // Check condition that we have a new grouping
160  if (int(hitCandidate.startTick) - lastTick > 1) {
161  mergedHitsVec.emplace_back(groupedHitVec);
162 
163  groupedHitVec.clear();
164  }
165 
166  // Add the current hit to the current group
167  groupedHitVec.emplace_back(hitCandidate);
168 
169  lastTick = hitCandidate.stopTick;
170  }
171 
172  // Check end condition
173  if (!groupedHitVec.empty()) mergedHitsVec.emplace_back(groupedHitVec);
174 
175  return;
176  }
reco_tool::ICandidateHitFinder::HitCandidateVec
std::vector< HitCandidate > HitCandidateVec
Definition: ICandidateHitFinder.h:36

Member Data Documentation

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

Definition at line 39 of file CandHitStandard_tool.cc.

const float reco_tool::CandHitStandard::fRoiThreshold
private

minimum maximum to minimum peak distance

Definition at line 37 of file CandHitStandard_tool.cc.

The documentation for this class was generated from the following file:
Generated by   doxygen 1.8.11