pmtana::AlgoSiPM Class Reference

#include <AlgoSiPM.h>

Inheritance diagram for pmtana::AlgoSiPM:

Public Member Functions
	AlgoSiPM (const fhicl::ParameterSet &pset, 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...

Detailed Description

Definition at line 22 of file AlgoSiPM.h.

Constructor & Destructor Documentation

pmtana::AlgoSiPM::AlgoSiPM	(	const fhicl::ParameterSet &	pset,
		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"       );

//    _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 36 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))
        {
          _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))
      {
        _pulse_v.push_back(_pulse);
        record_hit = false;
      }
      _pulse.reset_param();

    }

    return true;

  }

void pmtana::AlgoSiPM::Reset ( void  )

virtual

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

Reimplemented from pmtana::PMTPulseRecoBase.

Definition at line 28 of file AlgoSiPM.cxx.

  {

    PMTPulseRecoBase::Reset();

  }

Member Data Documentation

double pmtana::AlgoSiPM::_2nd_thres

protected

Definition at line 51 of file AlgoSiPM.h.

double pmtana::AlgoSiPM::_adc_thres

protected

Definition at line 45 of file AlgoSiPM.h.

int pmtana::AlgoSiPM::_min_width

protected

Definition at line 48 of file AlgoSiPM.h.

double pmtana::AlgoSiPM::_pedestal

protected

Definition at line 54 of file AlgoSiPM.h.

---

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

• larana/larana/OpticalDetector/OpHitFinder/AlgoSiPM.h
• larana/larana/OpticalDetector/OpHitFinder/AlgoSiPM.cxx