#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");


}