#include <AlgoThreshold.h>

Inheritance diagram for pmtana::AlgoThreshold:

Public Member Functions
	AlgoThreshold (const std::string name="AlgoThreshold")
	Default constructor. More...

	AlgoThreshold (const fhicl::ParameterSet &pset, std::unique_ptr< pmtana::RiseTimeCalculatorBase > risetimecalculator=nullptr, const std::string name="AlgoThreshold")
	Alternative constructor. More...

void	Reset ()
	Implementation of AlgoThreshold::reset() method. More...

Public Member Functions inherited from pmtana::PMTPulseRecoBase
	PMTPulseRecoBase (const std::string name="noname")
	Default constructor with fhicl parameters. More...

virtual	~PMTPulseRecoBase ()=default
	Default destructor. More...

const std::string &	Name () const
	Name getter. More...

bool	Status () const
	Status getter. More...

bool	Reconstruct (const pmtana::Waveform_t &, const pmtana::PedestalMean_t &, const pmtana::PedestalSigma_t &)

const pulse_param &	GetPulse (size_t index=0) const

const pulse_param_array &	GetPulses () const
	A getter for the whole array of pulse_param struct object. More...

size_t	GetNPulse () const
	A getter for the number of reconstructed pulses from the input waveform. More...

Protected Member Functions
bool	RecoPulse (const pmtana::Waveform_t &wf, const pmtana::PedestalMean_t &mean_v, const pmtana::PedestalSigma_t &sigma_v)
	Implementation of AlgoThreshold::reco() method. More...

Protected Member Functions inherited from pmtana::PMTPulseRecoBase
bool	Integral (const std::vector< short > &wf, double &result, size_t begin=0, size_t end=0) const

bool	Derivative (const std::vector< short > &wf, std::vector< int32_t > &diff, size_t begin=0, size_t end=0) const

size_t	Max (const std::vector< short > &wf, double &result, size_t begin=0, size_t end=0) const

size_t	Min (const std::vector< short > &wf, double &result, size_t begin=0, size_t end=0) const

Protected Attributes
double	_start_adc_thres
	A variable holder for a user-defined absolute ADC threshold value. More...

double	_end_adc_thres

double	_nsigma_start
	A variable holder for a multiplicative factor for the pedestal standard deviation to define the threshold. More...

double	_nsigma_end

Protected Attributes inherited from pmtana::PMTPulseRecoBase
pulse_param_array	_pulse_v
	A container array of pulse_param struct objects to store (possibly multiple) reconstructed pulse(s). More...

pulse_param	_pulse
	A subject pulse_param object to be filled with the last reconstructed pulse parameters. More...

std::unique_ptr < pmtana::RiseTimeCalculatorBase >	_risetime_calc_ptr = nullptr
	Tool for rise time calculation. More...

Detailed Description

This class implements threshold algorithm to AlgoThreshold class. The algorithm defines a pulse in user-specified time window. A typical usage is to set the beginning of the window to be 0 (= start of the waveform) and integrate over the time of interest. By default, the ending is set to index=0, in which case it uses the ending index of the input waveform (i.e. full integration).

Definition at line 35 of file AlgoThreshold.h.

Constructor & Destructor Documentation

pmtana::AlgoThreshold::AlgoThreshold ( const std::string name = "AlgoThreshold" )

Default constructor.

Definition at line 14 of file AlgoThreshold.cxx.

                                                    : PMTPulseRecoBase(name)
   //***************************************************************************
   {
     //_adc_thres = 3;
     //_nsigma = 5;
     Reset();
   }

pmtana::AlgoThreshold::AlgoThreshold	(	const fhicl::ParameterSet &	pset,
		std::unique_ptr< pmtana::RiseTimeCalculatorBase >	risetimecalculator = `nullptr`,
		const std::string	name = `"AlgoThreshold"`
	)

Alternative constructor.

Definition at line 23 of file AlgoThreshold.cxx.

     : PMTPulseRecoBase(name)
   //*******************************************************
   {
 
     _start_adc_thres = pset.get<double>("StartADCThreshold");
     _end_adc_thres = pset.get<double>("EndADCThreshold");
 
     //_nsigma = pset.get<double>("NSigmaThreshold");
 
     _nsigma_start = pset.get<double>("NSigmaThresholdStart");
     _nsigma_end   = pset.get<double>("NSigmaThresholdEnd");
 
     _risetime_calc_ptr = std::move(risetimecalculator);
 
     Reset();
 
   }

Member Function Documentation

bool pmtana::AlgoThreshold::RecoPulse	(	const pmtana::Waveform_t &	wf,
		const pmtana::PedestalMean_t &	mean_v,
		const pmtana::PedestalSigma_t &	sigma_v
	)

protectedvirtual

Implementation of AlgoThreshold::reco() method.

Implements pmtana::PMTPulseRecoBase.

Definition at line 53 of file AlgoThreshold.cxx.

   {
     bool fire = false;
 
     double counter=0;
 
     double ped_mean = mean_v.front();
     double ped_rms  = sigma_v.front();
 
     //double threshold = ( _adc_thres > (_nsigma * ped_rms) ? _adc_thres : (_nsigma * ped_rms) );
     auto start_threshold = ( _start_adc_thres > (_nsigma_start * ped_rms) ? _start_adc_thres : (_nsigma_start * ped_rms) );
     auto end_threshold   = ( _end_adc_thres   > (_nsigma_end   * ped_rms) ? _end_adc_thres   : (_nsigma_end * ped_rms) );
 
     //    threshold += ped_mean
 
     start_threshold += ped_mean;
     end_threshold   += ped_mean;
 
     Reset();
 
     for(auto const &value : wf){
 
       if( !fire && ((double)value) >= start_threshold ){
 
         // Found a new pulse
 
         fire = true;
 
         _pulse.ped_mean  = ped_mean;
         _pulse.ped_sigma = ped_rms;
 
         //vic: i move t_start back one, this helps with porch
 
         _pulse.t_start = counter - 1 > 0 ? counter - 1 : counter;
         //std::cout << "counter: " << counter << " tstart : " << _pulse.t_start << "\n";
 
       }
 
       if( fire && ((double)value) < end_threshold ){
 
         // Found the end of a pulse
 
         fire = false;
 
         //vic: i move t_start forward one, this helps with tail
         _pulse.t_end = counter < wf.size()  ? counter : counter - 1;
 
   if(_risetime_calc_ptr)
     _pulse.t_rise = _risetime_calc_ptr->RiseTime(
                       {wf.begin()+_pulse.t_start, wf.begin()+_pulse.t_end},
                       {mean_v.begin()+_pulse.t_start, mean_v.begin()+_pulse.t_end},
                       true);
 
   _pulse_v.push_back(_pulse);
 
         _pulse.reset_param();
 
       }
 
       //std::cout << "\tFire=" << fire << std::endl;
 
       if(fire){
 
         // Add this adc count to the integral
 
         _pulse.area += ((double)value - (double)ped_mean);
 
         if(_pulse.peak < ((double)value - (double)ped_mean)) {
 
           // Found a new maximum
 
           _pulse.peak = ((double)value - (double)ped_mean);
 
           _pulse.t_max = counter;
 
         }
 
       }
 
       counter++;
     }
 
     if(fire){
 
       // Take care of a pulse that did not finish within the readout window.
 
       fire = false;
 
       _pulse.t_end = counter - 1;
 
       if(_risetime_calc_ptr)
         _pulse.t_rise = _risetime_calc_ptr->RiseTime(
                           {wf.begin()+_pulse.t_start, wf.begin()+_pulse.t_end},
                           {mean_v.begin()+_pulse.t_start, mean_v.begin()+_pulse.t_end},
                           true);
 
       _pulse_v.push_back(_pulse);
 
       _pulse.reset_param();
 
     }
 
     return true;
 
   }

void pmtana::AlgoThreshold::Reset ( )

virtual

Implementation of AlgoThreshold::reset() method.

Reimplemented from pmtana::PMTPulseRecoBase.

Definition at line 46 of file AlgoThreshold.cxx.

   {
     PMTPulseRecoBase::Reset();
   }

Member Data Documentation

double pmtana::AlgoThreshold::_end_adc_thres

protected

Definition at line 61 of file AlgoThreshold.h.

double pmtana::AlgoThreshold::_nsigma_end

protected

Definition at line 66 of file AlgoThreshold.h.

double pmtana::AlgoThreshold::_nsigma_start

protected

A variable holder for a multiplicative factor for the pedestal standard deviation to define the threshold.

Definition at line 65 of file AlgoThreshold.h.

double pmtana::AlgoThreshold::_start_adc_thres

protected

A variable holder for a user-defined absolute ADC threshold value.

Definition at line 60 of file AlgoThreshold.h.

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation