#include <unistd.h>
#include <TMath.h>
#include <stdlib.h>
#include <vector>
#include "TRobustEstimator.h"
#include "TRandom3.h"

void testTRE(){

	const int nvectors=1000;
	const int nvariables=7;

	TRandom3 *rand = new TRandom3(0);

	TRobustEstimator *RE= new TRobustEstimator(nvectors,nvariables);
	TMatrixDSym *CovMat= new TMatrixDSym(nvariables);

	Double_t var;

	for (int i=0; i<nvectors; i++){
		printf("Row n. %i\n",i);
		Double_t Data[nvariables];
		for (int j=0; j< nvariables; j++){
			var=rand->Gaus(0,0.5);
			Data[j]=var;
			printf("%f\n", Data[j]);
		}
		RE->AddRow(Data);
	}

	printf("here\n");
	RE->Evaluate();
	RE->GetCovariance(*CovMat);
	CovMat->Print();

}