TString macroName = gROOT->GetTutorialDir();
macroName.Append("/foam/foam_demo.C");
gROOT->ProcessLine(TString::Format(".L %s+",macroName.Data()));

cout<<"====================== TestVector ================================"<<endl;
TFile fileA("foam_demo.root");
fileA.cd();
cout<<"------------------------------------------------------------------"<<endl;
fileA.ls();
cout<<"------------------------------------------------------------------"<<endl;
fileA.Map();
cout<<"------------------------------------------------------------------"<<endl;
fileA.ShowStreamerInfo();
cout<<"------------------------------------------------------------------"<<endl;
fileA.GetListOfKeys()->Print();
cout<<"------------------------------------------------------------------"<<endl;

TFoam  *FoamX = (TFoam*)fileA.Get("FoamX");

FoamX->CheckAll(1);

TFoamIntegrand * rho = (TFoamIntegrand*) gROOT->ProcessLine("return new TFDISTR();");
FoamX->SetRho(rho);

Double_t MCvect[2]; // 2-dim vector generated in the MC run

for(long loop=0; loop<50000; loop++){
   FoamX->MakeEvent();            // generate MC event
   FoamX->GetMCvect( MCvect);     // get generated vector (x,y)
   Double_t x=MCvect[0];
   Double_t y=MCvect[1];
   if(loop<10) cout<<"(x,y) =  ( "<< x <<", "<< y <<" )"<<endl;
   }// loop

Double_t IntNorm, Errel;
FoamX->Finalize(   IntNorm, Errel);     // final printout
Double_t MCresult, MCerror;
FoamX->GetIntegMC( MCresult, MCerror);  // get MC integral, should be one
cout << " MCresult= " << MCresult << " +- " << MCerror <<endl;
cout<<"===================== TestPers FINISHED ======================="<<endl;
return 0;