void makeclassification() {
Float_t *vars = new Float_t[10];
//initialize TMVA Reader (example here is diphoton mva from higgs->gamma gamma mva analysis)
TMVA::Reader* tmva = new TMVA::Reader();
tmva->AddVariable("masserrsmeared/mass", &vars[0]);
tmva->AddVariable("masserrsmearedwrongvtx/mass", &vars[1]);
tmva->AddVariable("vtxprob", &vars[2]);
tmva->AddVariable("ph1.pt/mass", &vars[3]);
tmva->AddVariable("ph2.pt/mass", &vars[4]);
tmva->AddVariable("ph1.eta", &vars[5]);
tmva->AddVariable("ph2.eta", &vars[6]);
tmva->AddVariable("TMath::Cos(ph1.phi-ph2.phi)" , &vars[7]);
tmva->AddVariable("ph1.idmva", &vars[8]);
tmva->AddVariable("ph2.idmva", &vars[9]);
tmva->BookMVA("BDTG","/afs/cern.ch/user/b/bendavid/cmspublic/diphotonmvaApr1/HggBambu_SMDipho_Jan16_BDTG.weights.xml");
//tmva->BookMVA("BDTG","/scratch/bendavid/root/HggBambu_SMDipho_Jan16_BDTG.weights.xml");
TMVA::MethodBDT *bdt = dynamic_cast<TMVA::MethodBDT*>(tmva->FindMVA("BDTG"));
//enable root i/o for objects with reflex dictionaries in standalone root mode
ROOT::Cintex::Cintex::Enable();
//open output root file
TFile *fout = new TFile("gbrtest.root","RECREATE");
//create GBRForest from tmva object
GBRForest *gbr = new GBRForest(bdt);
//write to file
fout->WriteObject(gbr,"gbrtest");
fout->Close();
}
