50 const recob::Wire::RegionsOfInterest_t::datarange_t& dataRange,
51 const size_t roiStartTick,
57 const Waveform& waveform = dataRange.data();
61 const size_t plane = wids[0].Plane;
64 findHitCandidates(waveform.begin(), waveform.end(), roiStartTick, plane, hitCandidateVec);
71 const size_t roiStartTick,
72 const size_t planeIdx,
76 if (std::distance(startItr, stopItr) > 4) {
78 auto maxItr = std::max_element(startItr, stopItr);
80 float maxValue = *maxItr;
81 int maxTime = std::distance(startItr, maxItr);
85 auto firstItr = std::distance(startItr, maxItr) > 2 ? maxItr - 1 : startItr;
87 while (firstItr != startItr) {
92 if (*firstItr < *(firstItr + 1) && *firstItr <= *(firstItr - 1))
break;
97 int firstTime = std::distance(startItr, firstItr);
100 findHitCandidates(startItr, firstItr + 1, roiStartTick, planeIdx, hitCandidateVec);
103 auto lastItr = std::distance(maxItr, stopItr) > 2 ? maxItr + 1 : stopItr - 1;
105 while (lastItr != stopItr - 1) {
110 if (*lastItr <= *(lastItr + 1) && *lastItr < *(lastItr - 1))
break;
115 int lastTime = std::distance(startItr, lastItr);
119 hitCandidate.
startTick = roiStartTick + firstTime;
120 hitCandidate.
stopTick = roiStartTick + lastTime;
121 hitCandidate.
maxTick = roiStartTick + firstTime;
122 hitCandidate.
minTick = roiStartTick + lastTime;
125 hitCandidate.
hitCenter = roiStartTick + maxTime;
126 hitCandidate.
hitSigma = std::max(2.,
float(lastTime - firstTime) / 6.);
129 hitCandidateVec.push_back(hitCandidate);
134 roiStartTick + std::distance(startItr, lastItr + 1),
148 if (hitCandidateVec.empty())
return;
153 int lastTick = hitCandidateVec.front().stopTick;
156 for (
const auto& hitCandidate : hitCandidateVec) {
158 if (
int(hitCandidate.startTick) - lastTick > 1) {
159 mergedHitsVec.emplace_back(groupedHitVec);
161 groupedHitVec.clear();
165 groupedHitVec.emplace_back(hitCandidate);
167 lastTick = hitCandidate.stopTick;
171 if (!groupedHitVec.empty()) mergedHitsVec.emplace_back(groupedHitVec);
Utilities related to art service access.
std::vector< WireID > ChannelToWire(raw::ChannelID_t const channel) const
Returns a list of wires connected to the specified TPC channel.
Description of geometry of one entire detector.
decltype(auto) get(T &&obj)
ADL-aware version of std::to_string.
This provides an interface for tools which are tasked with finding candidate hits on input waveforms...
art framework interface to geometry description