source: fact/tools/rootmacros/testfir.C@ 13408

#include "TRandom2.h"
#include "TH1F.h"
#include "TCanvas.h"

#include <stdio.h>
#include <iostream>


#define FAD_MAX_SAMPLES 1024

vector<float> Vmeas(FAD_MAX_SAMPLES);
vector <float> Vfir(FAD_MAX_SAMPLES);
vector <float> Vpulse(FAD_MAX_SAMPLES);
#include "factfir.C"


#define k 16
double a[k];
double b=k;

void testfir(){

        for (int j=0; j<k; j++){
                a[j]=1;
        }

        TRandom2 myrandom;
        myrandom.RndmArray(1024,Vmeas);
        
        int pulse_start=350;
        float pulse_height=0.6;
        int pulse_len = 60;
        for (int i = pulse_start; i<pulse_start + pulse_len; i++){
                Vpulse[i] += pulse_height/pulse_len * (pulse_len-(i-pulse_start));
                // in cas e you want the pulse in the data
                Vmeas[i] += Vpulse[i];   
        }


        factfir(b,a,k, Vfir);

        TH1F *meas = new TH1F("Vmeas","Vmeas test",FAD_MAX_SAMPLES, -0.5 , FAD_MAX_SAMPLES-0.5);
        TH1F *fir = new TH1F("Vfir","Vfir test",FAD_MAX_SAMPLES, -0.5 , FAD_MAX_SAMPLES-0.5);
        TH1F *pulse = new TH1F("Vpulse","Vpulse test",FAD_MAX_SAMPLES, -0.5 , FAD_MAX_SAMPLES-0.5);

        for (int i=0; i<FAD_MAX_SAMPLES; i++){
                meas->SetBinContent(i,Vmeas[i]);
                fir->SetBinContent(i,Vfir[i]);
                pulse->SetBinContent(i,Vpulse[i]);
        }

        TCanvas *c = new TCanvas();
        c->Divide(1,3);
        c->cd(1);
        meas->Draw("HIST");
        c->cd(2);
        fir->Draw("HIST");
        c->cd(3);
        pulse->Draw("HIST");


}