{
//////////////////////////////////////////////////////////
//
// SbottomAnalysis
//
// Author: Arnaud Duperrin & Eric Kajfasz      22-Dec-2003
// 
// Purpose: run Sbottom Analysis on TMBTree
//          
//          
//////////////////////////////////////////////////////////

// d0root_analysis and d0root_tmbtree
gInterpreter->ExecuteMacro ("d0root_analysis/compile.C");
gInterpreter->ExecuteMacro ("d0root_tmbtree/compile.C");

//b-tagging 
gSystem->Load("libd0root_btag_algorithms.so");
gInterpreter->ExecuteMacro ("d0root_csip/compile.C");
gInterpreter->ExecuteMacro ("d0root_jlip/compile.C");

//TRF 
gROOT->ProcessLine(".L sb_analysis/b-tag/pdgid.C+g");
gROOT->ProcessLine(".L sb_analysis/b-tag/JLIP_TRF/JLIP_TRF.C+g");

//JET_FLAVOR
gROOT->ProcessLine(".L sb_analysis/b-tag/JLIP/JET_FLAVOR.cpp+g");

//JLIP
gROOT->ProcessLine(".L sb_analysis/b-tag/JLIP/JLIP.cpp+g");


//load good and bad run list 
gROOT->ProcessLine(".L ./sb_analysis/runs/MapOfRuns.cpp+g");
gSystem->Exec("rm -f runs");
gSystem->Exec("ln -s sb_analysis/runs");

//load analysis code
gROOT->ProcessLine(".L ./sb_analysis/analysis/BranchDefs.cpp+g");
gROOT->ProcessLine(".L ./sb_analysis/analysis/SbottomAnalysis.cpp+g");


///----------------------------------------------------
/// TMBTree Input file
///----------------------------------------------------
TChain ch("TMBTree");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/");

// 121
ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/data_subskim03.root");

// Sb MC
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/TMBTree_sb125_n40_121sel_final.root");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/TMBTree_sb130_n65_121sel_final.root");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/TMBTree_sb160_n80_121sel_final.root");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/TMBTree_sb180_n60_121sel_final.root");

// SM MC
//ch.Add("/work/pariou-clued0/duperrin/MC/wenu/mc_wenuj_121sel.root");
//ch.Add("/work/pariou-clued0/duperrin/MC/wenu/mc_wenujj_121sel.root");

//ch.Add("/work/pariou-clued0/duperrin/MC/wmunu/mc_wmunuj_121sel.root"); // 0 evts
//ch.Add("/work/pariou-clued0/duperrin/MC/wmunu/mc_wmunujj_121sel.root");

//ch.Add("/work/pariou-clued0/duperrin/MC/wtaunu/mc_wtaunu_121sel.root");
//ch.Add("/work/pariou-clued0/duperrin/MC/wtaunu/mc_wtaunuj_121sel.root");
//ch.Add("/work/pariou-clued0/duperrin/MC/wtaunu/mc_wtaunujj_121sel.root");

//ch.Add("/work/pariou-clued0/duperrin/MC/zmumu/mc_zmumuj_121sel.root");
//ch.Add("/work/pariou-clued0/duperrin/MC/zmumu/mc_zmumujj_121sel.root");

//ch.Add("/work/pariou-clued0/duperrin/MC/znunu/mc_znunuj_121sel.root"); // bug
//ch.Add("/work/pariou-clued0/duperrin/MC/znunu/mc_znunujj_121sel.root");

//ch.Add("/work/pariou-clued0/duperrin/MC/ztautau/mc_ztautauj_121sel.root");

//ch.Add("/work/pariou-clued0/duperrin/MC/tt_wq_wq_3q_taunuq/TMBTree/mc_tt_wq_wq_3q_taunuq_121sel.root");

// SM MC b
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/mc_wenu_bb_jj_121sel_final.root");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/mc_wmunu_bb_121sel_final.root");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/mc_wtaunu_bb_121sel_final.root");

//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/mc_zee_bb_121sel_final.root");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/mc_zmumu_bb_121sel_final.root");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/mc_znunu_bb_121sel_final.root");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/mc_ztautau_bb_121sel_final.root");

//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/mc_wc_j_121sel_final.root");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/mc_wc_jj_121sel_final.root");
//ch.Add("/work/pariou-clued0/duperrin/LAL/selection121events/LQselection121events/TMBTreeFinal/mc_wc_jjj_121sel_final.root");

//ch.Add("/rooms/in2p3/projects/lebihan/ztautau_isomu_filter/TMBAnalyze_x-p14.05.02-Linux-ztautau_8358_list-847020/ztautau_8358_p14-d0correctv06.root");


cout << TMBTree->GetEntries() << " events to process in the TMBTree" << endl << endl; 


///----------------------------------------------------
/// Instanciate analysis class
///----------------------------------------------------
SbottomAnalysis* sb = new SbottomAnalysis(TMBTree, 14.0);

sb->SetDebug(false);
sb->SetNstarting(0);
//sb->SetNentries(2);

//TRF SM or Sbottom :
bool IsMC_Sbottom = false;
sb->SetIsMC_Sbottom(IsMC_Sbottom);
gSystem->Exec("rm -f p14_JLIP_TRF.root");
if (IsMC_Sbottom) {
  gSystem->Exec("ln -s sb_analysis/b-tag/sbottom/p14_JLIP_TRF_sbottom.root p14_JLIP_TRF.root");
} else {
  gSystem->Exec("ln -s sb_analysis/b-tag/JLIP_TRF/p14_JLIP_TRF.root p14_JLIP_TRF.root");
}

int fileNumber = 1;
char* jetAlg = "corrJCCB";
bool beamConstraint = false;


sb->Process("DATA",fileNumber,jetAlg,beamConstraint);
//sb->Process("MC",fileNumber,jetAlg,false);

}