TriggerPrimitiveFinderPass2 Class Reference

Inheritance diagram for TriggerPrimitiveFinderPass2:

Public Member Functions
	TriggerPrimitiveFinderPass2 (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 TriggerPrimitiveFinderPass1
	TriggerPrimitiveFinderPass1 (fhicl::ParameterSet const &p)
Public Member Functions inherited from TriggerPrimitiveFinderTool
virtual	~TriggerPrimitiveFinderTool ()=default

Protected Member Functions
void	hitFinding (const std::vector< short > &waveform, const std::vector< short > &iqr, std::vector< TriggerPrimitiveFinderTool::Hit > &hits, int channel)
Protected Member Functions inherited from TriggerPrimitiveFinderPass1
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)

Private Attributes
float	m_sigmaThreshold

Additional Inherited Members
Protected Attributes inherited from TriggerPrimitiveFinderPass1
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 15 of file TriggerPrimitiveFinderPass2_tool.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 37 of file TriggerPrimitiveFinderPass2_tool.cc.

    : TriggerPrimitiveFinderPass1(p),
      m_sigmaThreshold(p.get<float>("ThresholdInSigma", 500))
// Initialize member data here.
{
  std::cout << "Threshold in sigma is " << m_sigmaThreshold  << " (ignore the Threshold = 10 on previous line)\n";
}

Member Function Documentation

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

virtual

Reimplemented from TriggerPrimitiveFinderPass1.

Definition at line 83 of file TriggerPrimitiveFinderPass2_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> 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;
}

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

protected

Definition at line 46 of file TriggerPrimitiveFinderPass2_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=(float)adc>m_sigmaThreshold*iqr[isample];
        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;
    }
}

Member Data Documentation

float TriggerPrimitiveFinderPass2::m_sigmaThreshold

private

Definition at line 33 of file TriggerPrimitiveFinderPass2_tool.cc.

---

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

• duneana/duneana/DAQSimAna/TriggerPrimitiveFinderPass2_tool.cc