//TCanvas *c1 = new TCanvas("cn","cn",50,50,800,600);
#include <vector>
#include <iostream>
#include <sstream>
#include <iomanip>
#include <math.h>
#include <TH1D.h>
#include <TPave.h>
#include <TBox.h>
//#include "tdrstyle.C"
//#include "CMS_lumi.C"
#include "TH1.h"
#include "TH1F.h"

#include <TFile.h>
#include <TNtuple.h>
#include <TH2.h>
#include <TProfile.h>
#include <TCanvas.h>
#include <TFrame.h>
#include <TROOT.h>
#include <TSystem.h>
#include <TRandom3.h>
#include <TBenchmark.h>
#include <TInterpreter.h>

double calculateLumi(std::string fname, double xsec);


void plot(string obs,TString title,TString dir, TString yaxis, TString channel, int r){

  bool doBkgdScale = false; //true; //false;       //WZ is scaled if true
  bool doFitting = false; //true; //false;
  //  bool doFitting = false; //true; //false; //true; //false;

  //  setTDRStyle();

  int W = 800;
  int H = 600;
  int H_ref = 600;
  int W_ref = 800;

  float T = 0.08*H_ref;
  float B = 0.12*H_ref;
  float L = 0.12*W_ref;
  float R = 0.04*W_ref;

  TCanvas *c1 = new TCanvas("cn","cn",50,50,W,H);
  c1->SetFillColor(0);
  c1->SetBorderMode(0);
  c1->SetFrameFillStyle(0);
  c1->SetFrameBorderMode(0);
  c1->SetLeftMargin( L/W );
  c1->SetRightMargin( R/W );
  c1->SetTopMargin( T/H );
  c1->SetBottomMargin( B/H );
  c1->SetTickx(0);
  c1->SetTicky(0);  

 
  TFile *deb        = new TFile("Data13TeV_DoubleEG_2016B.root");
  TFile *dec        = new TFile("Data13TeV_DoubleEG_2016C.root");
  TFile *ded        = new TFile("Data13TeV_DoubleEG_2016D.root");
  TFile *dee        = new TFile("Data13TeV_DoubleEG_2016E.root");
  TFile *def        = new TFile("Data13TeV_DoubleEG_2016F.root");
  TFile *deg        = new TFile("Data13TeV_DoubleEG_2016G.root");
  TFile *deh2        = new TFile("Data13TeV_DoubleEG_2016Hv2.root");
  TFile *deh3        = new TFile("Data13TeV_DoubleEG_2016Hv3.root");

  TFile *dmub        = new TFile("Data13TeV_DoubleMuon_2016B_v1.root");
  TFile *dmuc        = new TFile("Data13TeV_DoubleMuon_2016C.root");
  TFile *dmud        = new TFile("Data13TeV_DoubleMuon_2016D.root");
  TFile *dmue        = new TFile("Data13TeV_DoubleMuon_2016E.root");
  TFile *dmuf        = new TFile("Data13TeV_DoubleMuon_2016F.root");
  TFile *dmug        = new TFile("Data13TeV_DoubleMuon_2016G.root");
  TFile *dmuh2        = new TFile("Data13TeV_DoubleMuon_2016Hv2.root");
  TFile *dmuh3        = new TFile("Data13TeV_DoubleMuon_2016Hv3.root");

  TFile *dmuegb        = new TFile("Data13TeV_MuonEG_2016B.root");
  TFile *dmuegc        = new TFile("Data13TeV_MuonEG_2016C.root");
  TFile *dmuegd        = new TFile("Data13TeV_MuonEG_2016D.root");
  TFile *dmuege        = new TFile("Data13TeV_MuonEG_2016E.root");
  TFile *dmuegf        = new TFile("Data13TeV_MuonEG_2016F.root");
  TFile *dmuegg        = new TFile("Data13TeV_MuonEG_2016G.root");
  TFile *dmuegh2        = new TFile("Data13TeV_MuonEG_2016Hv2.root");
  TFile *dmuegh3        = new TFile("Data13TeV_MuonEG_2016Hv3.root");

  // MC 
  TFile *tzq       = new TFile("MC13TeV_tZq_ll.root");	     
  TFile *ttbar     = new TFile("MC13TeV_TTJets.root");  	     
//I cannot find this
//TFile *ttw       = new TFile("MC13TeV_TTWToLNu.root");	     
  TFile *ttz       = new TFile("MC13TeV_TTZToLLNuNu.root");     
  TFile *zz        = new TFile("MC13TeV_ZZTo4L.root");	     
  TFile *zjets     = new TFile("MC13TeV_DY50toInf_nlo.root");   
  TFile *wzjets    = new TFile("MC13TeV_WZTo3LNu.root");	     
  TFile *ww        = new TFile("MC13TeV_WWTo2L2Nu.root");
  TFile *dy        = new TFile("MC13TeV_DY10to50.root");

  std::string f_tzq       = "MC13TeV_tZq_ll.root";	     
  std::string f_ttbar     = "MC13TeV_TTJets.root";  	     
//  std::string f_ttw       = "MC13TeV_TTWToLNu.root";	     
  std::string f_ttz       = "MC13TeV_TTZToLLNuNu.root";     
  std::string f_zz        = "MC13TeV_ZZTo4L.root";	     
  std::string f_zjets     = "MC13TeV_DY50toInf_nlo.root";   
  std::string f_wzjets    = "MC13TeV_WZTo3LNu.root";	     
  std::string f_ww        = "MC13TeV_WWTo2L2Nu.root";
  std::string f_dy        = "MC13TeV_DY10to50.root";


  //  gStyle->SetOptStat(0);
  /*
  TFile *seb        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleElectron_2016B.root");
  TFile *sec        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleElectron_2016C.root");
  TFile *sed        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleElectron_2016D.root");
  TFile *see        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleElectron_2016E.root");
  TFile *sef        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleElectron_2016F.root");
  TFile *seg        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleElectron_2016G.root");
  TFile *seh2        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleElectron_2016Hv2.root");
  TFile *seh3        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleElectron_2016Hv3.root");

  TFile *smub        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleMuon_2016B.root");
  TFile *smuc        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleMuon_2016C.root");
  TFile *smud        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleMuon_2016D.root");
  TFile *smue        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleMuon_2016E.root");
  TFile *smuf        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleMuon_2016F.root");
  TFile *smug        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleMuon_2016G.root");
  TFile *smuh2        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleMuon_2016Hv2.root");
  TFile *smuh3        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_SingleMuon_2016Hv3.root");
 
  TFile *deb        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleEG_2016B.root");
  TFile *dec        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleEG_2016C.root");
  TFile *ded        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleEG_2016D.root");
  TFile *dee        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleEG_2016E.root");
  TFile *def        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleEG_2016F.root");
  TFile *deg        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleEG_2016G.root");
  TFile *deh2        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleEG_2016Hv2.root");
  TFile *deh3        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleEG_2016Hv3.root");

  TFile *dmub        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleMuon_2016B_v1.root");
  TFile *dmuc        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleMuon_2016C.root");
  TFile *dmud        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleMuon_2016D.root");
  TFile *dmue        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleMuon_2016E.root");
  TFile *dmuf        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleMuon_2016F.root");
  TFile *dmug        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleMuon_2016G.root");
  TFile *dmuh2        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleMuon_2016Hv2.root");
  TFile *dmuh3        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_DoubleMuon_2016Hv3.root");

  TFile *dmuegb        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_MuonEG_2016B.root");
  TFile *dmuegc        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_MuonEG_2016C.root");
  TFile *dmuegd        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_MuonEG_2016D.root");
  TFile *dmuege        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_MuonEG_2016E.root");
  TFile *dmuegf        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_MuonEG_2016F.root");
  TFile *dmuegg        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_MuonEG_2016G.root");
  TFile *dmuegh2        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_MuonEG_2016Hv2.root");
  TFile *dmuegh3        = new TFile("../data/phy_2017/Files_Sep24/Data13TeV_MuonEG_2016Hv3.root");

  // MC 
  TFile *tzq       = new TFile("../data/phy_2017/Files_Sep24/MC13TeV_tZq_ll.root");	     
  TFile *ttbar     = new TFile("../data/phy_2017/Files_Sep24/MC13TeV_TTJets.root");  	     
  TFile *ttw       = new TFile("../data/phy_2017/Files_Sep24/MC13TeV_TTWToLNu.root");	     
  TFile *ttz       = new TFile("../data/phy_2017/Files_Sep24/MC13TeV_TTZToLLNuNu.root");     
  TFile *zz        = new TFile("../data/phy_2017/Files_Sep24/MC13TeV_ZZTo4L.root");	     
  TFile *zjets     = new TFile("../data/phy_2017/Fakes_Sep23/fakes.root");   
  TFile *wzjets    = new TFile("../data/phy_2017/Files_Sep24/MC13TeV_WZTo3LNu.root");	     

  std::string f_tzq       = "../data/phy_2017/Files_Sep24/MC13TeV_tZq_ll.root";	     
  std::string f_ttbar     = "../data/phy_2017/Files_Sep24/MC13TeV_TTJets.root";  	     
  std::string f_ttw       = "../data/phy_2017/Files_Sep24/MC13TeV_TTWToLNu.root";	     
  std::string f_ttz       = "../data/phy_2017/Files_Sep24/MC13TeV_TTZToLLNuNu.root";     
  std::string f_zz        = "../data/phy_2017/Files_Sep24/MC13TeV_ZZTo4L.root";	     
  std::string f_zjets     = "../data/phy_2017/Fakes_Sep23/fakes.root";   
  std::string f_wzjets    = "../data/phy_2017/Files_Sep24/MC13TeV_WZTo3LNu.root";	     
*/

  /*
  TFile *seb        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleElectron_2016B.root");
  TFile *sec        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleElectron_2016C.root");
  TFile *sed        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleElectron_2016D.root");
  TFile *see        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleElectron_2016E.root");
  TFile *sef        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleElectron_2016F.root");
  TFile *seg        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleElectron_2016G.root");
  TFile *seh2        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleElectron_2016Hv2.root");
  TFile *seh3        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleElectron_2016Hv3.root");

  TFile *smub        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleMuon_2016B.root");
  TFile *smuc        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleMuon_2016C.root");
  TFile *smud        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleMuon_2016D.root");
  TFile *smue        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleMuon_2016E.root");
  TFile *smuf        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleMuon_2016F.root");
  TFile *smug        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleMuon_2016G.root");
  TFile *smuh2        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleMuon_2016Hv2.root");
  TFile *smuh3        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_SingleMuon_2016Hv3.root");
 
  TFile *deb        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleEG_2016B.root");
  TFile *dec        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleEG_2016C.root");
  TFile *ded        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleEG_2016D.root");
  TFile *dee        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleEG_2016E.root");
  TFile *def        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleEG_2016F.root");
  TFile *deg        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleEG_2016G.root");
  TFile *deh2        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleEG_2016Hv2.root");
  TFile *deh3        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleEG_2016Hv3.root");

  TFile *dmub        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleMuon_2016B_v1.root");
  TFile *dmuc        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleMuon_2016C.root");
  TFile *dmud        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleMuon_2016D.root");
  TFile *dmue        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleMuon_2016E.root");
  TFile *dmuf        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleMuon_2016F.root");
  TFile *dmug        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleMuon_2016G.root");
  TFile *dmuh2        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleMuon_2016Hv2.root");
  TFile *dmuh3        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_DoubleMuon_2016Hv3.root");

  TFile *dmuegb        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_MuonEG_2016B.root");
  TFile *dmuegc        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_MuonEG_2016C.root");
  TFile *dmuegd        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_MuonEG_2016D.root");
  TFile *dmuege        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_MuonEG_2016E.root");
  TFile *dmuegf        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_MuonEG_2016F.root");
  TFile *dmuegg        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_MuonEG_2016G.root");
  TFile *dmuegh2        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_MuonEG_2016Hv2.root");
  TFile *dmuegh3        = new TFile("../data/phy_2017/Files_Sep23/Data13TeV_MuonEG_2016Hv3.root");

  // MC 
  TFile *tzq       = new TFile("../data/phy_2017/Files_Sep23/MC13TeV_tZq_ll.root");	     
  TFile *ttbar     = new TFile("../data/phy_2017/Files_Sep23/MC13TeV_TTJets.root");  	     
  TFile *ttw       = new TFile("../data/phy_2017/Files_Sep23/MC13TeV_TTWToLNu.root");	     
  TFile *ttz       = new TFile("../data/phy_2017/Files_Sep23/MC13TeV_TTZToLLNuNu.root");     
  TFile *zz        = new TFile("../data/phy_2017/Files_Sep23/MC13TeV_ZZTo4L.root");	     
  TFile *zjets     = new TFile("../data/phy_2017/Fakes_Sep23/fakes.root");   
  TFile *wzjets    = new TFile("../data/phy_2017/Files_Sep23/MC13TeV_WZTo3LNu.root");	     

  std::string f_tzq       = "../data/phy_2017/Files_Sep23/MC13TeV_tZq_ll.root";	     
  std::string f_ttbar     = "../data/phy_2017/Files_Sep23/MC13TeV_TTJets.root";  	     
  std::string f_ttw       = "../data/phy_2017/Files_Sep23/MC13TeV_TTWToLNu.root";	     
  std::string f_ttz       = "../data/phy_2017/Files_Sep23/MC13TeV_TTZToLLNuNu.root";     
  std::string f_zz        = "../data/phy_2017/Files_Sep23/MC13TeV_ZZTo4L.root";	     
  std::string f_zjets     = "../data/phy_2017/Fakes_Sep23/fakes.root";   
  std::string f_wzjets    = "../data/phy_2017/Files_Sep23/MC13TeV_WZTo3LNu.root";	     
*/

  /*
  TFile *seb        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleElectron_2016B.root");
  TFile *sec        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleElectron_2016C.root");
  TFile *sed        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleElectron_2016D.root");
  TFile *see        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleElectron_2016E.root");
  TFile *sef        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleElectron_2016F.root");
  TFile *seg        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleElectron_2016G.root");
  TFile *seh2        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleElectron_2016Hv2.root");
  TFile *seh3        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleElectron_2016Hv3.root");

  TFile *smub        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleMuon_2016B.root");
  TFile *smuc        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleMuon_2016C.root");
  TFile *smud        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleMuon_2016D.root");
  TFile *smue        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleMuon_2016E.root");
  TFile *smuf        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleMuon_2016F.root");
  TFile *smug        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleMuon_2016G.root");
  TFile *smuh2        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleMuon_2016Hv2.root");
  TFile *smuh3        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_SingleMuon_2016Hv3.root");
 
  TFile *deb        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleEG_2016B.root");
  TFile *dec        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleEG_2016C.root");
  TFile *ded        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleEG_2016D.root");
  TFile *dee        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleEG_2016E.root");
  TFile *def        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleEG_2016F.root");
  TFile *deg        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleEG_2016G.root");
  TFile *deh2        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleEG_2016Hv2.root");
  TFile *deh3        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleEG_2016Hv3.root");

  TFile *dmub        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleMuon_2016B_v1.root");
  TFile *dmuc        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleMuon_2016C.root");
  TFile *dmud        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleMuon_2016D.root");
  TFile *dmue        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleMuon_2016E.root");
  TFile *dmuf        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleMuon_2016F.root");
  TFile *dmug        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleMuon_2016G.root");
  TFile *dmuh2        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleMuon_2016Hv2.root");
  TFile *dmuh3        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_DoubleMuon_2016Hv3.root");

  TFile *dmuegb        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_MuonEG_2016B.root");
  TFile *dmuegc        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_MuonEG_2016C.root");
  TFile *dmuegd        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_MuonEG_2016D.root");
  TFile *dmuege        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_MuonEG_2016E.root");
  TFile *dmuegf        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_MuonEG_2016F.root");
  TFile *dmuegg        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_MuonEG_2016G.root");
  TFile *dmuegh2        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_MuonEG_2016Hv2.root");
  TFile *dmuegh3        = new TFile("../data/phy_2017/Files_Sep22/Data13TeV_MuonEG_2016Hv3.root");

  // MC 
  TFile *tzq       = new TFile("../data/phy_2017/Files_Sep22/MC13TeV_tZq_ll.root");	     
  TFile *ttbar     = new TFile("../data/phy_2017/Files_Sep22/MC13TeV_TTJets.root");  	     
  TFile *ttw       = new TFile("../data/phy_2017/Files_Sep22/MC13TeV_TTWToLNu.root");	     
  TFile *ttz       = new TFile("../data/phy_2017/Files_Sep22/MC13TeV_TTZToLLNuNu.root");     
  TFile *zz        = new TFile("../data/phy_2017/Files_Sep22/MC13TeV_ZZTo4L.root");	     
  TFile *zjets     = new TFile("../data/phy_2017/Fakes_Sep22/fakes.root");   
  TFile *wzjets    = new TFile("../data/phy_2017/Files_Sep22/MC13TeV_WZTo3LNu.root");	     

  std::string f_tzq       = "../data/phy_2017/Files_Sep22/MC13TeV_tZq_ll.root";	     
  std::string f_ttbar     = "../data/phy_2017/Files_Sep22/MC13TeV_TTJets.root";  	     
  std::string f_ttw       = "../data/phy_2017/Files_Sep22/MC13TeV_TTWToLNu.root";	     
  std::string f_ttz       = "../data/phy_2017/Files_Sep22/MC13TeV_TTZToLLNuNu.root";     
  std::string f_zz        = "../data/phy_2017/Files_Sep22/MC13TeV_ZZTo4L.root";	     
  std::string f_zjets     = "../data/phy_2017/Fakes_Sep22/fakes.root";   
  std::string f_wzjets    = "../data/phy_2017/Files_Sep22/MC13TeV_WZTo3LNu.root";	     
*/


  /*
  TFile *seb        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleElectron_2016B.root");
  TFile *sec        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleElectron_2016C.root");
  TFile *sed        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleElectron_2016D.root");
  TFile *see        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleElectron_2016E.root");
  TFile *sef        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleElectron_2016F.root");
  TFile *seg        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleElectron_2016G.root");
  TFile *seh2        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleElectron_2016Hv2.root");
  TFile *seh3        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleElectron_2016Hv3.root");

  TFile *smub        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleMuon_2016B.root");
  TFile *smuc        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleMuon_2016C.root");
  TFile *smud        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleMuon_2016D.root");
  TFile *smue        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleMuon_2016E.root");
  TFile *smuf        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleMuon_2016F.root");
  TFile *smug        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleMuon_2016G.root");
  TFile *smuh2        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleMuon_2016Hv2.root");
  TFile *smuh3        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_SingleMuon_2016Hv3.root");
 
  TFile *deb        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleEG_2016B.root");
  TFile *dec        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleEG_2016C.root");
  TFile *ded        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleEG_2016D.root");
  TFile *dee        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleEG_2016E.root");
  TFile *def        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleEG_2016F.root");
  TFile *deg        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleEG_2016G.root");
  TFile *deh2        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleEG_2016Hv2.root");
  TFile *deh3        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleEG_2016Hv3.root");

  TFile *dmub        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleMuon_2016B_v1.root");
  TFile *dmuc        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleMuon_2016C.root");
  TFile *dmud        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleMuon_2016D.root");
  TFile *dmue        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleMuon_2016E.root");
  TFile *dmuf        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleMuon_2016F.root");
  TFile *dmug        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleMuon_2016G.root");
  TFile *dmuh2        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleMuon_2016Hv2.root");
  TFile *dmuh3        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_DoubleMuon_2016Hv3.root");

  TFile *dmuegb        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_MuonEG_2016B.root");
  TFile *dmuegc        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_MuonEG_2016C.root");
  TFile *dmuegd        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_MuonEG_2016D.root");
  TFile *dmuege        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_MuonEG_2016E.root");
  TFile *dmuegf        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_MuonEG_2016F.root");
  TFile *dmuegg        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_MuonEG_2016G.root");
  TFile *dmuegh2        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_MuonEG_2016Hv2.root");
  TFile *dmuegh3        = new TFile("../data/phy_2017/Files_Sep21/Data13TeV_MuonEG_2016Hv3.root");

  // MC 
  TFile *tzq       = new TFile("../data/phy_2017/Files_Sep21/MC13TeV_tZq_ll.root");	     
  TFile *ttbar     = new TFile("../data/phy_2017/Files_Sep21/MC13TeV_TTJets.root");  	     
  TFile *ttw       = new TFile("../data/phy_2017/Files_Sep21/MC13TeV_TTWToLNu.root");	     
  TFile *ttz       = new TFile("../data/phy_2017/Files_Sep21/MC13TeV_TTZToLLNuNu.root");     
  TFile *zz        = new TFile("../data/phy_2017/Files_Sep21/MC13TeV_ZZTo4L.root");	     
  TFile *zjets     = new TFile("../data/phy_2017/Fakes_Aug05/fakes.root");   
  TFile *wzjets    = new TFile("../data/phy_2017/Files_Sep21/MC13TeV_WZTo3LNu.root");	     

  std::string f_tzq       = "../data/phy_2017/Files_Sep21/MC13TeV_tZq_ll.root";	     
  std::string f_ttbar     = "../data/phy_2017/Files_Sep21/MC13TeV_TTJets.root";  	     
  std::string f_ttw       = "../data/phy_2017/Files_Sep21/MC13TeV_TTWToLNu.root";	     
  std::string f_ttz       = "../data/phy_2017/Files_Sep21/MC13TeV_TTZToLLNuNu.root";     
  std::string f_zz        = "../data/phy_2017/Files_Sep21/MC13TeV_ZZTo4L.root";	     
  std::string f_zjets     = "../data/phy_2017/Fakes_Aug05/fakes.root";   
  std::string f_wzjets    = "../data/phy_2017/Files_Sep21/MC13TeV_WZTo3LNu.root";	     

*/


 /*
  TFile *seb        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleElectron_2016B.root");
  TFile *sec        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleElectron_2016C.root");
  TFile *sed        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleElectron_2016D.root");
  TFile *see        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleElectron_2016E.root");
  TFile *sef        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleElectron_2016F.root");
  TFile *seg        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleElectron_2016G.root");
  TFile *seh2        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleElectron_2016Hv2.root");
  TFile *seh3        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleElectron_2016Hv3.root");

  TFile *smub        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleMuon_2016B.root");
  TFile *smuc        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleMuon_2016C.root");
  TFile *smud        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleMuon_2016D.root");
  TFile *smue        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleMuon_2016E.root");
  TFile *smuf        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleMuon_2016F.root");
  TFile *smug        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleMuon_2016G.root");
  TFile *smuh2        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleMuon_2016Hv2.root");
  TFile *smuh3        = new TFile("../data/phy_2017/Files_Jul16/Data13TeV_SingleMuon_2016Hv3.root");
 

  TFile *deb        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleEG_2016B.root");
  TFile *dec        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleEG_2016C.root");
  TFile *ded        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleEG_2016D.root");
  TFile *dee        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleEG_2016E.root");
  TFile *def        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleEG_2016F.root");
  TFile *deg        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleEG_2016G.root");
  TFile *deh2        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleEG_2016Hv2.root");
  TFile *deh3        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleEG_2016Hv3.root");

  TFile *dmub        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleMuon_2016B_v1.root");
  TFile *dmuc        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleMuon_2016C.root");
  TFile *dmud        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleMuon_2016D.root");
  TFile *dmue        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleMuon_2016E.root");
  TFile *dmuf        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleMuon_2016F.root");
  TFile *dmug        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleMuon_2016G.root");
  TFile *dmuh2        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleMuon_2016Hv2.root");
  TFile *dmuh3        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_DoubleMuon_2016Hv3.root");

  TFile *dmuegb        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_MuonEG_2016B.root");
  TFile *dmuegc        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_MuonEG_2016C.root");
  TFile *dmuegd        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_MuonEG_2016D.root");
  TFile *dmuege        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_MuonEG_2016E.root");
  TFile *dmuegf        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_MuonEG_2016F.root");
  TFile *dmuegg        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_MuonEG_2016G.root");
  TFile *dmuegh2        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_MuonEG_2016Hv2.root");
  TFile *dmuegh3        = new TFile("../data/phy_2017/Files_Aug05/Data13TeV_MuonEG_2016Hv3.root");

  // MC 
  TFile *tzq       = new TFile("../data/phy_2017/Files_Aug05/MC13TeV_tZq_ll.root");	     
  TFile *ttbar     = new TFile("../data/phy_2017/Files_Aug05/MC13TeV_TTJets.root");  	     
  TFile *ttw       = new TFile("../data/phy_2017/Files_Aug05/MC13TeV_TTWToLNu.root");	     
  TFile *ttz       = new TFile("../data/phy_2017/Files_Aug05/MC13TeV_TTZToLLNuNu.root");     
  TFile *zz        = new TFile("../data/phy_2017/Files_Aug05/MC13TeV_ZZTo4L.root");	     
  TFile *zjets     = new TFile("../data/phy_2017/Fakes_Aug05/fakes.root");   
  TFile *wzjets    = new TFile("../data/phy_2017/Files_Aug05/MC13TeV_WZTo3LNu.root");	     

  std::string f_tzq       = "../data/phy_2017/Files_Aug05/MC13TeV_tZq_ll.root";	     
  std::string f_ttbar     = "../data/phy_2017/Files_Aug05/MC13TeV_TTJets.root";  	     
  std::string f_ttw       = "../data/phy_2017/Files_Aug05/MC13TeV_TTWToLNu.root";	     
  std::string f_ttz       = "../data/phy_2017/Files_Aug05/MC13TeV_TTZToLLNuNu.root";     
  std::string f_zz        = "../data/phy_2017/Files_Aug05/MC13TeV_ZZTo4L.root";	     
  std::string f_zjets     = "../data/phy_2017/Fakes_Aug05/fakes.root";   
  std::string f_wzjets    = "../data/phy_2017/Files_Aug05/MC13TeV_WZTo3LNu.root";	     
*/

  TH1D *DEB     = (TH1D*)deb->Get(dir+obs);
  TH1D *DEC     = (TH1D*)dec->Get(dir+obs);
  TH1D *DED     = (TH1D*)ded->Get(dir+obs);
  TH1D *DEE     = (TH1D*)dee->Get(dir+obs);
  TH1D *DEF     = (TH1D*)def->Get(dir+obs);
  TH1D *DEG     = (TH1D*)deg->Get(dir+obs);
  TH1D *DEH2    = (TH1D*)deh2->Get(dir+obs);
  TH1D *DEH3    = (TH1D*)deh3->Get(dir+obs);

  TH1D *DMuB     = (TH1D*)dmub->Get(dir+obs);
  TH1D *DMuC     = (TH1D*)dmuc->Get(dir+obs);
  TH1D *DMuD     = (TH1D*)dmud->Get(dir+obs);
  TH1D *DMuE     = (TH1D*)dmue->Get(dir+obs);
  TH1D *DMuF     = (TH1D*)dmuf->Get(dir+obs);
  TH1D *DMuG     = (TH1D*)dmug->Get(dir+obs);
  TH1D *DMuH2    = (TH1D*)dmuh2->Get(dir+obs);
  TH1D *DMuH3    = (TH1D*)dmuh3->Get(dir+obs);

  TH1D *DMuEGB     = (TH1D*)dmuegb->Get(dir+obs);
  TH1D *DMuEGC     = (TH1D*)dmuegc->Get(dir+obs);
  TH1D *DMuEGD     = (TH1D*)dmuegd->Get(dir+obs);
  TH1D *DMuEGE     = (TH1D*)dmuege->Get(dir+obs);
  TH1D *DMuEGF     = (TH1D*)dmuegf->Get(dir+obs);
  TH1D *DMuEGG     = (TH1D*)dmuegg->Get(dir+obs);
  TH1D *DMuEGH2    = (TH1D*)dmuegh2->Get(dir+obs);
  TH1D *DMuEGH3    = (TH1D*)dmuegh3->Get(dir+obs);

  TH1D *tZq    = (TH1D*)tzq->Get(dir+obs);
  TH1D *TTbar  = (TH1D*)ttbar->Get(dir+obs);
  TH1D *TTZ    = (TH1D*)ttz->Get(dir+obs);
//  TH1D *TTW    = (TH1D*)ttw->Get(dir+obs);
  TH1D *ZZ     = (TH1D*)zz->Get(dir+obs);
  TH1D *WZJets = (TH1D*)wzjets->Get(dir+obs);
  TH1D *ZJets  = (TH1D*)zjets->Get(dir+obs);
  TH1D *DY     = (TH1D*)dy->Get(dir+obs);
  TH1D *WW     = (TH1D*)ww->Get(dir+obs);


  double global_lumi = 35800; //2268;//166; //45; //2.942;//2.942; //1;//19700.;  //pb-1
  double data_lumi   = 35800;//2268;//166; //45; //2.942;//2.942;   //1;//19700.;  //pb-1

  double tzq_xsec    = 0.09418; //0.0758;   //pb-1
  double ttz_xsec    = 0.2529;   //pb-1
//  double ttw_xsec    = 0.2043;   //pb-1
  double zz_xsec     = 1.212 ; //0.0374; //1.212 ; //0.0374; //1.256; //0.0374; //1.212  ;    //pb-1
  double ttbar_xsec  = 831.76; //831.76*(0.1086*3)*(0.1086*3);//831.76; //831.76*(0.1086*3)*(0.1086*3);  //pb-1
  double wz_xsec     = 4.42965;     //pb-1
  double wzjets_xsec = 4.42965;     //5.26; // 5.26;     //pb-1
  double zjets_xsec  = 1921.8*3.; //6025.2; //1921.4; //6025.2;   //pb-1
  double ww_xsec     = 1.36 ; //12.178;
  double dy_xsec     = 18610. ; //pb-1

  double tzq_lumi    = calculateLumi(f_tzq,    tzq_xsec);
  double ttz_lumi    = calculateLumi(f_ttz,    ttz_xsec);
//  double ttw_lumi    = calculateLumi(f_ttw,    ttw_xsec);
  double zz_lumi     = calculateLumi(f_zz,     zz_xsec);
  double ttbar_lumi  = calculateLumi(f_ttbar,  ttbar_xsec);
  double wzjets_lumi = calculateLumi(f_wzjets, wzjets_xsec);
  double zjets_lumi  = calculateLumi(f_zjets,  zjets_xsec);
  double ww_lumi     = calculateLumi(f_ww,     ww_xsec);
  double dy_lumi     = calculateLumi(f_dy,     dy_xsec);

  tZq->Scale(global_lumi/tzq_lumi);
//  TTW->Scale(global_lumi/ttw_lumi);
  TTZ->Scale(global_lumi/ttz_lumi);
  ZZ->Scale(global_lumi/zz_lumi);
  TTbar->Scale(global_lumi/ttbar_lumi);
  WZJets->Scale(global_lumi/wzjets_lumi);
  ZJets->Scale(global_lumi/zjets_lumi);
  DY->Scale(global_lumi/dy_lumi);
  WW->Scale(global_lumi/ww_lumi);
  //  ZJets->Scale(0.041531); // uue
  //  ZJets->Scale(0.010147); // eeu
  //  ZJets->Scale(0.043770); // eee
  //  ZJets->Scale(0.044304); // eee
  //  ZJets->Scale(0.0062927); // uuu
  
  //  std::cout<< " Channel : " << channel << " hisot name :" << obs << std::endl;
  //  ZJets->Scale(0.067771 * 0.013749 * 0.006801 * 0.075466); // eee
      
  //  if(channel=="uue")  ZJets->Scale(0.067771); // uue
  //  if(channel=="eeu")  ZJets->Scale(0.013749); // eeu
  //  if(channel=="eee")  ZJets->Scale(0.075466); // eee
  //  if(channel=="uuu")  ZJets->Scale(0.006801); // uuu
  
  /*
  TH1D *ZJets1 = (TH1D*)ZJets->Clone();
  double bin1 =  ZJets->GetBinContent(1);
  ZJets1->Sumw2();
  double bin2 =  ZJets->GetBinContent(2);
  double bin3 =  ZJets->GetBinContent(3);
  double bin4 =  ZJets->GetBinContent(4);
  double bin5 =  ZJets->GetBinContent(5);
  bin2=bin2*0.008556;
  bin3=bin3*0.074635;
  bin4=bin4*0.068209;
  bin5=bin5*0.013773;
  //  bin1=bin1*0.008556*0.074635*0.068209*0.008556;
  bin1=bin2+bin3+bin4+bin5;
  ZJets1->SetBinContent(1,bin1);
  ZJets1->SetBinContent(2,bin2);
  ZJets1->SetBinContent(3,bin3);
  ZJets1->SetBinContent(4,bin4);
  ZJets1->SetBinContent(5,bin5);
*/

  tZq->Rebin(r);
//  TTW->Rebin(r);
  TTZ->Rebin(r);
  ZZ->Rebin(r);
  TTbar->Rebin(r);
  WZJets->Rebin(r);
  ZJets->Rebin(r);
  DY->Rebin(r);
  WW->Rebin(r);
  
  ZJets->Add(TTbar);
  ZJets->Add(DY);
  ZJets->Add(WW);

  //--- Adding Data ----
  TH1D *data = (TH1D*)DEB->Clone();
  data->Add(DEC);
  data->Add(DED);
  data->Add(DEE);
  data->Add(DEF);
  data->Add(DEG);
  data->Add(DEH2);
  data->Add(DEH3);
  
  data->Add(DMuB);
  data->Add(DMuC);
  data->Add(DMuD);
  data->Add(DMuE);
  data->Add(DMuF);
  data->Add(DMuG);
  data->Add(DMuH2);
  data->Add(DMuH3);
  
  data->Add(DMuEGB);
  data->Add(DMuEGC);
  data->Add(DMuEGD);
  data->Add(DMuEGE);
  data->Add(DMuEGF);
  data->Add(DMuEGG);
  data->Add(DMuEGH2);
  data->Add(DMuEGH3);
  
  data->Rebin(r);
  data->SetMarkerStyle(20);
  data->SetMarkerSize(1.2);

  tZq->SetLineColor(kBlack);
  tZq->SetFillColor(kRed+2);

  TTZ->SetLineColor(kBlack);
  TTZ->SetFillColor(kYellow-7);

//  TTW->SetLineColor(kBlack);
//  TTW->SetFillColor(kCyan+3);

  ZZ->SetLineColor(kBlack);
  ZZ->SetFillColor(kRed-10);

  WZJets->SetLineColor(kBlack);
  WZJets->SetFillColor(kOrange);

  ZJets->SetLineColor(kBlack);
  ZJets->SetFillColor(kOrange+7);
    
  TPad*    upperPad = new TPad("upperPad", "upperPad",0, 0.3,  1, 1.);
  upperPad->SetTopMargin(0.07);
  upperPad->SetBottomMargin(0);
  upperPad->SetTickx(1);
  upperPad->SetTicky(1);
  upperPad->Draw();
  upperPad->cd();

  THStack *hs = new THStack("hs","");
  hs->Add(ZJets);
  hs->Add(WZJets);
  hs->Add(TTZ);
//  hs->Add(TTW);
  hs->Add(ZZ);
  hs->Add(tZq);


  double ymax = hs->GetMaximum();
  //  cout<<"max y bin "<<ymax<<endl;
  data->GetYaxis()->SetRangeUser(0., ymax* 1.8);

  //  ZJets->Draw("HIST");
  data->Draw("E1P");
  hs->Draw("HIST sames");
  //  hs->Draw("HIST");
  data->Draw("E1P sames");
  //  hs->Draw("SAME hist");
  /*  
  TString step[]={"Total ","(#mu^{+}#mu^{-})#mu^{#pm}","(e^{+}e^{-})e^{#pm}","(#mu^{+}#mu^{-})e^{#pm}","(e^{+}e^{-})#mu^{#pm}"};
  for(int ii=0; ii<sizeof(step)/sizeof(TString); ii++)
    {
      data->GetXaxis()->SetBinLabel(ii+1,step[ii]);
    }
  */

  TH1D *totMC = (TH1D*)tZq->Clone();
  TH1D *totData = (TH1D*)data->Clone();
  totMC->Add(TTZ);
//  totMC->Add(TTW);
  totMC->Add(ZJets);
  totMC->Add(ZZ);
  totMC->Add(WZJets);  
  totMC->SetFillStyle(3006);
  totMC->SetFillColor(kBlack);
  //totMC->Draw("SAME HIST");
  totMC->Draw("same, E2");
  /*
  hs->GetXaxis()->SetTitle(title);
  hs->GetXaxis()->SetLabelFont(42);
  hs->GetXaxis()->SetLabelSize(0.05);
  hs->GetXaxis()->SetTitleSize(0.06);
  hs->GetXaxis()->SetTitleOffset(0.9);
  hs->GetXaxis()->SetTitleFont(42);
  hs->GetXaxis()->SetNdivisions(6,5,0);

  hs->GetYaxis()->SetTitle(yaxis);
  hs->GetYaxis()->SetLabelFont(42);
  hs->GetYaxis()->SetLabelOffset(0.004);
  hs->GetYaxis()->SetLabelSize(0.05);
  hs->GetYaxis()->SetTitleSize(0.06);
  hs->GetYaxis()->SetTitleOffset(0.9);
  hs->GetYaxis()->SetTitleFont(42);
  hs->GetYaxis()->SetNdivisions(6,5,0);
  */

      
 data->GetXaxis()->SetTitle(title);
 // data->GetXaxis()->SetLabelFont(42);
 // data->GetXaxis()->SetLabelSize(0.05);
 // data->GetXaxis()->SetTitleSize(0.06);
 // data->GetXaxis()->SetTitleOffset(0.9);
 // data->GetXaxis()->SetTitleFont(42);
 data->GetXaxis()->SetNdivisions(6,5,0);

 data->GetYaxis()->SetTitle(yaxis);
 // data->GetYaxis()->SetLabelFont(42);
 // data->GetYaxis()->SetLabelOffset(0.004);
 // data->GetYaxis()->SetLabelSize(0.05);
 // data->GetYaxis()->SetTitleSize(0.07);
 // data->GetYaxis()->SetTitleOffset(0.6);
 data->GetYaxis()->SetTitleOffset(1);
 // data->GetYaxis()->SetTitleFont(42);
 data->GetYaxis()->SetNdivisions(6,5,0);
  
 c1->SetTickx(1);
 c1->SetTicky(1);

 TLegend *leg = new TLegend( 0.6449843,0.65038,0.8377743,0.8838219,NULL,"brNDC");
 leg->SetFillColor(0);
 leg->SetBorderSize(0);
 leg->SetTextFont(42);
 leg->SetTextSize(0.04);
 leg->SetLineColor(1);
 leg->SetLineStyle(1);
 leg->SetLineWidth(1);
 leg->SetMargin(0.55);
 leg->SetNColumns(2);
 leg->SetColumnSeparation(.6);

 leg->AddEntry(data,"Data","lpe");//"PL");//"lpe");
 leg->AddEntry(tZq,"tZq","f");//"lpe");
 leg->AddEntry(ZZ,"ZZ","F");
// leg->AddEntry(TTW,"t#bar{t}W","F");
 leg->AddEntry(TTZ,"t#bar{t}Z","F");
 leg->AddEntry(ZJets,"Fakes","F");
 leg->AddEntry(WZJets,"WZJets","F");
 leg->AddEntry(totMC,"Stat.","F");
 leg->Draw();

 /*
 TLatex *   tex1 = new TLatex(0.86,0.936,"#sqrt[]{s}=13TeV, L_{int}= 36 fb^{-1}");
 tex1->SetNDC();
 tex1->SetTextAlign(31);
 tex1->SetTextFont(42);
 tex1->SetTextSize(0.042);
 tex1->SetLineWidth(2);
 tex1->Draw();

 TLatex * tex2 = new TLatex(0.12,0.936,channel);
 tex2->SetNDC();
 tex2->SetTextAlign(61);
 tex2->SetTextFont(42);
 tex2->SetTextSize(0.058);
 tex2->SetLineWidth(2);
 tex2->Draw();

 TLatex * tex3 = new TLatex(0.25,0.936," CMS Preliminary");
 tex3->SetNDC();
 tex3->SetTextAlign(10);
 tex3->SetTextFont(52);
 tex3->SetTextSize(0.05);
 tex3->SetLineWidth(2);
 tex3->Draw();
 */
 c1->cd();
 
 TPad*    lowerPad = new TPad("lowerPad", "lowerPad",0., 0.05, 1, 0.30);
 lowerPad->SetTopMargin(0);
 lowerPad->SetBottomMargin(0.3);
 lowerPad->SetGridy();
 lowerPad->SetTickx(1);
 lowerPad->SetTicky(1);
 lowerPad->Draw();
 lowerPad->cd();

 // totData->Sumw2();
 totData->Divide(totMC);
 totData->GetYaxis()->SetTitle("Data/MC");
 totData->GetYaxis()->CenterTitle(true);
 totData->GetXaxis()->SetTitle("");

 totData->GetXaxis()->SetTitle(title);
 // totData->GetXaxis()->SetLabelFont(42);
 // totData->GetXaxis()->SetLabelSize(0.12);
 // totData->GetXaxis()->SetTitleSize(0.2);
 // totData->GetXaxis()->SetTitleOffset(0.55);
 totData->GetXaxis()->SetTitleOffset(1);
 // totData->GetXaxis()->SetTitleFont(42);
 totData->GetXaxis()->SetNdivisions(6,5,0);

 // totData->GetYaxis()->SetLabelSize(0.1);
 // totData->GetYaxis()->SetTitleSize(0.15);
 // totData->GetYaxis()->SetTitleOffset(0.27);
 totData->GetYaxis()->SetTitleOffset(1);
 totData->GetYaxis()->SetNdivisions(6,5,0);
 totData->GetYaxis()->SetRangeUser(0.5,1.5);
 totData->SetStats(0);
 gStyle->SetOptStat(0);
 totData->Draw("E1P");
 totData->SetTitle("");
 c1->cd();
 
 // string s1 = "Results/2017/Files_Sep24/Unscaled_"+obs+".png";
 string s1 = obs+".png";

 c1->SaveAs(s1.c_str());
 
 deb     -> Close();
 dec     -> Close();
 ded     -> Close();
 dee     -> Close();
 def     -> Close();
 deg     -> Close();
 deh2     -> Close();
 deh3     -> Close();

 dmub     -> Close();
 dmuc     -> Close();
 dmud     -> Close();
 dmue     -> Close();
 dmuf     -> Close();
 dmug     -> Close();
 dmuh2     -> Close();
 dmuh3     -> Close();

 dmuegb     -> Close();
 dmuegc     -> Close();
 dmuegd     -> Close();
 dmuege     -> Close();
 dmuegf     -> Close();
 dmuegg     -> Close();
 dmuegh2     -> Close();
 dmuegh3     -> Close();

 tzq     -> Close();
 ttbar   -> Close();
// ttw     -> Close();
 ttz     -> Close();
 wzjets  -> Close();
 zz      -> Close();
 zjets   -> Close();
 ww   -> Close();
 dy   -> Close();

  }


void Nabila(){

//  plot("Z_Inv_mass","Muon Multi","", "Entries / GeV", "eee",10);
  std::cout<< " ===== : " << std::endl;
//    plot("w_tm_uuu_cr1","W_{tMass}","", "Entries / GeV", "uuu",10);
  //  plot("w_tm_eee_cr1","W_{tMass}","", "Entries / GeV", "eee",10);
  /*  
  plot("muon_multi_eee_cr1","Muon Multi","", "Entries / GeV", "eee",1);
  plot("muon_multi_uuu_cr1","Muon Multi","", "Entries / GeV", "uuu",1);
  plot("muon_multi_eeu_cr1","Muon Multi","", "Entries / GeV", "eeu",1);
  plot("muon_multi_uue_cr1","Muon Multi","", "Entries / GeV", "uue",1);

  plot("elec_multi_eee_cr1","Electron Multi","", "Entries / GeV", "eee",1);
  plot("elec_multi_uuu_cr1","Electron Multi","", "Entries / GeV", "uuu",1);
  plot("elec_multi_eeu_cr1","Electron Multi","", "Entries / GeV", "eeu",1);
  plot("elec_multi_uue_cr1","Electron Multi","", "Entries / GeV", "uue",1);

  plot("lmuon_pt_uuu_cr1","Muon p_{T}","", "Entries / GeV", "uuu",10);
  plot("lmuon_pt_uue_cr1","Muon p_{T}","", "Entries / GeV", "uue",10);
  plot("lmuon_pt_eeu_cr1","Muon p_{T}","", "Entries / GeV", "eeu",10);
  plot("lmuon_pt_eee_cr1","Muon p_{T}","", "Entries / GeV", "eee",10);
  
  plot("lelectron_pt_uuu_cr1","Electron p_{T}","", "Entries / GeV", "uuu",10);
  plot("lelectron_pt_uue_cr1","Electron p_{T}","", "Entries / GeV", "uue",10);
  plot("lelectron_pt_eeu_cr1","Electron p_{T}","", "Entries / GeV", "eeu",10);
  plot("lelectron_pt_eee_cr1","Electron p_{T}","", "Entries / GeV", "eee",10);

  plot("lmuon_eta_uuu_cr1","Muon #eta","", "Entries / GeV", "uuu",10);
  plot("lmuon_eta_uue_cr1","Muon #eta","", "Entries / GeV", "uue",10);
  plot("lmuon_eta_eeu_cr1","Muon #eta","", "Entries / GeV", "eeu",10);
  plot("lmuon_eta_eee_cr1","Muon #eta","", "Entries / GeV", "eee",10);
  
  plot("lelectron_eta_uuu_cr1","Electron #eta","", "Entries / GeV", "uuu",10);
  plot("lelectron_eta_uue_cr1","Electron #eta","", "Entries / GeV", "uue",10);
  plot("lelectron_eta_eeu_cr1","Electron #eta","", "Entries / GeV", "eeu",10);
  plot("lelectron_eta_eee_cr1","Electron #eta","", "Entries / GeV", "eee",10);

  plot("lbjet_eta_uuu_cr1","bJet #eta","", "Entries / GeV", "uuu",10);
  plot("lbjet_eta_uue_cr1","bJet #eta","", "Entries / GeV", "uue",10);
  plot("lbjet_eta_eeu_cr1","bJet #eta","", "Entries / GeV", "eeu",10);
  plot("lbjet_eta_eee_cr1","bJet #eta","", "Entries / GeV", "eee",10);

  plot("lbjet_pt_uuu_cr1","bJet p_{T}","", "Entries / GeV", "uuu",10);
  plot("lbjet_pt_uue_cr1","bJet p_{T}","", "Entries / GeV", "uue",10);
  plot("lbjet_pt_eeu_cr1","bJet p_{T}","", "Entries / GeV", "eeu",10);
  plot("lbjet_pt_eee_cr1","bJet p_{T}","", "Entries / GeV", "eee",10);

  //  plot("tzq_lightjets_multi_","light jet multi","", "Entries / GeV", "",1);
  //  plot("tzq_bjets_multi","bjet multi","", "Entries / GeV", "",1);
  //  plot("tzq_totjets_multi_","jet multi","", "Entries / GeV", "",1);

  //  plot("ttz_lightjets_multi_","light jet multi","", "Entries / GeV", "",1);
  //  plot("ttz_bjets_multi_","bjet multi","", "Entries / GeV", "",1);
  //  plot("ttz_totjets_multi_","jet multi","", "Entries / GeV", "",1);

  
  plot("lmuon_pt_uuu_cr1","Muon p_{T}","", "Entries / GeV", "uuu",10);
  plot("lmuon_pt_uue_cr1","Muon p_{T}","", "Entries / GeV", "uue",10);
  plot("lmuon_pt_eeu_cr1","Muon p_{T}","", "Entries / GeV", "eeu",10);
  plot("lmuon_pt_eee_cr1","Muon p_{T}","", "Entries / GeV", "eee",10);
  
  plot("lelectron_pt_uuu_cr1","Electron p_{T}","", "Entries / GeV", "uuu",10);
  plot("lelectron_pt_uue_cr1","Electron p_{T}","", "Entries / GeV", "uue",10);
  plot("lelectron_pt_eeu_cr1","Electron p_{T}","", "Entries / GeV", "eeu",10);
  plot("lelectron_pt_eee_cr1","Electron p_{T}","", "Entries / GeV", "eee",10);

  plot("lmuon_eta_uuu_cr1","Muon #eta","", "Entries / GeV", "uuu",10);
  plot("lmuon_eta_uue_cr1","Muon #eta","", "Entries / GeV", "uue",10);
  plot("lmuon_eta_eeu_cr1","Muon #eta","", "Entries / GeV", "eeu",10);
  plot("lmuon_eta_eee_cr1","Muon #eta","", "Entries / GeV", "eee",10);
  
  plot("lelectron_eta_uuu_cr1","Electron #eta","", "Entries / GeV", "uuu",10);
  plot("lelectron_eta_uue_cr1","Electron #eta","", "Entries / GeV", "uue",10);
  plot("lelectron_eta_eeu_cr1","Electron #eta","", "Entries / GeV", "eeu",10);
  plot("lelectron_eta_eee_cr1","Electron #eta","", "Entries / GeV", "eee",10);

  plot("z_mass_eee_cr1","Z_{e^{+}e^{-}}","", "Entries / GeV", "eee",2);
  plot("z_mass_uuu_cr1","Z_{#mu^{+}#mu^{-}}","", "Entries / GeV", "uuu",2);
  plot("z_mass_eeu_cr1","Z_{e^{+}e^{-}}","", "Entries / GeV", "eeu",2);
  plot("z_mass_uue_cr1","Z_{#mu^{+}#mu^{-}}","", "Entries / GeV", "uue",2);
  
  plot("w_tm_eee_cr1","m_{T}^{W}[GeV]","", "Entries / GeV", "eee",10);
  plot("w_tm_uuu_cr1","m_{T}^{W}[GeV]","", "Entries / GeV", "uuu",10);
  plot("w_tm_eeu_cr1","m_{T}^{W}[GeV]","", "Entries / GeV", "eeu",10);
  plot("w_tm_uue_cr1","m_{T}^{W}[GeV]","", "Entries / GeV", "uue",10);
  
  plot("met_uuu_cr1","E_{T}^{miss}","", "Entries / GeV", "uuu",10);
  plot("met_uue_cr1","E_{T}^{miss}","", "Entries / GeV", "uue",10);
  plot("met_eeu_cr1","E_{T}^{miss}","", "Entries / GeV", "eeu",10);
  plot("met_eee_cr1","E_{T}^{miss}","", "Entries / GeV", "eee",10);

  plot("ljet_eta_uuu_cr1","Jet #eta","", "Entries / GeV", "uuu",10);
  plot("ljet_eta_uue_cr1","Jet #eta","", "Entries / GeV", "uue",10);
  plot("ljet_eta_eeu_cr1","Jet #eta","", "Entries / GeV", "eeu",10);
  plot("ljet_eta_eee_cr1","Jet #eta","", "Entries / GeV", "eee",10);

  plot("jets_multi_uuu_cr1","Jet #eta","", "Entries / GeV", "uuu",1);
  plot("jets_multi_uue_cr1","Jet #eta","", "Entries / GeV", "uue",1);
  plot("jets_multi_eeu_cr1","Jet #eta","", "Entries / GeV", "eeu",1);
  plot("jets_multi_eee_cr1","Jet #eta","", "Entries / GeV", "eee",1);


  plot("ljet_pt_uuu_cr1","Jet p_{T}","", "Entries / GeV", "uuu",10);
  plot("ljet_pt_uue_cr1","Jet p_{T}","", "Entries / GeV", "uue",10);
  plot("ljet_pt_eeu_cr1","Jet p_{T}","", "Entries / GeV", "eeu",10);
  plot("ljet_pt_eee_cr1","Jet p_{T}","", "Entries / GeV", "eee",10);

  plot("dphi_ElMet_uue_cr1","#Delta#Phi W_{(e,MET)}","", "Entries / GeV", "uue",4);
  plot("dphi_ElMet_eee_cr1","#Delta#Phi W_{(e,MET)}","", "Entries / GeV", "eee",4);

  plot("dphi_MuMet_uuu_cr1","#Delta#Phi W_{(#mu,MET)}","", "Entries / GeV", "uuu",4);
  plot("dphi_MuMet_eeu_cr1","#Delta#Phi W_{(#mu,MET)}","", "Entries / GeV", "eeu",4);

  plot("deltaPhi_mu_uuu_cr1","#Delta#Phi (#mu,#mu)","", "Entries / GeV", "uuu",2);
  plot("deltaPhi_mu_uue_cr1","#Delta#Phi (#mu,#mu)","", "Entries / GeV", "uue",2);

  plot("deltaPhi_el_eeu_cr1","#Delta#Phi (#el,#el)","", "Entries / GeV", "eeu",2);
  plot("deltaPhi_el_eee_cr1","#Delta#Phi (#el,#el)","", "Entries / GeV", "eee",2);
  */

  

  //  plot("finalCutflow_ttz","","", "Entries", "",1); 
  //  plot("cr1_cutflow","","", "Entries", "",1); 
  // plot("w_tm_eee_cr1","W_{tMass}","", "Entries / GeV", "eee",10);
  

  ///eee
//  plot("z_mass_eee_cr1","M_ll","", "Entries / GeV", "eee",2);
  plot("w_tm_eee_cr1","W_{tMass}","", "Entries / GeV", "eee",10);
/*  plot("met_eee_cr1","E_{T}^{miss}","", "Entries / GeV", "eee",10);
  
  ///uuu
  plot("z_mass_uuu_cr1","M_ll","", "Entries / GeV", "uuu",2);
  plot("w_tm_uuu_cr1","W_{tMass}","", "Entries / GeV", "uuu",10);
  plot("met_uuu_cr1","E_{T}^{miss}","", "Entries / GeV", "uuu",10);
  ///eeu
  plot("z_mass_eeu_cr1","M_ll","", "Entries / GeV", "eeu",2);
  plot("w_tm_eeu_cr1","W_{tMass}","", "Entries / GeV", "eeu",10);
  plot("met_uue_cr1","E_{T}^{miss}","", "Entries / GeV", "uue",10);
  ///uue
  plot("z_mass_uue_cr1","M_ll","", "Entries / GeV", "uue",2);
  plot("w_tm_uue_cr1","W_{tMass}","", "Entries / GeV", "uue",10);  
  plot("met_eeu_cr1","E_{T}^{miss}","", "Entries / GeV", "eeu",10);
  */
  /*
  
  ///eee
  plot("f_z_mass_eee","Z_{ll}","", "Entries / GeV", "eee",2);
  plot("f_w_tm_eee","m_{T}^{W}[GeV]","", "Entries / GeV", "eee",10);
  plot("f_met_eee","E_{T}^{miss}","", "Entries / GeV", "eee",10);
  plot("f_top_invM_eee","M_{top}","", "Entries / GeV", "eee",25);

  
  ///uuu
  plot("f_z_mass_uuu","M_{ll}","", "Entries / GeV", "uuu",2);
  plot("f_w_tm_uuu","m_{T}^{W}[GeV]","", "Entries / GeV", "uuu",10);
  plot("f_met_uuu","E_{T}^{miss}","", "Entries / GeV", "uuu",10);
  plot("f_top_invM_uuu","M_{top}","", "Entries / GeV", "uuu",10);
  ///eeu
  plot("f_z_mass_eeu","M_{ll}","", "Entries / GeV", "eeu",2);
  plot("f_w_tm_eeu","m_{T}^{W}[GeV]","", "Entries / GeV", "eeu",10);
  plot("f_met_eeu","E_{T}^{miss}","", "Entries / GeV", "eeu",10);
  plot("f_top_invM_eeu","M_{top}","", "Entries / GeV", "eeu",10);
  ///uue
  plot("f_z_mass_uue","M_{ll}","", "Entries / GeV", "uue",2);
  plot("f_w_tm_uue","m_{T}^{W}[GeV]","", "Entries / GeV", "uue",10);  
  plot("f_met_uue","E_{T}^{miss}","", "Entries / GeV", "uue",10);
  plot("f_top_invM_uue","M_{top}","", "Entries / GeV", "uue",10);
*/
 
  /*
  ///eee
  plot("ttz_z_mass_eee","M_ll","", "Entries / GeV", "eee",2);
  plot("ttz_w_tm_eee","W_{tMass}","", "Entries / GeV", "eee",10);
  plot("ttz_met_eee","E_{T}^{miss}","", "Entries / GeV", "eee",10);
  ///uuu
  plot("ttz_z_mass_uuu","M_ll","", "Entries / GeV", "uuu",2);
  plot("ttz_w_tm_uuu","W_{tMass}","", "Entries / GeV", "uuu",10);
  plot("ttz_met_uuu","E_{T}^{miss}","", "Entries / GeV", "uuu",10);
  ///eeu
  plot("ttz_z_mass_eeu","M_ll","", "Entries / GeV", "eeu",2);
  plot("ttz_w_tm_eeu","W_{tMass}","", "Entries / GeV", "eeu",10);
  plot("ttz_met_uue","E_{T}^{miss}","", "Entries / GeV", "uue",10);
  ///uue
  plot("ttz_z_mass_uue","M_ll","", "Entries / GeV", "uue",2);
  plot("ttz_w_tm_uue","W_{tMass}","", "Entries / GeV", "uue",10);  
  plot("ttz_met_eeu","E_{T}^{miss}","", "Entries / GeV", "eeu",10);
  */


  //  plot("w_tm_elnu","W_{e#nu}","", "Entries", "",10); 
  //  plot("w_tm_munu","W_{#mu#nu}","", "Entries", "",10);        
  //  plot("w_tm_munu_new","W_{#mu#nu}","", "Entries", "",10);        
  //  plot("w_tm_elnu_new","W_{#mu#nu}","", "Entries", "",10);        
  //  plot("bjets_multi","BJets Multiplicity","", "Entries", "",1);
  //  plot("w_tm_eee_cr1","W_{#mu#nu}","", "Entries", "",10);        
  //  plot("deltaPhi_el","#Delta #Phi","", "Entries / GeV", "",5);    
  //  plot("z_mass_eee_cr1","M_ll","", "Entries / GeV", "uu",2);

  //  plot("z_pt_uu","p_{T}","", "Entries / GeV", "uu",2);
  //  plot("z_mass_ee","M_ll","", "Entries / GeV", "ee",2);
  //  plot("z_mass_uuu_cr1","M_ll","", "Entries / GeV", "eee",2);
  //  plot("ttz_z_mass_eee","M_ll","", "Entries / GeV", "eee",2);
  //  plot("z_mass_uue_cr1","M_ll","", "Entries / GeV", "eee",2);
  //  plot("met_uue_cr1","M_ll","", "Entries / GeV", "eee",10);
  //  plot("met_eee","E_{T}^{miss}","", "Entries / GeV", "uuu",10);
  //  plot("tel_iso_bc","electron iso","", "Entries / GeV", "",5);    
  /*
  plot("dR_emu_","#Delta R(#mu,e)","", "Entries", "",4);
  plot("dphi_ElMet","","", "Entries", "",5);  
  plot("dphi_MuMet","","", "Entries", "",5);
    

  plot("tmu_iso","muon iso","", "Entries / GeV", "",5);    
    
  plot("lmuon_pt_","","", "Entries", "",5);
  
  plot("Mu_Pt","Muon p_{T}","", "Entries", "",5);
  plot("Mu_Phi","Muon p_{T}","", "Entries", "",5);
  plot("Mu_Eta","Muon p_{T}","", "Entries", "",5);
  plot("El_Pt", "Electroon p_{T}","", "Entries", "",5);
  plot("El_Phi","Electroon p_{T}","", "Entries", "",5);
  plot("El_Eta","Electroon p_{T}","", "Entries", "",5);
  plot("tel_iso_bc","electron iso","", "Entries / GeV", "",5);    
  plot("tmu_iso_bc","muon iso","", "Entries / GeV", "",5);    
  plot("vetoLep_multi_","","", "Entries", "",1);
  
  //  plot("w_tm_munu_new","W_{#mu#nu}","", "Entries", "",10);
  //  plot("w_tm_elnu_new","W_{e#nu}","", "Entries", "",10);         
  //  plot("cr1_cutflow","","", "Entries", "",1);
  //  plot("cr0_cutflow","","", "Entries", "",1);
  //  plot("presel_cutflow","","", "Entries", "",1); 
  
  //  plot("muon_multi_","","", "Entries", "",1);
  //  plot("electron_iso","electron iso","", "Entries / GeV", "",1);
  //  plot("met_","E_{T}^{miss}","", "Entries / GeV", "",10);
  //  plot("elec_multi_","","", "Entries", "",1);
  
   
  //  plot("w_tm_munu","W_{#mu#nu}","", "Entries", "",10);        
  //  plot("w_tm_munu_cr0","W_{#mu#nu}","", "Entries", "",10);        
  //  plot("w_tm_munu_cr1","W_{#mu#nu}","", "Entries", "",10);        
  
  plot("w_tm_elnu","W_{e#nu}","", "Entries", "",10); 
  plot("w_tm_elnu_cr0","W_{e#nu}","", "Entries", "",10); 
  plot("w_tm_elnu_cr1","W_{e#nu}","", "Entries", "",10); 
  
  //  plot("w_tm_eee_cr1","W_{tMass}","", "Entries / GeV", "eee",10);
  //  plot("met_eee_cr1","E_{T}^{miss}","", "Entries / GeV", "uue",5);
  //  plot("w_tm_uuu_cr1","W_{tMass}","", "Entries / GeV", "uuu",10);
  //  plot("w_tm_eeu_cr1","W_{tMass}","", "Entries / GeV", "eeu",10);
  //  plot("w_tm_uue_cr1","W_{tMass}","", "Entries / GeV", "uue",10);
  
  plot("finalCutflow_tzq2","","", "Entries", "",1);
  plot("w_tm_eee_cr1","W_{tMass}","", "Entries / GeV", "eee",10);
  plot("w_tm_uuu_cr1","W_{tMass}","", "Entries / GeV", "uuu",10);
  plot("w_tm_eeu_cr1","W_{tMass}","", "Entries / GeV", "eeu",10);
  plot("w_tm_uue_cr1","W_{tMass}","", "Entries / GeV", "uue",10); 
  
  plot("finalCutflow","","", "Entries", "",1);
  plot("HT","HT","", "Entries", "",20);
  plot("ST","ST","", "Entries", "",20);

  plot("cutflow","","", "Entries", "",1);
  plot("cr0_cutflow","","", "Entries", "",1);
  plot("muon_multi_","","", "Entries", "",1);
  plot("elec_multi_","","", "Entries", "",1);

  plot("top_invM_mu","M_{top}","", "Entries / GeV", "t->b+W_{#mu#nu_{#mu}}",20);
  plot("top_pt_mu","p_{T}_{top}","", "Entries / GeV", "t->b+W_{#mu#nu_{#mu}}",20);
  plot("top_tm_mu","tM_{top}","", "Entries / GeV", "t->b+W_{#mu#nu_{#mu}}",20);
  plot("top_invM_el","M_{top}","", "Entries / GeV", "t->b+W_{e#nu_{e}}",20);
  plot("top_pt_el","p_{T}_{top}","", "Entries / GeV", "t->b+W_{e#nu_{e}}",20);
  plot("top_tm_el","tM_{top}","", "Entries / GeV", "t->b+W_{e#nu_{e}}",10);

  plot("vetoLep_multi_","veto Lepton Multiplicity","", "Entries", "",1);
  plot("finalCutflow","","", "Entries", "",1);
  plot("finalCutflow_ttz","","", "Entries", "",1);
  plot("finalCutflow_tzq1","","", "Entries", "",1);
  
  plot("presel_cutflow","","", "Entries", "",1);
  plot("gen_wgt_","generator weight","", "Entries / GeV", "",1);
  plot("csv_","csv","", "Entries / GeV", "",5);
  
  plot("Nvtx_before","N_{vtx}","", "Entries", "",2);
  plot("Nvtx_after","N_{vtx}","", "Entries", "",2);
  plot("lmuon_pt_","muon p_{T}","", "Entries / GeV", "",10);
  plot("lmuon_eta_","muon #eta","", "Entries / GeV", "",10);
  plot("lelectron_pt_","electron p_{T}","", "Entries / GeV", "",10);
  plot("lelectron_eta_","electron #eta","", "Entries / GeV", "",10);
  plot("electron_iso","electron iso","", "Entries / GeV", "",10);
  plot("muon_iso","muon iso","", "Entries / GeV", "",10);
  
  plot("lepton_multi_bc","Lepton Multiplicity","", "Entries", "",1);
  plot("bjets_multi","BJets Multiplicity","", "Entries", "",1);
  plot("jets_multi","Jets Multiplicity","", "Entries", "",1);
  plot("bjets_multi_bc","BJets Multiplicity","", "Entries", "",1);
  plot("ljets_multi","Jets Multiplicity","", "Entries", "",1);
  
  plot("jet_eta_","bjet #eta","", "Entries / GeV", "",10);
  plot("jet_pt_","jet p_{T}","", "Entries / GeV", "",15);
  
  plot("lightjets_multi_","Jets Multiplicity","", "Entries", "",1);
  plot("totjets_multi_","Jets Multiplicity","", "Entries", "",1);
  plot("bjets_multi","BJets Multiplicity","", "Entries", "",1);
  
  plot("ljet_eta_","jet #eta","", "Entries / GeV", "",10);
  plot("ljet_pt_","jet p_{T}","", "Entries / GeV", "",15);
  
  plot("lbjet_eta_","bjet #eta","", "Entries / GeV", "",10);
  plot("lbjet_pt_","bjet p_{T}","", "Entries / GeV", "",15);
  plot("bjets_multi","bJets Multiplicity","", "Entries", "",1);
    
  plot("met_uuu","E_{T}^{miss}","", "Entries / GeV", "uuu",10);
  plot("met_uue","E_{T}^{miss}","", "Entries / GeV", "uue",10);
  plot("met_eeu","E_{T}^{miss}","", "Entries / GeV", "eeu",10);
  plot("met_eee","E_{T}^{miss}","", "Entries / GeV", "eee",10);

  plot("w_tm_eee","W_{tMass}","", "Entries / GeV", "eee",10); 
  plot("w_tm_uuu","W_{tMass}","", "Entries / GeV", "uuu",10); 
  plot("w_tm_eeu","W_{tMass}","", "Entries / GeV", "eeu",10); 
  plot("w_tm_uue","W_{tMass}","", "Entries / GeV", "uue",10); 


  plot("z_mass_eee","M_ll","", "Entries / GeV", "eee",2);
  plot("z_mass_uuu","M_ll","", "Entries / GeV", "uuu",2);
  plot("z_mass_eeu","M_ll","", "Entries / GeV", "eeu",2);
  plot("z_mass_uue","M_ll","", "Entries / GeV", "uue",2);
  */
  /*
  plot("z_pt_eee","p_{T} of Z","", "Entries / GeV", "eee",10); 
  plot("z_pt_uuu","p_{T} of Z","", "Entries / GeV", "uuu",10); 
  plot("z_pt_eeu","p_{T} of Z","", "Entries / GeV", "eeu",10); 
  plot("z_pt_uue","p_{T} of Z","", "Entries / GeV", "uue",10); 

  plot("z_mass_eee_cr0","M_ll","", "Entries / GeV", "eee",2);
  plot("z_mass_uuu_cr0","M_ll","", "Entries / GeV", "uuu",2);
  plot("z_mass_eeu_cr0","M_ll","", "Entries / GeV", "eeu",2);
  plot("z_mass_uue_cr0","M_ll","", "Entries / GeV", "uue",2);
  
  plot("w_tm_eee_cr0","W_{tMass}","", "Entries / GeV", "eee",10);
  plot("w_tm_uuu_cr0","W_{tMass}","", "Entries / GeV", "uuu",10);
  plot("w_tm_eeu_cr0","W_{tMass}","", "Entries / GeV", "eeu",10);
  plot("w_tm_uue_cr0","W_{tMass}","", "Entries / GeV", "uue",10);

  plot("met_uuu_cr0","E_{T}^{miss}","", "Entries / GeV", "uuu",10);
  plot("met_uue_cr0","E_{T}^{miss}","", "Entries / GeV", "uue",10);
  plot("met_eeu_cr0","E_{T}^{miss}","", "Entries / GeV", "eeu",10);
  plot("met_eee_cr0","E_{T}^{miss}","", "Entries / GeV", "eee",10);
  

  plot("cr0_cutflow","","", "Entries", "",1);
  plot("cr1_cutflow","","", "Entries", "",1);
  
  plot("z_mass_eee_cr1","M_ll","", "Entries / GeV", "eee",2);
  plot("z_mass_uuu_cr1","M_ll","", "Entries / GeV", "uuu",2);
  plot("z_mass_eeu_cr1","M_ll","", "Entries / GeV", "eeu",2);
  plot("z_mass_uue_cr1","M_ll","", "Entries / GeV", "uue",2);
  
  plot("w_tm_eee_cr1","W_{tMass}","", "Entries / GeV", "eee",10);
  plot("w_tm_uuu_cr1","W_{tMass}","", "Entries / GeV", "uuu",10);
  plot("w_tm_eeu_cr1","W_{tMass}","", "Entries / GeV", "eeu",10);
  plot("w_tm_uue_cr1","W_{tMass}","", "Entries / GeV", "uue",10);
  
  plot("met_uuu_cr1","E_{T}^{miss}","", "Entries / GeV", "uuu",10);
  plot("met_uue_cr1","E_{T}^{miss}","", "Entries / GeV", "uue",10);
  plot("met_eeu_cr1","E_{T}^{miss}","", "Entries / GeV", "eeu",10);
  plot("met_eee_cr1","E_{T}^{miss}","", "Entries / GeV", "eee",10);
*/  
  
}


double calculateLumi(std::string fname, double xsec){
  TFile *file = new TFile ( fname.c_str() );
  //  std::cout<<"pointer : "<<file<<"  fname "<< file->GetName() <<std::endl;

  if(!file) {
    std::cout<<"FILE CAN'T BE OPENED"<<std::endl;
    return 0;
  }

  //  TH1D * h =  (TH1D*)gDirectory->Get("cutflow");
  TH1D * h =  (TH1D*)gDirectory->Get("counter_");
   //   cout<<"h name "<<h->GetName()<<endl;
   double all_events =  h->GetBinContent(2);
   double lumi = all_events/xsec;
   //double lumi = all_events;                                                                                           
   //   std::cout<<"lumi : "<<lumi<<"  evts : "<<all_events<<"  xsec : "<<xsec<<std::endl;
  file->Close();
  return lumi;
}