RunningSumTPFinderPass3 Class Reference

#include <RunningSumTPFinderPass3.h>

Inheritance diagram for RunningSumTPFinderPass3:

Public Member Functions
	RunningSumTPFinderPass3 (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
short	m_runningSumAlpha
float	m_sigmaThreshold

Detailed Description

Definition at line 10 of file RunningSumTPFinderPass3.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 18 of file RunningSumTPFinderPass3_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_runningSumAlpha    (p.get<short>             ("RunningSumAlpha"      ,                          97)),
    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.
{

}

Member Function Documentation

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

protected

Definition at line 37 of file RunningSumTPFinderPass3_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 > RunningSumTPFinderPass3::filter ( const std::vector< short > &  orig )

protected

Definition at line 68 of file RunningSumTPFinderPass3_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 > RunningSumTPFinderPass3::findHits ( const std::vector< unsigned int > &  channel_numbers, const std::vector< std::vector< short >> &  collection_samples  )

virtual

Implements RunningSumTPFinderTool.

Definition at line 137 of file RunningSumTPFinderPass3_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);
    std::vector<short> runingSum = runSum(filtered);
    hitFinding(runingSum, 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 > RunningSumTPFinderPass3::findPedestal ( const std::vector< short > &  orig )

protected

Definition at line 54 of file RunningSumTPFinderPass3_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 RunningSumTPFinderPass3::hitFinding ( const std::vector< short > &  waveform, const std::vector< short > &  iqr, std::vector< RunningSumTPFinderTool::Hit > &  hits, int  channel  )

protected

Definition at line 104 of file RunningSumTPFinderPass3_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 >  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 > RunningSumTPFinderPass3::runSum ( const std::vector< short > &  orig )

protected

Definition at line 87 of file RunningSumTPFinderPass3_tool.cc.

                                                                                 {

  //---------------------------------------------
  // Running Sum
  //---------------------------------------------
  std::vector<short> runningSum(filtered.size(), 0); runningSum[0] = 0;
  for (size_t i=0; i<filtered.size(); ++i) {
    short sumValue = (filtered[i]/10) + ((runningSum[i-1]/10)*m_runningSumAlpha);
    sumValue /= 10;
    if (sumValue < 0) runningSum[i] = 0;
    else              runningSum[i] = sumValue;
  }
  return runningSum;
}

Member Data Documentation

bool RunningSumTPFinderPass3::m_doFiltering

protected

Definition at line 36 of file RunningSumTPFinderPass3.h.

unsigned int RunningSumTPFinderPass3::m_downsampleFactor

protected

Definition at line 37 of file RunningSumTPFinderPass3.h.

std::vector<short> RunningSumTPFinderPass3::m_filterTaps

protected

Definition at line 38 of file RunningSumTPFinderPass3.h.

short RunningSumTPFinderPass3::m_frugalNContig

protected

Definition at line 35 of file RunningSumTPFinderPass3.h.

int RunningSumTPFinderPass3::m_multiplier

protected

Definition at line 39 of file RunningSumTPFinderPass3.h.

short RunningSumTPFinderPass3::m_runningSumAlpha

protected

Definition at line 40 of file RunningSumTPFinderPass3.h.

float RunningSumTPFinderPass3::m_sigmaThreshold

protected

Definition at line 41 of file RunningSumTPFinderPass3.h.

short RunningSumTPFinderPass3::m_signalKillLookahead

protected

Definition at line 32 of file RunningSumTPFinderPass3.h.

short RunningSumTPFinderPass3::m_signalKillNContig

protected

Definition at line 34 of file RunningSumTPFinderPass3.h.

short RunningSumTPFinderPass3::m_signalKillThreshold

protected

Definition at line 33 of file RunningSumTPFinderPass3.h.

bool RunningSumTPFinderPass3::m_useSignalKill

protected

Definition at line 31 of file RunningSumTPFinderPass3.h.

---

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

• duneana/duneana/DAQSimAna/RunningSumHitFinder/RunningSumTPFinderPass3.h
• duneana/duneana/DAQSimAna/RunningSumHitFinder/RunningSumTPFinderPass3_tool.cc