PndEmcPSAFPGAFilterMWD.cxx

Go to the documentation of this file.

#include "PndEmcPSAFPGAFilterMWD.h"
#include <math.h>

PndEmcPSAFPGAFilterMWD::PndEmcPSAFPGAFilterMWD() :
        buffer(0)
        ,bpointer(0)
{
}

PndEmcPSAFPGAFilterMWD::~PndEmcPSAFPGAFilterMWD()
{
        buffer.resize(0);
}

void PndEmcPSAFPGAFilterMWD::set(float tau, unsigned int newBufferSize)
{
        buffer.resize(newBufferSize);
        for(unsigned int i=0; i<newBufferSize; i++)
                buffer[i]=0;
        bpointer=0;
        sum = 0.0;
        lambda = 1.0/tau;
}

//#include <iostream>
float PndEmcPSAFPGAFilterMWD::put(float valueToStore){
        //std::cout<<"buffer# "<<buffer.size()<<std::endl;
        bpointer %= buffer.size();
        float diff = valueToStore - buffer[bpointer];
        float mwd = diff + lambda*sum;
        sum += diff;
        buffer[bpointer] = valueToStore;
        bpointer++;
        return mwd;
}

void PndEmcPSAFPGAFilterMWD::resetToZero() {
  for(unsigned int i=0; i<buffer.size(); i++) {
    buffer[i]=0;
  }
    bpointer = 0.;
    sum = 0.;
}