//-----------------------------------------------------------
//
// Description:
//      Pulseshape analysis for ADC waveforms
//      (see header file for more details)
//
//      Dima Melnychuk
//
//-----------------------------------------------------------


// This Class' Header ------------------
#include "PndEmcPSAParabolic.h"
#include "PndEmcWaveform.h"

// Class Member definitions -----------

void PndEmcPSAParabolic::Process(const PndEmcWaveform *waveform, Double_t &amplitude, Double_t &time)
{
	const std::vector<Double_t> signal=waveform->GetSignal();
	FitPeak(signal, amplitude, time);
}

void PndEmcPSAParabolic::FitPeak(const std::vector<Double_t> & signal, Double_t& ampl, Double_t& pos, Int_t peakBin) const
{
	Int_t waveformLength=signal.size();
	
	ampl=pos=-1.;
	if (peakBin>0 && peakBin<waveformLength-1)
	{
		long theBin(peakBin);
	
		Double_t pValue = signal[peakBin];
		Double_t pPosition = Double_t(peakBin);
	
		Double_t leftValue = signal[theBin-1];
		Double_t rightValue = signal[theBin+1];
		if (leftValue<pValue && rightValue<pValue) {
			Double_t d = 0.25*(rightValue-leftValue);
			Double_t b = pValue-0.5*(leftValue+rightValue);
			pValue += d*d/b;
			pPosition += d/b;
		}
		ampl=pValue;
		pos=pPosition;
  }
}

void PndEmcPSAParabolic::FitPeak(const std::vector<Double_t> & signal, Double_t& ampl, Double_t& pos, Int_t start, Int_t end) const
{
	std::vector<Double_t>::const_iterator p;
	p=max_element(signal.begin()+start,signal.begin()+end);
	Int_t pPosition = distance(signal.begin(),p);
	FitPeak(signal,ampl,pos,pPosition);
}
  
void PndEmcPSAParabolic::FitPeak(const std::vector<Double_t> & signal, Double_t& ampl, Double_t& pos) const
{
	std::vector<Double_t>::const_iterator p;
	p=max_element(signal.begin(),signal.end());
	Int_t pPosition = distance(signal.begin(),p);
	FitPeak(signal,ampl,pos,pPosition);
}