void QGAnalysis::InitJJ(){
Log(__FUNCTION__,"INFO","Booking Histos for DiJet");
for ( string l : AllLabel() ) {
Book ("QGAnalysis/CutFlow/CutFlowDiJet_"+ l ,"CutFlow;Events",10,-.5,10-.5);
//
Book ("QGAnalysis/DiJet/Npv_"+ l ,"Npv jj", 50,0,50);
//
for (const string& t : jetTypes ) {
for (const string& v : jetVars ){
if ( l == "Data" and t != "U") continue; // Data have no MC info
float xmin=-1, xmax=1.0;
int nbins=200;
if ( v == "axis2") { xmin=0; xmax=10;}
if ( v == "mult") { xmin=0; xmax=100;nbins=100;}
Book ("QGAnalysis/DiJet/"+v +"_"+ t+"_"+ l ,v+"_" + t, 200,-1,1.0);
for(int ptb=0;ptb < ptBins.size()-1; ++ptb)
for(int aetab=0;aetab < aetaBins.size()-1; ++aetab)
{
Book ("QGAnalysis/DiJet/"+v+"_"+t+Form("_pt%.0f_%.0f",ptBins[ptb],ptBins[ptb+1]) + Form("_eta%.1f_%.1f",aetaBins[aetab],aetaBins[aetab+1])+"_"+ l ,v +"_" + t, nbins,xmin,xmax);
} // pt-eta
} // vars
} // types
} // labels
}
void TauOptimization::Init(){
cout<<"[TauOptimization]::[Init]::[INFO] Booking Histo "<<endl;
for ( string& l : AllLabel() )
{
Book2D("TauOpt/Iso_Rho_" + l , "tau Iso delta Beta corrections",1000,0,200,1000,0,20);
Book2D("TauOpt/Iso2_Rho_" + l , "tau Iso2 delta Beta corrections",1000,0,200,1000,0,20);
Book("TauOpt/Iso_Match_"+l,"tau Iso",1000,0.,20.);
Book("TauOpt/Iso_Lead_"+l,"tau Iso",1000,0.,20.);
Book("TauOpt/IsLeadMatched_"+l,"Leading Matched",10,-4.5,5.5);
Book("TauOpt/Pt_Match_"+l,"tau Pt",1000,0.,1000.);
Book("TauOpt/Pt_Lead_"+l,"tau Pt",1000,0.,1000.);
Book("TauOpt/Iso2_Match_"+l,"tau Iso Delta Beta",1500,-10,20.);
Book("TauOpt/Iso2_Lead_"+l,"tau Iso Delta Beta",1500,-10,20.);
}
string treename = "TauOptTree";
InitTree(treename);
Branch(treename, "tau1Pt" ,'F');
Branch(treename, "tau2Pt" ,'F');
Branch(treename, "tau1Eta" ,'F');
Branch(treename, "tau2Eta" ,'F');
Branch(treename, "tau1Iso" ,'F');
Branch(treename, "tau2Iso" ,'F');
Branch(treename, "ttM" ,'F');
}
void DYAnalysis::Init(){
for ( string l : AllLabel() ) {
Log(__FUNCTION__,"INFO","Booking Histo Mass");
Book ("DYAnalysis/Vars/Mee_"+ l ,"Mee;m^{ee} [GeV];Events" , 100,50,200);
Book ("DYAnalysis/Vars/Mmm_"+ l ,"Mmm;m^{#mu#mu} [GeV];Events", 100,50,200);
Book ("DYAnalysis/Vars/Mem_"+ l ,"Mem;m^{e#mu} [GeV];Events" , 100,50,200);
//
Book ("DYAnalysis/Vars/MuonIso_"+ l ,"Muon Isolation;Iso^{#mu} [GeV];Events", 1000,0,0.1);
Book ("DYAnalysis/Vars/Ptee_"+ l ,"Ptee;p_{T}^{ee} [GeV];Events" , 1000,0,1000);
Book ("DYAnalysis/Vars/Ptmm_"+ l ,"Ptmm;p_{T}^{#mu#mu} [GeV];Events", 1000,0,1000);
Book ("DYAnalysis/Vars/Ptem_"+ l ,"Ptem;p_{T}^{e#mu} [GeV];Events" , 1000,0,1000);
Book ("DYAnalysis/Vars/PtemNoMCut_"+ l ,"Ptem;p_{T}^{e#mu} [GeV];Events" , 1000,0,1000);
// counting unweighted events
Book ("DYAnalysis/Vars/PtmmUW_"+ l ,"Ptmm;p_{T}^{#mu#mu} [GeV];Events", 1000,0,1000);
// trilepton
Book ("DYAnalysis/Vars/Ptlll_"+ l ,"Ptlll;p_{T}^{ll} [GeV];Events" , 100,0,1000);
// eta
Book ("DYAnalysis/Vars/EtaMu1_"+ l ,"EtaMu1;#eta^{#mu1} ;Events", 1000,-5,5);
//
Book ("DYAnalysis/Vars/NJee_"+ l ,"NJee;N_{jets}^{ee};Events" , 10,0,10);
Book ("DYAnalysis/Vars/NJmm_"+ l ,"NJmm;N_{jets}^{#mu#mu};Events", 10,0,10);
Book ("DYAnalysis/Vars/NJem_"+ l ,"NJem;N_{jets}^{e#mu};Events" , 10,0,10);
//
Book ("DYAnalysis/Vars/Npvee_"+ l ,"Npvee", 50,0,50);
Book ("DYAnalysis/Vars/Npvmm_"+ l ,"Npvmm", 50,0,50);
Book ("DYAnalysis/Vars/Npvem_"+ l ,"Npvem", 50,0,50);
//
Book ("DYAnalysis/Vars/MHighPtee_"+ l ,"MHighPtee", 1000,0,1000);
Book ("DYAnalysis/Vars/MHighPtmm_"+ l ,"MHighPtmm", 1000,0,1000);
Book ("DYAnalysis/Vars/MJHighPtee_"+ l ,"MJHighPtee", 1000,0,1000);
Book ("DYAnalysis/Vars/MJHighPtmm_"+ l ,"MJHighPtmm", 1000,0,1000);
//
Book ("DYAnalysis/Vars/MCloseMt_"+ l ,"MCloseMt", 1000,0,1000);
Book ("DYAnalysis/Vars/MCloseJZ_"+ l ,"MCloseJZ", 1000,0,1000);
//
Book ("DYAnalysis/Vars/Pt20mm_"+ l ,"Ptmm20", 1000,0,1000);
Book ("DYAnalysis/Vars/Accmm_"+ l ,";Events", 10,0,10);
{
GetHisto("DYAnalysis/Vars/Accmm_"+ l,"") ->GetXaxis() ->SetBinLabel(1,"Gen-Reco");
GetHisto("DYAnalysis/Vars/Accmm_"+ l,"") ->GetXaxis() ->SetBinLabel(2,"Gen-!Reco");
GetHisto("DYAnalysis/Vars/Accmm_"+ l,"") ->GetXaxis() ->SetBinLabel(3,"!Gen-Reco");
}
Book ("DYAnalysis/Vars/Accee_"+ l ,";Events", 10,0,10);
{
GetHisto("DYAnalysis/Vars/Accee_"+ l,"") ->GetXaxis() ->SetBinLabel(1,"Gen-Reco");
GetHisto("DYAnalysis/Vars/Accee_"+ l,"") ->GetXaxis() ->SetBinLabel(2,"Gen-!Reco");
GetHisto("DYAnalysis/Vars/Accee_"+ l,"") ->GetXaxis() ->SetBinLabel(3,"!Gen-Reco");
}
}
}
void ChargedHiggsTauNu::initGen(){
for ( string& l : AllLabel() ) {
Book( "ChargedHiggsTauNu/Gen/CutFlow_"+l,"CutFlow",100,0,100);
GetHisto("ChargedHiggsTauNu/Gen/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(1, "All");
GetHisto("ChargedHiggsTauNu/Gen/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(2, "H#pm");
GetHisto("ChargedHiggsTauNu/Gen/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(3, "H#pm#rightarrow#tau");
GetHisto("ChargedHiggsTauNu/Gen/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(4, "H#pm#rightarrow#tau^{h} (B)");
GetHisto("ChargedHiggsTauNu/Gen/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(5, "(B) + Lep Veto Gen Level");
GetHisto("ChargedHiggsTauNu/Gen/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(6, "(B) + LooseIso Veto");
GetHisto("ChargedHiggsTauNu/Gen/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(7, "(B) + MiniIso Veto");
}
}
void ChargedHiggsHW::Init()
{
for ( string l : AllLabel() ) {
cout <<"[ChargedHiggsHW]::[Init]::[INFO] Boking Histo CutFlow_" <<l<<endl;
Book( ("ChargedHiggsHW/CutFlow/CutFlow_"+ l ) .c_str(), ("CutFlow "+ l).c_str(),100,-.5,100-.5);
GetHisto(("ChargedHiggsHW/CutFlow/CutFlow_"+l).c_str(),"")->GetXaxis()->SetBinLabel(1,"Total");
//Histograms should be booked here for details see src/AnalysisChargedHiggsTauNu.cpp
}
}
void ChargedHiggsQCDPurity::Init()
{
// -- simplified selection
for ( string& l : AllLabel() )
for (size_t iBin = -1 ; iBin + 1 < PtBins.size() ; ++iBin )
{
float pt = -1;
if (iBin>=0 ) pt= PtBins[iBin];
// direct, fullSel
Book( dir + HistName(pt, true , false)+"_"+ l , ("EtMiss "+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, false)+"_"+ l , ("EtMissIsoInv "+ l).c_str(),250,0.,500.);
// ----
Book( dir + HistName(pt, true , false, "Upar")+"_"+ l , ("EtMissParallel "+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, false, "Upar")+"_"+ l , ("EtMissParallelIsoInv "+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, true , false, "Uperp")+"_"+ l , ("EtMissPerp "+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, false, "Uperp")+"_"+ l , ("EtMissPerpIsoInv "+ l).c_str(),250,0.,500.);
// Used in case of spline reweighting iterations
Book( dir + HistName(pt, true, false, "TauPt")+"_"+ l , ("PtTau "+ l).c_str(),1000,0.,1000.);
Book( dir + HistName(pt, false, false, "TauPt")+"_"+ l , ("PtTauIsoInv "+ l).c_str(),1000,0.,1000.);
// QG ? Tau
Book( dir + HistName(pt, true , false)+"_Q_" + l , ("EtMiss Q"+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, false)+"_Q_" + l , ("EtMissIsoInv Q"+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, true , false)+"_G_" + l , ("EtMiss G"+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, false)+"_G_" + l , ("EtMissIsoInv G"+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, true , false)+"_T_" + l , ("EtMiss G"+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, false)+"_T_" + l , ("EtMissIsoInv G"+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, true , false)+"_U_" + l , ("EtMiss G"+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, false)+"_U_" + l , ("EtMissIsoInv G"+ l).c_str(),250,0.,500.);
// Prong
Book( dir + HistName(pt, true , false)+"_1p_"+ l , ("EtMiss "+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, false)+"_1p_"+ l , ("EtMissIsoInv "+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, true , false)+"_3p_"+ l , ("EtMiss "+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, false)+"_3p_"+ l , ("EtMissIsoInv "+ l).c_str(),250,0.,500.);
}
// --- for event NOT BINNED! CONTROL PLOTS
for (string &l : AllLabel() )
{
// Cut flow with QCD
Book( "ChargedHiggsTauNu/CutFlowQCD/CutFlowQCD_"+l , ("CutFlowQCD "+ l).c_str(), 10, -.5, 10.0-.5);
GetHisto("ChargedHiggsTauNu/CutFlowQCD/CutFlowQCD_"+l,"")->GetXaxis()->SetBinLabel(1, "N_{taus} #geq 1 \\ N_{leps} = 0 \\ N_{jets} #geq 3 \\ Trigger");
GetHisto("ChargedHiggsTauNu/CutFlowQCD/CutFlowQCD_"+l,"")->GetXaxis()->SetBinLabel(2,"N_{Bjets} #geq 1");
GetHisto("ChargedHiggsTauNu/CutFlowQCD/CutFlowQCD_"+l,"")->GetXaxis()->SetBinLabel(3,"MET > 60 GeV");
GetHisto("ChargedHiggsTauNu/CutFlowQCD/CutFlowQCD_"+l,"")->GetXaxis()->SetBinLabel(4,"R_{bb}^{Min}>40");
// Used for control plots
string hist= "TauPt_IsoInv_Control";
Book( dir + hist+"_"+ l , ("PtTauIsoInv "+ l).c_str(),1000,0.,1000.); // this guys has the R factor applied in the loos selection
hist= "TauPt_Control";
Book( dir + hist+"_"+ l , ("PtTau "+ l).c_str(),1000,0.,1000.); // this guys has the R factor applied in the loos selection
hist= "EtMiss_IsoInv_Control";
Book( dir + hist+"_"+ l , ("EtMissTauIsoInv "+ l).c_str(),250,0,500);
hist= "EtMiss_Control";
Book( dir + hist+"_"+ l , ("EtMiss "+ l).c_str(),250,0,500);
}
// -- full selection
for ( string& l : AllLabel() )
{
for (size_t iBin = -1 ; iBin + 1 < PtBins.size() ; ++iBin )
{
float pt = -1;
if (iBin>=0 ) pt= PtBins[iBin];
// direct full
Book( dir + HistName(pt, true , true)+"_"+ l , ("EtMiss "+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, true)+"_"+ l , ("EtMissIsoInv "+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, true , true)+"_1p_"+ l , ("EtMiss "+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, true)+"_1p_"+ l , ("EtMissIsoInv "+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, true , true)+"_3p_"+ l , ("EtMiss "+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, true)+"_3p_"+ l , ("EtMissIsoInv "+ l).c_str(),250,0.,500.);
// --- NoR
Book( dir + HistName(pt, false, true)+"_NoR_"+ l , ("EtMissIsoInv "+ l).c_str(),250,0.,500.);
//
Book( dir + HistName(pt, true , true, "Upar")+"_"+ l , ("EtMissParallel "+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, true, true, "Uperp")+"_"+ l , ("EtMissPerp "+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, false , true, "Upar")+"_"+ l , ("EtMissParallelIsoInv "+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, true, "Uperp")+"_"+ l , ("EtMissPerpIsoInv "+ l).c_str(),250,0.,500.);
// QG ?
Book( dir + HistName(pt, true , true )+"_Q_" + l , ("EtMiss Q"+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, true )+"_Q_" + l , ("EtMissIsoInv Q"+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, true , true )+"_G_" + l , ("EtMiss G"+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, true )+"_G_" + l , ("EtMissIsoInv G"+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, true , true )+"_T_" + l , ("EtMiss G"+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, true )+"_T_" + l , ("EtMissIsoInv G"+ l).c_str(),250,0.,500.);
Book( dir + HistName(pt, true , true )+"_U_" + l , ("EtMiss G"+ l).c_str(),250,0.,500);
Book( dir + HistName(pt, false, true )+"_U_" + l , ("EtMissIsoInv G"+ l).c_str(),250,0.,500.);
}
// I don't need to split it by pt
Book( dir + "Mt"+"_"+ l , ("Mt "+ l).c_str(),8000,0.,8000); // same binning in TauNu
Book( dir + "MtIsoInv"+"_"+ l , ("MtIsoInv "+ l).c_str(),8000,0.,8000.);
Book( none + "EtMissIsoInv"+"_"+ l , ("EtMissIsoInv "+ l).c_str(),1000,0.,1000.);
Book( none + "EtMiss"+"_"+ l , ("EtMiss "+ l).c_str(),1000,0.,1000.); // copy of the Tau Nu ?
Book( none+"RbbMin_"+l,"RbbMin "+l+";R_{bb}^{min}",100,0,2*TMath::Pi());
Book( none+"RCollMin_"+l,"RCollMin "+l+";R_{coll}^{min}",100,0,2*TMath::Pi());
Book( none+"RbbMinIsoInv_"+l,"RbbMin IsoInv "+l+";R_{bb}^{min}",100,0,2*TMath::Pi());
Book( none+"RCollMinIsoInv_"+l,"RCollMin IsoInv "+l+";R_{coll}^{min}",100,0,2*TMath::Pi());
Book( none+"NBjets_"+l,"NBjets "+l+";NBjets",1000,0,1000);
Book( none+"NBjetsIsoInv_"+l,"NBjets IsoInv "+l+";NBjets",1000,0,1000);
Book( none+"Tau1Pt_"+l,"Tau1Pt "+l+";Tau1Pt",1000,0,1000);
Book( none+"Tau1PtIsoInv_"+l,"Tau1Pt IsoInv "+l+";Tau1Pt",1000,0,1000);
}
}
void ChargedHiggsTauNu::Init()
{
Log(__FUNCTION__,"INFO",Form("nProngs selected=%d",nprongs) );
if ( doGen ) initGen();
for ( string& l : AllLabel() ) {
Log(__FUNCTION__,"INFO", "Booking Histo CutFlow_" + l);
Book( "ChargedHiggsTauNu/CutFlow/CutFlow_"+ l , ("CutFlow "+ l).c_str(),100,-.5,100-.5);
{
GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(Total+1,"Total"); // bin shift
GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(OneTau+1,"N_{taus} #geq 1");
GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(NoLep+1,"N_{leps} = 0");
GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(ThreeJets,"N_{jets} #geq 3");
GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(OneBjet,"N_{Bjets} #geq 1");
GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(Met,"MET > 60 GeV");
//GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(7,"R_{sr}^{Max} < 1.96");
//GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(7,"R_{sr}^{Max} < 140");
GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(AngColl,"R_{coll}^{Min}>40");
GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(AngRbb,"R_{bb}^{Min}>40");
//GetHisto("ChargedHiggsTauNu/CutFlow/CutFlow_"+l,"")->GetXaxis()->SetBinLabel(7,"#Delta #phi (Tau,Jet1)<2.09");
}
// Cut flow with QCD
Book( "ChargedHiggsTauNu/CutFlowQCD/CutFlowQCD_"+l , ("CutFlowQCD "+ l).c_str(), 100, -.5, 100-.5);
{
GetHisto("ChargedHiggsTauNu/CutFlowQCD/CutFlowQCD_"+l,"")->GetXaxis()->SetBinLabel(1, "N_{taus} #geq 1 \\ N_{leps} = 0 \\ N_{jets} #geq 3 \\ Trigger");
GetHisto("ChargedHiggsTauNu/CutFlowQCD/CutFlowQCD_"+l,"")->GetXaxis()->SetBinLabel(2,"N_{Bjets} #geq 1");
GetHisto("ChargedHiggsTauNu/CutFlowQCD/CutFlowQCD_"+l,"")->GetXaxis()->SetBinLabel(3,"MET > 60 GeV");
GetHisto("ChargedHiggsTauNu/CutFlowQCD/CutFlowQCD_"+l,"")->GetXaxis()->SetBinLabel(4,"R_{bb}^{Min}>40");
}
// -- cut flow 2 is the one used to synchronized with HIP
Book( "ChargedHiggsTauNu/CutFlow/CutFlow2_"+ l , ("CutFlow "+ l).c_str(),100,-.5,100-.5);
// study on base selection
Book( "ChargedHiggsTauNu/Base/Tau1Pt_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_0Pi_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_1Pi_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_2Pi_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_0Pi_TrSingleTau_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_1Pi_TrSingleTau_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_0Pi_TrMet110_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_1Pi_TrMet110_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_TrMet110_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_Met150_TrMet110_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_TrSingleTau_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_TrBoth_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
Book( "ChargedHiggsTauNu/Base/Tau1Pt_Met150_TrBoth_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
//
Book( "ChargedHiggsTauNu/NOne/NTaus_"+l,"NTaus "+l +";Number of selected hadronic tau" ,1000,0,1000);
Book( "ChargedHiggsTauNu/NOne/Tau1Pt_"+l,"Tau1Pt "+l +";p_{T}^{#tau}(leading #tau) [GeV]" ,1000,0,1000);
//Book( "ChargedHiggsTauNu/NOne/Tau1Eta_"+l,"Tau1Eta "+l +";#eta^{#tau}(leading #tau) [GeV]",100,-5,5);
Book( "ChargedHiggsTauNu/NOne/NJets_"+l,"NJets "+l + "; Number of selected jets",1000,0,1000);
Book( "ChargedHiggsTauNu/NOne/Jet1Pt_"+l,"Jet1Pt "+l+";p_{T}^{jet}(leading jet) [GeV]",1000,0,1000);
Book( "ChargedHiggsTauNu/NOne/Jet1Eta_"+l,"Jet1Eta "+l+";#eta^{jet}(leading jet) [GeV]",100,-5,5);
//Book( "ChargedHiggsTauNu/NOne/Jet2Pt_"+l,"Jet2Pt "+l+";p_{T}^{jet}(sub-leading jet) [GeV]",1000,0,1000);
//Book( "ChargedHiggsTauNu/NOne/Jet2Eta_"+l,"Jet2Eta "+l+";#eta^{jet}(sub-leading jet) [GeV]",100,-5,5);
//Book( "ChargedHiggsTauNu/NOne/Jet3Pt_"+l,"Jet3Pt "+l+";p_{T}^{jet}(subsub-leading jet) [GeV]",1000,0,1000);
//Book( "ChargedHiggsTauNu/NOne/Jet3Eta_"+l,"Jet3Eta "+l + ";#eta^{jet}(subsub-leading jet) [GeV]",100,-5,5);
Book( "ChargedHiggsTauNu/NOne/NBjets_"+l,"NBjets "+l + ";Number of selected b jets",1000,0,1000);
Book( "ChargedHiggsTauNu/NOne/Bdiscr_"+l,"BDiscr "+l + ";Btag Discr",1000,-2,2);
Book( "ChargedHiggsTauNu/NOne/Bjet1Pt_"+l,"Bjet1Pt "+l+"p_{T}^{b-tagged jet} [GeV]",1000,0,1000);
//Book( "ChargedHiggsTauNu/NOne/Bjet1Eta_"+l,"Bjet1Eta "+l+";#eta^{b-tagged jet} [GeV]",100,-5,5);
//Book( "ChargedHiggsTauNu/NOne/MaxDEtaBjetJets_"+l,"MaxDEtaBjetJets "+l+";Max #Delta#eta (b-tagged jet,jets)",100,0,10);
//Book( "ChargedHiggsTauNu/NOne/MaxInvMassBjetJets_"+l,"MaxInvMassBjetJets "+l+";Max invariant mass (b-tagged jet,jets) [GeV]",2000,0,2000);
Book( "ChargedHiggsTauNu/NOne/EtMiss_"+l,"EtMiss "+l+";E_{T}^{miss} [GeV]",1000,0,1000);
//Book2D( "ChargedHiggsTauNu/NOne/EtMissVsMt_"+l,"EtMissVsMt "+l+";m_{T} [GeV];E_{T}^{miss} [GeV]",1000,0,1000,1000,0,1000);
Book2D( "ChargedHiggsTauNu/NOne/DPhiJet1MetVsDPhiTauMet_"+l,"DPhiJetMetVsDPhiTauMet "+l + ";#Delta #phi(#tau,MET);#Delta #phi(Jet1,MET)",50,0,TMath::Pi(),50,0,TMath::Pi());
Book2D( "ChargedHiggsTauNu/NOne/DPhiJet2MetVsDPhiTauMet_"+l,"DPhiJetMetVsDPhiTauMet "+l+";Delta #phi(#tau,MET);#Delta #phi(Jet2,MET)",50,0,TMath::Pi(),50,0,TMath::Pi());
Book2D( "ChargedHiggsTauNu/NOne/DPhiJet3MetVsDPhiTauMet_"+l,"DPhiJetMetVsDPhiTauMet "+l+";#Delta #phi(#tau,MET);#Delta #phi(Jet3,MET)",50,0,TMath::Pi(),50,0,TMath::Pi());
Book( "ChargedHiggsTauNu/NOne/RCollMin_"+l,"RCollMin "+l+";R_{coll}^{min}",100,0,2*TMath::Pi());
Book( "ChargedHiggsTauNu/NOne/RbbMin_"+l,"RbbMin "+l+";R_{bb}^{min}",100,0,2*TMath::Pi());
//Book( "ChargedHiggsTauNu/NOne/RsrMax_"+l,"RsrMax "+l+";R_{sr}^{max}",100,0,2*TMath::Pi());
//Book( "ChargedHiggsTauNu/NOne/DPhiTauJet1_"+l,"DPhiTauJet1 "+l+";#Delta #phi (#tau, leading jet)",50,0,TMath::Pi());
//Book2D( "ChargedHiggsTauNu/NOne/RCollMinVsMt_"+l,"RCollMinVsMt "+l + ";m_{T} [GeV];R_{Coll}^{min}",1000,0,1000,100,0,2*TMath::Pi());
//Book2D( "ChargedHiggsTauNu/NOne/RbbMinVsMt_"+l,"RbbMinVsMt "+l+";m_{T} [GeV];R_{bb}^{min}",1000,0,1000,100,0,2*TMath::Pi());
//Book2D( "ChargedHiggsTauNu/NOne/RsrMaxVsMt_"+l,"RsrMaxVsMt "+l+";m_{T} [GeV];R_{sr}^{Max}",1000,0,1000,100,0,2*TMath::Pi());
//Book2D( "ChargedHiggsTauNu/NOne/DPhiTauJet1VsMt_"+l,"DPhiTauJet1VsMt "+l + ";m_{T} [GeV];#Delta #phi (#tau, leading jet)",1000,0,1000,50,0,TMath::Pi());
/*********************************************
* FULL SELECTION *
*********************************************/
// for the TTJets background
Book( "ChargedHiggsTauNu/Vars/JetInvMass_" + l ,"Jet Invariant mass;m_{jj} [GeV]",1000,0,1000) ;
//Book( "ChargedHiggsTauNu/Vars/Jet13InvMass_" + l ,"Jet Invariant mass;m_{jj} [GeV]",1000,0,1000) ;
//Book( "ChargedHiggsTauNu/Vars/Jet23InvMass_" + l ,"Jet Invariant mass;m_{jj} [GeV]",1000,0,1000) ;
//Book( "ChargedHiggsTauNu/Vars/Jet1QGL_" + l ,"QGL of the leading jet;QGL",300,-1.5,1.5) ;
//Book( "ChargedHiggsTauNu/Vars/Jet2QGL_" + l ,"QGL of the subleading jet;QGL",300,-1.5,1.5) ;
//Book( "ChargedHiggsTauNu/Vars/MtDecoQ_" + l ,"Mt deco;MtQ",1000,0,1000) ;
//Book( "ChargedHiggsTauNu/Vars/MtDecoCosPhi_" + l ,"Mt deco;MtCosPhi",1000,-1,1) ;
Book( "ChargedHiggsTauNu/Vars/Mt_matchTau_"+l,"Mt "+l + ";m_{T} [GeV]",8000,0,8000); // the Vars directory contains the full selection
Book( "ChargedHiggsTauNu/Vars/Mt_matchEle_"+l,"Mt "+l + ";m_{T} [GeV]",8000,0,8000); // the Vars directory contains the full selection
Book( "ChargedHiggsTauNu/Vars/Mt_matchMu_"+l,"Mt "+l + ";m_{T} [GeV]",8000,0,8000); // the Vars directory contains the full selection
Book( "ChargedHiggsTauNu/Vars/Mt_matchJet_"+l,"Mt "+l + ";m_{T} [GeV]",8000,0,8000); // the Vars directory contains the full selection
Book( "ChargedHiggsTauNu/Vars/Mt_matchOther_"+l,"Mt "+l + ";m_{T} [GeV]",8000,0,8000); // the Vars directory contains the full selection
/**********************************************
* MT *
**********************************************/
Book( "ChargedHiggsTauNu/Vars/Mt_"+l,"Mt "+l + ";m_{T} [GeV]",8000,0,8000); // the Vars directory contains the full selection
AddFinalHisto("ChargedHiggsTauNu/Vars/Mt_"+l);
// study categorization
for(int i=0;i<=1;++i)
{
Book( string("ChargedHiggsTauNu/Vars/Mt_")+Form("cat%d",i)+"_"+l,"Mt "+l + ";m_{T} [GeV]",8000,0,8000); // the Vars directory contains the full selection
}
}
}