#include <TriggerPrimitiveFinderPass1.h>

Inheritance diagram for TriggerPrimitiveFinderPass1:

Public Member Functions
	TriggerPrimitiveFinderPass1 (fhicl::ParameterSet const &p)

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

Public Member Functions inherited from TriggerPrimitiveFinderTool
virtual	~TriggerPrimitiveFinderTool ()=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< TriggerPrimitiveFinderTool::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 TriggerPrimitiveFinderPass1.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 18 of file TriggerPrimitiveFinderPass1_tool.cc.

     : m_threshold(p.get<unsigned int>("Threshold", 10)),
       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)),
       // Default filter taps calculated by:
       //  np.round(scipy.signal.firwin(7, 0.1)*100)
       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))
 
 // Initialize member data here.
 {
     std::cout << "Threshold is " << m_threshold << std::endl;
 }

Member Function Documentation

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

protected

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

protected

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

virtual

Implements TriggerPrimitiveFinderTool.

Reimplemented in TriggerPrimitiveFinderPass2.

Definition at line 125 of file TriggerPrimitiveFinderPass1_tool.cc.

 {
     auto hits=std::vector<TriggerPrimitiveFinderTool::Hit>();
     std::cout << "findHits called with " << collection_samples.size()
               << " channels. First chan has " << collection_samples[0].size() << " samples" << std::endl;
     // std::cout << "First few samples: ";
     // for(int i=0; i<10; ++i) std::cout << collection_samples[0][i] << " ";
     // std::cout << 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> filtered=filter(pedsub);
         hitFinding(filtered, 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 > TriggerPrimitiveFinderPass1::findPedestal ( const std::vector< short > & orig )

protected

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

protected

Definition at line 89 of file TriggerPrimitiveFinderPass1_tool.cc.

 {
     //---------------------------------------------
     // Hit finding
     //---------------------------------------------
     bool is_hit=false;
     bool was_hit=false;
     TriggerPrimitiveFinderTool::Hit hit(channel, 0, 0, 0);
     for(size_t isample=0; isample<waveform.size()-1; ++isample){
         // if(ich>11510) std::cout << isample << " " << std::flush;
         int sample_time=isample*m_downsampleFactor;
         short adc=waveform[isample];
         is_hit=adc>(short)m_threshold;
         if(is_hit && !was_hit){
             // We just started a hit. Set the start time
             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){
             // The hit is over. Push it to the output vector
             hit.charge/=m_multiplier;
             hits.push_back(hit);
         }
         was_hit=is_hit;
     }
 }