#define histoLeptonId_cxx
#include "histoLeptonId.h"

const Float_t Pi = 3.14159; 

void histoLeptonId::BookLeptonId(void){

  char title[50];
  char* var;

  var="#Delta R";
  char* xlabel="#Delta R";
  sprintf(title, "e (HEPG-REC) %s", var);
  fHist.hele0_delR = new TH1F("DeltaR0_Ele",title,200,0.,0.1);
  fHist.hele0_delR ->GetXaxis()->SetTitle(xlabel);

  var="im35";
  xlabel="Mll(GeV)";
  sprintf(title, " %s", var);
  fHist.hepg_im35 = new TH1F("im35",title,800,155.,165.);
  fHist.hepg_im35 ->GetXaxis()->SetTitle(xlabel);

  var="pt35";
  xlabel="pt25(GeV)";
  sprintf(title, " %s", var);
  fHist.hepg_pt35 = new TH1F("pt35",title,200,0.,200.);
  fHist.hepg_pt35 ->GetXaxis()->SetTitle(xlabel);

  var="#Delta R";
  xlabel="#Delta R";
  sprintf(title, "#mu (HEPG-REC) %s", var);
  fHist.hmu0_delR = new TH1F("DeltaR0_Mu",title,200,0.,0.1);
  fHist.hmu0_delR ->GetXaxis()->SetTitle(xlabel);


  var="#Delta R";
  xlabel="#Delta R";
  sprintf(title, "e (HEPG-REC) %s", var);
  fHist.hele_delR = new TH1F("DeltaR_Ele",title,200,0.,0.1);
  fHist.hele_delR ->GetXaxis()->SetTitle(xlabel);

  var="#Delta R";
  xlabel="#Delta R";
  sprintf(title, "#mu (HEPG-REC) %s", var);
  fHist.hmu_delR = new TH1F("DeltaR_Mu",title,200,0.,0.1);
  fHist.hmu_delR ->GetXaxis()->SetTitle(xlabel);

  var="Nele";
  xlabel="Nele";
  sprintf(title, " %s", var);
  fHist.he_n = new TH1F("Ne",title,5,0.,5.);
  fHist.he_n ->GetXaxis()->SetTitle(xlabel);


  var="RawTotalE";
  xlabel="RawTotalE";
  sprintf(title, " %s", var);
  fHist.he_rawtotale = new TH1F("RawTotalE",title,100,0.,500.);
  fHist.he_rawtotale ->GetXaxis()->SetTitle(xlabel);

  var="RawTotalEt";
  xlabel="RawTotalEt";
  sprintf(title, " %s", var);
  fHist.he_rawtotalet = new TH1F("RawTotalEt",title,100,0.,500.);
  fHist.he_rawtotalet ->GetXaxis()->SetTitle(xlabel);

  var="CorrTotalEt";
  xlabel="CorrTotalEt";
  sprintf(title, " %s", var);
  fHist.he_corrtotalet = new TH1F("CorrTotalEt",title,100,0.,500.);
  fHist.he_corrtotalet ->GetXaxis()->SetTitle(xlabel);


  var="Nelec";
  xlabel="Nelec";
  sprintf(title, " %s", var);
  fHist.hnelec = new TH1F("Nelec",title,4,0.,4.);
  fHist.hnelec ->GetXaxis()->SetTitle(xlabel);

  var="Nmuon";
  xlabel="Nmuon";
  sprintf(title, " %s", var);
  fHist.hnmuon = new TH1F("Nmuon",title,4,0.,4.);
  fHist.hnmuon ->GetXaxis()->SetTitle(xlabel);


  /////ID VARIABLES BASELINE ELECTRON- Before Base Line Cuts
  
  var="Et(two-twer Em sum * sin(theta) of track)";
  xlabel="Et(GeV)";
  sprintf(title, " %s", var);
  fHist.hbe_et = new TH1F("Et_BLEle",title,50,0.,200.);
  fHist.hbe_et ->GetXaxis()->SetTitle(xlabel);

  var="Track Pt (COT-only beam-constrained track)";
  xlabel="Pt(GeV/c)";
  sprintf(title, " %s", var);
  fHist.hbe_pt = new TH1F("Pt_BLEle",title,100,0.,200.);
  fHist.hbe_pt ->GetXaxis()->SetTitle(xlabel);  

  var="EoP";
  xlabel="EoP";
  sprintf(title, " %s", var);
  fHist.hbe_eop = new TH1F("EoP_BLEle",title,30,0.,3.);
  fHist.hbe_eop ->GetXaxis()->SetTitle(xlabel);  

  var="z0";
  xlabel="z0 cm";
  sprintf(title, " %s", var);
  fHist.hbe_z0 = new TH1F("Z0_BLEle",title,100,-100.,100.);
  fHist.hbe_z0 ->GetXaxis()->SetTitle(xlabel);  

  var="Had/Em";
  xlabel="Had/Em";
  sprintf(title, " %s", var);
  fHist.hbe_hadoem = new TH1F("HadoEm_BLEle",title,40,0.,0.3);
  fHist.hbe_hadoem ->GetXaxis()->SetTitle(xlabel);  

  var="Lshr";
  xlabel="Lshr";
  sprintf(title, " %s", var);
  fHist.hbe_lshr = new TH1F("Lshr_BLEle",title,60,-1.5,1.5);
  fHist.hbe_lshr ->GetXaxis()->SetTitle(xlabel);  

  var="CSDeltaX";
  xlabel="CSDeltaX";
  sprintf(title, " %s", var);
  fHist.hbe_csdeltax = new TH1F("CSdeltaX_BLEle",title,50,-5.,5.);
  fHist.hbe_csdeltax ->GetXaxis()->SetTitle(xlabel); 

  var="DeltaZ";
  xlabel="DeltaZ";
  sprintf(title, " %s", var);
  fHist.hbe_deltaz = new TH1F("DeltaZ_BLEle",title,50,-5.,5.);
  fHist.hbe_deltaz ->GetXaxis()->SetTitle(xlabel); 

  var="Chi2";
  xlabel="Chi2";
  sprintf(title, " %s", var);
  fHist.hbe_chi2 = new TH1F("Chi2_BLEle",title,50,-1.,20.);
  fHist.hbe_chi2 ->GetXaxis()->SetTitle(xlabel); 

  var="FCalIsoPJW";
  xlabel="FCalIsoPJW";
  sprintf(title, " %s", var);
  fHist.hbe_fcaliso = new TH1F("FCalIsoPJW_BLEle",title,50,-0.05,0.15);
  fHist.hbe_fcaliso ->GetXaxis()->SetTitle(xlabel); 

  var="TrkPhi";
  xlabel="TrkPhi";
  sprintf(title, " %s", var);
  fHist.hbe_trkphi = new TH1F("TrkPhi_BLEle",title,35,0.,7.);
  fHist.hbe_trkphi ->GetXaxis()->SetTitle(xlabel); 

  var="TrkEta";
  xlabel="TrkEta";
  sprintf(title, " %s", var);
  fHist.hbe_trketa = new TH1F("TrkEta_BLEle",title,60,-2.,2.);
  fHist.hbe_trketa ->GetXaxis()->SetTitle(xlabel); 


  /////ID VARIABLES BASELINE ELECTRON- After Base Line Cuts
  
  var="Et(two-twer Em sum * sin(theta) of track)";
  xlabel="Et(GeV)";
  sprintf(title, " %s", var);
  fHist.hbec_et = new TH1F("Et_BLEleC",title,50,0.,200.);
  fHist.hbec_et ->GetXaxis()->SetTitle(xlabel);

  var="Track Pt (COT-only beam-constrained track)";
  xlabel="Pt(GeV/c)";
  sprintf(title, " %s", var);
  fHist.hbec_pt = new TH1F("Pt_BLEleC",title,100,0.,200.);
  fHist.hbec_pt ->GetXaxis()->SetTitle(xlabel);  

  var="EoP";
  xlabel="EoP";
  sprintf(title, " %s", var);
  fHist.hbec_eop = new TH1F("EoP_BLEleC",title,30,0.,3.);
  fHist.hbec_eop ->GetXaxis()->SetTitle(xlabel);  

  var="z0";
  xlabel="z0 cm";
  sprintf(title, " %s", var);
  fHist.hbec_z0 = new TH1F("Z0_BLEleC",title,100,-100.,100.);
  fHist.hbec_z0 ->GetXaxis()->SetTitle(xlabel);  

  var="Had/Em";
  xlabel="Had/Em";
  sprintf(title, " %s", var);
  fHist.hbec_hadoem = new TH1F("HadoEm_BLEleC",title,40,0.,0.3);
  fHist.hbec_hadoem ->GetXaxis()->SetTitle(xlabel);  

  var="Lshr";
  xlabel="Lshr";
  sprintf(title, " %s", var);
  fHist.hbec_lshr = new TH1F("Lshr_BLEleC",title,60,-1.5,1.5);
  fHist.hbec_lshr ->GetXaxis()->SetTitle(xlabel);  

  var="CSDeltaX";
  xlabel="CSDeltaX";
  sprintf(title, " %s", var);
  fHist.hbec_csdeltax = new TH1F("CSdeltaX_BLEleC",title,50,-5.,5.);
  fHist.hbec_csdeltax ->GetXaxis()->SetTitle(xlabel); 

  var="DeltaZ";
  xlabel="DeltaZ";
  sprintf(title, " %s", var);
  fHist.hbec_deltaz = new TH1F("DeltaZ_BLEleC",title,50,-5.,5.);
  fHist.hbec_deltaz ->GetXaxis()->SetTitle(xlabel); 

  var="Chi2";
  xlabel="Chi2";
  sprintf(title, " %s", var);
  fHist.hbec_chi2 = new TH1F("Chi2_BLEleC",title,50,-1.,20.);
  fHist.hbec_chi2 ->GetXaxis()->SetTitle(xlabel); 

  var="FCalIsoPJW";
  xlabel="FCalIsoPJW";
  sprintf(title, " %s", var);
  fHist.hbec_fcaliso = new TH1F("FCalIsoPJW_BLEleC",title,50,-0.05,0.15);
  fHist.hbec_fcaliso ->GetXaxis()->SetTitle(xlabel); 

  var="TrkPhi";
  xlabel="TrkPhi";
  sprintf(title, " %s", var);
  fHist.hbec_trkphi = new TH1F("TrkPhi_BLEleC",title,35,0.,7.);
  fHist.hbec_trkphi ->GetXaxis()->SetTitle(xlabel); 

  var="TrkEta";
  xlabel="TrkEta";
  sprintf(title, " %s", var);
  fHist.hbec_trketa = new TH1F("TrkEta_BLEleC",title,60,-2.,2.);
  fHist.hbec_trketa ->GetXaxis()->SetTitle(xlabel); 

  /////ID VARIABLES BASELINE MUON- Before Base Line Cuts

  var="Track Pt (COT-only beam-constrained track)";
  xlabel="Pt(GeV/c)";
  sprintf(title, " %s", var);
  fHist.hbm_pt = new TH1F("Pt_BLMu",title,100,0.,200.);
  fHist.hbm_pt ->GetXaxis()->SetTitle(xlabel);  

  var="z0";
  xlabel="z0 cm";
  sprintf(title, " %s", var);
  fHist.hbm_z0 = new TH1F("Z0_BLMu",title,100,-100.,100.);
  fHist.hbm_z0 ->GetXaxis()->SetTitle(xlabel);  

  var="Em";
  xlabel="Em";
  sprintf(title, " %s", var);
  fHist.hbm_em = new TH1F("Em_BLMu",title,60,0.,6.);
  fHist.hbm_em ->GetXaxis()->SetTitle(xlabel);  


  var="Had";
  xlabel="Had";
  sprintf(title, " %s", var);
  fHist.hbm_had = new TH1F("Had_BLMu",title,50,0.,10.);
  fHist.hbm_had ->GetXaxis()->SetTitle(xlabel);  

  var="d0WS";
  xlabel="d0 (cm)";
  sprintf(title, " %s", var);
  fHist.hbm_d0ws = new TH1F("D0WS_BLMu",title,50,-0.05,0.05);
  fHist.hbm_d0ws ->GetXaxis()->SetTitle(xlabel);  

  var="d0WoS";
  xlabel="d0 (cm)";
  sprintf(title, " %s", var);
  fHist.hbm_d0wos = new TH1F("D0WOS_BLMu",title,50,-0.5,0.5);
  fHist.hbm_d0wos ->GetXaxis()->SetTitle(xlabel);  

  var="FCIso";
  xlabel="FCIso";
  sprintf(title, " %s", var);
  fHist.hbm_fciso = new TH1F("FCIso_BLMu",title,25,-0.1,0.4);
  fHist.hbm_fciso ->GetXaxis()->SetTitle(xlabel);  

  var="DeltaX";
  xlabel="DeltaX";
  sprintf(title, " %s", var);
  fHist.hbm_deltax_cmu = new TH1F("DeltaX_CMU_BLMu",title,60,-8.,8.);
  fHist.hbm_deltax_cmu ->GetXaxis()->SetTitle(xlabel);  

  var="DeltaX";
  xlabel="DeltaX";
  sprintf(title, " %s", var);
  fHist.hbm_deltax_cmp = new TH1F("DeltaX_CMP_BLMu",title,60,-8.,8.);
  fHist.hbm_deltax_cmp ->GetXaxis()->SetTitle(xlabel); 

  var="DeltaX";
  xlabel="DeltaX";
  sprintf(title, " %s", var);
  fHist.hbm_deltax_cmx = new TH1F("DeltaX_CMX_BLMu",title,60,-8.,8.);
  fHist.hbm_deltax_cmx ->GetXaxis()->SetTitle(xlabel); 


  var="TrkPhi-Cmup";
  xlabel="TrkPhi-Cmup";
  sprintf(title, " %s", var);
  fHist.hbm_trkphi_cmup = new TH1F("TrkPhi_CMUP_BLMu",title,70,0.,7.);
  fHist.hbm_trkphi_cmup ->GetXaxis()->SetTitle(xlabel); 

  var="TrkEta-Cmup";
  xlabel="TrkEta-Cmup";
  sprintf(title, " %s", var);
  fHist.hbm_trketa_cmup = new TH1F("TrkEta_CMUP_BLMu",title,80,-2.,2.);
  fHist.hbm_trketa_cmup ->GetXaxis()->SetTitle(xlabel); 



  var="TrkPhi-Cmu";
  xlabel="TrkPhi-Cmu";
  sprintf(title, " %s", var);
  fHist.hbm_trkphi_cmu = new TH1F("TrkPhi_CMU_BLMu",title,70,0.,7.);
  fHist.hbm_trkphi_cmu ->GetXaxis()->SetTitle(xlabel); 

  var="TrkEta-Cmu";
  xlabel="TrkEta-Cmu";
  sprintf(title, " %s", var);
  fHist.hbm_trketa_cmu = new TH1F("TrkEta_CMU_BLMu",title,80,-2.,2.);
  fHist.hbm_trketa_cmu ->GetXaxis()->SetTitle(xlabel); 



  var="TrkPhi-Cmp";
  xlabel="TrkPhi-Cmp";
  sprintf(title, " %s", var);
  fHist.hbm_trkphi_cmp = new TH1F("TrkPhi_CMP_BLMu",title,70,0.,7.);
  fHist.hbm_trkphi_cmp ->GetXaxis()->SetTitle(xlabel); 

  var="TrkEta-Cmp";
  xlabel="TrkEta-Cmp";
  sprintf(title, " %s", var);
  fHist.hbm_trketa_cmp = new TH1F("TrkEta_CMP_BLMu",title,80,-2.,2.);
  fHist.hbm_trketa_cmp ->GetXaxis()->SetTitle(xlabel); 



  var="TrkPhi-Cmx";
  xlabel="TrkPhi-Cmx";
  sprintf(title, " %s", var);
  fHist.hbm_trkphi_cmx = new TH1F("TrkPhi_CMX_BLMu",title,70,0.,7.);
  fHist.hbm_trkphi_cmx ->GetXaxis()->SetTitle(xlabel); 

  var="TrkEta-Cmx";
  xlabel="TrkEta-Cmx";
  sprintf(title, " %s", var);
  fHist.hbm_trketa_cmx = new TH1F("TrkEta_CMX_BLMu",title,80,-2.,2.);
  fHist.hbm_trketa_cmx ->GetXaxis()->SetTitle(xlabel); 



  /////ID VARIABLES BASELINE MUON- After Base Line Cuts

  var="Track Pt (COT-only beam-constrained track)";
  xlabel="Pt(GeV/c)";
  sprintf(title, " %s", var);
  fHist.hbmc_pt = new TH1F("Pt_BLMuC",title,100,0.,200.);
  fHist.hbmc_pt ->GetXaxis()->SetTitle(xlabel);  

  var="z0";
  xlabel="z0 cm";
  sprintf(title, " %s", var);
  fHist.hbmc_z0 = new TH1F("Z0_BLMuC",title,100,-100.,100.);
  fHist.hbmc_z0 ->GetXaxis()->SetTitle(xlabel);  

  var="Em";
  xlabel="Em";
  sprintf(title, " %s", var);
  fHist.hbmc_em = new TH1F("Em_BLMuC",title,60,0.,6.);
  fHist.hbmc_em ->GetXaxis()->SetTitle(xlabel);  


  var="Had";
  xlabel="Had";
  sprintf(title, " %s", var);
  fHist.hbmc_had = new TH1F("Had_BLMuC",title,50,0.,10.);
  fHist.hbmc_had ->GetXaxis()->SetTitle(xlabel);  

  var="d0 WS";
  xlabel="d0 (cm)";
  sprintf(title, " %s", var);
  fHist.hbmc_d0ws = new TH1F("D0WS_BLMuC",title,50,-0.05,0.05);
  fHist.hbmc_d0ws ->GetXaxis()->SetTitle(xlabel);  

  var="d0 WoS";
  xlabel="d0 (cm)";
  sprintf(title, " %s", var);
  fHist.hbmc_d0wos = new TH1F("D0WOS_BLMuC",title,50,-0.5,0.5);
  fHist.hbmc_d0wos ->GetXaxis()->SetTitle(xlabel);  

  var="FCIso";
  xlabel="FCIso";
  sprintf(title, " %s", var);
  fHist.hbmc_fciso = new TH1F("FCIso_BLMuC",title,25,-0.1,0.4);
  fHist.hbmc_fciso ->GetXaxis()->SetTitle(xlabel);  

  var="DeltaX";
  xlabel="DeltaX";
  sprintf(title, " %s", var);
  fHist.hbmc_deltax_cmu = new TH1F("DeltaX_CMU_BLMuC",title,60,-8.,8.);
  fHist.hbmc_deltax_cmu ->GetXaxis()->SetTitle(xlabel);  

  var="DeltaX";
  xlabel="DeltaX";
  sprintf(title, " %s", var);
  fHist.hbmc_deltax_cmp = new TH1F("DeltaX_CMP_BLMuC",title,60,-8.,8.);
  fHist.hbmc_deltax_cmp ->GetXaxis()->SetTitle(xlabel); 

  var="DeltaX";
  xlabel="DeltaX";
  sprintf(title, " %s", var);
  fHist.hbmc_deltax_cmx = new TH1F("DeltaX_CMX_BLMuC",title,60,-8.,8.);
  fHist.hbmc_deltax_cmx ->GetXaxis()->SetTitle(xlabel); 



  var="TrkPhi-Cmup";
  xlabel="TrkPhi-Cmup";
  sprintf(title, " %s", var);
  fHist.hbmc_trkphi_cmup = new TH1F("TrkPhi_CMUP_BLMuC",title,70,0.,7.);
  fHist.hbmc_trkphi_cmup ->GetXaxis()->SetTitle(xlabel); 

  var="TrkEta-Cmup";
  xlabel="TrkEta-Cmup";
  sprintf(title, " %s", var);
  fHist.hbmc_trketa_cmup = new TH1F("TrkEta_CMUP_BLMuC",title,80,-2.,2.);
  fHist.hbmc_trketa_cmup ->GetXaxis()->SetTitle(xlabel); 



  var="TrkPhi-Cmu";
  xlabel="TrkPhi-Cmu";
  sprintf(title, " %s", var);
  fHist.hbmc_trkphi_cmu = new TH1F("TrkPhi_CMU_BLMuC",title,70,0.,7.);
  fHist.hbmc_trkphi_cmu ->GetXaxis()->SetTitle(xlabel); 

  var="TrkEta-Cmu";
  xlabel="TrkEta-Cmu";
  sprintf(title, " %s", var);
  fHist.hbmc_trketa_cmu = new TH1F("TrkEta_CMU_BLMuC",title,80,-2.,2.);
  fHist.hbmc_trketa_cmu ->GetXaxis()->SetTitle(xlabel); 



  var="TrkPhi-Cmp";
  xlabel="TrkPhi-Cmp";
  sprintf(title, " %s", var);
  fHist.hbmc_trkphi_cmp = new TH1F("TrkPhi_CMP_BLMuC",title,70,0.,7.);
  fHist.hbmc_trkphi_cmp ->GetXaxis()->SetTitle(xlabel); 

  var="TrkEta-Cmp";
  xlabel="TrkEta-Cmp";
  sprintf(title, " %s", var);
  fHist.hbmc_trketa_cmp = new TH1F("TrkEta_CMP_BLMuC",title,80,-2.,2.);
  fHist.hbmc_trketa_cmp ->GetXaxis()->SetTitle(xlabel); 



  var="TrkPhi-Cmx";
  xlabel="TrkPhi-Cmx";
  sprintf(title, " %s", var);
  fHist.hbmc_trkphi_cmx = new TH1F("TrkPhi_CMX_BLMuC",title,70,0.,7.);
  fHist.hbmc_trkphi_cmx ->GetXaxis()->SetTitle(xlabel); 

  var="TrkEta-Cmx";
  xlabel="TrkEta-Cmx";
  sprintf(title, " %s", var);
  fHist.hbmc_trketa_cmx = new TH1F("TrkEta_CMX_BLMuC",title,80,-2.,2.);
  fHist.hbmc_trketa_cmx ->GetXaxis()->SetTitle(xlabel); 

  var="Flag-Dilep";
  xlabel="DileptonFlag";
  sprintf(title, " %s", var);
  fHist.hdilepflag = new TH1F("DileptonFlag",title,25,0.,25.);
  fHist.hdilepflag ->GetXaxis()->SetTitle(xlabel); 

  var="Njets";
  xlabel="Njets";
  sprintf(title, " %s", var);
  fHist.hjets_njets = new TH1F("Njets",title,5,0.,5.);
  fHist.hjets_njets ->GetXaxis()->SetTitle(xlabel); 

  var="Et-jet-1";
  xlabel="Et-jet-1";
  sprintf(title, " %s", var);
  fHist.hjets_et1 = new TH1F("Etj1",title,50,0.,100.);
  fHist.hjets_et1 ->GetXaxis()->SetTitle(xlabel); 


  var="Et-jet-2";
  xlabel="Et-jet-2";
  sprintf(title, " %s", var);
  fHist.hjets_et2 = new TH1F("Etj2",title,50,0.,100.);
  fHist.hjets_et2 ->GetXaxis()->SetTitle(xlabel); 


  var="Et-jet-3";
  xlabel="Et-jet-3";
  sprintf(title, " %s", var);
  fHist.hjets_et3 = new TH1F("Etj3",title,50,0.,100.);
  fHist.hjets_et3 ->GetXaxis()->SetTitle(xlabel); 

  var="Eta-jet-1";
  xlabel="Eta-jet-1";
  sprintf(title, " %s", var);
  fHist.hjets_eta1 = new TH1F("Etaj1",title,60,-3.,3.);
  fHist.hjets_eta1 ->GetXaxis()->SetTitle(xlabel);

  var="Eta-jet-2";
  xlabel="Eta-jet-2";
  sprintf(title, " %s", var);
  fHist.hjets_eta2 = new TH1F("Etaj2",title,60,-3.,3.);
  fHist.hjets_eta2 ->GetXaxis()->SetTitle(xlabel);

  var="Eta-jet-3";
  xlabel="Eta-jet-3";
  sprintf(title, " %s", var);
  fHist.hjets_eta3 = new TH1F("Etaj3",title,60,-3.,3.);
  fHist.hjets_eta3 ->GetXaxis()->SetTitle(xlabel);

  var="Phi-jet-1";
  xlabel="Phi-jet-1";
  sprintf(title, " %s", var);
  fHist.hjets_phi1 = new TH1F("Phij1",title,70,0.,7.);
  fHist.hjets_phi1 ->GetXaxis()->SetTitle(xlabel);
  
  var="Phi-jet-2";
  xlabel="Phi-jet-2";
  sprintf(title, " %s", var);
  fHist.hjets_phi2 = new TH1F("Phij2",title,70,0.,7.);
  fHist.hjets_phi2 ->GetXaxis()->SetTitle(xlabel);
  
  var="Phi-jet-3";
  xlabel="Phi-jet-3";
  sprintf(title, " %s", var);
  fHist.hjets_phi3 = new TH1F("Phij3",title,70,0.,7.);
  fHist.hjets_phi3 ->GetXaxis()->SetTitle(xlabel);
  

  //MET vs JET multiplicity plots


  var="Met";
  xlabel="Met";
  sprintf(title, "Met Nj==0 ee %s", var);
  fHist.hmet_j0_ee = new TH1F("Met_j0_ee",title,70,0.,140);
  fHist.hmet_j0_ee ->GetXaxis()->SetTitle(xlabel);

  var="Met";
  xlabel="Met";
  sprintf(title, "Met Nj==1 ee %s", var);
  fHist.hmet_j1_ee = new TH1F("Met_j1_ee",title,70,0.,140);
  fHist.hmet_j1_ee ->GetXaxis()->SetTitle(xlabel);

  var="Met";
  xlabel="Met";
  sprintf(title, "Met Nj>=2 ee %s", var);
  fHist.hmet_j2_ee = new TH1F("Met_j2_ee",title,70,0.,140);
  fHist.hmet_j2_ee ->GetXaxis()->SetTitle(xlabel);

  
  var="Met";
  xlabel="Met";
  sprintf(title, "Met Nj==0 emu %s", var);
  fHist.hmet_j0_emu = new TH1F("Met_j0_e#mu",title,70,0.,140);
  fHist.hmet_j0_emu ->GetXaxis()->SetTitle(xlabel);

  var="Met";
  xlabel="Met";
  sprintf(title, "Met Nj==1 emu %s", var);
  fHist.hmet_j1_emu = new TH1F("Met_j1_e#mu",title,70,0.,140);
  fHist.hmet_j1_emu ->GetXaxis()->SetTitle(xlabel);

  var="Met";
  xlabel="Met";
  sprintf(title, "Met Nj>=2 emu %s", var);
  fHist.hmet_j2_emu = new TH1F("Met_j2_e#mu",title,70,0.,140);
  fHist.hmet_j2_emu ->GetXaxis()->SetTitle(xlabel);

  
  var="Met";
  xlabel="Met";
  sprintf(title, "Met Nj==0 #mu#mu %s", var);
  fHist.hmet_j0_mumu = new TH1F("Met_j0_mumu",title,70,0.,140);
  fHist.hmet_j0_mumu ->GetXaxis()->SetTitle(xlabel);

  var="Met";
  xlabel="Met";
  sprintf(title, "Met Nj==1 #mu#mu %s", var);
  fHist.hmet_j1_mumu = new TH1F("Met_j1_mumu",title,70,0.,140);
  fHist.hmet_j1_mumu ->GetXaxis()->SetTitle(xlabel);

  var="Met";
  xlabel="Met";
  sprintf(title, "Met Nj>=2 #mu#mu %s", var);
  fHist.hmet_j2_mumu = new TH1F("Met_j2_mumu",title,70,0.,140);
  fHist.hmet_j2_mumu ->GetXaxis()->SetTitle(xlabel);


  //DeltaPhi vs MET plots

  //HEE

  var="Met";
  xlabel="Met";
  char* ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: all-jet", var);
  fHist.heeDphimljMetj = new TH2F("heeDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.heeDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.heeDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: 0-jet", var);
  fHist.heeDphimljMetj0 = new TH2F("heeDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.heeDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.heeDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: 1-jet", var);
  fHist.heeDphimljMetj1 = new TH2F("heeDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.heeDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.heeDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: >=2-jet", var);
  fHist.heeDphimljMetj2 = new TH2F("heeDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.heeDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.heeDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  //HMUMU

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: all-jet", var);
  fHist.hmumuDphimljMetj = new TH2F("hmumuDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.hmumuDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.hmumuDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: 0-jet", var);
  fHist.hmumuDphimljMetj0 = new TH2F("hmumuDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.hmumuDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hmumuDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: 1-jet", var);
  fHist.hmumuDphimljMetj1 = new TH2F("hmumuDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.hmumuDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.hmumuDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet", var);
  fHist.hmumuDphimljMetj2 = new TH2F("hmumuDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.hmumuDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hmumuDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  //HEMU

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: all-jet", var);
  fHist.hemuDphimljMetj = new TH2F("hemuDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.hemuDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.hemuDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: 0-jet", var);
  fHist.hemuDphimljMetj0 = new TH2F("hemuDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.hemuDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hemuDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: 1-jet", var);
  fHist.hemuDphimljMetj1 = new TH2F("hemuDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.hemuDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.hemuDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet", var);
  fHist.hemuDphimljMetj2 = new TH2F("hemuDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.hemuDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hemuDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);


  //HTEE

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: all-jet Ht>200", var);
  fHist.hhteeDphimljMetj = new TH2F("hhteeDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.hhteeDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.hhteeDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: 0-jet Ht>200", var);
  fHist.hhteeDphimljMetj0 = new TH2F("hhteeDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.hhteeDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hhteeDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: 1-jet Ht>200", var);
  fHist.hhteeDphimljMetj1 = new TH2F("hhteeDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.hhteeDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.hhteeDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: >=2-jet Ht>200", var);
  fHist.hhteeDphimljMetj2 = new TH2F("hhteeDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.hhteeDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hhteeDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: all-jet Ht<200", var);
  fHist.nothhteeDphimljMetj = new TH2F("nothhteeDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.nothhteeDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhteeDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: 0-jet Ht<200", var);
  fHist.nothhteeDphimljMetj0 = new TH2F("nothhteeDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.nothhteeDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhteeDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: 1-jet Ht<200", var);
  fHist.nothhteeDphimljMetj1 = new TH2F("nothhteeDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.nothhteeDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhteeDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: >=2-jet Ht<200", var);
  fHist.nothhteeDphimljMetj2 = new TH2F("nothhteeDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.nothhteeDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhteeDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  //HQEE

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: all-jet OS and Ht>200", var);
  fHist.hqeeDphimljMetj = new TH2F("hqeeDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.hqeeDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.hqeeDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: 0-jet OS and Ht>200", var);
  fHist.hqeeDphimljMetj0 = new TH2F("hqeeDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.hqeeDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqeeDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: 1-jet OS and Ht>200", var);
  fHist.hqeeDphimljMetj1 = new TH2F("hqeeDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.hqeeDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqeeDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: >=2-jet OS and Ht>200", var);
  fHist.hqeeDphimljMetj2 = new TH2F("hqeeDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.hqeeDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqeeDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:l-l) vs Met: >=2-jet OS and Ht>200", var);
  fHist.hqeeDphillj2 = new TH2F("hqeeDphillj2",title,70,0.,140.,140,0.,3.5);
  fHist.hqeeDphillj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqeeDphillj2 ->GetYaxis()->SetTitle(ylabel); 

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:l-l) vs Met: >=2-jet OS and Ht>200", var);
  fHist.hqeeDphillj0 = new TH2F("hqeeDphillj0",title,70,0.,140.,140,0.,3.5);
  fHist.hqeeDphillj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqeeDphillj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: all-jet SSor Ht<200", var);
  fHist.nothqeeDphimljMetj = new TH2F("nothqeeDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.nothqeeDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqeeDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: 0-jet SS or Ht<200", var);
  fHist.nothqeeDphimljMetj0 = new TH2F("nothqeeDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.nothqeeDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqeeDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: 1-jet SS or Ht<200", var);
  fHist.nothqeeDphimljMetj1 = new TH2F("nothqeeDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.nothqeeDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqeeDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:Met-l,j) vs Met: >=2-jet SS or Ht<200", var);
  fHist.nothqeeDphimljMetj2 = new TH2F("nothqeeDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.nothqeeDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqeeDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);
 
  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:l-l) vs Met: >=2-jet SS or Ht<200", var);
  fHist.nothqeeDphillj2 = new TH2F("nothqeeDphillj2",title,70,0.,140.,140,0.,3.5);
  fHist.nothqeeDphillj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqeeDphillj2 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "ee #Delta(#phi:l-l) vs Met: >=2-jet SS or Ht<200", var);
  fHist.nothqeeDphillj0 = new TH2F("nothqeeDphillj0",title,70,0.,140.,140,0.,3.5);
  fHist.nothqeeDphillj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqeeDphillj0 ->GetYaxis()->SetTitle(ylabel);

  //HTMUMU

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: all-jet Ht>200", var);
  fHist.hhtmumuDphimljMetj = new TH2F("hhtmumuDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.hhtmumuDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.hhtmumuDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: 0-jet Ht>200", var);
  fHist.hhtmumuDphimljMetj0 = new TH2F("hhtmumuDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.hhtmumuDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hhtmumuDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: 1-jet Ht>200", var);
  fHist.hhtmumuDphimljMetj1 = new TH2F("hhtmumuDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.hhtmumuDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.hhtmumuDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet Ht>200", var);
  fHist.hhtmumuDphimljMetj2 = new TH2F("hhtmumuDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.hhtmumuDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hhtmumuDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

   var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: all-jet Ht<200", var);
  fHist.nothhtmumuDphimljMetj = new TH2F("nothhtmumuDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.nothhtmumuDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhtmumuDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: 0-jet Ht<200", var);
  fHist.nothhtmumuDphimljMetj0 = new TH2F("nothhtmumuDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.nothhtmumuDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhtmumuDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: 1-jet Ht<200", var);
  fHist.nothhtmumuDphimljMetj1 = new TH2F("nothhtmumuDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.nothhtmumuDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhtmumuDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet Ht<200", var);
  fHist.nothhtmumuDphimljMetj2 = new TH2F("nothhtmumuDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.nothhtmumuDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhtmumuDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  // HQMUMU

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: all-jet OS and Ht>200", var);
  fHist.hqmumuDphimljMetj = new TH2F("hqmumuDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.hqmumuDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.hqmumuDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: 0-jet OS and Ht>200", var);
  fHist.hqmumuDphimljMetj0 = new TH2F("hqmumuDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.hqmumuDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqmumuDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: 1-jet OS and Ht>200", var);
  fHist.hqmumuDphimljMetj1 = new TH2F("hqmumuDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.hqmumuDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqmumuDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet OS and Ht>200", var);
  fHist.hqmumuDphimljMetj2 = new TH2F("hqmumuDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.hqmumuDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqmumuDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:l-l) vs Met: >=2-jet OS and Ht>200", var);
  fHist.hqmumuDphillj2 = new TH2F("hqmumuDphillj2",title,70,0.,140.,140,0.,3.5);
  fHist.hqmumuDphillj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqmumuDphillj2 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:l-l) vs Met: >=2-jet OS and Ht>200", var);
  fHist.hqmumuDphillj0 = new TH2F("hqmumuDphillj0",title,70,0.,140.,140,0.,3.5);
  fHist.hqmumuDphillj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqmumuDphillj0 ->GetYaxis()->SetTitle(ylabel);

var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: all-jet SS and Ht<200", var);
  fHist.nothqmumuDphimljMetj = new TH2F("nothqmumuDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.nothqmumuDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqmumuDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: 0-jet SS and Ht<200", var);
  fHist.nothqmumuDphimljMetj0 = new TH2F("nothqmumuDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.nothqmumuDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqmumuDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: 1-jet SS and Ht<200", var);
  fHist.nothqmumuDphimljMetj1 = new TH2F("nothqmumuDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.nothqmumuDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqmumuDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet SS and Ht<200", var);
  fHist.nothqmumuDphimljMetj2 = new TH2F("nothqmumuDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.nothqmumuDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqmumuDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:l-l) vs Met: >=2-jet SS and Ht<200", var);
  fHist.nothqmumuDphillj2 = new TH2F("nothqmumuDphillj2",title,70,0.,140.,140,0.,3.5);
  fHist.nothqmumuDphillj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqmumuDphillj2 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "#mu#mu #Delta(#phi:l-l) vs Met: >=2-jet SS and Ht<200", var);
  fHist.nothqmumuDphillj0 = new TH2F("nothqmumuDphillj0",title,70,0.,140.,140,0.,3.5);
  fHist.nothqmumuDphillj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqmumuDphillj0 ->GetYaxis()->SetTitle(ylabel);

  // HTEMU


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: all-jet Ht>200", var);
  fHist.hhtemuDphimljMetj = new TH2F("hhtemuDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.hhtemuDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.hhtemuDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: 0-jet Ht>200", var);
  fHist.hhtemuDphimljMetj0 = new TH2F("hhtemuDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.hhtemuDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hhtemuDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: 1-jet Ht>200", var);
  fHist.hhtemuDphimljMetj1 = new TH2F("hhtemuDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.hhtemuDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.hhtemuDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet Ht>200", var);
  fHist.hhtemuDphimljMetj2 = new TH2F("hhtemuDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.hhtemuDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hhtemuDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);


  
  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: all-jet Ht<200", var);
  fHist.nothhtemuDphimljMetj = new TH2F("nothhtemuDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.nothhtemuDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhtemuDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: 0-jet Ht<200", var);
  fHist.nothhtemuDphimljMetj0 = new TH2F("nothhtemuDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.nothhtemuDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhtemuDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: 1-jet Ht<200", var);
  fHist.nothhtemuDphimljMetj1 = new TH2F("nothhtemuDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.nothhtemuDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhtemuDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet Ht<200", var);
  fHist.nothhtemuDphimljMetj2 = new TH2F("nothhtemuDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.nothhtemuDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothhtemuDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  //HQEMU

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: all-jet SS and Ht>200", var);
  fHist.hqemuDphimljMetj = new TH2F("hqemuDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.hqemuDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.hqemuDphimljMetj ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: 0-jet SS and Ht>200", var);
  fHist.hqemuDphimljMetj0 = new TH2F("hqemuDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.hqemuDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqemuDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: 1-jet SS and Ht>200", var);
  fHist.hqemuDphimljMetj1 = new TH2F("hqemuDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.hqemuDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqemuDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet SS and Ht>200", var);
  fHist.hqemuDphimljMetj2 = new TH2F("hqemuDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.hqemuDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqemuDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:l-l) vs Met: >=2-jet SS and Ht>200", var);
  fHist.hqemuDphillj2 = new TH2F("hqemuDphimllj2",title,70,0.,140.,140,0.,3.5);
  fHist.hqemuDphillj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqemuDphillj2 ->GetYaxis()->SetTitle(ylabel);

   var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:l-l) vs Met: >=2-jet SS and Ht>200", var);
  fHist.hqemuDphillj0 = new TH2F("hqemuDphimllj0",title,70,0.,140.,140,0.,3.5);
  fHist.hqemuDphillj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.hqemuDphillj0 ->GetYaxis()->SetTitle(ylabel);


  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: all-jet OS and Ht<200", var);
  fHist.nothqemuDphimljMetj = new TH2F("nothqemuDphimljMetj",title,70,0.,140.,140,0.,3.5);
  fHist.nothqemuDphimljMetj ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqemuDphimljMetj ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: 0-jet OS and Ht<200", var);
  fHist.nothqemuDphimljMetj0 = new TH2F("nothqemuDphimljMetj0",title,70,0.,140.,140,0.,3.5);
  fHist.nothqemuDphimljMetj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqemuDphimljMetj0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: 1-jet OS and Ht<200", var);
  fHist.nothqemuDphimljMetj1 = new TH2F("nothqemuDphimljMetj1",title,70,0.,140.,140,0.,3.5);
  fHist.nothqemuDphimljMetj1 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqemuDphimljMetj1 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:Met-l,j)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet OS and Ht<200", var);
  fHist.nothqemuDphimljMetj2 = new TH2F("nothqemuDphimljMetj2",title,70,0.,140.,140,0.,3.5);
  fHist.nothqemuDphimljMetj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqemuDphimljMetj2 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:l-l)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet OS and Ht<200", var);
  fHist.nothqemuDphillj2 = new TH2F("nothqemuDphillj2",title,70,0.,140.,140,0.,3.5);
  fHist.nothqemuDphillj2 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqemuDphillj2 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="#Delta(#phi:l-l)";
  sprintf(title, "e#mu #Delta(#phi:Met-l,j) vs Met: >=2-jet OS and Ht<200", var);
  fHist.nothqemuDphillj0 = new TH2F("nothqemuDphillj0",title,70,0.,140.,140,0.,3.5);
  fHist.nothqemuDphillj0 ->GetXaxis()->SetTitle(xlabel);
  fHist.nothqemuDphillj0 ->GetYaxis()->SetTitle(ylabel);

  //Zmass plots


  var="M_{ee}";
  xlabel="[GeV/c]";
  sprintf(title, "WW->dil  %s", var);
  fHist.hmee = new TH1F("ZeeMass",title,100,0.,200.);
  fHist.hmee ->GetXaxis()->SetTitle(xlabel);

  var="M_{mumu}";
  xlabel="[GeV/c]";
  sprintf(title, "WW->dil %s", var);
  fHist.hmmumu = new TH1F("ZmumMass",title,100,0.,200.);
  fHist.hmmumu ->GetXaxis()->SetTitle(xlabel);

  var="M_{ee}";
  xlabel="[GeV/c]";
  sprintf(title, "WW->dil  %s", var);
  fHist.hmee0 = new TH1F("ZeeMass0",title,100,0.,200.);
  fHist.hmee0 ->GetXaxis()->SetTitle(xlabel);

  var="M_{mumu}";
  xlabel="[GeV/c]";
  sprintf(title, "WW->dil %s", var);
  fHist.hmmumu0 = new TH1F("ZmumMass0",title,100,0.,200.);
  fHist.hmmumu0 ->GetXaxis()->SetTitle(xlabel);

  var="DPhiJEle";
  xlabel="[rad]";
  sprintf(title, "WW->dil %s", var);
  fHist.hdphijele = new TH1F("DPhiJEle",title,70,0.,7.);
  fHist.hdphijele ->GetXaxis()->SetTitle(xlabel);

  var="DRJEle";
  xlabel="[adim]";
  sprintf(title, "WW->dil %s", var);
  fHist.hRjele = new TH1F("DRJEle",title,40,0.,2.);
  fHist.hRjele ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->all %s", var);
  fHist.hallHtj = new TH1F("Htall",title,50,0.,500.);
  fHist.hallHtj ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->all %s", var);
  fHist.hallHtj0 = new TH1F("Htall0",title,50,0.,500.);
  fHist.hallHtj0 ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->all %s", var);
  fHist.hallHtj1inc = new TH1F("Htall1inc",title,50,0.,500.);
  fHist.hallHtj1inc ->GetXaxis()->SetTitle(xlabel);


  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "TT->all %s", var);
  fHist.hallHtj2 = new TH1F("Htall2",title,50,0.,500.);
  fHist.hallHtj2 ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->ee %s", var);
  fHist.heeHt = new TH1F("HtEE",title,50,0.,500.);
  fHist.heeHt ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->ee %s", var);
  fHist.heeHtTow = new TH1F("HtEETow",title,75,0.,750.);
  fHist.heeHtTow ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->ee 0-jets %s", var);
  fHist.heeHtj0 = new TH1F("HtEE0",title,50,0.,500.);
  fHist.heeHtj0 ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->ee 1-jet  %s", var);
  fHist.heeHtj1 = new TH1F("HtEE1",title,50,0.,500.);
  fHist.heeHtj1 ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->ee >=1-jet  %s", var);
  fHist.heeHtj1inc = new TH1F("HtEE1inc",title,50,0.,500.);
  fHist.heeHtj1inc ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->ee 2-jets %s", var);
  fHist.heeHtj2 = new TH1F("HtEE2",title,50,0.,500.);
  fHist.heeHtj2 ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->mumu %s", var);
  fHist.hmumuHtTow = new TH1F("HtMUMUTow",title,75,0.,750.);
  fHist.hmumuHtTow ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->mumu %s", var);
  fHist.hmumuHt = new TH1F("HtMUMU",title,50,0.,500.);
  fHist.hmumuHt ->GetXaxis()->SetTitle(xlabel);
  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->mumu 0-jets %s", var);
  fHist.hmumuHtj0 = new TH1F("HtMUMU0",title,50,0.,500.);
  fHist.hmumuHtj0 ->GetXaxis()->SetTitle(xlabel);
  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->mumu 1-jet  %s", var);
  fHist.hmumuHtj1 = new TH1F("HtMUMU1",title,50,0.,500.);
  fHist.hmumuHtj1 ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->mumu >=1-jet  %s", var);
  fHist.hmumuHtj1inc = new TH1F("HtMUMU1inc",title,50,0.,500.);
  fHist.hmumuHtj1inc ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->mumu 2-jets %s", var);
  fHist.hmumuHtj2 = new TH1F("HtMUMU2",title,50,0.,500.);
  fHist.hmumuHtj2 ->GetXaxis()->SetTitle(xlabel);


  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->emu %s", var);
  fHist.hemuHt = new TH1F("HtEMU",title,50,0.,500.);
  fHist.hemuHt ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->emu %s", var);
  fHist.hemuHtTow = new TH1F("HtEMUTow",title,75,0.,750.);
  fHist.hemuHtTow ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->emu 0-jets %s", var);
  fHist.hemuHtj0 = new TH1F("HtEMU0",title,50,0.,500.);
  fHist.hemuHtj0 ->GetXaxis()->SetTitle(xlabel);
  var="Ht";

  xlabel="[GeV]";
  sprintf(title, "WW->emu 1-jet  %s", var);
  fHist.hemuHtj1 = new TH1F("HtEMU1",title,50,0.,500.);
  fHist.hemuHtj1 ->GetXaxis()->SetTitle(xlabel);

  xlabel="[GeV]";
  sprintf(title, "WW->emu >= 1-jet  %s", var);
  fHist.hemuHtj1inc = new TH1F("HtEMU1inc",title,50,0.,500.);
  fHist.hemuHtj1inc ->GetXaxis()->SetTitle(xlabel);

  var="Ht";
  xlabel="[GeV]";
  sprintf(title, "WW->emu 2-jets %s", var);
  fHist.hemuHtj2 = new TH1F("HtEMU2",title,50,0.,500.);
  fHist.hemuHtj2 ->GetXaxis()->SetTitle(xlabel);


  var="Et Ele (Met<20GeV AND Ht>200 GeV)";
  xlabel="Et(GeV)";
  sprintf(title, " %s", var);
  fHist.h_Ete_j0met20 = new TH1F("h_Ete_j0met20",title,60,0.,360.);
  fHist.h_Ete_j0met20 ->GetXaxis()->SetTitle(xlabel);

  var="Ht (Met<20GeV AND Ht>200 GeV)";
  xlabel="Et(GeV)";
  sprintf(title, " %s", var);
  fHist.h_Ht_j0met20 = new TH1F("h_Ht_j0met20",title,60,0.,360.);
  fHist.h_Ht_j0met20 ->GetXaxis()->SetTitle(xlabel);

  var="Ht-bis (Met<20GeV AND Ht>200 GeV)";
  xlabel="Et(GeV)";
  sprintf(title, " %s", var);
  fHist.h_Htbis_j0met20 = new TH1F("h_Htbis_j0met20",title,60,0.,360.);
  fHist.h_Htbis_j0met20 ->GetXaxis()->SetTitle(xlabel);

  var="DrLepJet";
  xlabel="[adim]";
  sprintf(title, "WW->dil %s", var);
  fHist.h_dRlj = new TH1F("DrLepJet",title,40,0.,2.);
  fHist.h_dRlj ->GetXaxis()->SetTitle(xlabel);

  var="Met";
  xlabel="Met";
  ylabel="Ht";
  sprintf(title, "ee Ht vs Met: 0-jets ", var);
  fHist.h_Ht_Met_j0 = new TH2F("Ht_Met_j0",title,120,0.,360.,120,0.,360.);
  fHist.h_Ht_Met_j0 ->GetXaxis()->SetTitle(xlabel);
  fHist.h_Ht_Met_j0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="Ete1+Ete2";
  sprintf(title, "ee Ete1+Ete2 vs Met: 0-jets ", var);
  fHist.h_Etee_Met_j0 = new TH2F("Ht_Etee_j0",title,120,0.,360.,120,0.,360.);
  fHist.h_Etee_Met_j0 ->GetXaxis()->SetTitle(xlabel);
  fHist.h_Etee_Met_j0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="Ete1+Ete2";
  sprintf(title, "ee Ete1+Ete2 >100 vs Met: 0-jets ", var);
  fHist.h1_Etee_PhiMet_j0 = new TH2F("h1_Etee_PhiMet_j0",title,70,0.,7.,120,0.,360.);
  fHist.h1_Etee_PhiMet_j0 ->GetXaxis()->SetTitle(xlabel);
  fHist.h1_Etee_PhiMet_j0 ->GetYaxis()->SetTitle(ylabel);

  var="Met";
  xlabel="Met";
  ylabel="Ete1+Ete2";
  sprintf(title, "ee Ete1+Ete2 <100 vs Met: 0-jets ", var);
  fHist.h2_Etee_PhiMet_j0 = new TH2F("h2_Etee_PhiMet_j0",title,70,0.,7.,120,0.,360.);
  fHist.h2_Etee_PhiMet_j0 ->GetXaxis()->SetTitle(xlabel);
  fHist.h2_Etee_PhiMet_j0 ->GetYaxis()->SetTitle(ylabel);


  var="Ht";
  xlabel="Ht";
  ylabel="Ete1+Ete2+met";
  sprintf(title, "ee Ht vs Ete1+Ete2+met: 0-jets ", var);
  fHist.h_Ht_EteeMet_j0 = new TH2F("Ht_EteeMet_j0",title,120,0.,360.,120,0.,360.);
  fHist.h_Ht_EteeMet_j0 ->GetXaxis()->SetTitle(xlabel);
  fHist.h_Ht_EteeMet_j0 ->GetYaxis()->SetTitle(ylabel);

  var="Ht";
  xlabel="Etmin";
  ylabel="Etmax";
  sprintf(title, "ee Et max vs Et min: 0-jets ", var);
  fHist.h_Ete1_Ete2 = new TH2F("Ete1_Ete2",title,120,0.,360.,120,0.,360.);
  fHist.h_Ete1_Ete2 ->GetXaxis()->SetTitle(xlabel);
  fHist.h_Ete1_Ete2 ->GetYaxis()->SetTitle(ylabel);

  //WW->all events after all cuts.
 
  var="zmass";
  xlabel="zmass of ll";
  ylabel="no. of evnts";
  sprintf(title, "z mass of ww candidates ", var);
  fHist.h_zmass_wwall = new TH1F("zmass_wwall",title,40,0.,200.);
  fHist.h_zmass_wwall->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_wwall->GetYaxis()->SetTitle(ylabel);
 
  var="zmass";
  xlabel="zmass of ll";
  ylabel="no. of evnts";
  sprintf(title, "z mass of ww candidates BHT", var);
  fHist.h_zmass_wwallbht = new TH1F("zmass_wwallbht",title,40,0.,200.);
  fHist.h_zmass_wwallbht->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_wwallbht->GetYaxis()->SetTitle(ylabel);

  var="ptllPair";
  xlabel="pt of ll Pair";
  ylabel="no. of evnts";
  sprintf(title, "Pt(ll) of ww ee candidates ", var);
  fHist.h_ptllPair_wwall = new TH1F("ptllPair_wwall",title,120,0.,300.); 
  fHist.h_ptllPair_wwall->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptllPair_wwall->GetYaxis()->SetTitle(ylabel);
 
  var="pte1pte2_wwee";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of ww candidates ", var);
  fHist.h_pte1pte2_wwee = new TH2F("pte1pte2_wwee",title,200,0.,200.,150,0.,150.); 
  fHist.h_pte1pte2_wwee->GetXaxis()->SetTitle(xlabel);
  fHist.h_pte1pte2_wwee->GetYaxis()->SetTitle(ylabel);
 
  var="ptmu1ptmu2_wwmumu";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of ww mumu candidates ", var);
  fHist.h_ptmu1ptmu2_wwmumu = new TH2F("ptmu1ptmu2_wwmumu",title,200,0.,200.,150,0.,150.); 
  fHist.h_ptmu1ptmu2_wwmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptmu1ptmu2_wwmumu->GetYaxis()->SetTitle(ylabel);  

  var="pteptmu_wwemu";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of ww emu candidates ", var);
  fHist.h_pteptmu_wwemu = new TH2F("pteptmu_wwemu",title,200,0.,200.,150,0.,150.); 
  fHist.h_pteptmu_wwemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_pteptmu_wwemu->GetYaxis()->SetTitle(ylabel);

  
  var="pte1pte2_1jee";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of 1j candidates ", var);
  fHist.h_pte1pte2_1jee = new TH2F("pte1pte2_1jee",title,200,0.,200.,150,0.,150.); 
  fHist.h_pte1pte2_1jee->GetXaxis()->SetTitle(xlabel);
  fHist.h_pte1pte2_1jee->GetYaxis()->SetTitle(ylabel);
 
  var="ptmu1ptmu2_1jmumu";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of 1j mumu candidates ", var);
  fHist.h_ptmu1ptmu2_1jmumu = new TH2F("ptmu1ptmu2_1jmumu",title,200,0.,200.,150,0.,150.); 
  fHist.h_ptmu1ptmu2_1jmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptmu1ptmu2_1jmumu->GetYaxis()->SetTitle(ylabel);  

  var="pteptmu_1jemu";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of 1j emu candidates ", var);
  fHist.h_pteptmu_1jemu = new TH2F("pteptmu_1jemu",title,200,0.,200.,150,0.,150.); 
  fHist.h_pteptmu_1jemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_pteptmu_1jemu->GetYaxis()->SetTitle(ylabel);


  
  var="pte1pte2_ttee";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of tt candidates ", var);
  fHist.h_pte1pte2_ttee = new TH2F("pte1pte2_ttee",title,200,0.,200.,150,0.,150.); 
  fHist.h_pte1pte2_ttee->GetXaxis()->SetTitle(xlabel);
  fHist.h_pte1pte2_ttee->GetYaxis()->SetTitle(ylabel);
 
  var="ptmu1ptmu2_ttmumu";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of tt mumu candidates ", var);
  fHist.h_ptmu1ptmu2_ttmumu = new TH2F("ptmu1ptmu2_ttmumu",title,200,0.,200.,150,0.,150.); 
  fHist.h_ptmu1ptmu2_ttmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptmu1ptmu2_ttmumu->GetYaxis()->SetTitle(ylabel);  

  var="pteptmu_ttemu";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of tt emu candidates ", var);
  fHist.h_pteptmu_ttemu = new TH2F("pteptmu_ttemu",title,200,0.,200.,150,0.,150.); 
  fHist.h_pteptmu_ttemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_pteptmu_ttemu->GetYaxis()->SetTitle(ylabel);

  var="pte1pte2_ttee";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of tt candidates BHT", var);
  fHist.h_pte1pte2_tteebht = new TH2F("pte1pte2_tteebht",title,200,0.,200.,150,0.,150.); 
  fHist.h_pte1pte2_tteebht->GetXaxis()->SetTitle(xlabel);
  fHist.h_pte1pte2_tteebht->GetYaxis()->SetTitle(ylabel);
 
  var="ptmu1ptmu2_ttmumu";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of tt mumu candidates BHT", var);
  fHist.h_ptmu1ptmu2_ttmumubht = new TH2F("ptmu1ptmu2_ttmumubht",title,200,0.,200.,150,0.,150.); 
  fHist.h_ptmu1ptmu2_ttmumubht->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptmu1ptmu2_ttmumubht->GetYaxis()->SetTitle(ylabel);  

  var="pteptmu_ttemu";
  xlabel="pt of l1";
  ylabel="pt of l2";
  sprintf(title, "Pt(l1)xPt(l2) of tt emu candidates BHT", var);
  fHist.h_pteptmu_ttemubht = new TH2F("pteptmu_ttemubht",title,200,0.,200.,150,0.,150.); 
  fHist.h_pteptmu_ttemubht->GetXaxis()->SetTitle(xlabel);
  fHist.h_pteptmu_ttemubht->GetYaxis()->SetTitle(ylabel);


  var="ptl";
  xlabel="pt of lepton";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each lepton of ww candidates ", var);
  fHist.h_ptl_wwall = new TH1F("ptl_wwall",title,40,0.,200.);
  fHist.h_ptl_wwall->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptl_wwall->GetYaxis()->SetTitle(ylabel);
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of ww cands after all cuts ", var);
  fHist.h_met_wwall = new TH1F("met_wwall",title,30,0.,150.);
  fHist.h_met_wwall->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_wwall->GetYaxis()->SetTitle(ylabel);

  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of ww cands after all cuts BHT", var);
  fHist.h_met_wwallbht = new TH1F("met_wwallbht",title,30,0.,150.);
  fHist.h_met_wwallbht->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_wwallbht->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdee";
  xlabel="dphimetdee";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest l/j ww cands after all cuts ",
var);
  fHist.h_dphimetdee_wwall = new
TH1F("dphimetdee_wwall",title,72,0.,Pi);
  fHist.h_dphimetdee_wwall->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_wwall->GetYaxis()->SetTitle(ylabel);
 
  var="dphillmet";
  xlabel="dphillmet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met ll ww cands after all cuts ", var);
  fHist.h_dphillmet_wwall = new TH1F("dphillmet_wwall",title,72,0.,Pi);
  fHist.h_dphillmet_wwall->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphillmet_wwall->GetYaxis()->SetTitle(ylabel);

  //WW->ee only events after all cuts
 
 
  var="zmass";
  xlabel="zmass of ee";
  ylabel="no. of evnts";
  sprintf(title, "z mass of ww->evev candidates ", var);
  fHist.h_zmass_wwee = new TH1F("zmass_wwee",title,40,0.,200.);
  fHist.h_zmass_wwee->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_wwee->GetYaxis()->SetTitle(ylabel);

  var="zmass";
  xlabel="zmass of ee";
  ylabel="no. of evnts";
  sprintf(title, "z mass of ww->evev candidates BHT", var);
  fHist.h_zmass_wweebht = new TH1F("zmass_wweebht",title,40,0.,200.);
  fHist.h_zmass_wweebht->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_wweebht->GetYaxis()->SetTitle(ylabel);
 
  var="pteePair";
  xlabel="pteePair of ee";
  ylabel="no. of evnts";
  sprintf(title, "z mass of ww->evev candidates ", var);
  fHist.h_pteePair_wwee = new TH1F("pteePair_wwee",title,120,0.,300.);
  fHist.h_pteePair_wwee->GetXaxis()->SetTitle(xlabel);
  fHist.h_pteePair_wwee->GetYaxis()->SetTitle(ylabel);
 
var="pte";
  xlabel="pt of electron";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each electron of ww candidates ", var);
  fHist.h_pte_wwee = new TH1F("pte_wwee",title,100,0.,200.);
  fHist.h_pte_wwee->GetXaxis()->SetTitle(xlabel);
  fHist.h_pte_wwee->GetYaxis()->SetTitle(ylabel);
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of ww->evev cands after all cuts ", var);
  fHist.h_met_wwee = new TH1F("met_wwee",title,30,0.,150.);
  fHist.h_met_wwee->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_wwee->GetYaxis()->SetTitle(ylabel);
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of ww->evev cands after all cuts ", var);
  fHist.h_met_wweebht = new TH1F("met_wweebht",title,30,0.,150.);
  fHist.h_met_wweebht->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_wweebht->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdee";
  xlabel="dphimetdee";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest l/j ww->evev cands after all cuts ", var);
  fHist.h_dphimetdee_wwee = new TH1F("dphimetdee_wwee",title,72,0.,Pi);
  fHist.h_dphimetdee_wwee->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_wwee->GetYaxis()->SetTitle(ylabel);
 
  var="dphieemet";
  xlabel="dphieemet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met ee ww->evev cands after all cuts ",
var);
  fHist.h_dphieemet_wwee = new TH1F("dphieemet_wwee",title,72,0.,Pi);
  fHist.h_dphieemet_wwee->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphieemet_wwee->GetYaxis()->SetTitle(ylabel);

  //ww_mumu
 
  var="zmass";
  xlabel="zmass of mumu";
  ylabel="no. of evnts";
  sprintf(title, "z mass of ww->muvmuv candidates ", var);
  fHist.h_zmass_wwmumu = new TH1F("zmass_wwmumu",title,40,0.,200.);
  fHist.h_zmass_wwmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_wwmumu->GetYaxis()->SetTitle(ylabel);

  var="zmass";
  xlabel="zmass of mumu";
  ylabel="no. of evnts";
  sprintf(title, "z mass of ww->muvmuv candidates BHT", var);
  fHist.h_zmass_wwmumubht = new TH1F("zmass_wwmumubht",title,40,0.,200.);
  fHist.h_zmass_wwmumubht->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_wwmumubht->GetYaxis()->SetTitle(ylabel);
 
  var="ptmumuPair";
  xlabel="pt of mumu Pair";
  ylabel="no. of evnts";
  sprintf(title, "pt(mumu) of ww->muvmuv candidates ", var);
  fHist.h_ptmumuPair_wwmumu = new
TH1F("ptmumuPair_wwmumu",title,120,0.,300.);
  fHist.h_ptmumuPair_wwmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptmumuPair_wwmumu->GetYaxis()->SetTitle(ylabel);
 
  var="ptmu";
  xlabel="pt of muon";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each muon of ww candidates ", var);
  fHist.h_ptmu_wwmumu = new TH1F("ptmu_wwmumu",title,100,0.,200.);
  fHist.h_ptmu_wwmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptmu_wwmumu->GetYaxis()->SetTitle(ylabel);
 
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of ww->muvmuv cands after all cuts ", var);
  fHist.h_met_wwmumu = new TH1F("met_wwmumu",title,30,0.,150.);
  fHist.h_met_wwmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_wwmumu->GetYaxis()->SetTitle(ylabel);

  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of ww->muvmuv cands after all cuts BHT ", var);
  fHist.h_met_wwmumubht = new TH1F("met_wwmumubht",title,30,0.,150.);
  fHist.h_met_wwmumubht->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_wwmumubht->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdmumu";
  xlabel="dphimetdmumu";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest mu/j ww->muvmuv cands after all cuts ", var);
  fHist.h_dphimetdee_wwmumu = new
TH1F("dphimetdee_wwmumu",title,72,0.,Pi);
  fHist.h_dphimetdee_wwmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_wwmumu->GetYaxis()->SetTitle(ylabel);
 
var="dphimumumet";
  xlabel="dphimumumet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met mumu ww->muvmuv cands after all cuts ", var);
  fHist.h_dphimumumet_wwmumu = new
TH1F("dphimumumet_wwmumu",title,72,0.,Pi);
  fHist.h_dphimumumet_wwmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimumumet_wwmumu->GetYaxis()->SetTitle(ylabel);
  
  //WW emu only after all cuts:
 
  var="zmass";
  xlabel="zmass of emu";
  ylabel="no. of evnts";
  sprintf(title, "z mass of ww->evmuv candidates ", var);
  fHist.h_zmass_wwemu = new TH1F("zmass_wwemu",title,40,0.,200.);
  fHist.h_zmass_wwemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_wwemu->GetYaxis()->SetTitle(ylabel);

  var="zmass";
  xlabel="zmass of emu";
  ylabel="no. of evnts";
  sprintf(title, "z mass of ww->evmuv candidates BHT", var);
  fHist.h_zmass_wwemubht = new TH1F("zmass_wwemubht",title,40,0.,200.);
  fHist.h_zmass_wwemubht->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_wwemubht->GetYaxis()->SetTitle(ylabel);
 
  var="ptemuPair";
  xlabel="ptemuPair of ee";
  ylabel="no. of evnts";
  sprintf(title, "Pt(emu) of ww->evmuv candidates ", var);
  fHist.h_ptemuPair_wwemu = new
TH1F("ptemuPair_wwemu",title,120,0.,300.);
  fHist.h_ptemuPair_wwemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptemuPair_wwemu->GetYaxis()->SetTitle(ylabel);
 
  var="ptl";
  xlabel="pt of lepton";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each lepton of ww candidates ", var);
  fHist.h_ptl_wwemu = new TH1F("ptl_wwemu",title,100,0.,200.);
  fHist.h_ptl_wwemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptl_wwemu->GetYaxis()->SetTitle(ylabel);
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of ww->evmuv cands after all cuts ", var);
  fHist.h_met_wwemu = new TH1F("met_wwemu",title,30,0.,150.);
  fHist.h_met_wwemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_wwemu->GetYaxis()->SetTitle(ylabel);

  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of ww->evmuv cands after all cuts BHT", var);
  fHist.h_met_wwemubht = new TH1F("met_wwemubht",title,30,0.,150.);
  fHist.h_met_wwemubht->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_wwemubht->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdemu";
  xlabel="dphimetdemu";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest l/j ww->evmuv cands after all cuts ", var);
  fHist.h_dphimetdee_wwemu = new
TH1F("dphimetdee_wwemu",title,72,0.,Pi);
  fHist.h_dphimetdee_wwemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_wwemu->GetYaxis()->SetTitle(ylabel);
 
  var="dphiemumet";
  xlabel="dphiemumet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met emu ww->evmuv cands after all cuts ",
var);
  fHist.h_dphiemumet_wwemu = new
TH1F("dphiemumet_wwemu",title,72,0.,Pi);
  fHist.h_dphiemumet_wwemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphiemumet_wwemu->GetYaxis()->SetTitle(ylabel);

  //ttbar->all events after all cuts.
 
  var="zmass";
  xlabel="zmass of ll";
  ylabel="no. of evnts";
  sprintf(title, "z mass of tt candidates ", var);
  fHist.h_zmass_ttall = new TH1F("zmass_ttall",title,40,0.,200.);
  fHist.h_zmass_ttall->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_ttall->GetYaxis()->SetTitle(ylabel);

  var="zmass";
  xlabel="zmass of ll";
  ylabel="no. of evnts";
  sprintf(title, "z mass of tt candidates BHT", var);
  fHist.h_zmass_ttallbht = new TH1F("zmass_ttallbht",title,40,0.,200.);
  fHist.h_zmass_ttallbht->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_ttallbht->GetYaxis()->SetTitle(ylabel); 

  var="ptllPair";
  xlabel="pt of ll Pair";
  ylabel="no. of evnts";
  sprintf(title, "Pt(ll) of tt candidates ", var);
  fHist.h_ptllPair_ttall = new TH1F("ptllPair_ttall",title,120,0.,300.);
 
  fHist.h_ptllPair_ttall->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptllPair_ttall->GetYaxis()->SetTitle(ylabel);
 
  var="ptl";
  xlabel="pt of lepton";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each lepton of tt candidates ", var);
  fHist.h_ptl_ttall = new TH1F("ptl_ttall",title,100,0.,200.);
  fHist.h_ptl_ttall->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptl_ttall->GetYaxis()->SetTitle(ylabel);
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of tt cands after all cuts ", var);
  fHist.h_met_ttall = new TH1F("met_ttall",title,30,0.,150.);
  fHist.h_met_ttall->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_ttall->GetYaxis()->SetTitle(ylabel);

  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of tt cands after all cuts BHT", var);
  fHist.h_met_ttallbht = new TH1F("met_ttallbht",title,30,0.,150.);
  fHist.h_met_ttallbht->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_ttallbht->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdee";
  xlabel="dphimetdee";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest l/j tt cands after all cuts ",
var);
  fHist.h_dphimetdee_ttall = new
TH1F("dphimetdee_ttall",title,72,0.,Pi);
  fHist.h_dphimetdee_ttall->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_ttall->GetYaxis()->SetTitle(ylabel);
 
  var="dphillmet";
  xlabel="dphillmet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met ll tt cands after all cuts ", var);
  fHist.h_dphillmet_ttall = new TH1F("dphillmet_ttall",title,72,0.,Pi);
  fHist.h_dphillmet_ttall->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphillmet_ttall->GetYaxis()->SetTitle(ylabel);

  //TT->ee only events after all cuts
 
 
  var="zmass";
  xlabel="zmass of ee";
  ylabel="no. of evnts";
  sprintf(title, "z mass of tt->evev candidates ", var);
  fHist.h_zmass_ttee = new TH1F("zmass_ttee",title,40,0.,200.);
  fHist.h_zmass_ttee->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_ttee->GetYaxis()->SetTitle(ylabel);

   var="zmass";
  xlabel="zmass of ee";
  ylabel="no. of evnts";
  sprintf(title, "z mass of tt->evev candidates BHT", var);
  fHist.h_zmass_tteebht = new TH1F("zmass_tteebht",title,40,0.,200.);
  fHist.h_zmass_tteebht->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_tteebht->GetYaxis()->SetTitle(ylabel);
 
  var="pteePair";
  xlabel="pteePair of ee";
  ylabel="no. of evnts";
  sprintf(title, "z mass of tt->evev candidates ", var);
  fHist.h_pteePair_ttee = new TH1F("pteePair_ttee",title,120,0.,300.);
  fHist.h_pteePair_ttee->GetXaxis()->SetTitle(xlabel);
  fHist.h_pteePair_ttee->GetYaxis()->SetTitle(ylabel);
 
  var="pte";
  xlabel="pt of electron";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each electron of tt candidates ", var);
  fHist.h_pte_ttee = new TH1F("pte_ttee",title,100,0.,200.);
  fHist.h_pte_ttee->GetXaxis()->SetTitle(xlabel);
  fHist.h_pte_ttee->GetYaxis()->SetTitle(ylabel);
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of tt->evev cands after all cuts ", var);
  fHist.h_met_ttee = new TH1F("met_ttee",title,30,0.,150.);
  fHist.h_met_ttee->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_ttee->GetYaxis()->SetTitle(ylabel);

  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of tt->evev cands after all cuts BHT", var);
  fHist.h_met_tteebht = new TH1F("met_tteebht",title,30,0.,150.);
  fHist.h_met_tteebht->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_tteebht->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdee";
  xlabel="dphimetdee";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest l/j tt->evev cands after all cuts ", var);
  fHist.h_dphimetdee_ttee = new TH1F("dphimetdee_ttee",title,72,0.,Pi);
  fHist.h_dphimetdee_ttee->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_ttee->GetYaxis()->SetTitle(ylabel);
 
  var="dphieemet";
  xlabel="dphieemet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met ee tt->evev cands after all cuts ",
var);
  fHist.h_dphieemet_ttee = new TH1F("dphieemet_ttee",title,72,0.,Pi);
  fHist.h_dphieemet_ttee->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphieemet_ttee->GetYaxis()->SetTitle(ylabel);

  //tt_mumu
 
  var="zmass";
  xlabel="zmass of mumu";
  ylabel="no. of evnts";
  sprintf(title, "z mass of tt->muvmuv candidates ", var);
  fHist.h_zmass_ttmumu = new TH1F("zmass_ttmumu",title,40,0.,200.);
  fHist.h_zmass_ttmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_ttmumu->GetYaxis()->SetTitle(ylabel);
 
  var="zmass";
  xlabel="zmass of mumu";
  ylabel="no. of evnts";
  sprintf(title, "z mass of tt->muvmuv candidates BHT", var);
  fHist.h_zmass_ttmumubht = new TH1F("zmass_ttmumubht",title,40,0.,200.);
  fHist.h_zmass_ttmumubht->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_ttmumubht->GetYaxis()->SetTitle(ylabel);
 
  var="ptmumuPair";
  xlabel="pt of mumu Pair";
  ylabel="no. of evnts";
  sprintf(title, "pt(mumu) of tt->muvmuv candidates ", var);
  fHist.h_ptmumuPair_ttmumu = new
TH1F("ptmumuPair_ttmumu",title,120,0.,300.);
  fHist.h_ptmumuPair_ttmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptmumuPair_ttmumu->GetYaxis()->SetTitle(ylabel);
 
  var="ptmu";
  xlabel="pt of muon";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each muon of tt candidates ", var);
  fHist.h_ptmu_ttmumu = new TH1F("ptmu_ttmumu",title,100,0.,200.);
  fHist.h_ptmu_ttmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptmu_ttmumu->GetYaxis()->SetTitle(ylabel);
 
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of tt->muvmuv cands after all cuts ", var);
  fHist.h_met_ttmumu = new TH1F("met_ttmumu",title,30,0.,150.);
  fHist.h_met_ttmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_ttmumu->GetYaxis()->SetTitle(ylabel);

   var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of tt->muvmuv cands after all cuts BHT", var);
  fHist.h_met_ttmumubht = new TH1F("met_ttmumubht",title,30,0.,150.);
  fHist.h_met_ttmumubht->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_ttmumubht->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdmumu";
  xlabel="dphimetdmumu";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest mu/j tt->muvmuv cands after all cuts ", var);
  fHist.h_dphimetdee_ttmumu = new
TH1F("dphimetdee_ttmumu",title,72,0.,Pi);
  fHist.h_dphimetdee_ttmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_ttmumu->GetYaxis()->SetTitle(ylabel);
 
var="dphimumumet";
  xlabel="dphimumumet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met mumu tt->muvmuv cands after all cuts ", var);
  fHist.h_dphimumumet_ttmumu = new
TH1F("dphimumumet_ttmumu",title,72,0.,Pi);
  fHist.h_dphimumumet_ttmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimumumet_ttmumu->GetYaxis()->SetTitle(ylabel);
  
  //TT emu only after all cuts:
 
  var="zmass";
  xlabel="zmass of emu";
  ylabel="no. of evnts";
  sprintf(title, "z mass of tt->evmuv candidates ", var);
  fHist.h_zmass_ttemu = new TH1F("zmass_ttemu",title,40,0.,200.);
  fHist.h_zmass_ttemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_ttemu->GetYaxis()->SetTitle(ylabel);
 
  var="zmass";
  xlabel="zmass of emu";
  ylabel="no. of evnts";
  sprintf(title, "z mass of tt->evmuv candidates BHT", var);
  fHist.h_zmass_ttemubht = new TH1F("zmass_ttemubht",title,40,0.,200.);
  fHist.h_zmass_ttemubht->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_ttemubht->GetYaxis()->SetTitle(ylabel);

  var="ptemuPair";
  xlabel="ptemuPair of ee";
  ylabel="no. of evnts";
  sprintf(title, "Pt(emu) of tt->evmuv candidates ", var);
  fHist.h_ptemuPair_ttemu = new
TH1F("ptemuPair_ttemu",title,120,0.,300.);
  fHist.h_ptemuPair_ttemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptemuPair_ttemu->GetYaxis()->SetTitle(ylabel);
 
  var="ptl";
  xlabel="pt of lepton";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each lepton of tt candidates ", var);
  fHist.h_ptl_ttemu = new TH1F("ptl_ttemu",title,100,0.,200.);
  fHist.h_ptl_ttemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptl_ttemu->GetYaxis()->SetTitle(ylabel);
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of tt->evmuv cands after all cuts ", var);
  fHist.h_met_ttemu = new TH1F("met_ttemu",title,30,0.,150.);
  fHist.h_met_ttemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_ttemu->GetYaxis()->SetTitle(ylabel);

  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of tt->evmuv cands after all cuts BHT", var);
  fHist.h_met_ttemubht = new TH1F("met_ttemubht",title,30,0.,150.);
  fHist.h_met_ttemubht->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_ttemubht->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdemu";
  xlabel="dphimetdemu";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest l/j tt->evmuv cands after all cuts ", var);
  fHist.h_dphimetdee_ttemu = new
TH1F("dphimetdee_ttemu",title,72,0.,Pi);
  fHist.h_dphimetdee_ttemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_ttemu->GetYaxis()->SetTitle(ylabel);
 
  var="dphiemumet";
  xlabel="dphiemumet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met emu tt->evmuv cands after all cuts ",var);
  fHist.h_dphiemumet_ttemu = new
TH1F("dphiemumet_ttemu",title,72,0.,Pi);
  fHist.h_dphiemumet_ttemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphiemumet_ttemu->GetYaxis()->SetTitle(ylabel);

  //1J->all events after all cuts.
 
  var="zmass";
  xlabel="zmass of ll";
  ylabel="no. of evnts";
  sprintf(title, "z mass of 1j candidates ", var);
  fHist.h_zmass_1jall = new TH1F("zmass_1jall",title,40,0.,200.);
  fHist.h_zmass_1jall->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_1jall->GetYaxis()->SetTitle(ylabel);
 
  var="ptllPair";
  xlabel="pt of ll Pair";
  ylabel="no. of evnts";
  sprintf(title, "Pt(ll) of 1j candidates ", var);
  fHist.h_ptllPair_1jall = new TH1F("ptllPair_1jall",title,120,0.,300.);
 
  fHist.h_ptllPair_1jall->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptllPair_1jall->GetYaxis()->SetTitle(ylabel);
 
  var="ptl";
  xlabel="pt of lepton";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each lepton of 1j candidates ", var);
  fHist.h_ptl_1jall = new TH1F("ptl_1jall",title,100,0.,200.);
  fHist.h_ptl_1jall->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptl_1jall->GetYaxis()->SetTitle(ylabel);
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of 1j cands after all cuts ", var);
  fHist.h_met_1jall = new TH1F("met_1jall",title,120,0.,300.);
  fHist.h_met_1jall->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_1jall->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdee";
  xlabel="dphimetdee";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest l/j 1j cands after all cuts ",
var);
  fHist.h_dphimetdee_1jall = new
TH1F("dphimetdee_1jall",title,72,0.,Pi);
  fHist.h_dphimetdee_1jall->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_1jall->GetYaxis()->SetTitle(ylabel);
 
  var="dphillmet";
  xlabel="dphillmet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met ll 1j cands after all cuts ", var);
  fHist.h_dphillmet_1jall = new TH1F("dphillmet_1jall",title,72,0.,Pi);
  fHist.h_dphillmet_1jall->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphillmet_1jall->GetYaxis()->SetTitle(ylabel);

  //1J->ee only events after all cuts
 
 
  var="zmass";
  xlabel="zmass of ee";
  ylabel="no. of evnts";
  sprintf(title, "z mass of 1j->evev candidates ", var);
  fHist.h_zmass_1jee = new TH1F("zmass_1jee",title,40,0.,200.);
  fHist.h_zmass_1jee->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_1jee->GetYaxis()->SetTitle(ylabel);
 
  var="pteePair";
  xlabel="pteePair of ee";
  ylabel="no. of evnts";
  sprintf(title, "z mass of 1j->evev candidates ", var);
  fHist.h_pteePair_1jee = new TH1F("pteePair_1jee",title,120,0.,300.);
  fHist.h_pteePair_1jee->GetXaxis()->SetTitle(xlabel);
  fHist.h_pteePair_1jee->GetYaxis()->SetTitle(ylabel);
 
var="pte";
  xlabel="pt of electron";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each electron of 1j candidates ", var);
  fHist.h_pte_1jee = new TH1F("pte_1jee",title,100,0.,200.);
  fHist.h_pte_1jee->GetXaxis()->SetTitle(xlabel);
  fHist.h_pte_1jee->GetYaxis()->SetTitle(ylabel);
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of 1j->evev cands after all cuts ", var);
  fHist.h_met_1jee = new TH1F("met_1jee",title,120,0.,300.);
  fHist.h_met_1jee->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_1jee->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdee";
  xlabel="dphimetdee";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest l/j 1j->evev cands after all cuts ", var);
  fHist.h_dphimetdee_1jee = new TH1F("dphimetdee_1jee",title,72,0.,Pi);
  fHist.h_dphimetdee_1jee->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_1jee->GetYaxis()->SetTitle(ylabel);
 
  var="dphieemet";
  xlabel="dphieemet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met ee 1j->evev cands after all cuts ",
var);
  fHist.h_dphieemet_1jee = new TH1F("dphieemet_1jee",title,72,0.,Pi);
  fHist.h_dphieemet_1jee->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphieemet_1jee->GetYaxis()->SetTitle(ylabel);

  //1j_mumu
 
  var="zmass";
  xlabel="zmass of mumu";
  ylabel="no. of evnts";
  sprintf(title, "z mass of 1j->muvmuv candidates ", var);
  fHist.h_zmass_1jmumu = new TH1F("zmass_1jmumu",title,40,0.,200.);
  fHist.h_zmass_1jmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_1jmumu->GetYaxis()->SetTitle(ylabel);
 
  var="ptmumuPair";
  xlabel="pt of mumu Pair";
  ylabel="no. of evnts";
  sprintf(title, "pt(mumu) of 1j->muvmuv candidates ", var);
  fHist.h_ptmumuPair_1jmumu = new
TH1F("ptmumuPair_1jmumu",title,120,0.,300.);
  fHist.h_ptmumuPair_1jmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptmumuPair_1jmumu->GetYaxis()->SetTitle(ylabel);
 
  var="ptmu";
  xlabel="pt of muon";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each muon of 1j candidates ", var);
  fHist.h_ptmu_1jmumu = new TH1F("ptmu_1jmumu",title,100,0.,200.);
  fHist.h_ptmu_1jmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptmu_1jmumu->GetYaxis()->SetTitle(ylabel);
 
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of 1j->muvmuv cands after all cuts ", var);
  fHist.h_met_1jmumu = new TH1F("met_1jmumu",title,120,0.,300.);
  fHist.h_met_1jmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_1jmumu->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdmumu";
  xlabel="dphimetdmumu";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest mu/j 1j->muvmuv cands after all cuts ", var);
  fHist.h_dphimetdee_1jmumu = new
TH1F("dphimetdee_1jmumu",title,72,0.,Pi);
  fHist.h_dphimetdee_1jmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_1jmumu->GetYaxis()->SetTitle(ylabel);
 
var="dphimumumet";
  xlabel="dphimumumet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met mumu 1j->muvmuv cands after all cuts ", var);
  fHist.h_dphimumumet_1jmumu = new
TH1F("dphimumumet_1jmumu",title,72,0.,Pi);
  fHist.h_dphimumumet_1jmumu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimumumet_1jmumu->GetYaxis()->SetTitle(ylabel);
  
  //1J emu only after all cuts:
 
  var="zmass";
  xlabel="zmass of emu";
  ylabel="no. of evnts";
  sprintf(title, "z mass of 1j->evmuv candidates ", var);
  fHist.h_zmass_1jemu = new TH1F("zmass_1jemu",title,40,0.,200.);
  fHist.h_zmass_1jemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_zmass_1jemu->GetYaxis()->SetTitle(ylabel);
 
  var="ptemuPair";
  xlabel="ptemuPair of ee";
  ylabel="no. of evnts";
  sprintf(title, "Pt(emu) of 1j->evmuv candidates ", var);
  fHist.h_ptemuPair_1jemu = new
TH1F("ptemuPair_1jemu",title,120,0.,300.);
  fHist.h_ptemuPair_1jemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptemuPair_1jemu->GetYaxis()->SetTitle(ylabel);
 
  var="ptl";
  xlabel="pt of lepton";
  ylabel="no. of evnts";
  sprintf(title, "Pt of each lepton of 1j candidates ", var);
  fHist.h_ptl_1jemu = new TH1F("ptl_1jemu",title,100,0.,200.);
  fHist.h_ptl_1jemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_ptl_1jemu->GetYaxis()->SetTitle(ylabel);
 
  var="met";
  xlabel="met";
  ylabel="no. of evnts";
  sprintf(title, "met of 1j->evmuv cands after all cuts ", var);
  fHist.h_met_1jemu = new TH1F("met_1jemu",title,120,0.,300.);
  fHist.h_met_1jemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_met_1jemu->GetYaxis()->SetTitle(ylabel);
 
  var="dphimetdemu";
  xlabel="dphimetdemu";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met nearest l/j 1j->evmuv cands after all cuts ", var);
  fHist.h_dphimetdee_1jemu = new
TH1F("dphimetdee_1jemu",title,72,0.,Pi);
  fHist.h_dphimetdee_1jemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphimetdee_1jemu->GetYaxis()->SetTitle(ylabel);
 
  var="dphiemumet";
  xlabel="dphiemumet";
  ylabel="no. of evnts";
  sprintf(title, "delta-phi met emu 1j->evmuv cands after all cuts ",
var);
  fHist.h_dphiemumet_1jemu = new
TH1F("dphiemumet_1jemu",title,72,0.,Pi);
  fHist.h_dphiemumet_1jemu->GetXaxis()->SetTitle(xlabel);
  fHist.h_dphiemumet_1jemu->GetYaxis()->SetTitle(ylabel);


   fHist.hepg_im34    = new TH1F("im34","im34",100,0.,200.);  
   fHist.hepg_im56    = new TH1F("im56","im56",100,0.,200.);  
   fHist.hepg_pt34    = new TH1F("pt34","pt34",50,0.,200.);  
   fHist.hepg_pt56    = new TH1F("pt56","pt56",50,0.,200.);
  //
   fHist.hepg_imw   = new TH1F("imw","imw",100,0.,200.); 
   fHist.hepg_ptw   = new TH1F("ptw","ptw",100,0.,200.); 
  //
   fHist.hepg_im3456   = new TH1F("im3456","im3456",250,0.,1000.);  
   fHist.hepg_pt3456   = new TH1F("pt3456","pt3456",100,0.,200.);
  //
   fHist.hepg_Pt3    = new TH1F("hPt3","pt3",50,0.,200.);  
   fHist.hepg_Pt4    = new TH1F("hPt4","pt4",50,0.,200.);  
   fHist.hepg_Pt5    = new TH1F("hPt5","pt5",50,0.,200.);  
   fHist.hepg_Pt6    = new TH1F("hPt6","pt6",50,0.,200.);
   fHist.hepg_Ptmin   = new TH1F("hPtmin","ptmin",50,0.,200.);
   fHist.hepg_Ptmax   = new TH1F("hPtmax","ptmax",50,0.,200.);
  //
   fHist.hepg_E3    = new TH1F("hE3","E3",50,0.,200.);  
   fHist.hepg_E4    = new TH1F("hE4","E4",50,0.,200.);  
   fHist.hepg_E5    = new TH1F("hE5","E5",50,0.,200.);  
   fHist.hepg_E6    = new TH1F("hE6","E6",50,0.,200.);
  //
   fHist.hepg_Ptl    = new TH1F("hPtl","ptl",50,0.,200.);
   fHist.hepg_Ptnunu = new TH1F("hPtnunu","ptnunu",50,0.,200.);
  //
  //
   fHist.hepg_dR35 = new TH1F("dR35","dR35",50,0.,4.);
   fHist.hepg_dphi35 = new TH1F("dphi35","dphi35",80,0.,4.);
   fHist.hepg_dphi565 = new TH1F("dphi565","dphi565",30,0.,4.);
   fHist.hepg_eta56arest = new TH1F("eta56arest","eta56arest",40,-4.,4.);
}


void histoLeptonId::WriteHistoFile(){

   TFile *fout = new TFile("wwmadgraph_unki.root","RECREATE");
   fout->cd();
   //
   fHist.he_n->Write();
   fHist.he_rawtotale->Write();
   fHist.he_rawtotalet->Write();
   fHist.he_corrtotalet->Write();
   //
   fHist.hnelec->Write();
   fHist.hnmuon->Write();
   //
   cout<<"111111"<<endl;
   //
    fHist.hbe_et->Write();
    fHist.hbe_pt->Write();
    fHist.hbe_eop->Write();
    fHist.hbe_z0->Write();
    fHist.hbe_hadoem->Write();
    fHist.hbe_lshr->Write();
    fHist.hbe_csdeltax->Write();
    fHist.hbe_deltaz->Write();
    fHist.hbe_chi2->Write();
    fHist.hbe_fcaliso->Write();
    fHist.hbe_trkphi->Write();
    fHist.hbe_trketa->Write();

    fHist.hbec_et->Write();
    fHist.hbec_pt->Write();
    fHist.hbec_eop->Write();
    fHist.hbec_z0->Write();
    fHist.hbec_hadoem->Write();
    fHist.hbec_lshr->Write();
    fHist.hbec_csdeltax->Write();
    fHist.hbec_deltaz->Write();
    fHist.hbec_chi2->Write();
    fHist.hbec_fcaliso->Write();
    fHist.hbec_trkphi->Write();
    fHist.hbec_trketa->Write();
    //
    cout<<"2222222"<<endl;
    //
    fHist.hbm_pt->Write();
    fHist.hbm_z0->Write();
    fHist.hbm_em->Write();
    fHist.hbm_had->Write();
    fHist.hbm_d0ws->Write();
    fHist.hbm_d0wos->Write();
    fHist.hbm_fciso->Write();
    fHist.hbm_deltax_cmu->Write();
    fHist.hbm_deltax_cmp->Write();
    fHist.hbm_deltax_cmx->Write();
    fHist.hbm_trketa_cmup->Write();
    fHist.hbm_trketa_cmu->Write();
    fHist.hbm_trketa_cmp->Write();
    fHist.hbm_trketa_cmx->Write();
    fHist.hbm_trkphi_cmup->Write();
    fHist.hbm_trkphi_cmu->Write();
    fHist.hbm_trkphi_cmp->Write();
    fHist.hbm_trkphi_cmx->Write();

    fHist.hbmc_pt->Write();
    fHist.hbmc_z0->Write();
    fHist.hbmc_em->Write();
    fHist.hbmc_had->Write();
    fHist.hbmc_d0ws->Write();
    fHist.hbmc_d0wos->Write();
    fHist.hbmc_fciso->Write();
    fHist.hbmc_deltax_cmu->Write();
    fHist.hbmc_deltax_cmp->Write();
    fHist.hbmc_deltax_cmx->Write();
    fHist.hbmc_trketa_cmup->Write();
    fHist.hbmc_trketa_cmu->Write();
    fHist.hbmc_trketa_cmp->Write();
    fHist.hbmc_trketa_cmx->Write();
    fHist.hbmc_trkphi_cmup->Write();
    fHist.hbmc_trkphi_cmu->Write();
    fHist.hbmc_trkphi_cmp->Write();
    fHist.hbmc_trkphi_cmx->Write();
   //
    cout<<"333333"<<endl;
    //
    fHist.hjets_njets->Write();
    fHist.hjets_et1->Write();
    fHist.hjets_eta1->Write();
    fHist.hjets_phi1->Write();
    fHist.hjets_et2->Write();
    fHist.hjets_eta2->Write();
    fHist.hjets_phi2->Write();
    fHist.hjets_et3->Write();
    fHist.hjets_eta3->Write();
    fHist.hjets_phi3->Write();
    //
    fHist.hdilepflag->Write();
    //
    fHist.hmee->Write();
    fHist.hmmumu->Write();
    fHist.hmee0->Write();
    fHist.hmmumu0->Write();
    //
    cout<<"444444"<<endl;
    //
    fHist.hmet_j0_ee->Write();
    fHist.hmet_j1_ee->Write();
    fHist.hmet_j2_ee->Write();
    fHist.hmet_j0_emu->Write();
    fHist.hmet_j1_emu->Write();
    fHist.hmet_j2_emu->Write();
    fHist.hmet_j0_mumu->Write();
    fHist.hmet_j1_mumu->Write();
    fHist.hmet_j2_mumu->Write();
    //
    fHist.heeDphimljMetj->Write();
    fHist.heeDphimljMetj0->Write();
    fHist.heeDphimljMetj1->Write();
    fHist.heeDphimljMetj2->Write();
    //
    fHist.hemuDphimljMetj->Write();
    fHist.hemuDphimljMetj0->Write();
    fHist.hemuDphimljMetj1->Write();
    fHist.hemuDphimljMetj2->Write();
    //
    fHist.hmumuDphimljMetj->Write();
    fHist.hmumuDphimljMetj0->Write();
    fHist.hmumuDphimljMetj1->Write();
    fHist.hmumuDphimljMetj2->Write();
    //
    fHist.hhteeDphimljMetj->Write();
    fHist.hhteeDphimljMetj0->Write();
    fHist.hhteeDphimljMetj1->Write();
    fHist.hhteeDphimljMetj2->Write();
    //
    fHist.hhtemuDphimljMetj->Write();
    fHist.hhtemuDphimljMetj0->Write();
    fHist.hhtemuDphimljMetj1->Write();
    fHist.hhtemuDphimljMetj2->Write();
    //
    fHist.hhtmumuDphimljMetj->Write();
    fHist.hhtmumuDphimljMetj0->Write();
    fHist.hhtmumuDphimljMetj1->Write();
    fHist.hhtmumuDphimljMetj2->Write();
    //
    fHist.hqeeDphimljMetj->Write();
    fHist.hqeeDphimljMetj0->Write();
    fHist.hqeeDphimljMetj1->Write();
    fHist.hqeeDphimljMetj2->Write();
    fHist.hqeeDphillj2->Write();
    fHist.hqeeDphillj0->Write();
    //
    cout<<"55555555"<<endl;
    
    fHist.hqemuDphimljMetj->Write();
    fHist.hqemuDphimljMetj0->Write();
    fHist.hqemuDphimljMetj1->Write();
    fHist.hqemuDphimljMetj2->Write();
    fHist.hqemuDphillj2->Write();
    fHist.hqemuDphillj0->Write();
    //
    fHist.hqmumuDphimljMetj->Write();
    fHist.hqmumuDphimljMetj0->Write();
    fHist.hqmumuDphimljMetj1->Write();
    fHist.hqmumuDphimljMetj2->Write();
    fHist.hqmumuDphillj2->Write();
    fHist.hqmumuDphillj0->Write();
    //
    fHist.nothhteeDphimljMetj->Write();
    fHist.nothhteeDphimljMetj0->Write();
    fHist.nothhteeDphimljMetj1->Write();
    fHist.nothhteeDphimljMetj2->Write();
    //
    fHist.nothhtemuDphimljMetj->Write();
    fHist.nothhtemuDphimljMetj0->Write();
    fHist.nothhtemuDphimljMetj1->Write();
    //
    fHist.nothhtmumuDphimljMetj->Write();
    fHist.nothhtmumuDphimljMetj0->Write();
    fHist.nothhtmumuDphimljMetj1->Write();
    fHist.nothhtmumuDphimljMetj2->Write();
    //
    fHist.nothqeeDphimljMetj->Write();
    fHist.nothqeeDphimljMetj0->Write();
    fHist.nothqeeDphimljMetj1->Write();
    fHist.nothqeeDphimljMetj2->Write();
    fHist.nothqeeDphillj2->Write();
    fHist.nothqeeDphillj0->Write();
    //
    fHist.nothqemuDphimljMetj->Write();
    fHist.nothqemuDphimljMetj0->Write();
    fHist.nothqemuDphimljMetj1->Write();
    fHist.nothqemuDphimljMetj2->Write();
    fHist.nothqemuDphillj2->Write();
    fHist.nothqemuDphillj0->Write();
    //
    fHist.nothqmumuDphimljMetj->Write();
    fHist.nothqmumuDphimljMetj0->Write();
    fHist.nothqmumuDphimljMetj1->Write();
    fHist.nothqmumuDphimljMetj2->Write();
    fHist.nothqmumuDphillj2->Write();
    fHist.nothqmumuDphillj0->Write();
    //
    fHist.hdphijele->Write();
    fHist.hRjele->Write();
    //
    cout<<"6666666"<<endl;
    //
    fHist.hallHtj->Write();
    fHist.hallHtj0->Write();
    fHist.hallHtj1inc->Write();
    fHist.hallHtj2->Write();
    fHist.heeHt->Write();
    fHist.heeHtTow->Write();
    fHist.heeHtj0->Write();
    fHist.heeHtj1->Write();
    fHist.heeHtj1inc->Write();
    fHist.heeHtj2->Write();
    fHist.hmumuHt->Write();
    fHist.hmumuHtTow->Write();
    fHist.hmumuHtj0->Write();
    fHist.hmumuHtj1->Write();
    fHist.hmumuHtj1inc->Write();
    fHist.hmumuHtj2->Write();
    fHist.hemuHt->Write();
    fHist.hemuHtTow->Write();
    fHist.hemuHtj0->Write();
    fHist.hemuHtj1->Write();
    fHist.hemuHtj1inc->Write();
    fHist.hemuHtj2->Write();
    //
    fHist.h_Ete_j0met20->Write();
    fHist.h_Ht_j0met20->Write();
    fHist.h_Htbis_j0met20->Write();
    fHist.h_dRlj->Write();
    //
    fHist.h_Ht_Met_j0->Write();
    fHist.h_Etee_Met_j0->Write();
    fHist.h1_Etee_PhiMet_j0->Write();
    fHist.h2_Etee_PhiMet_j0->Write();
    fHist.h_Ht_EteeMet_j0->Write();
    //
    fHist.h_Ete1_Ete2->Write();
    //
    fHist.h_dphimetdee_wwall->Write();
    fHist.h_dphillmet_wwall->Write();
    fHist.h_met_wwall->Write();
    fHist.h_zmass_wwall->Write();
    fHist.h_met_wwallbht->Write();
    fHist.h_zmass_wwallbht->Write();
    fHist.h_ptllPair_wwall->Write();
    fHist.h_ptl_wwall->Write();
    //ee
    fHist.h_dphimetdee_wwee->Write();
    fHist.h_dphieemet_wwee->Write();
    fHist.h_met_wwee->Write();
    fHist.h_zmass_wweebht->Write();
    fHist.h_met_wweebht->Write();
    fHist.h_zmass_wwee->Write();
    fHist.h_pteePair_wwee->Write();
    fHist.h_pte_wwee->Write();
    //mumu
    fHist.h_dphimetdee_wwmumu->Write();
    fHist.h_dphimumumet_wwmumu->Write();
    fHist.h_met_wwmumu->Write();
    fHist.h_zmass_wwmumu->Write();
    fHist.h_met_wwmumubht->Write();
    fHist.h_zmass_wwmumubht->Write();
    fHist.h_ptmumuPair_wwmumu->Write();
    fHist.h_ptmu_wwmumu->Write();
    //emu
    fHist.h_dphimetdee_wwemu->Write();
    fHist.h_dphiemumet_wwemu->Write();
    fHist.h_met_wwemu->Write();
    fHist.h_zmass_wwemu->Write();
    fHist.h_met_wwemubht->Write();
    fHist.h_zmass_wwemubht->Write();
    fHist.h_ptemuPair_wwemu->Write();
    fHist.h_ptl_wwemu->Write();
        
    cout<<"777777777"<<endl;   
 
    fHist.h_dphimetdee_ttall->Write();
    fHist.h_dphillmet_ttall->Write();
    fHist.h_met_ttall->Write();
    fHist.h_met_ttallbht->Write();
    fHist.h_zmass_ttall->Write();
    fHist.h_zmass_ttallbht->Write();
    fHist.h_ptllPair_ttall->Write();
    fHist.h_ptl_ttall->Write();
    //ee
    fHist.h_dphimetdee_ttee->Write();
    fHist.h_dphieemet_ttee->Write();
    fHist.h_met_ttee->Write();
    fHist.h_met_tteebht->Write();
    fHist.h_zmass_ttee->Write();
    fHist.h_zmass_tteebht->Write();
    fHist.h_pteePair_ttee->Write();
    fHist.h_pte_ttee->Write();
    //mumu
    fHist.h_dphimetdee_ttmumu->Write();
    fHist.h_dphimumumet_ttmumu->Write();
    fHist.h_met_ttmumu->Write();
    fHist.h_met_ttmumubht->Write();
    fHist.h_zmass_ttmumu->Write();
    fHist.h_zmass_ttmumubht->Write();
    fHist.h_ptmumuPair_ttmumu->Write();
    fHist.h_ptmu_ttmumu->Write();
    //emu
    fHist.h_dphimetdee_ttemu->Write();
    fHist.h_dphiemumet_ttemu->Write();
    fHist.h_met_ttemu->Write();
    fHist.h_met_ttemubht->Write();
    fHist.h_zmass_ttemu->Write();
    fHist.h_zmass_ttemubht->Write();
    fHist.h_ptemuPair_ttemu->Write();
    fHist.h_ptl_ttemu->Write();

    fHist.h_dphimetdee_1jall->Write();
    fHist.h_dphillmet_1jall->Write();
    fHist.h_met_1jall->Write();
    fHist.h_zmass_1jall->Write();
    fHist.h_ptllPair_1jall->Write();
    fHist.h_ptl_1jall->Write();
    //ee
    fHist.h_dphimetdee_1jee->Write();
    fHist.h_dphieemet_1jee->Write();
    fHist.h_met_1jee->Write();
    fHist.h_zmass_1jee->Write();
    fHist.h_pteePair_1jee->Write();
    fHist.h_pte_1jee->Write();
    //mumu
    fHist.h_dphimetdee_1jmumu->Write();
    fHist.h_dphimumumet_1jmumu->Write(); 
    fHist.h_met_1jmumu->Write();
    fHist.h_zmass_1jmumu->Write();
    fHist.h_ptmumuPair_1jmumu->Write();
    fHist.h_ptmu_1jmumu->Write();
    //emu
    cout<<"88888888"<<endl;

    fHist.h_dphimetdee_1jemu->Write();
    fHist.h_dphiemumet_1jemu->Write();
    fHist.h_met_1jemu->Write();
    fHist.h_zmass_1jemu->Write();
    fHist.h_ptemuPair_1jemu->Write();
    fHist.h_ptl_1jemu->Write();

    fHist.h_pte1pte2_wwee->Write();
    fHist.h_ptmu1ptmu2_wwmumu->Write();
    fHist.h_pteptmu_wwemu->Write();
    fHist.h_pte1pte2_1jee->Write();
    fHist.h_ptmu1ptmu2_1jmumu->Write();
    fHist.h_pteptmu_1jemu->Write();
    cout<<"9999999"<<endl; 
    fHist.h_pte1pte2_ttee->Write();
    fHist.h_ptmu1ptmu2_ttmumu->Write();
    fHist.h_pteptmu_ttemu->Write();
    fHist.h_pte1pte2_tteebht->Write();
    fHist.h_ptmu1ptmu2_ttmumubht->Write();
    fHist.h_pteptmu_ttemubht->Write();

    fHist.hepg_im34->Write();
    fHist.hepg_im56->Write();
    fHist.hepg_pt34->Write();
    fHist.hepg_pt56->Write();
    fHist.hepg_im3456->Write();
    fHist.hepg_pt3456->Write();
    fHist.hepg_imw->Write();
    fHist.hepg_ptw->Write();
    fHist.hepg_E3->Write();
    fHist.hepg_E4->Write();
    fHist.hepg_E5->Write();
    fHist.hepg_E6->Write();
    fHist.hepg_Pt3->Write();
    fHist.hepg_Pt4->Write();
    fHist.hepg_Pt5->Write();
    fHist.hepg_Pt6->Write();
    fHist.hepg_Ptl->Write();
    fHist.hepg_Ptnunu->Write();
    fHist.hepg_Ptl->Write();
    fHist.hepg_Ptmin->Write();
    fHist.hepg_Ptmax->Write();
    fHist.hepg_dR35->Write();
    fHist.hepg_dphi35->Write();   

    fHist.hepg_dphi565->Write(); 
    fHist.hepg_eta56arest->Write();     

    //Dean adds:
    fHist.hepg_im35->Write();	
    fHist.hepg_pt35->Write();

	
    fout->Write();

}