OptimumFilter.hh

Go to the documentation of this file.

#ifndef __QOPTIMUMFILTER_HH__
#define __QOPTIMUMFILTER_HH__
#include <vector>
//#include "Sampler.hh"

#include "QVector.hh"
#include "QVectorC.hh"
#include "QError.hh"
#include "cmath"
using namespace std; 
using namespace Cuore; 


struct PulseInfo
{
int max_index;
float amplitude,frac_delay,delay,second_amp,second_delay;
};



class OptimumFilter {
  
public:

  OptimumFilter();


  virtual ~OptimumFilter();
  
  QError Compute();

  QError Compute_d();
  QError SetFilters(const QVector &avg, const QVector &noise);
  QError SetFilters_d(const QVector &avg, const QVector &noise);
  float Amplitude(QVector &pSamples, float &amp) ;
  float PulseAmplitude(QVector &pSamples) ;
  float GetTVR();
  float GetTVL();
  float GetRT(){return RT;};

  float GetDT(){return DT;}
 
  float GetAdc2mV(){return pow(2.,-16.)*pow(10.,4);}; 
  float GetDelay(){return data.delay; };//fDelay;}
  float GetFracDelay(){return data.frac_delay; };//fDelay;}

  void SetTriggerPos(int trigger) { fTriggerPos = trigger; }; 




  void  ShapeParameters(QVector &SamplesFiltered);
  QVector PolFit(QVector &x,QVector& y);
  QError PowerSpectrum(const QVector &pSamples, QVector &samplesFiltered);
  QError Filter(const QVector &pSamples, QVector &samplesFiltered);
  QError Filter_d(const QVector &pSamples, QVector &samplesFiltered);
  QError FAmplitude(const QVector &pSamples, double &amp);

  float MultiplePulseCorrection(QVector &, QVector & , QVector &,bool);
  float GetAmplitude(){return data.amplitude;};
  float GetSecondAmplitude(){return data.second_amp;}; 
  float GetSecondDelay(){return data.second_delay;}; 
  void CallPileUp(){ fMultiPulse = true; } ; 
  void BiPol(QVector &, QVector &, QVector &); 
  const QVector& averagePulse(){return AvgPulse;}
  const QVector& averagePulse_f(){return AvgPulse_f;}
  const QVector& noiseSpectrum(){return fPS;}

  QVector ReconstructionFit(QVector & , QVector &, QVector &);

  QVector GetReconstructedFit(){return fReconstructedPulse;}; 


  void Decay_and_Rise(const QVector &);


  QVector fShift(QVector &, float &); 
  float GetBLS(){return fBS;};
void SetPulseInfo(const QVector &); 
//float Correlate(QVector & pSamples);
private:
// real and imaginary part of the average spectrum
//
PulseInfo data; 
  QVector AvgPulse;             
  int nSigPulses;
  
  // noise power spectrum
  QVector fPS;
  bool fMultiPulse; 
  int nNoisePS;
  unsigned int avg_nmax,fTriggerPos;
  //number of points used for TVL&TVR
  int rms_left, rms_right,fMaxRange,fSize, fBaseLineStart,fBaseLineNpt; 

  //conversion factors and delay 
  double conv,conv_d,fDelay;

  //shape parameters
  float tvl,tvr,RT,DT,fBS;

  // filter transfer function 
  QVectorC H;// real part
  // filtered average pulse
  QVector AvgPulse_f;
  //  Sampler * sam; 
  QVector fReconstructedPulse;

  };

#endif

---