Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
RunningSumTPFinderPass2 Class Reference

#include <RunningSumTPFinderPass2.h>

Inheritance diagram for RunningSumTPFinderPass2:
RunningSumTPFinderTool

Public Member Functions

 RunningSumTPFinderPass2 (fhicl::ParameterSet const &p)
 
virtual std::vector< RunningSumTPFinderTool::Hit > findHits (const std::vector< unsigned int > &channel_numbers, const std::vector< std::vector< short >> &collection_samples)
 
- Public Member Functions inherited from RunningSumTPFinderTool
virtual ~RunningSumTPFinderTool ()=default
 

Protected Member Functions

std::vector< short > downSample (const std::vector< short > &orig)
 
std::vector< short > findPedestal (const std::vector< short > &orig)
 
std::vector< short > filter (const std::vector< short > &orig)
 
void hitFinding (const std::vector< short > &waveform, std::vector< RunningSumTPFinderTool::Hit > &hits, int channel)
 

Protected Attributes

unsigned int m_threshold
 
bool m_useSignalKill
 
short m_signalKillLookahead
 
short m_signalKillThreshold
 
short m_signalKillNContig
 
short m_frugalNContig
 
bool m_doFiltering
 
unsigned int m_downsampleFactor
 
std::vector< short > m_filterTaps
 
int m_multiplier
 

Detailed Description

Definition at line 10 of file RunningSumTPFinderPass2.h.

Constructor & Destructor Documentation

RunningSumTPFinderPass2::RunningSumTPFinderPass2 ( fhicl::ParameterSet const &  p)
explicit

Definition at line 18 of file RunningSumTPFinderPass2_tool.cc.

19  : m_threshold (p.get<unsigned int> ("Threshold" , 10)),
20  m_useSignalKill (p.get<bool> ("UseSignalKill" , true)),
21  m_signalKillLookahead(p.get<short> ("SignalKillLookahead" , 5)),
22  m_signalKillThreshold(p.get<short> ("SignalKillThreshold" , 15)),
23  m_signalKillNContig (p.get<short> ("SignalKillNContig" , 1)),
24  m_frugalNContig (p.get<short> ("FrugalPedestalNContig", 10)),
25  m_doFiltering (p.get<bool> ("DoFiltering" , true)),
26  m_downsampleFactor (p.get<unsigned int> ("DownsampleFactor" , 1)),
27  m_filterTaps (p.get<std::vector<short>>("FilterCoeffs" , {2, 9, 23, 31, 23, 9, 2})),
28  m_multiplier (std::accumulate(m_filterTaps.begin(), m_filterTaps.end(), 0))
29  // Default filter taps calculated by:
30  // np.round(scipy.signal.firwin(7, 0.1)*100)
31  // Initialize member data here.
32 {
33  std::cout << "Threshold is " << m_threshold << std::endl;
34 }
test.p
p
Definition: test.py:223
RunningSumTPFinderPass2::m_filterTaps
std::vector< short > m_filterTaps
Definition: RunningSumTPFinderPass2.h:35
endl
QTextStream & endl(QTextStream &s)
Definition: qtextstream.cpp:2030

Member Function Documentation

std::vector< short > RunningSumTPFinderPass2::downSample ( const std::vector< short > &  orig)
protected

Definition at line 36 of file RunningSumTPFinderPass2_tool.cc.

36  {
37 
38  //---------------------------------------------
39  // Do the downsampling
40  //---------------------------------------------
41  if (m_downsampleFactor==1) {
42  return orig;
43  }
44  else {
45  std::vector<short> waveform;
46  for(size_t i=0; i<orig.size(); i+=m_downsampleFactor) {
47  waveform.push_back(orig[i]);
48  }
49  return waveform;
50  }
51 }
std::vector< short > RunningSumTPFinderPass2::filter ( const std::vector< short > &  orig)
protected

Definition at line 67 of file RunningSumTPFinderPass2_tool.cc.

67  {
68 
69  //---------------------------------------------
70  // Filtering
71  //---------------------------------------------
72  const size_t ntaps = m_filterTaps.size();
73  const short* taps = m_filterTaps.data();
74 
75  std::vector<short> filtered(m_doFiltering ?
76  apply_fir_filter(pedsub, ntaps, taps) :
77  pedsub);
78  if (!m_doFiltering) {
79  std::transform(filtered.begin(), filtered.end(),
80  filtered.begin(),
81  [=](short a) { return a*m_multiplier; });
82  }
83  return filtered;
84 }
a
const double a
Definition: gUpMuFluxGen.cxx:164
apply_fir_filter
std::vector< short > apply_fir_filter(const std::vector< short > &input, const size_t ntaps, const short *taps)
Definition: AlgParts.h:104
RunningSumTPFinderPass2::m_filterTaps
std::vector< short > m_filterTaps
Definition: RunningSumTPFinderPass2.h:35
std::vector< RunningSumTPFinderTool::Hit > RunningSumTPFinderPass2::findHits ( const std::vector< unsigned int > &  channel_numbers,
const std::vector< std::vector< short >> &  collection_samples 
)
virtual

Implements RunningSumTPFinderTool.

Definition at line 119 of file RunningSumTPFinderPass2_tool.cc.

120  {
121 
122  auto hits = std::vector<RunningSumTPFinderTool::Hit>();
123  std::cout << "findHits called with " << collection_samples.size()
124  << " channels. First chan has " << collection_samples[0].size() << " samples" << std::endl;
125 
126  for(size_t ich=0; ich<collection_samples.size(); ++ich){
127  const std::vector<short>& waveformOrig = collection_samples[ich];
128 
129  std::vector<short> waveform = downSample (waveformOrig);
130  std::vector<short> pedestal = findPedestal(waveform );
131  std::vector<short> pedsub(waveform.size(), 0);
132  for(size_t i=0; i<pedsub.size(); ++i)
133  pedsub[i]=waveform[i]-pedestal[i];
134  std::vector<short> filtered = filter(pedsub);
135  hitFinding(filtered, hits, channel_numbers[ich]);
136  }
137  std::cout << "Returning " << hits.size() << " hits" << std::endl;
138  std::cout << "hits/channel=" << float(hits.size())/collection_samples .size() << std::endl;
139  std::cout << "hits/tick=" << float(hits.size())/collection_samples[0].size() << std::endl;
140  return hits;
141 }
reco_momentum_tuples.float
float
Definition: reco_momentum_tuples.py:12
RunningSumTPFinderPass2::filter
std::vector< short > filter(const std::vector< short > &orig)
Definition: RunningSumTPFinderPass2_tool.cc:67
RunningSumTPFinderPass2::hitFinding
void hitFinding(const std::vector< short > &waveform, std::vector< RunningSumTPFinderTool::Hit > &hits, int channel)
Definition: RunningSumTPFinderPass2_tool.cc:87
util::size
decltype(auto) constexpr size(T &&obj)
ADL-aware version of std::size.
Definition: StdUtils.h:92
RunningSumTPFinderPass2::downSample
std::vector< short > downSample(const std::vector< short > &orig)
Definition: RunningSumTPFinderPass2_tool.cc:36
RunningSumTPFinderPass2::findPedestal
std::vector< short > findPedestal(const std::vector< short > &orig)
Definition: RunningSumTPFinderPass2_tool.cc:53
endl
QTextStream & endl(QTextStream &s)
Definition: qtextstream.cpp:2030
std::vector< short > RunningSumTPFinderPass2::findPedestal ( const std::vector< short > &  orig)
protected

Definition at line 53 of file RunningSumTPFinderPass2_tool.cc.

54 {
55  //---------------------------------------------
56  // Pedestal subtraction
57  //---------------------------------------------
58  const std::vector<short>& pedestal=m_useSignalKill ?
59  frugal_pedestal_sigkill(waveform,
60  m_signalKillLookahead,
61  m_signalKillThreshold,
62  m_signalKillNContig) :
63  frugal_pedestal(waveform, m_frugalNContig);
64  return pedestal;
65 }
frugal_pedestal_sigkill
std::vector< short > frugal_pedestal_sigkill(const std::vector< short > &raw_in, const int lookahead, const int threshold, const int ncontig)
Definition: AlgParts.h:21
frugal_pedestal
std::vector< short > frugal_pedestal(const std::vector< short > &raw_in, const int ncontig)
Definition: AlgParts.h:86
void RunningSumTPFinderPass2::hitFinding ( const std::vector< short > &  waveform,
std::vector< RunningSumTPFinderTool::Hit > &  hits,
int  channel 
)
protected

Definition at line 87 of file RunningSumTPFinderPass2_tool.cc.

89  {
90 
91  //---------------------------------------------
92  // Hit finding
93  //---------------------------------------------
94  bool is_hit = false;
95  bool was_hit = false;
96  RunningSumTPFinderTool::Hit hit(channel, 0, 0, 0);
97  for(size_t isample=0; isample<waveform.size()-1; ++isample){
98  int sample_time = isample * m_downsampleFactor;
99  short adc = waveform[isample];
100  is_hit = adc > (short)m_threshold*10;
101  if(is_hit && !was_hit) {
102  hit.startTime = sample_time;
103  hit.charge = adc;
104  hit.timeOverThreshold = m_downsampleFactor;
105  }
106  if(is_hit && was_hit) {
107  hit.charge += adc*m_downsampleFactor;
108  hit.timeOverThreshold += m_downsampleFactor;
109  }
110  if(!is_hit && was_hit) {
111  hit.charge /= m_multiplier;
112  hits.push_back(hit);
113  }
114  was_hit = is_hit;
115  }
116 }
adc
int16_t adc
Definition: CRTFragment.hh:202
channel
uint8_t channel
Definition: CRTFragment.hh:201
hit
Detector simulation of raw signals on wires.
Definition: NumberOfHitsFilter_module.cc:25

Member Data Documentation

bool RunningSumTPFinderPass2::m_doFiltering
protected

Definition at line 33 of file RunningSumTPFinderPass2.h.

unsigned int RunningSumTPFinderPass2::m_downsampleFactor
protected

Definition at line 34 of file RunningSumTPFinderPass2.h.

std::vector<short> RunningSumTPFinderPass2::m_filterTaps
protected

Definition at line 35 of file RunningSumTPFinderPass2.h.

short RunningSumTPFinderPass2::m_frugalNContig
protected

Definition at line 32 of file RunningSumTPFinderPass2.h.

int RunningSumTPFinderPass2::m_multiplier
protected

Definition at line 36 of file RunningSumTPFinderPass2.h.

short RunningSumTPFinderPass2::m_signalKillLookahead
protected

Definition at line 29 of file RunningSumTPFinderPass2.h.

short RunningSumTPFinderPass2::m_signalKillNContig
protected

Definition at line 31 of file RunningSumTPFinderPass2.h.

short RunningSumTPFinderPass2::m_signalKillThreshold
protected

Definition at line 30 of file RunningSumTPFinderPass2.h.

unsigned int RunningSumTPFinderPass2::m_threshold
protected

Definition at line 27 of file RunningSumTPFinderPass2.h.

bool RunningSumTPFinderPass2::m_useSignalKill
protected

Definition at line 28 of file RunningSumTPFinderPass2.h.

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