{
//compile and load libraries for d0root_analysis, em_util and jetcorr
gInterpreter->ExecuteMacro ("d0root_analysis/compile.C");
gInterpreter->ExecuteMacro ("d0root_tmbtree/compile.C");

// load JES stuff
//gSystem->Load("jetcorr/src/JetResponse_cpp.so");
//gSystem->Load("jetcorr/src/JetOffset_cpp.so");
//gSystem->Load("jetcorr/src/JetOOC_cpp.so");
//gSystem->Load("jetcorr/src/JetDecay_cpp.so");
//gSystem->Load("jetcorr/src/ParticleJetCorr_cpp.so");

//load d0root_btag algorithms
gSystem->Load("libd0root_btag_algorithms.so");

//compile and load libraries for V0 removal
gInterpreter->ExecuteMacro ("d0root_csip/compile.C");

//compile and load libraries for JLIP tagger
 gInterpreter->ExecuteMacro ("d0root_jlip/compile.C");

//load analysis code
gROOT->ProcessLine(".L ./d0root_btag/analysis/bbAnalysisI.cpp+g");


///----------------------------------------------------
/// Input file
///----------------------------------------------------


TChain ch("TMBTree");

//Input file
ch.Add("<filename1>");



///----------------------------------------------------
/// Instanciate analysis class
///----------------------------------------------------
cout << TMBTree->GetEntries() << " events to process" << endl << endl; 
bbAnalysisI* b = new bbAnalysisI(TMBTree, 14.00);
 b->Analyze("<filename2>.root",15,3.0);
}