67 <<
"IndWidth must be 1 at minimum. Resetting width to one time tick";
72 <<
"ColWidth must be 1 at minimum. Resetting width to one time tick";
96 std::vector<recob::Wire>
const& wireVec(*wireVecHandle);
104 float threshold_peak = 0;
105 float threshold_tail = -99;
109 for (
unsigned int wireIter = 0; wireIter < wireVec.size(); wireIter++) {
115 std::vector<float> signal(wire->
Signal());
134 float half_width = ((float)width - 1) / 2.;
139 for (timeIter = signal.begin(); timeIter < signal.end(); timeIter++) {
144 peak_val = *timeIter;
145 else if (width % 2 == 0)
146 peak_val = 0.5 * (*timeIter + *(timeIter + 1));
149 if (peak_val < threshold_peak)
continue;
152 if (time_bin - half_width < 0)
continue;
153 if (time_bin + half_width > signal.size())
continue;
156 int begin_tail_tick = std::floor(time_bin - half_width);
157 float totalCharge =
getTotalCharge(&signal[begin_tail_tick], width, threshold_tail);
158 if (totalCharge == -999) {
160 <<
"Rejecting would be hit at (plane,wire,time_bin,first_bin,last_bin)=(" 161 << wire_id.
Plane <<
"," << wire_id.
Wire <<
"," << time_bin <<
"," << begin_tail_tick
162 <<
"," << begin_tail_tick + width - 1 <<
"): " << signal.at(time_bin - 1) <<
" " 163 << signal.at(time_bin) <<
" " << signal.at(time_bin + 1);
169 float hit_time = time_bin;
170 if (width % 2 == 0) hit_time = time_bin + 0.5;
173 const raw::TDCtick_t start_tick = hit_time - width, end_tick = hit_time + width;
193 if (wire_id.
Plane == 0)
195 else if (wire_id.
Plane == 1)
197 else if (wire_id.
Plane == 2)
204 <<
"\tTotal hits (U,V,Y)= (" << hitCollection_U.
size() <<
"," 205 << hitCollection_V.
size() <<
"," << hitCollection_Y.
size() <<
")";
210 mf::LogInfo(
"TTHitFinder") <<
"Total TTHitFinder hits (U,V,Y)=(" << hitCollection_U.
size()
211 <<
"," << hitCollection_V.
size() <<
"," << hitCollection_Y.
size()
222 float threshold = -99)
225 float totalCharge = 0;
227 if (signal_vector[
tick] < threshold) {
231 totalCharge += signal_vector[
tick];
int fColWidth
Induction wire hit width (in time ticks)
float fMinSigTailInd
Collection wire signal height threshold at peak.
std::vector< WireID > ChannelToWire(raw::ChannelID_t const channel) const
Returns a list of wires connected to the specified TPC channel.
float fMinSigPeakCol
Induction wire signal height threshold at peak.
MaybeLogger_< ELseverityLevel::ELsev_info, false > LogInfo
EDProducer(fhicl::ParameterSet const &pset)
Definition of basic raw digits.
WireID_t Wire
Index of the wire within its plane.
static void declare_products(art::ProducesCollector &collector, std::string instance_name="", bool doWireAssns=true, bool doRawDigitAssns=true)
Declares the hit products we are going to fill.
int TDCtick_t
Type representing a TDC tick.
Class managing the creation of a new recob::Hit object.
Helper functions to create a hit.
A class handling a collection of hits and its associations.
#define DEFINE_ART_MODULE(klass)
Signal from induction planes.
raw::ChannelID_t Channel() const
Returns the ID of the channel (or InvalidChannelID)
enum geo::_plane_sigtype SigType_t
Enumerate the possible plane projections.
tick_as<> tick
Tick number, represented by std::ptrdiff_t.
void emplace_back(recob::Hit &&hit, art::Ptr< recob::Wire > const &wire=art::Ptr< recob::Wire >(), art::Ptr< raw::RawDigit > const &digits=art::Ptr< raw::RawDigit >())
Adds the specified hit to the data collection.
PlaneID_t Plane
Index of the plane within its TPC.
float fMinSigTailCol
Induction wire signal height threshold outside peak.
Definition of data types for geometry description.
void put_into(art::Event &)
Moves the data into an event.
Detector simulation of raw signals on wires.
std::vector< float > Signal() const
Return a zero-padded full length vector filled with RoI signal.
ProducesCollector & producesCollector() noexcept
std::string fCalDataModuleLabel
size_t size() const
Returns the number of hits currently in the collection.
SigType_t SignalType(PlaneID const &pid) const
Returns the type of signal on the channels of specified TPC plane.
bool getByLabel(std::string const &label, std::string const &instance, Handle< PROD > &result) const
float fMinSigPeakInd
Input caldata module name.
MaybeLogger_< ELseverityLevel::ELsev_warning, false > LogWarning
Declaration of basic channel signal object.
void produce(art::Event &evt) override
int fIndWidth
Collection wire signal height threshold outside peak.
TTHitFinder(fhicl::ParameterSet const &pset)
recob::Hit && move()
Prepares the constructed hit to be moved away.
art framework interface to geometry description
float getTotalCharge(const float *, int, float)
Collection wire hit width (in time ticks)
Signal from collection planes.