#include <AlgoSiPM.h>

Inheritance diagram for pmtana::AlgoSiPM:

Public Member Functions
	AlgoSiPM (const fhicl::ParameterSet &pset, std::unique_ptr< pmtana::RiseTimeCalculatorBase > risetimecalculator=nullptr, const std::string name="AlgoSiPM")

void	Reset ()
	A method to be called event-wise to reset parameters. 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 &, const pmtana::PedestalMean_t &, const pmtana::PedestalSigma_t &)

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	_adc_thres

int	_min_width

double	_2nd_thres

double	_pedestal

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

Definition at line 22 of file AlgoSiPM.h.

Constructor & Destructor Documentation

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

Definition at line 12 of file AlgoSiPM.cxx.

     : PMTPulseRecoBase(name)
   {
 
     _adc_thres = pset.get< float >("ADCThreshold"   );
     _min_width = pset.get< float >("MinWidth"       );
     _2nd_thres = pset.get< float >("SecondThreshold");
     _pedestal  = pset.get< float >("Pedestal"       );
 
     _risetime_calc_ptr = std::move(risetimecalculator);
 
 //    _nsigma = 5;
 
     Reset();
 
   }

Member Function Documentation

bool pmtana::AlgoSiPM::RecoPulse	(	const pmtana::Waveform_t &	wf,
		const pmtana::PedestalMean_t &	ped_mean,
		const pmtana::PedestalSigma_t &	ped_rms
	)

protectedvirtual

Implements pmtana::PMTPulseRecoBase.

Definition at line 40 of file AlgoSiPM.cxx.

   {
 
     bool   fire          = false;
     bool   first_found   = false;
     bool   record_hit    = false;
     int    counter       = 0;
     //    double threshold   = (_2nd_thres > (_nsigma*_ped_rms) ? _2nd_thres
     //                                                    : (_nsigma*_ped_rms));
     //double pedestal      = _pedestal;
     double pedestal = ped_mean.front();  //Switch pedestal definition to incoroprate pedestal finder - K.S. 04/18/2019
 
     double threshold     = _adc_thres;
     threshold           += pedestal;
     double pre_threshold = _2nd_thres;
     pre_threshold       += pedestal;
 
     Reset();
 
     for (short const &value : wf) {
 
       if (!fire && (double(value) >= pre_threshold)) {
 
         // Found a new pulse
         fire           = true;
         first_found    = false;
         record_hit     = false;
         _pulse.t_start = counter;
 
       }
 
       if (fire && (double(value) < pre_threshold)) {
 
         // Found the end of a pulse
         fire = false;
         _pulse.t_end = counter - 1;
         if (record_hit && ((_pulse.t_end - _pulse.t_start) >= _min_width))
         {
           if(_risetime_calc_ptr)
             _pulse.t_rise = _risetime_calc_ptr->RiseTime(
                               {wf.begin()+_pulse.t_start, wf.begin()+_pulse.t_end},
                               {ped_mean.begin()+_pulse.t_start, ped_mean.begin()+_pulse.t_end},
                               true);
 
           _pulse_v.push_back(_pulse);
           record_hit = false;
         }
         _pulse.reset_param();
 
       }
 
       if (fire) {
 
         // We want to record the hit only if _adc_thres is reached
         if (!record_hit && (double(value) >= threshold)) record_hit = true;
 
         // Add this ADC count to the integral
         _pulse.area += (double(value) - double(pedestal));
 
         if (!first_found &&
             (_pulse.peak < (double(value) - double(pedestal)))) {
 
           // Found a new maximum
           _pulse.peak  = (double(value) - double(pedestal));
           _pulse.t_max = counter;
 
         }
         else if (!first_found)
           // Found the first peak
           first_found = true;
 
       }
 
       counter++;
 
     }
 
     if (fire) {
 
       // Take care of a pulse that did not finish within the readout window
       fire = false;
       _pulse.t_end = counter - 1;
       if (record_hit && ((_pulse.t_end - _pulse.t_start) >= _min_width))
       {
         if(_risetime_calc_ptr)
           _pulse.t_rise = _risetime_calc_ptr->RiseTime(
                             {wf.begin()+_pulse.t_start, wf.begin()+_pulse.t_end},
                             {ped_mean.begin()+_pulse.t_start, ped_mean.begin()+_pulse.t_end},
                             true);
 
         _pulse_v.push_back(_pulse);
         record_hit = false;
       }
       _pulse.reset_param();
 
     }
 
     return true;
 
   }

void pmtana::AlgoSiPM::Reset ( )

virtual

A method to be called event-wise to reset parameters.

Reimplemented from pmtana::PMTPulseRecoBase.

Definition at line 32 of file AlgoSiPM.cxx.

   {
 
     PMTPulseRecoBase::Reset();
 
   }

Member Data Documentation

double pmtana::AlgoSiPM::_2nd_thres

protected

Definition at line 53 of file AlgoSiPM.h.

double pmtana::AlgoSiPM::_adc_thres

protected

Definition at line 47 of file AlgoSiPM.h.

int pmtana::AlgoSiPM::_min_width

protected

Definition at line 50 of file AlgoSiPM.h.

double pmtana::AlgoSiPM::_pedestal

protected

Definition at line 56 of file AlgoSiPM.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