PndEmcPSAParabolicBaseline.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

//-----------------------------------------------------------
// Description:
//      Pulseshape analysis for ADC waveforms
//      Amplitude and time are determined from the simple parabolic fit
//      by three points (maximum bin and its left and right neigbours)
//
//-----------------------------------------------------------

#ifndef EMCPSAPARABOLICBASELINE_HH
#define EMCPSAPARABOLICBASELINE_HH

// Base Class Headers ----------------
#include <vector>
#include "PndEmcAbsPSA.h"
#include "TObject.h"

class PndEmcWaveform;

class PndEmcPSAParabolicBaseline : public PndEmcAbsPSA {
 public:
  // Constructors/Destructors ---------
  PndEmcPSAParabolicBaseline(Int_t baselength = 0) : lBaselength(baselength) {}
  virtual ~PndEmcPSAParabolicBaseline() { ; }

  virtual void Process(const PndEmcWaveform *waveform, Double_t &amplitude, Double_t &time);

  virtual Int_t Process(const PndEmcWaveform *waveform);
  virtual void GetHit(Int_t i, Double_t &energy, Double_t &time);

  void FitPeak(const std::vector<Double_t> &signal, Double_t &ampl, Double_t &pos, Int_t peakBin) const;
  void FitPeak(const std::vector<Double_t> &signal, Double_t &ampl, Double_t &pos, Int_t start, Int_t end) const;
  void FitPeak(const std::vector<Double_t> &signal, Double_t &ampl, Double_t &pos) const;
  void GetBaseline(const std::vector<Double_t> &signal, Double_t &baseline) const;

 private:
  Int_t lBaselength;
  Double_t lAmplitude;
  Double_t lTime;

  ClassDef(PndEmcPSAParabolicBaseline, 1);
};

#endif