{

TFile * outputfile = TFile::Open ("outputTMVA.root","RECREATE");

TMVA::Factory * TMVAtest = new TMVA::Factory ("TMVAtest", outputfile, "S") ;

//PG get the signal and deliver it to the TMVA factory

TChain signalTree ("sample") ;

signalTree.Add ("data/sig_0.root") ;

std::cout << "READ " << signalTree.GetEntries () << " signal events\n" ;

TMVAtest->AddSignalTree (&signalTree, 1) ;

//PG get the bkg and deliver it to the TMVA factory

TChain bkgTree ("sample") ;

bkgTree.Add ("data/bkg_0.root") ;

std::cout << "READ " << bkgTree.GetEntries () << " bkg events\n" ;

TMVAtest->AddBackgroundTree (&bkgTree, 1) ;

//PG get the training and test samples and deliver them to the TMVA factory

TChain signalTrainTree ("sample") ;

signalTrainTree.Add ("data/sig_1.root") ;

std::cout << "READ " << signalTrainTree.GetEntries () << " signal train events\n" ;

TChain bkgTrainTree ("sample") ;

bkgTrainTree.Add ("data/bkg_1.root") ;

std::cout << "READ " << bkgTrainTree.GetEntries () << " bkg train events\n" ;

TMVAtest->SetInputTrees (signalTrainTree.GetTree (), bkgTrainTree.GetTree (), 1., 1.) ;

//PG variables to be used for the selection

//PG must be defined in the TTrees

TMVAtest->AddVariable ("vars.x", 'F') ;

TMVAtest->AddVariable ("vars.y" , 'F') ;

int signalNumTrain = signalTrainTree.GetEntries () * 4 / 5 ;

int bkgNumTrain = bkgTrainTree.GetEntries () * 4 / 5 ;

int signalNumTest = signalTrainTree.GetEntries () - signalNumTrain ;

int bkgNumTest = bkgTrainTree.GetEntries () - bkgNumTrain ;

char trainOptions[120] ;

sprintf (trainOptions,"NSigTrain=%d:NBkgTrain=%d:NSigTest=%d:NBkgTest=%d",

signalNumTrain, bkgNumTrain,

signalNumTest, bkgNumTest) ;

sprintf (trainOptions,"NSigTrain=%d:NBkgTrain=%d:NSigTest=%d:NBkgTest=%d",

0,0,0,0) ;

std::cout << "TRAINING CONFIGURATION : " << trainOptions << "\n" ;

TMVAtest->PrepareTrainingAndTestTree ("",trainOptions) ;

//PG prepare the classifier

//PG cut-based, default params

TMVAtest->BookMethod (TMVA::Types::kCuts, "Cuts") ;

TMVAtest->TrainAllMethods () ;

TMVAtest->TestAllMethods () ;

TMVAtest->EvaluateAllMethods () ;

delete TMVAtest ;

delete outputfile ;

}