#include <AlgoSSPLeadingEdge.h>

Inheritance diagram for pmtana::AlgoSSPLeadingEdge:

Public Member Functions
	AlgoSSPLeadingEdge (const fhicl::ParameterSet &pset, const std::string name="AlgoSSPLeadingEdge")

	~AlgoSSPLeadingEdge ()

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...

Public Attributes
std::vector< int >	trg_wvf

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

double	_pedestal

int	_dwindow

size_t	_rdwindow

size_t	_pretrg

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 21 of file AlgoSSPLeadingEdge.h.

Constructor & Destructor Documentation

pmtana::AlgoSSPLeadingEdge::AlgoSSPLeadingEdge	(	const fhicl::ParameterSet &	pset,
		const std::string	name = `"AlgoSSPLeadingEdge"`
	)

Definition at line 17 of file AlgoSSPLeadingEdge.cc.

     : 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"       );
     _dwindow   = pset.get< int    >("DWindow"        );
     _rdwindow  = pset.get< size_t >("ReadoutWd"      );
     _pretrg    = pset.get< size_t >("PreTrg"         );
 
  //  std::cout << "VITOR DEBUG" << std::endl;
  //  std::cout << "ADCThreshold = " << _adc_thres << std::endl;
  //  std::cout << "DWindow = " << _dwindow << std::endl;
  //  std::cout << "Readout Window = " << _rdwindow << std::endl;
  //  std::cout << "Pre Trigger = " << _pretrg << std::endl;
  //  std::cout << "VITOR DEBUG" << std::endl;
 
 
 
 
 //    _nsigma = 5;
 
     Reset();
 
   }

pmtana::AlgoSSPLeadingEdge::~AlgoSSPLeadingEdge ( )

Definition at line 46 of file AlgoSSPLeadingEdge.cc.

47 {}

Member Function Documentation

bool pmtana::AlgoSSPLeadingEdge::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 58 of file AlgoSSPLeadingEdge.cc.

   {
 
     bool   fire          = false;
     int    counter       = 0;
     double pedestal      = _pedestal;
     double threshold     = _adc_thres;
     double threshold2    = _adc_thres;
     threshold           += pedestal;
  //   double pre_threshold = _2nd_thres;
  //   pre_threshold       += pedestal;
 
     int    d_window       = _dwindow; 
     int    timer          = 0;  //time to delay the wf and compare amplitudes
     int    readout_window = _rdwindow; //from 0 to 2046 in SSP
     int    readout_pretrigger = _pretrg; //(from 0 to 2047 in SSP)
     double threshold_cmp;
     int    tamanho, cnt2;
 
     Reset();
     tamanho = wf.size();
 
 
     for (short const &value : wf) {
 
          cnt2 = value;
          if(cnt2 != value){std::cout << "teste VITOR" << std::endl; }
 
          if (counter < (tamanho - d_window )) {
 
                 threshold_cmp = wf[counter+d_window]-wf[counter];       
                         
 
                 if (threshold_cmp >= threshold2){
                         
                         if(timer==0){trg_wvf.emplace_back(counter+d_window);}
 
                         if(!fire && timer==0){
                                 fire = true;
                                 if(readout_pretrigger < (counter + d_window)){
                                         _pulse.t_start = counter + d_window - readout_pretrigger;
                                 }else{
                                         _pulse.t_start = 1;
                                 }
                                 timer++;
                         }else if(timer>0 && timer < d_window){
                                 timer++;
                         }else if(timer == d_window){
                                 timer=0;
                         }
 
                 }else {
                         if(timer>0 && timer < d_window){timer++;}
                         else if(timer==d_window || timer==0 ){timer=0;}
 
                 }
                 if (fire==true && counter == (_pulse.t_start + readout_window)){ //we need to configure the _min_widht to be 700 bins
                          fire = false;
                          _pulse.t_end = counter - 1;
                          _pulse_v.push_back(_pulse);
                          _pulse.reset_param();
 
                 }
         }
 
         counter++;
 
     }
 
     if (fire) { //It's not mine, but is a good strategy to handle the pulses at the end of waveform
 
       // Take care of a pulse that did not finish within the readout window
    
       fire = false;
       _pulse.t_end = tamanho;
       _pulse_v.push_back(_pulse);
       _pulse.reset_param();
 
     }
 
     return true;
 
   }

void pmtana::AlgoSSPLeadingEdge::Reset ( )

virtual

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

Reimplemented from pmtana::PMTPulseRecoBase.

Definition at line 50 of file AlgoSSPLeadingEdge.cc.

   {
 
     PMTPulseRecoBase::Reset();
 
   }

Member Data Documentation

double pmtana::AlgoSSPLeadingEdge::_adc_thres

protected

Definition at line 42 of file AlgoSSPLeadingEdge.h.

int pmtana::AlgoSSPLeadingEdge::_dwindow

protected

Definition at line 48 of file AlgoSSPLeadingEdge.h.

double pmtana::AlgoSSPLeadingEdge::_pedestal

protected

Definition at line 45 of file AlgoSSPLeadingEdge.h.

size_t pmtana::AlgoSSPLeadingEdge::_pretrg

protected

Definition at line 52 of file AlgoSSPLeadingEdge.h.

size_t pmtana::AlgoSSPLeadingEdge::_rdwindow

protected

Definition at line 51 of file AlgoSSPLeadingEdge.h.

std::vector<int> pmtana::AlgoSSPLeadingEdge::trg_wvf

Definition at line 25 of file AlgoSSPLeadingEdge.h.

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

duneopdet/duneopdet/OpticalDetector/AlgoSSPLeadingEdge.h
duneopdet/duneopdet/OpticalDetector/AlgoSSPLeadingEdge.cc

Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation