RunningSumTPFinderPass4 Class Reference

#include <RunningSumTPFinderPass4.h>

Inheritance diagram for RunningSumTPFinderPass4:

Public Member Functions
	RunningSumTPFinderPass4 (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)
std::vector< short >	runSum (const std::vector< short > &orig)
void	hitFinding (const std::vector< short > &waveform, const std::vector< short > &iqr, std::vector< RunningSumTPFinderTool::Hit > &hits, int channel)

Protected Attributes
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
float	m_sigmaThreshold

Detailed Description

Definition at line 10 of file RunningSumTPFinderPass4.h.

Constructor & Destructor Documentation

RunningSumTPFinderPass4::RunningSumTPFinderPass4 ( fhicl::ParameterSet const &  p )

explicit

Definition at line 18 of file RunningSumTPFinderPass4_tool.cc.

  : m_useSignalKill      (p.get<bool>              ("UseSignalKill"        ,                        true)),
    m_signalKillLookahead(p.get<short>             ("SignalKillLookahead"  ,                           5)),
    m_signalKillThreshold(p.get<short>             ("SignalKillThreshold"  ,                          15)),
    m_signalKillNContig  (p.get<short>             ("SignalKillNContig"    ,                           1)),
    m_frugalNContig      (p.get<short>             ("FrugalPedestalNContig",                          10)),
    m_doFiltering        (p.get<bool>              ("DoFiltering"          ,                        true)),
    m_downsampleFactor   (p.get<unsigned int>      ("DownsampleFactor"     ,                           1)),
    m_filterTaps         (p.get<std::vector<short>>("FilterCoeffs"         , {2,  9, 23, 31, 23,  9,  2})),
    m_multiplier         (std::accumulate(m_filterTaps.begin(), m_filterTaps.end(), 0)),
    m_sigmaThreshold     (p.get<float>             ("SigmaThreshold"       ,                           5))
    // Default filter taps calculated by:
    // np.round(scipy.signal.firwin(7, 0.1)*100)
    // Initialize member data here.
{
  std::cout << "Sigma Threshold = " << m_sigmaThreshold << std::endl;
}

Member Function Documentation

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

protected

Definition at line 36 of file RunningSumTPFinderPass4_tool.cc.

                                                                                 {

  //---------------------------------------------
  // Do the downsampling
  //---------------------------------------------
  if (m_downsampleFactor==1) {
    return orig;
  }
  else {
    std::vector<short> waveform;
    for(size_t i=0; i<orig.size(); i+=m_downsampleFactor) {
      waveform.push_back(orig[i]);
    }
    return waveform;
  }
}

std::vector< short > RunningSumTPFinderPass4::filter ( const std::vector< short > &  orig )

protected

Definition at line 67 of file RunningSumTPFinderPass4_tool.cc.

                                                                               {

  //---------------------------------------------
  // Filtering
  //---------------------------------------------
  const size_t ntaps = m_filterTaps.size();
  const short*  taps = m_filterTaps.data();

  std::vector<short> filtered(m_doFiltering ? 
                              apply_fir_filter(pedsub, ntaps, taps) :
                              pedsub);
  if (!m_doFiltering) {
    std::transform(filtered.begin(), filtered.end(),
                   filtered.begin(), 
                   [=](short a) { return a*m_multiplier; });
  }
  return filtered;
}

std::vector< RunningSumTPFinderTool::Hit > RunningSumTPFinderPass4::findHits ( const std::vector< unsigned int > &  channel_numbers, const std::vector< std::vector< short >> &  collection_samples  )

virtual

Implements RunningSumTPFinderTool.

Definition at line 120 of file RunningSumTPFinderPass4_tool.cc.

                                                                                         {

  auto hits = std::vector<RunningSumTPFinderTool::Hit>();
  std::cout << "findHits called with "      << collection_samples.size()
            << " channels. First chan has " << collection_samples[0].size() << " samples" << std::endl;

  for(size_t ich=0; ich<collection_samples.size(); ++ich){
    const std::vector<short>& waveformOrig = collection_samples[ich];

    std::vector<short> waveform  = downSample  (waveformOrig);
    std::vector<short> pedestal  = findPedestal(waveform    );
    std::vector<short> pedsub(waveform.size(), 0);
    for(size_t i=0; i<pedsub.size(); ++i)
      pedsub[i]=waveform[i]-pedestal[i];
    std::vector<short> iqr=frugal_iqr(waveform, pedestal, m_frugalNContig);
    std::vector<short> filtered  = filter(pedsub);
    hitFinding(filtered, iqr, hits, channel_numbers[ich]);
  }
  std::cout << "Returning " << hits.size() << " hits" << std::endl;
  std::cout << "hits/channel=" << float(hits.size())/collection_samples   .size() << std::endl;
  std::cout << "hits/tick="    << float(hits.size())/collection_samples[0].size() << std::endl;
  return hits;
}

std::vector< short > RunningSumTPFinderPass4::findPedestal ( const std::vector< short > &  orig )

protected

Definition at line 53 of file RunningSumTPFinderPass4_tool.cc.

{
  //---------------------------------------------
  // Pedestal subtraction
  //---------------------------------------------
  const std::vector<short>& pedestal=m_useSignalKill ?
    frugal_pedestal_sigkill(waveform,
                            m_signalKillLookahead,
                            m_signalKillThreshold,
                            m_signalKillNContig) :
    frugal_pedestal(waveform, m_frugalNContig);
  return pedestal;
}

void RunningSumTPFinderPass4::hitFinding ( const std::vector< short > &  waveform, const std::vector< short > &  iqr, std::vector< RunningSumTPFinderTool::Hit > &  hits, int  channel  )

protected

Definition at line 87 of file RunningSumTPFinderPass4_tool.cc.

                                                 {

  //---------------------------------------------
  // Hit finding
  //---------------------------------------------
  bool is_hit  = false;
  bool was_hit = false;
  RunningSumTPFinderTool::Hit hit(channel, 0, 0, 0);
  for(size_t isample=0; isample<waveform.size()-1; ++isample){
    int   sample_time = isample * m_downsampleFactor;
    short adc         = waveform[isample];
    is_hit = (float)adc >  100*m_sigmaThreshold*iqr[isample];
    if(is_hit && !was_hit) {
      hit.startTime         = sample_time;
      hit.charge            = adc;
      hit.timeOverThreshold = m_downsampleFactor;
    }
    if(is_hit && was_hit) {
      hit.charge            += adc*m_downsampleFactor;
      hit.timeOverThreshold += m_downsampleFactor;
    }
    if(!is_hit && was_hit) {
      hit.charge /= m_multiplier;
      hits.push_back(hit);
    }
    was_hit = is_hit;
  }
}

std::vector<short> RunningSumTPFinderPass4::runSum ( const std::vector< short > &  orig )

protected

Member Data Documentation

bool RunningSumTPFinderPass4::m_doFiltering

protected

Definition at line 36 of file RunningSumTPFinderPass4.h.

unsigned int RunningSumTPFinderPass4::m_downsampleFactor

protected

Definition at line 37 of file RunningSumTPFinderPass4.h.

std::vector<short> RunningSumTPFinderPass4::m_filterTaps

protected

Definition at line 38 of file RunningSumTPFinderPass4.h.

short RunningSumTPFinderPass4::m_frugalNContig

protected

Definition at line 35 of file RunningSumTPFinderPass4.h.

int RunningSumTPFinderPass4::m_multiplier

protected

Definition at line 39 of file RunningSumTPFinderPass4.h.

float RunningSumTPFinderPass4::m_sigmaThreshold

protected

Definition at line 40 of file RunningSumTPFinderPass4.h.

short RunningSumTPFinderPass4::m_signalKillLookahead

protected

Definition at line 32 of file RunningSumTPFinderPass4.h.

short RunningSumTPFinderPass4::m_signalKillNContig

protected

Definition at line 34 of file RunningSumTPFinderPass4.h.

short RunningSumTPFinderPass4::m_signalKillThreshold

protected

Definition at line 33 of file RunningSumTPFinderPass4.h.

bool RunningSumTPFinderPass4::m_useSignalKill

protected

Definition at line 31 of file RunningSumTPFinderPass4.h.

---

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

• duneana/duneana/DAQSimAna/RunningSumHitFinder/RunningSumTPFinderPass4.h
• duneana/duneana/DAQSimAna/RunningSumHitFinder/RunningSumTPFinderPass4_tool.cc