#include <fstream>
using namespace std;

void VLQvalid(){
  
  gStyle->SetOptStat(0);
  //gStyle->SetOptStat(111111);
  gStyle->SetOptTitle(0);


  // xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx

  // -----
  // Files
  // -----

  TFile *f_mT700  = new TFile("histos_mT700.root");
  TFile *f_mT800  = new TFile("histos_mT800.root");
  TFile *f_mT900  = new TFile("histos_mT900.root");

  TFile *f_mB700  = new TFile("histos_mB700.root");
  TFile *f_mB800  = new TFile("histos_mB800.root");
  TFile *f_mB900  = new TFile("histos_mB900.root");

  // Histos  
  TH1F* hT700_Q_mass  = (TH1F*)f_mT700->Get("h_Q_mass"); hT700_Q_mass->Scale(1./hT700_Q_mass->Integral());
  TH1F* hT800_Q_mass  = (TH1F*)f_mT800->Get("h_Q_mass"); hT800_Q_mass->Scale(1./hT800_Q_mass->Integral());
  TH1F* hT900_Q_mass  = (TH1F*)f_mT900->Get("h_Q_mass"); hT900_Q_mass->Scale(1./hT900_Q_mass->Integral());

  TH1F* hB700_Q_mass  = (TH1F*)f_mB700->Get("h_Q_mass"); hB700_Q_mass->Scale(1./hB700_Q_mass->Integral());
  TH1F* hB800_Q_mass  = (TH1F*)f_mB800->Get("h_Q_mass"); hB800_Q_mass->Scale(1./hB800_Q_mass->Integral());
  TH1F* hB900_Q_mass  = (TH1F*)f_mB900->Get("h_Q_mass"); hB900_Q_mass->Scale(1./hB900_Q_mass->Integral());

  TH1F* hT700_Q_dID   = (TH1F*)f_mT700->Get("h_Q_dID"); hT700_Q_dID->Scale(1./hT700_Q_dID->Integral());
  TH1F* hT800_Q_dID   = (TH1F*)f_mT800->Get("h_Q_dID"); hT800_Q_dID->Scale(1./hT800_Q_dID->Integral());
  TH1F* hT900_Q_dID   = (TH1F*)f_mT900->Get("h_Q_dID"); hT900_Q_dID->Scale(1./hT900_Q_dID->Integral());

  TH1F* hB700_Q_dID   = (TH1F*)f_mB700->Get("h_Q_dID"); hB700_Q_dID->Scale(1./hB700_Q_dID->Integral());
  TH1F* hB800_Q_dID   = (TH1F*)f_mB800->Get("h_Q_dID"); hB800_Q_dID->Scale(1./hB800_Q_dID->Integral());
  TH1F* hB900_Q_dID   = (TH1F*)f_mB900->Get("h_Q_dID"); hB900_Q_dID->Scale(1./hB900_Q_dID->Integral());

  //

  TH1F* hT700_Q_pt    = (TH1F*)f_mT700->Get("h_Q_pT"); hT700_Q_pt->Scale(1./hT700_Q_pt->Integral());
  TH1F* hT800_Q_pt    = (TH1F*)f_mT800->Get("h_Q_pT"); hT800_Q_pt->Scale(1./hT800_Q_pt->Integral());
  TH1F* hT900_Q_pt    = (TH1F*)f_mT900->Get("h_Q_pT"); hT900_Q_pt->Scale(1./hT900_Q_pt->Integral());

  TH1F* hB700_Q_pt    = (TH1F*)f_mB700->Get("h_Q_pT"); hB700_Q_pt->Scale(1./hB700_Q_pt->Integral());
  TH1F* hB800_Q_pt    = (TH1F*)f_mB800->Get("h_Q_pT"); hB800_Q_pt->Scale(1./hB800_Q_pt->Integral());
  TH1F* hB900_Q_pt    = (TH1F*)f_mB900->Get("h_Q_pT"); hB900_Q_pt->Scale(1./hB900_Q_pt->Integral());

  TH1F* hT700_Q_eta   = (TH1F*)f_mT700->Get("h_Q_eta"); hT700_Q_eta->Scale(1./hT700_Q_eta->Integral());
  TH1F* hT800_Q_eta   = (TH1F*)f_mT800->Get("h_Q_eta"); hT800_Q_eta->Scale(1./hT800_Q_eta->Integral());
  TH1F* hT900_Q_eta   = (TH1F*)f_mT900->Get("h_Q_eta"); hT900_Q_eta->Scale(1./hT900_Q_eta->Integral());

  TH1F* hB700_Q_eta   = (TH1F*)f_mB700->Get("h_Q_eta"); hB700_Q_eta->Scale(1./hB700_Q_eta->Integral());
  TH1F* hB800_Q_eta   = (TH1F*)f_mB800->Get("h_Q_eta"); hB800_Q_eta->Scale(1./hB800_Q_eta->Integral());
  TH1F* hB900_Q_eta   = (TH1F*)f_mB900->Get("h_Q_eta"); hB900_Q_eta->Scale(1./hB900_Q_eta->Integral());
  
  //

  TH1F* hT700_W_pt    = (TH1F*)f_mT700->Get("h_W_pT"); hT700_W_pt->Scale(1./hT700_W_pt->Integral());
  TH1F* hT800_W_pt    = (TH1F*)f_mT800->Get("h_W_pT"); hT800_W_pt->Scale(1./hT800_W_pt->Integral());
  TH1F* hT900_W_pt    = (TH1F*)f_mT900->Get("h_W_pT"); hT900_W_pt->Scale(1./hT900_W_pt->Integral());

  TH1F* hT800_Z_pt    = (TH1F*)f_mT800->Get("h_Z_pT"); hT800_Z_pt->Scale(1./hT800_Z_pt->Integral());
  TH1F* hT800_H_pt    = (TH1F*)f_mT800->Get("h_H_pT"); hT800_H_pt->Scale(1./hT800_H_pt->Integral());

  TH1F* hT700_W_eta   = (TH1F*)f_mT700->Get("h_W_eta"); hT700_W_eta->Scale(1./hT700_W_eta->Integral());
  TH1F* hT800_W_eta   = (TH1F*)f_mT800->Get("h_W_eta"); hT800_W_eta->Scale(1./hT800_W_eta->Integral());
  TH1F* hT900_W_eta   = (TH1F*)f_mT900->Get("h_W_eta"); hT900_W_eta->Scale(1./hT900_W_eta->Integral());

  TH1F* hT800_Z_eta   = (TH1F*)f_mT800->Get("h_Z_eta"); hT800_Z_eta->Scale(1./hT800_Z_eta->Integral());
  TH1F* hT800_H_eta   = (TH1F*)f_mT800->Get("h_H_eta"); hT800_H_eta->Scale(1./hT800_H_eta->Integral());

  //

  TH1F* hB700_W_pt    = (TH1F*)f_mB700->Get("h_W_pT"); hB700_W_pt->Scale(1./hB700_W_pt->Integral());
  TH1F* hB800_W_pt    = (TH1F*)f_mB800->Get("h_W_pT"); hB800_W_pt->Scale(1./hB800_W_pt->Integral());
  TH1F* hB900_W_pt    = (TH1F*)f_mB900->Get("h_W_pT"); hB900_W_pt->Scale(1./hB900_W_pt->Integral());

  TH1F* hB800_Z_pt    = (TH1F*)f_mB800->Get("h_Z_pT"); hB800_Z_pt->Scale(1./hB800_Z_pt->Integral());
  TH1F* hB800_H_pt    = (TH1F*)f_mB800->Get("h_H_pT"); hB800_H_pt->Scale(1./hB800_H_pt->Integral());

  TH1F* hB700_W_eta   = (TH1F*)f_mB700->Get("h_W_eta"); hB700_W_eta->Scale(1./hB700_W_eta->Integral());
  TH1F* hB800_W_eta   = (TH1F*)f_mB800->Get("h_W_eta"); hB800_W_eta->Scale(1./hB800_W_eta->Integral());
  TH1F* hB900_W_eta   = (TH1F*)f_mB900->Get("h_W_eta"); hB900_W_eta->Scale(1./hB900_W_eta->Integral());

  TH1F* hB800_Z_eta   = (TH1F*)f_mB800->Get("h_Z_eta"); hB800_Z_eta->Scale(1./hB800_Z_eta->Integral());
  TH1F* hB800_H_eta   = (TH1F*)f_mB800->Get("h_H_eta"); hB800_H_eta->Scale(1./hB800_H_eta->Integral());

  //

  TH1F* hT800_qW_pt   = (TH1F*)f_mT800->Get("h_qW_pT"); hT800_qW_pt->Scale(1./hT800_qW_pt->Integral());
  TH1F* hT800_qZ_pt   = (TH1F*)f_mT800->Get("h_qZ_pT"); hT800_qZ_pt->Scale(1./hT800_qZ_pt->Integral());
  TH1F* hT800_qH_pt   = (TH1F*)f_mT800->Get("h_qH_pT"); hT800_qH_pt->Scale(1./hT800_qH_pt->Integral());

  TH1F* hB800_qW_pt   = (TH1F*)f_mB800->Get("h_qW_pT"); hB800_qW_pt->Scale(1./hB800_qW_pt->Integral());
  TH1F* hB800_qZ_pt   = (TH1F*)f_mB800->Get("h_qZ_pT"); hB800_qZ_pt->Scale(1./hB800_qZ_pt->Integral());
  TH1F* hB800_qH_pt   = (TH1F*)f_mB800->Get("h_qH_pT"); hB800_qH_pt->Scale(1./hB800_qH_pt->Integral());

  TH1F* hT800_qW_eta  = (TH1F*)f_mT800->Get("h_qW_eta"); hT800_qW_eta->Scale(1./hT800_qW_eta->Integral());
  TH1F* hT800_qZ_eta  = (TH1F*)f_mT800->Get("h_qZ_eta"); hT800_qZ_eta->Scale(1./hT800_qZ_eta->Integral());
  TH1F* hT800_qH_eta  = (TH1F*)f_mT800->Get("h_qH_eta"); hT800_qH_eta->Scale(1./hT800_qH_eta->Integral());

  TH1F* hB800_qW_eta  = (TH1F*)f_mB800->Get("h_qW_eta"); hB800_qW_eta->Scale(1./hB800_qW_eta->Integral());
  TH1F* hB800_qZ_eta  = (TH1F*)f_mB800->Get("h_qZ_eta"); hB800_qZ_eta->Scale(1./hB800_qZ_eta->Integral());
  TH1F* hB800_qH_eta  = (TH1F*)f_mB800->Get("h_qH_eta"); hB800_qH_eta->Scale(1./hB800_qH_eta->Integral());

  //

  TH1F* hT700_evtype  = (TH1F*)f_mT700->Get("h_evtype"); hT700_evtype->Scale(1./hT700_evtype->Integral());
  TH1F* hT800_evtype  = (TH1F*)f_mT800->Get("h_evtype"); hT800_evtype->Scale(1./hT800_evtype->Integral());
  TH1F* hT900_evtype  = (TH1F*)f_mT900->Get("h_evtype"); hT900_evtype->Scale(1./hT900_evtype->Integral());

  TH1F* hB700_evtype  = (TH1F*)f_mB700->Get("h_evtype"); hB700_evtype->Scale(1./hB700_evtype->Integral());
  TH1F* hB800_evtype  = (TH1F*)f_mB800->Get("h_evtype"); hB800_evtype->Scale(1./hB800_evtype->Integral());
  TH1F* hB900_evtype  = (TH1F*)f_mB900->Get("h_evtype"); hB900_evtype->Scale(1./hB900_evtype->Integral());

  //

  TH1F* hT700_qW_dR   = (TH1F*)f_mT700->Get("h_qW_dR"); hT700_qW_dR->Scale(1./hT700_qW_dR->Integral());
  TH1F* hT800_qW_dR   = (TH1F*)f_mT800->Get("h_qW_dR"); hT800_qW_dR->Scale(1./hT800_qW_dR->Integral());
  TH1F* hT900_qW_dR   = (TH1F*)f_mT900->Get("h_qW_dR"); hT900_qW_dR->Scale(1./hT900_qW_dR->Integral());

  TH1F* hT800_qZ_dR   = (TH1F*)f_mT800->Get("h_qZ_dR"); hT800_qZ_dR->Scale(1./hT800_qZ_dR->Integral());
  TH1F* hT800_qH_dR   = (TH1F*)f_mT800->Get("h_qH_dR"); hT800_qH_dR->Scale(1./hT800_qH_dR->Integral());

  TH1F* hB700_qW_dR   = (TH1F*)f_mB700->Get("h_qW_dR"); hB700_qW_dR->Scale(1./hB700_qW_dR->Integral());
  TH1F* hB800_qW_dR   = (TH1F*)f_mB800->Get("h_qW_dR"); hB800_qW_dR->Scale(1./hB800_qW_dR->Integral());
  TH1F* hB900_qW_dR   = (TH1F*)f_mB900->Get("h_qW_dR"); hB900_qW_dR->Scale(1./hB900_qW_dR->Integral());

  TH1F* hB800_qZ_dR   = (TH1F*)f_mB800->Get("h_qZ_dR"); hB800_qZ_dR->Scale(1./hB800_qZ_dR->Integral());
  TH1F* hB800_qH_dR   = (TH1F*)f_mB800->Get("h_qH_dR"); hB800_qH_dR->Scale(1./hB800_qH_dR->Integral());

  //

  TH1F* hT700_lep_n   = (TH1F*)f_mT700->Get("h_lep_n"); hT700_lep_n->Scale(1./hT700_lep_n->Integral());
  TH1F* hT800_lep_n   = (TH1F*)f_mT800->Get("h_lep_n"); hT800_lep_n->Scale(1./hT800_lep_n->Integral());
  TH1F* hT900_lep_n   = (TH1F*)f_mT900->Get("h_lep_n"); hT900_lep_n->Scale(1./hT900_lep_n->Integral());

  TH1F* hB700_lep_n   = (TH1F*)f_mB700->Get("h_lep_n"); hB700_lep_n->Scale(1./hB700_lep_n->Integral());
  TH1F* hB800_lep_n   = (TH1F*)f_mB800->Get("h_lep_n"); hB800_lep_n->Scale(1./hB800_lep_n->Integral());
  TH1F* hB900_lep_n   = (TH1F*)f_mB900->Get("h_lep_n"); hB900_lep_n->Scale(1./hB900_lep_n->Integral());

  TH1F* hT800_lep_n_WW = (TH1F*)f_mT800->Get("h_lep_n_WW"); hT800_lep_n_WW->Scale(1./hT800_lep_n_WW->Integral());
  TH1F* hT800_lep_n_ZZ = (TH1F*)f_mT800->Get("h_lep_n_ZZ"); hT800_lep_n_ZZ->Scale(1./hT800_lep_n_ZZ->Integral());
  TH1F* hT800_lep_n_HH = (TH1F*)f_mT800->Get("h_lep_n_HH"); hT800_lep_n_HH->Scale(1./hT800_lep_n_HH->Integral());
  TH1F* hT800_lep_n_WZ = (TH1F*)f_mT800->Get("h_lep_n_WZ"); hT800_lep_n_WZ->Scale(1./hT800_lep_n_WZ->Integral());
  TH1F* hT800_lep_n_WH = (TH1F*)f_mT800->Get("h_lep_n_WH"); hT800_lep_n_WH->Scale(1./hT800_lep_n_WH->Integral());
  TH1F* hT800_lep_n_ZH = (TH1F*)f_mT800->Get("h_lep_n_ZH"); hT800_lep_n_ZH->Scale(1./hT800_lep_n_ZH->Integral());

  TH1F* hB800_lep_n_WW = (TH1F*)f_mB800->Get("h_lep_n_WW"); hB800_lep_n_WW->Scale(1./hB800_lep_n_WW->Integral());
  TH1F* hB800_lep_n_ZZ = (TH1F*)f_mB800->Get("h_lep_n_ZZ"); hB800_lep_n_ZZ->Scale(1./hB800_lep_n_ZZ->Integral());
  TH1F* hB800_lep_n_HH = (TH1F*)f_mB800->Get("h_lep_n_HH"); hB800_lep_n_HH->Scale(1./hB800_lep_n_HH->Integral());
  TH1F* hB800_lep_n_WZ = (TH1F*)f_mB800->Get("h_lep_n_WZ"); hB800_lep_n_WZ->Scale(1./hB800_lep_n_WZ->Integral());
  TH1F* hB800_lep_n_WH = (TH1F*)f_mB800->Get("h_lep_n_WH"); hB800_lep_n_WH->Scale(1./hB800_lep_n_WH->Integral());
  TH1F* hB800_lep_n_ZH = (TH1F*)f_mB800->Get("h_lep_n_ZH"); hB800_lep_n_ZH->Scale(1./hB800_lep_n_ZH->Integral());

  //

  TH1F* hT700_el_pt   = (TH1F*)f_mT700->Get("h_el_pt"); hT700_el_pt->Scale(1./hT700_el_pt->Integral());      
  TH1F* hT800_el_pt   = (TH1F*)f_mT800->Get("h_el_pt"); hT800_el_pt->Scale(1./hT800_el_pt->Integral());
  TH1F* hT900_el_pt   = (TH1F*)f_mT900->Get("h_el_pt"); hT900_el_pt->Scale(1./hT900_el_pt->Integral());

  TH1F* hT700_el_eta  = (TH1F*)f_mT700->Get("h_el_eta"); hT700_el_eta->Scale(1./hT700_el_eta->Integral());
  TH1F* hT800_el_eta  = (TH1F*)f_mT800->Get("h_el_eta"); hT800_el_eta->Scale(1./hT800_el_eta->Integral());
  TH1F* hT900_el_eta  = (TH1F*)f_mT900->Get("h_el_eta"); hT900_el_eta->Scale(1./hT900_el_eta->Integral());

  TH1F* hB700_el_pt   = (TH1F*)f_mB700->Get("h_el_pt"); hB700_el_pt->Scale(1./hB700_el_pt->Integral());
  TH1F* hB800_el_pt   = (TH1F*)f_mB800->Get("h_el_pt"); hB800_el_pt->Scale(1./hB800_el_pt->Integral());
  TH1F* hB900_el_pt   = (TH1F*)f_mB900->Get("h_el_pt"); hB900_el_pt->Scale(1./hB900_el_pt->Integral());

  TH1F* hB700_el_eta  = (TH1F*)f_mB700->Get("h_el_eta"); hB700_el_eta->Scale(1./hB700_el_eta->Integral());
  TH1F* hB800_el_eta  = (TH1F*)f_mB800->Get("h_el_eta"); hB800_el_eta->Scale(1./hB800_el_eta->Integral());
  TH1F* hB900_el_eta  = (TH1F*)f_mB900->Get("h_el_eta"); hB900_el_eta->Scale(1./hB900_el_eta->Integral());

  //

  TH1F* hT700_mu_pt   = (TH1F*)f_mT700->Get("h_mu_pt"); hT700_mu_pt->Scale(1./hT700_mu_pt->Integral());
  TH1F* hT800_mu_pt   = (TH1F*)f_mT800->Get("h_mu_pt"); hT800_mu_pt->Scale(1./hT800_mu_pt->Integral());
  TH1F* hT900_mu_pt   = (TH1F*)f_mT900->Get("h_mu_pt"); hT900_mu_pt->Scale(1./hT900_mu_pt->Integral());

  TH1F* hT700_mu_eta  = (TH1F*)f_mT700->Get("h_mu_eta"); hT700_mu_eta->Scale(1./hT700_mu_eta->Integral());
  TH1F* hT800_mu_eta  = (TH1F*)f_mT800->Get("h_mu_eta"); hT800_mu_eta->Scale(1./hT800_mu_eta->Integral());
  TH1F* hT900_mu_eta  = (TH1F*)f_mT900->Get("h_mu_eta"); hT900_mu_eta->Scale(1./hT900_mu_eta->Integral());

  TH1F* hB700_mu_pt   = (TH1F*)f_mB700->Get("h_mu_pt"); hB700_mu_pt->Scale(1./hB700_mu_pt->Integral());
  TH1F* hB800_mu_pt   = (TH1F*)f_mB800->Get("h_mu_pt"); hB800_mu_pt->Scale(1./hB800_mu_pt->Integral());
  TH1F* hB900_mu_pt   = (TH1F*)f_mB900->Get("h_mu_pt"); hB900_mu_pt->Scale(1./hB900_mu_pt->Integral());

  TH1F* hB700_mu_eta  = (TH1F*)f_mB700->Get("h_mu_eta"); hB700_mu_eta->Scale(1./hB700_mu_eta->Integral());
  TH1F* hB800_mu_eta  = (TH1F*)f_mB800->Get("h_mu_eta"); hB800_mu_eta->Scale(1./hB800_mu_eta->Integral());
  TH1F* hB900_mu_eta  = (TH1F*)f_mB900->Get("h_mu_eta"); hB900_mu_eta->Scale(1./hB900_mu_eta->Integral());

  //

  TH1F* hT700_MET  = (TH1F*)f_mT700->Get("h_MET"); hT700_MET->Scale(1./hT700_MET->Integral());
  TH1F* hT800_MET  = (TH1F*)f_mT800->Get("h_MET"); hT800_MET->Scale(1./hT800_MET->Integral());
  TH1F* hT900_MET  = (TH1F*)f_mT900->Get("h_MET"); hT900_MET->Scale(1./hT900_MET->Integral());

  TH1F* hT800_MET_WW = (TH1F*)f_mT800->Get("h_MET_WW"); hT800_MET_WW->Scale(1./hT800_MET_WW->Integral());
  TH1F* hT800_MET_ZZ = (TH1F*)f_mT800->Get("h_MET_ZZ"); hT800_MET_ZZ->Scale(1./hT800_MET_ZZ->Integral());
  TH1F* hT800_MET_HH = (TH1F*)f_mT800->Get("h_MET_HH"); hT800_MET_HH->Scale(1./hT800_MET_HH->Integral());
  TH1F* hT800_MET_WZ = (TH1F*)f_mT800->Get("h_MET_WZ"); hT800_MET_WZ->Scale(1./hT800_MET_WZ->Integral());
  TH1F* hT800_MET_WH = (TH1F*)f_mT800->Get("h_MET_WH"); hT800_MET_WH->Scale(1./hT800_MET_WH->Integral());
  TH1F* hT800_MET_ZH = (TH1F*)f_mT800->Get("h_MET_ZH"); hT800_MET_ZH->Scale(1./hT800_MET_ZH->Integral());

  TH1F* hB700_MET  = (TH1F*)f_mB700->Get("h_MET"); hB700_MET->Scale(1./hB700_MET->Integral());
  TH1F* hB800_MET  = (TH1F*)f_mB800->Get("h_MET"); hB800_MET->Scale(1./hB800_MET->Integral());
  TH1F* hB900_MET  = (TH1F*)f_mB900->Get("h_MET"); hB900_MET->Scale(1./hB900_MET->Integral());

  TH1F* hB800_MET_WW = (TH1F*)f_mB800->Get("h_MET_WW"); hB800_MET_WW->Scale(1./hB800_MET_WW->Integral());
  TH1F* hB800_MET_ZZ = (TH1F*)f_mB800->Get("h_MET_ZZ"); hB800_MET_ZZ->Scale(1./hB800_MET_ZZ->Integral());
  TH1F* hB800_MET_HH = (TH1F*)f_mB800->Get("h_MET_HH"); hB800_MET_HH->Scale(1./hB800_MET_HH->Integral());
  TH1F* hB800_MET_WZ = (TH1F*)f_mB800->Get("h_MET_WZ"); hB800_MET_WZ->Scale(1./hB800_MET_WZ->Integral());
  TH1F* hB800_MET_WH = (TH1F*)f_mB800->Get("h_MET_WH"); hB800_MET_WH->Scale(1./hB800_MET_WH->Integral());
  TH1F* hB800_MET_ZH = (TH1F*)f_mB800->Get("h_MET_ZH"); hB800_MET_ZH->Scale(1./hB800_MET_ZH->Integral());

  //

  TH1F* hT700_njet  = (TH1F*)f_mT700->Get("h_jet4pre_n"); hT700_njet->Scale(1./hT700_njet->Integral());
  TH1F* hT800_njet  = (TH1F*)f_mT800->Get("h_jet4pre_n"); hT800_njet->Scale(1./hT800_njet->Integral());
  TH1F* hT900_njet  = (TH1F*)f_mT900->Get("h_jet4pre_n"); hT900_njet->Scale(1./hT900_njet->Integral());

  TH1F* hB700_njet  = (TH1F*)f_mB700->Get("h_jet4pre_n"); hB700_njet->Scale(1./hB700_njet->Integral());
  TH1F* hB800_njet  = (TH1F*)f_mB800->Get("h_jet4pre_n"); hB800_njet->Scale(1./hB800_njet->Integral());
  TH1F* hB900_njet  = (TH1F*)f_mB900->Get("h_jet4pre_n"); hB900_njet->Scale(1./hB900_njet->Integral());

  TH1F* hT700_jmass  = (TH1F*)f_mT700->Get("h_jet4_mass"); hT700_jmass->Scale(1./hT700_jmass->Integral());
  TH1F* hT800_jmass  = (TH1F*)f_mT800->Get("h_jet4_mass"); hT800_jmass->Scale(1./hT800_jmass->Integral());
  TH1F* hT900_jmass  = (TH1F*)f_mT900->Get("h_jet4_mass"); hT900_jmass->Scale(1./hT900_jmass->Integral());

  TH1F* hB700_jmass  = (TH1F*)f_mB700->Get("h_jet4_mass"); hB700_jmass->Scale(1./hB700_jmass->Integral());
  TH1F* hB800_jmass  = (TH1F*)f_mB800->Get("h_jet4_mass"); hB800_jmass->Scale(1./hB800_jmass->Integral());
  TH1F* hB900_jmass  = (TH1F*)f_mB900->Get("h_jet4_mass"); hB900_jmass->Scale(1./hB900_jmass->Integral());

  //

  TH1F* hT700_jpt   = (TH1F*)f_mT700->Get("h_jet4_pt"); hT700_jpt->Scale(1./hT700_jpt->Integral());
  TH1F* hT800_jpt   = (TH1F*)f_mT800->Get("h_jet4_pt"); hT800_jpt->Scale(1./hT800_jpt->Integral());
  TH1F* hT900_jpt   = (TH1F*)f_mT900->Get("h_jet4_pt"); hT900_jpt->Scale(1./hT900_jpt->Integral());

  TH1F* hB700_jpt   = (TH1F*)f_mB700->Get("h_jet4_pt"); hB700_jpt->Scale(1./hB700_jpt->Integral());
  TH1F* hB800_jpt   = (TH1F*)f_mB800->Get("h_jet4_pt"); hB800_jpt->Scale(1./hB800_jpt->Integral());
  TH1F* hB900_jpt   = (TH1F*)f_mB900->Get("h_jet4_pt"); hB900_jpt->Scale(1./hB900_jpt->Integral());

  TH1F* hT700_jeta  = (TH1F*)f_mT700->Get("h_jet4_eta"); hT700_jeta->Scale(1./hT700_jeta->Integral());
  TH1F* hT800_jeta  = (TH1F*)f_mT800->Get("h_jet4_eta"); hT800_jeta->Scale(1./hT800_jeta->Integral());
  TH1F* hT900_jeta  = (TH1F*)f_mT900->Get("h_jet4_eta"); hT900_jeta->Scale(1./hT900_jeta->Integral());

  TH1F* hB700_jeta  = (TH1F*)f_mB700->Get("h_jet4_eta"); hB700_jeta->Scale(1./hB700_jeta->Integral());
  TH1F* hB800_jeta  = (TH1F*)f_mB800->Get("h_jet4_eta"); hB800_jeta->Scale(1./hB800_jeta->Integral());
  TH1F* hB900_jeta  = (TH1F*)f_mB900->Get("h_jet4_eta"); hB900_jeta->Scale(1./hB900_jeta->Integral());

  // **************

  TCanvas *c1 = new TCanvas("c1","1: Heavy Quark Mass and Daughters",1000,700);
  c1->SetFillColor(10);
  c1->Divide(2,2);

  TPad *pad11 = (TCanvas*)c1->GetPad(1);
  pad11->cd(); 
  pad11->SetTopMargin(0.06);
  pad11->SetBottomMargin(0.14);
  pad11->SetLeftMargin(0.12);
  pad11->SetRightMargin(0.04);

  hT700_Q_mass->GetXaxis()->SetLabelFont(42);
  hT700_Q_mass->GetXaxis()->SetLabelSize(0.05);
  hT700_Q_mass->GetXaxis()->SetTitleFont(42);
  hT700_Q_mass->GetXaxis()->SetTitleSize(0.06);
  hT700_Q_mass->GetXaxis()->SetTitle("m_{T} [GeV]");

  hT700_Q_mass->GetYaxis()->SetLabelFont(42);
  hT700_Q_mass->GetYaxis()->SetLabelSize(0.05);
  hT700_Q_mass->GetYaxis()->SetTitleFont(42);
  hT700_Q_mass->GetYaxis()->SetTitleSize(0.06);
  hT700_Q_mass->GetYaxis()->SetTitleOffset(1.0);
  hT700_Q_mass->GetYaxis()->SetTitle("Fraction of Entries / 5 GeV");

  hT700_Q_mass->SetLineColor(2);
  hT700_Q_mass->Draw();	

  hT800_Q_mass->SetLineColor(2);
  hT800_Q_mass->SetLineStyle(2);
  hT800_Q_mass->Draw("sames");	

  hT900_Q_mass->SetLineColor(2);
  hT900_Q_mass->SetLineStyle(5);
  hT900_Q_mass->Draw("sames");	

  TLegend* leg1 = new TLegend(0.60,0.55,0.90,0.85);
  leg1->SetFillColor(10);
  leg1->SetBorderSize(0);
  leg1->SetTextFont(62);
  leg1->SetTextSize(0.06);
  leg1->AddEntry(hT700_Q_mass,"#font[42]{m_{T} = 700 GeV}","L");
  leg1->AddEntry(hT800_Q_mass,"#font[42]{m_{T} = 800 GeV}","L");
  leg1->AddEntry(hT900_Q_mass,"#font[42]{m_{T} = 900 GeV}","L");
  leg1->Draw();

  TPad *pad12 = (TCanvas*)c1->GetPad(2);
  pad12->cd(); 
  pad12->SetTopMargin(0.06);
  pad12->SetBottomMargin(0.14);
  pad12->SetLeftMargin(0.12);
  pad12->SetRightMargin(0.04);

  hB700_Q_mass->GetXaxis()->SetLabelFont(42);
  hB700_Q_mass->GetXaxis()->SetLabelSize(0.05);
  hB700_Q_mass->GetXaxis()->SetTitleFont(42);
  hB700_Q_mass->GetXaxis()->SetTitleSize(0.06);
  hB700_Q_mass->GetXaxis()->SetTitle("m_{B} [GeV]");

  hB700_Q_mass->GetYaxis()->SetLabelFont(42);
  hB700_Q_mass->GetYaxis()->SetLabelSize(0.05);
  hB700_Q_mass->GetYaxis()->SetTitleFont(42);
  hB700_Q_mass->GetYaxis()->SetTitleSize(0.06);
  hB700_Q_mass->GetYaxis()->SetTitleOffset(1.0);
  hB700_Q_mass->GetYaxis()->SetTitle("Fraction of Entries / 5 GeV");

  hB700_Q_mass->SetLineColor(4);
  hB700_Q_mass->Draw();	

  hB800_Q_mass->SetLineColor(4);
  hB800_Q_mass->SetLineStyle(2);
  hB800_Q_mass->Draw("sames");	

  hB900_Q_mass->SetLineColor(4);
  hB900_Q_mass->SetLineStyle(5);
  hB900_Q_mass->Draw("sames");	

  TLegend* leg2 = new TLegend(0.60,0.55,0.90,0.85);
  leg2->SetFillColor(10);
  leg2->SetBorderSize(0);
  leg2->SetTextFont(62);
  leg2->SetTextSize(0.06);
  leg2->AddEntry(hB700_Q_mass,"#font[42]{m_{B} = 700 GeV}","L");
  leg2->AddEntry(hB800_Q_mass,"#font[42]{m_{B} = 800 GeV}","L");
  leg2->AddEntry(hB900_Q_mass,"#font[42]{m_{B} = 900 GeV}","L");
  leg2->Draw();

  TPad *pad13 = (TCanvas*)c1->GetPad(3);
  pad13->cd(); 
  pad13->SetTopMargin(0.06);
  pad13->SetBottomMargin(0.14);
  pad13->SetLeftMargin(0.12);
  pad13->SetRightMargin(0.04);

  hT700_Q_dID->GetXaxis()->SetLabelFont(42);
  hT700_Q_dID->GetXaxis()->SetLabelSize(0.05);
  hT700_Q_dID->GetXaxis()->SetTitleFont(42);
  hT700_Q_dID->GetXaxis()->SetTitleSize(0.06);
  hT700_Q_dID->GetXaxis()->SetTitle("T daughter PIDs");

  hT700_Q_dID->GetYaxis()->SetLabelFont(42);
  hT700_Q_dID->GetYaxis()->SetLabelSize(0.05);
  hT700_Q_dID->GetYaxis()->SetTitleFont(42);
  hT700_Q_dID->GetYaxis()->SetTitleSize(0.06);
  hT700_Q_dID->GetYaxis()->SetTitleOffset(1.0);
  hT700_Q_dID->GetYaxis()->SetTitle("Fraction of Entries / PID");

  hT700_Q_dID->SetLineColor(2);
  hT700_Q_dID->Draw();	

  hT800_Q_dID->SetLineColor(2);
  hT800_Q_dID->SetLineStyle(2);
  hT800_Q_dID->Draw("sames");	

  hT900_Q_dID->SetLineColor(2);
  hT900_Q_dID->SetLineStyle(5);
  hT900_Q_dID->Draw("sames");

  TPad *pad14 = (TCanvas*)c1->GetPad(4);
  pad14->cd(); 
  pad14->SetTopMargin(0.06);
  pad14->SetBottomMargin(0.14);
  pad14->SetLeftMargin(0.12);
  pad14->SetRightMargin(0.04);

  hB700_Q_dID->GetXaxis()->SetLabelFont(42);
  hB700_Q_dID->GetXaxis()->SetLabelSize(0.05);
  hB700_Q_dID->GetXaxis()->SetTitleFont(42);
  hB700_Q_dID->GetXaxis()->SetTitleSize(0.06);
  hB700_Q_dID->GetXaxis()->SetTitle("B daughter PIDs");

  hB700_Q_dID->GetYaxis()->SetLabelFont(42);
  hB700_Q_dID->GetYaxis()->SetLabelSize(0.05);
  hB700_Q_dID->GetYaxis()->SetTitleFont(42);
  hB700_Q_dID->GetYaxis()->SetTitleSize(0.06);
  hB700_Q_dID->GetYaxis()->SetTitleOffset(1.0);
  hB700_Q_dID->GetYaxis()->SetTitle("Fraction of Entries / PID");

  hB700_Q_dID->SetLineColor(4);
  hB700_Q_dID->Draw();	

  hB800_Q_dID->SetLineColor(4);
  hB800_Q_dID->SetLineStyle(2);
  hB800_Q_dID->Draw("sames");	

  hB900_Q_dID->SetLineColor(4);
  hB900_Q_dID->SetLineStyle(5);
  hB900_Q_dID->Draw("sames");

  // **************

  TCanvas *c2 = new TCanvas("c2","2: Heavy Quark pT and Eta",1000,700);
  c2->SetFillColor(10);
  c2->Divide(2,2);

  TPad *pad21 = (TCanvas*)c2->GetPad(1);
  pad21->cd(); 
  pad21->SetTopMargin(0.06);
  pad21->SetBottomMargin(0.14);
  pad21->SetLeftMargin(0.12);
  pad21->SetRightMargin(0.04);

  hT700_Q_pt->GetXaxis()->SetLabelFont(42);
  hT700_Q_pt->GetXaxis()->SetLabelSize(0.05);
  hT700_Q_pt->GetXaxis()->SetTitleFont(42);
  hT700_Q_pt->GetXaxis()->SetTitleSize(0.06);
  hT700_Q_pt->GetXaxis()->SetTitle("p_{T}(T) [GeV]");

  hT700_Q_pt->GetYaxis()->SetLabelFont(42);
  hT700_Q_pt->GetYaxis()->SetLabelSize(0.04);
  hT700_Q_pt->GetYaxis()->SetTitleFont(42);
  hT700_Q_pt->GetYaxis()->SetTitleSize(0.06);
  hT700_Q_pt->GetYaxis()->SetTitleOffset(1.0);
  hT700_Q_pt->GetYaxis()->SetTitle("Fraction of Entries / 20 GeV");

  hT700_Q_pt->SetLineColor(2);
  hT700_Q_pt->Draw();	

  hT800_Q_pt->SetLineColor(2);
  hT800_Q_pt->SetLineStyle(2);
  hT800_Q_pt->Draw("sames");	

  hT900_Q_pt->SetLineColor(2);
  hT900_Q_pt->SetLineStyle(5);
  hT900_Q_pt->Draw("sames");

  leg1->Draw();

  TPad *pad22 = (TCanvas*)c2->GetPad(2);
  pad22->cd(); 
  pad22->SetTopMargin(0.06);
  pad22->SetBottomMargin(0.14);
  pad22->SetLeftMargin(0.12);
  pad22->SetRightMargin(0.04);

  hB700_Q_pt->GetXaxis()->SetLabelFont(42);
  hB700_Q_pt->GetXaxis()->SetLabelSize(0.05);
  hB700_Q_pt->GetXaxis()->SetTitleFont(42);
  hB700_Q_pt->GetXaxis()->SetTitleSize(0.06);
  hB700_Q_pt->GetXaxis()->SetTitle("p_{T}(B) [GeV]");

  hB700_Q_pt->GetYaxis()->SetLabelFont(42);
  hB700_Q_pt->GetYaxis()->SetLabelSize(0.04);
  hB700_Q_pt->GetYaxis()->SetTitleFont(42);
  hB700_Q_pt->GetYaxis()->SetTitleSize(0.06);
  hB700_Q_pt->GetYaxis()->SetTitleOffset(1.0);
  hB700_Q_pt->GetYaxis()->SetTitle("Fraction of Entries / 20 GeV");

  hB700_Q_pt->SetLineColor(4);
  hB700_Q_pt->Draw();	

  hB800_Q_pt->SetLineColor(4);
  hB800_Q_pt->SetLineStyle(2);
  hB800_Q_pt->Draw("sames");	

  hB900_Q_pt->SetLineColor(4);
  hB900_Q_pt->SetLineStyle(5);
  hB900_Q_pt->Draw("sames");

  leg2->Draw();

  TPad *pad23 = (TCanvas*)c2->GetPad(3);
  pad23->cd(); 
  pad23->SetTopMargin(0.06);
  pad23->SetBottomMargin(0.14);
  pad23->SetLeftMargin(0.12);
  pad23->SetRightMargin(0.04);

  hT700_Q_eta->GetXaxis()->SetLabelFont(42);
  hT700_Q_eta->GetXaxis()->SetLabelSize(0.05);
  hT700_Q_eta->GetXaxis()->SetTitleFont(42);
  hT700_Q_eta->GetXaxis()->SetTitleSize(0.06);
  hT700_Q_eta->GetXaxis()->SetTitle("#eta(T)");

  hT700_Q_eta->GetYaxis()->SetLabelFont(42);
  hT700_Q_eta->GetYaxis()->SetLabelSize(0.04);
  hT700_Q_eta->GetYaxis()->SetTitleFont(42);
  hT700_Q_eta->GetYaxis()->SetTitleSize(0.06);
  hT700_Q_eta->GetYaxis()->SetTitleOffset(1.0);
  hT700_Q_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hT700_Q_eta->SetMaximum(0.05);
  hT700_Q_eta->SetLineColor(2);
  hT700_Q_eta->Draw();	

  hT800_Q_eta->SetLineColor(2);
  hT800_Q_eta->SetLineStyle(2);
  hT800_Q_eta->Draw("sames");	

  hT900_Q_eta->SetLineColor(2);
  hT900_Q_eta->SetLineStyle(5);
  hT900_Q_eta->Draw("sames");

  TPad *pad24 = (TCanvas*)c2->GetPad(4);
  pad24->cd(); 
  pad24->SetTopMargin(0.06);
  pad24->SetBottomMargin(0.14);
  pad24->SetLeftMargin(0.12);
  pad24->SetRightMargin(0.04);

  hB700_Q_eta->GetXaxis()->SetLabelFont(42);
  hB700_Q_eta->GetXaxis()->SetLabelSize(0.05);
  hB700_Q_eta->GetXaxis()->SetTitleFont(42);
  hB700_Q_eta->GetXaxis()->SetTitleSize(0.06);
  hB700_Q_eta->GetXaxis()->SetTitle("#eta(B)");

  hB700_Q_eta->GetYaxis()->SetLabelFont(42);
  hB700_Q_eta->GetYaxis()->SetLabelSize(0.04);
  hB700_Q_eta->GetYaxis()->SetTitleFont(42);
  hB700_Q_eta->GetYaxis()->SetTitleSize(0.06);
  hB700_Q_eta->GetYaxis()->SetTitleOffset(1.0);
  hB700_Q_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hB700_Q_eta->SetMaximum(0.05);
  hB700_Q_eta->SetLineColor(4);
  hB700_Q_eta->Draw();	

  hB800_Q_eta->SetLineColor(4);
  hB800_Q_eta->SetLineStyle(2);
  hB800_Q_eta->Draw("sames");	

  hB900_Q_eta->SetLineColor(4);
  hB900_Q_eta->SetLineStyle(5);
  hB900_Q_eta->Draw("sames");

  // **************

  TCanvas *c3 = new TCanvas("c3","3: W Daughter pT and Eta",1000,700);
  c3->SetFillColor(10);
  c3->Divide(2,2);

  TPad *pad31 = (TCanvas*)c3->GetPad(1);
  pad31->cd(); 
  pad31->SetTopMargin(0.06);
  pad31->SetBottomMargin(0.14);
  pad31->SetLeftMargin(0.12);
  pad31->SetRightMargin(0.04);

  hT700_W_pt->GetXaxis()->SetLabelFont(42);
  hT700_W_pt->GetXaxis()->SetLabelSize(0.05);
  hT700_W_pt->GetXaxis()->SetTitleFont(42);
  hT700_W_pt->GetXaxis()->SetTitleSize(0.06);
  hT700_W_pt->GetXaxis()->SetTitle("p_{T}(W) [GeV]");

  hT700_W_pt->GetYaxis()->SetLabelFont(42);
  hT700_W_pt->GetYaxis()->SetLabelSize(0.05);
  hT700_W_pt->GetYaxis()->SetTitleFont(42);
  hT700_W_pt->GetYaxis()->SetTitleSize(0.06);
  hT700_W_pt->GetYaxis()->SetTitleOffset(1.0);
  hT700_W_pt->GetYaxis()->SetTitle("Fraction of Entries / 20 GeV");

  hT700_W_pt->SetLineColor(2);
  hT700_W_pt->Draw();	

  hT800_W_pt->SetLineColor(2);
  hT800_W_pt->SetLineStyle(2);
  hT800_W_pt->Draw("sames");	

  hT900_W_pt->SetLineColor(2);
  hT900_W_pt->SetLineStyle(5);
  hT900_W_pt->Draw("sames");

  leg1->Draw();

  TLatex* txt=new TLatex(0.65, 0.35, "T #rightarrow Wb");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.06);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  TPad *pad32 = (TCanvas*)c3->GetPad(2);
  pad32->cd(); 
  pad32->SetTopMargin(0.06);
  pad32->SetBottomMargin(0.14);
  pad32->SetLeftMargin(0.12);
  pad32->SetRightMargin(0.04);

  hB700_W_pt->GetXaxis()->SetLabelFont(42);
  hB700_W_pt->GetXaxis()->SetLabelSize(0.05);
  hB700_W_pt->GetXaxis()->SetTitleFont(42);
  hB700_W_pt->GetXaxis()->SetTitleSize(0.06);
  hB700_W_pt->GetXaxis()->SetTitle("p_{T}(W) [GeV]");

  hB700_W_pt->GetYaxis()->SetLabelFont(42);
  hB700_W_pt->GetYaxis()->SetLabelSize(0.05);
  hB700_W_pt->GetYaxis()->SetTitleFont(42);
  hB700_W_pt->GetYaxis()->SetTitleSize(0.06);
  hB700_W_pt->GetYaxis()->SetTitleOffset(1.0);
  hB700_W_pt->GetYaxis()->SetTitle("Fraction of Entries / 20 GeV");

  hB700_W_pt->SetLineColor(4);
  hB700_W_pt->Draw();	

  hB800_W_pt->SetLineColor(4);
  hB800_W_pt->SetLineStyle(2);
  hB800_W_pt->Draw("sames");	

  hB900_W_pt->SetLineColor(4);
  hB900_W_pt->SetLineStyle(5);
  hB900_W_pt->Draw("sames");

  leg2->Draw();

  TLatex* txt=new TLatex(0.65, 0.35, "B #rightarrow Wt");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.06);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  TPad *pad33 = (TCanvas*)c3->GetPad(3);
  pad33->cd(); 
  pad33->SetTopMargin(0.06);
  pad33->SetBottomMargin(0.14);
  pad33->SetLeftMargin(0.12);
  pad33->SetRightMargin(0.04);

  hT700_W_eta->GetXaxis()->SetLabelFont(42);
  hT700_W_eta->GetXaxis()->SetLabelSize(0.05);
  hT700_W_eta->GetXaxis()->SetTitleFont(42);
  hT700_W_eta->GetXaxis()->SetTitleSize(0.06);
  hT700_W_eta->GetXaxis()->SetTitle("#eta(W)");

  hT700_W_eta->GetYaxis()->SetLabelFont(42);
  hT700_W_eta->GetYaxis()->SetLabelSize(0.05);
  hT700_W_eta->GetYaxis()->SetTitleFont(42);
  hT700_W_eta->GetYaxis()->SetTitleSize(0.06);
  hT700_W_eta->GetYaxis()->SetTitleOffset(1.0);
  hT700_W_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hT700_W_eta->SetMaximum(0.08);
  hT700_W_eta->SetLineColor(2);
  hT700_W_eta->Draw();	

  hT800_W_eta->SetLineColor(2);
  hT800_W_eta->SetLineStyle(2);
  hT800_W_eta->Draw("sames");	

  hT900_W_eta->SetLineColor(2);
  hT900_W_eta->SetLineStyle(5);
  hT900_W_eta->Draw("sames");

  TPad *pad34 = (TCanvas*)c3->GetPad(4);
  pad34->cd(); 
  pad34->SetTopMargin(0.06);
  pad34->SetBottomMargin(0.14);
  pad34->SetLeftMargin(0.12);
  pad34->SetRightMargin(0.04);

  hB700_W_eta->GetXaxis()->SetLabelFont(42);
  hB700_W_eta->GetXaxis()->SetLabelSize(0.05);
  hB700_W_eta->GetXaxis()->SetTitleFont(42);
  hB700_W_eta->GetXaxis()->SetTitleSize(0.06);
  hB700_W_eta->GetXaxis()->SetTitle("#eta(W)");

  hB700_W_eta->GetYaxis()->SetLabelFont(42);
  hB700_W_eta->GetYaxis()->SetLabelSize(0.05);
  hB700_W_eta->GetYaxis()->SetTitleFont(42);
  hB700_W_eta->GetYaxis()->SetTitleSize(0.06);
  hB700_W_eta->GetYaxis()->SetTitleOffset(1.0);
  hB700_W_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hB700_W_eta->SetMaximum(0.08);
  hB700_W_eta->SetLineColor(4);
  hB700_W_eta->Draw();	

  hB800_W_eta->SetLineColor(4);
  hB800_W_eta->SetLineStyle(2);
  hB800_W_eta->Draw("sames");	

  hB900_W_eta->SetLineColor(4);
  hB900_W_eta->SetLineStyle(5);
  hB900_W_eta->Draw("sames");

  // **************

  TCanvas *c4 = new TCanvas("c4","4: Boson Daughter pT and Eta",1000,700);
  c4->SetFillColor(10);
  c4->Divide(2,2);

  TPad *pad41 = (TCanvas*)c4->GetPad(1);
  pad41->cd(); 
  pad41->SetTopMargin(0.06);
  pad41->SetBottomMargin(0.14);
  pad41->SetLeftMargin(0.12);
  pad41->SetRightMargin(0.04);

  TH1F* hcT800_W_pt = hT800_W_pt->Clone("hcT800_W_pt");
  hcT800_W_pt->GetXaxis()->SetLabelFont(42);
  hcT800_W_pt->GetXaxis()->SetLabelSize(0.05);
  hcT800_W_pt->GetXaxis()->SetTitleFont(42);
  hcT800_W_pt->GetXaxis()->SetTitleSize(0.06);
  hcT800_W_pt->GetXaxis()->SetTitle("p_{T}(W,Z,H) [GeV]");

  hcT800_W_pt->GetYaxis()->SetLabelFont(42);
  hcT800_W_pt->GetYaxis()->SetLabelSize(0.04);
  hcT800_W_pt->GetYaxis()->SetTitleFont(42);
  hcT800_W_pt->GetYaxis()->SetTitleSize(0.06);
  hcT800_W_pt->GetYaxis()->SetTitleOffset(1.0);
  hcT800_W_pt->GetYaxis()->SetTitle("Fraction of Entries / 20 GeV");

  hcT800_W_pt->SetMaximum(0.05);
  hcT800_W_pt->SetLineStyle(1);
  hcT800_W_pt->SetLineColor(2);
  hcT800_W_pt->Draw();	

  hT800_Z_pt->SetLineColor(2);
  hT800_Z_pt->SetLineStyle(2);
  hT800_Z_pt->Draw("sames");	

  hT800_H_pt->SetLineColor(2);
  hT800_H_pt->SetLineStyle(5);
  hT800_H_pt->Draw("sames");

  TLatex* txt=new TLatex(0.60, 0.35, "m_{T} = 800 GeV");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.06);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  TLegend* leg3 = new TLegend(0.60,0.55,0.90,0.85);
  leg3->SetFillColor(10);
  leg3->SetBorderSize(0);
  leg3->SetTextFont(62);
  leg3->SetTextSize(0.06);
  leg3->AddEntry(hcT800_W_pt,"#font[42]{T #rightarrow Wb}","L");
  leg3->AddEntry(hT800_Z_pt,"#font[42]{T #rightarrow Zt}","L");
  leg3->AddEntry(hT800_H_pt,"#font[42]{T #rightarrow Ht}","L");
  leg3->Draw();

  TPad *pad42 = (TCanvas*)c4->GetPad(2);
  pad42->cd(); 
  pad42->SetTopMargin(0.06);
  pad42->SetBottomMargin(0.14);
  pad42->SetLeftMargin(0.12);
  pad42->SetRightMargin(0.04);

  TH1F* hcB800_W_pt = hB800_W_pt->Clone("hcB800_W_pt");
  hcB800_W_pt->GetXaxis()->SetLabelFont(42);
  hcB800_W_pt->GetXaxis()->SetLabelSize(0.05);
  hcB800_W_pt->GetXaxis()->SetTitleFont(42);
  hcB800_W_pt->GetXaxis()->SetTitleSize(0.06);
  hcB800_W_pt->GetXaxis()->SetTitle("p_{T}(W,Z,H) [GeV]");

  hcB800_W_pt->GetYaxis()->SetLabelFont(42);
  hcB800_W_pt->GetYaxis()->SetLabelSize(0.04);
  hcB800_W_pt->GetYaxis()->SetTitleFont(42);
  hcB800_W_pt->GetYaxis()->SetTitleSize(0.06);
  hcB800_W_pt->GetYaxis()->SetTitleOffset(1.0);
  hcB800_W_pt->GetYaxis()->SetTitle("Fraction of Entries / 20 GeV");

  hcB800_W_pt->SetMaximum(0.05);
  hcB800_W_pt->SetLineStyle(1);
  hcB800_W_pt->SetLineColor(4);
  hcB800_W_pt->Draw();	

  hB800_Z_pt->SetLineColor(4);
  hB800_Z_pt->SetLineStyle(2);
  hB800_Z_pt->Draw("sames");	

  hB800_H_pt->SetLineColor(4);
  hB800_H_pt->SetLineStyle(5);
  hB800_H_pt->Draw("sames");

  TLatex* txt=new TLatex(0.60, 0.35, "m_{B} = 800 GeV");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.06);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  TLegend* leg4 = new TLegend(0.60,0.55,0.90,0.85);
  leg4->SetFillColor(10);
  leg4->SetBorderSize(0);
  leg4->SetTextFont(62);
  leg4->SetTextSize(0.06);
  leg4->AddEntry(hcB800_W_pt,"#font[42]{B #rightarrow Wt}","L");
  leg4->AddEntry(hB800_Z_pt,"#font[42]{B #rightarrow Zb}","L");
  leg4->AddEntry(hB800_H_pt,"#font[42]{B #rightarrow Hb}","L");
  leg4->Draw();

  TPad *pad43 = (TCanvas*)c4->GetPad(3);
  pad43->cd(); 
  pad43->SetTopMargin(0.06);
  pad43->SetBottomMargin(0.14);
  pad43->SetLeftMargin(0.12);
  pad43->SetRightMargin(0.04);

  TH1F* hcT800_W_eta = hT800_W_eta->Clone("hcT800_W_eta");
  hcT800_W_eta->GetXaxis()->SetLabelFont(42);
  hcT800_W_eta->GetXaxis()->SetLabelSize(0.05);
  hcT800_W_eta->GetXaxis()->SetTitleFont(42);
  hcT800_W_eta->GetXaxis()->SetTitleSize(0.06);
  hcT800_W_eta->GetXaxis()->SetTitle("#eta(W,Z,H)");

  hcT800_W_eta->GetYaxis()->SetLabelFont(42);
  hcT800_W_eta->GetYaxis()->SetLabelSize(0.04);
  hcT800_W_eta->GetYaxis()->SetTitleFont(42);
  hcT800_W_eta->GetYaxis()->SetTitleSize(0.06);
  hcT800_W_eta->GetYaxis()->SetTitleOffset(1.0);
  hcT800_W_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hcT800_W_eta->SetLineColor(2);
  hcT800_W_eta->SetLineStyle(1);
  hcT800_W_eta->Draw();	

  hT800_Z_eta->SetLineColor(2);
  hT800_Z_eta->SetLineStyle(2);
  hT800_Z_eta->Draw("sames");	

  hT800_H_eta->SetLineColor(2);
  hT800_H_eta->SetLineStyle(5);
  hT800_H_eta->Draw("sames");

  TPad *pad44 = (TCanvas*)c4->GetPad(4);
  pad44->cd(); 
  pad44->SetTopMargin(0.06);
  pad44->SetBottomMargin(0.14);
  pad44->SetLeftMargin(0.12);
  pad44->SetRightMargin(0.04);

  TH1F* hcB800_W_eta = hB800_W_eta->Clone("hcB800_W_eta");
  hcB800_W_eta->GetXaxis()->SetLabelFont(42);
  hcB800_W_eta->GetXaxis()->SetLabelSize(0.05);
  hcB800_W_eta->GetXaxis()->SetTitleFont(42);
  hcB800_W_eta->GetXaxis()->SetTitleSize(0.06);
  hcB800_W_eta->GetXaxis()->SetTitle("#eta(W,Z,H)");

  hcB800_W_eta->GetYaxis()->SetLabelFont(42);
  hcB800_W_eta->GetYaxis()->SetLabelSize(0.04);
  hcB800_W_eta->GetYaxis()->SetTitleFont(42);
  hcB800_W_eta->GetYaxis()->SetTitleSize(0.06);
  hcB800_W_eta->GetYaxis()->SetTitleOffset(1.0);
  hcB800_W_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hcB800_W_eta->SetLineColor(4);
  hcB800_W_eta->SetLineStyle(1);
  hcB800_W_eta->Draw();	

  hB800_Z_eta->SetLineColor(4);
  hB800_Z_eta->SetLineStyle(2);
  hB800_Z_eta->Draw("sames");	

  hB800_H_eta->SetLineColor(4);
  hB800_H_eta->SetLineStyle(5);
  hB800_H_eta->Draw("sames");

  // **************

  TCanvas *c5 = new TCanvas("c5","5: Quark Daughter pT and Eta",1000,700);
  c5->SetFillColor(10);
  c5->Divide(2,2);

  TPad *pad51 = (TCanvas*)c5->GetPad(1);
  pad51->cd(); 
  pad51->SetTopMargin(0.06);
  pad51->SetBottomMargin(0.14);
  pad51->SetLeftMargin(0.12);
  pad51->SetRightMargin(0.04);

  hT800_qW_pt->GetXaxis()->SetLabelFont(42);
  hT800_qW_pt->GetXaxis()->SetLabelSize(0.05);
  hT800_qW_pt->GetXaxis()->SetTitleFont(42);
  hT800_qW_pt->GetXaxis()->SetTitleSize(0.06);
  hT800_qW_pt->GetXaxis()->SetTitle("p_{T}(t,b) [GeV]");

  hT800_qW_pt->GetYaxis()->SetLabelFont(42);
  hT800_qW_pt->GetYaxis()->SetLabelSize(0.04);
  hT800_qW_pt->GetYaxis()->SetTitleFont(42);
  hT800_qW_pt->GetYaxis()->SetTitleSize(0.06);
  hT800_qW_pt->GetYaxis()->SetTitleOffset(1.0);
  hT800_qW_pt->GetYaxis()->SetTitle("Fraction of Entries / 20 GeV");

  hT800_qW_pt->SetMaximum(0.05);
  hT800_qW_pt->SetLineStyle(1);
  hT800_qW_pt->SetLineColor(2);
  hT800_qW_pt->Draw();	

  hT800_qZ_pt->SetLineColor(2);
  hT800_qZ_pt->SetLineStyle(2);
  hT800_qZ_pt->Draw("sames");	

  hT800_qH_pt->SetLineColor(2);
  hT800_qH_pt->SetLineStyle(5);
  hT800_qH_pt->Draw("sames");

  TLatex* txt=new TLatex(0.60, 0.35, "m_{T} = 800 GeV");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.06);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  leg3->Draw();

  TPad *pad52 = (TCanvas*)c5->GetPad(2);
  pad52->cd(); 
  pad52->SetTopMargin(0.06);
  pad52->SetBottomMargin(0.14);
  pad52->SetLeftMargin(0.12);
  pad52->SetRightMargin(0.04);

  hB800_qW_pt->GetXaxis()->SetLabelFont(42);
  hB800_qW_pt->GetXaxis()->SetLabelSize(0.05);
  hB800_qW_pt->GetXaxis()->SetTitleFont(42);
  hB800_qW_pt->GetXaxis()->SetTitleSize(0.06);
  hB800_qW_pt->GetXaxis()->SetTitle("p_{T}(t,b) [GeV]");

  hB800_qW_pt->GetYaxis()->SetLabelFont(42);
  hB800_qW_pt->GetYaxis()->SetLabelSize(0.04);
  hB800_qW_pt->GetYaxis()->SetTitleFont(42);
  hB800_qW_pt->GetYaxis()->SetTitleSize(0.06);
  hB800_qW_pt->GetYaxis()->SetTitleOffset(1.0);
  hB800_qW_pt->GetYaxis()->SetTitle("Fraction of Entries / 20 GeV");

  hB800_qW_pt->SetMaximum(0.05);
  hB800_qW_pt->SetLineStyle(1);
  hB800_qW_pt->SetLineColor(4);
  hB800_qW_pt->Draw();	

  hB800_qZ_pt->SetLineColor(4);
  hB800_qZ_pt->SetLineStyle(2);
  hB800_qZ_pt->Draw("sames");	

  hB800_qH_pt->SetLineColor(4);
  hB800_qH_pt->SetLineStyle(5);
  hB800_qH_pt->Draw("sames");

  TLatex* txt=new TLatex(0.60, 0.35, "m_{B} = 800 GeV");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.06);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  leg4->Draw();

  TPad *pad53 = (TCanvas*)c5->GetPad(3);
  pad53->cd(); 
  pad53->SetTopMargin(0.06);
  pad53->SetBottomMargin(0.14);
  pad53->SetLeftMargin(0.12);
  pad53->SetRightMargin(0.04);

  hT800_qW_eta->GetXaxis()->SetLabelFont(42);
  hT800_qW_eta->GetXaxis()->SetLabelSize(0.05);
  hT800_qW_eta->GetXaxis()->SetTitleFont(42);
  hT800_qW_eta->GetXaxis()->SetTitleSize(0.06);
  hT800_qW_eta->GetXaxis()->SetTitle("#eta(t,b)");

  hT800_qW_eta->GetYaxis()->SetLabelFont(42);
  hT800_qW_eta->GetYaxis()->SetLabelSize(0.04);
  hT800_qW_eta->GetYaxis()->SetTitleFont(42);
  hT800_qW_eta->GetYaxis()->SetTitleSize(0.06);
  hT800_qW_eta->GetYaxis()->SetTitleOffset(1.0);
  hT800_qW_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hT800_qW_eta->SetMaximum(0.08);
  hT800_qW_eta->SetLineStyle(1);
  hT800_qW_eta->SetLineColor(2);
  hT800_qW_eta->Draw();	

  hT800_qZ_eta->SetLineColor(2);
  hT800_qZ_eta->SetLineStyle(2);
  hT800_qZ_eta->Draw("sames");	

  hT800_qH_eta->SetLineColor(2);
  hT800_qH_eta->SetLineStyle(5);
  hT800_qH_eta->Draw("sames");


  TPad *pad54 = (TCanvas*)c5->GetPad(4);
  pad54->cd(); 
  pad54->SetTopMargin(0.06);
  pad54->SetBottomMargin(0.14);
  pad54->SetLeftMargin(0.12);
  pad54->SetRightMargin(0.04);

  hB800_qW_eta->GetXaxis()->SetLabelFont(42);
  hB800_qW_eta->GetXaxis()->SetLabelSize(0.05);
  hB800_qW_eta->GetXaxis()->SetTitleFont(42);
  hB800_qW_eta->GetXaxis()->SetTitleSize(0.06);
  hB800_qW_eta->GetXaxis()->SetTitle("#eta(t,b)");

  hB800_qW_eta->GetYaxis()->SetLabelFont(42);
  hB800_qW_eta->GetYaxis()->SetLabelSize(0.04);
  hB800_qW_eta->GetYaxis()->SetTitleFont(42);
  hB800_qW_eta->GetYaxis()->SetTitleSize(0.06);
  hB800_qW_eta->GetYaxis()->SetTitleOffset(1.0);
  hB800_qW_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hB800_qW_eta->SetMaximum(0.08);
  hB800_qW_eta->SetLineStyle(1);
  hB800_qW_eta->SetLineColor(4);
  hB800_qW_eta->Draw();	

  hB800_qZ_eta->SetLineColor(4);
  hB800_qZ_eta->SetLineStyle(2);
  hB800_qZ_eta->Draw("sames");	

  hB800_qH_eta->SetLineColor(4);
  hB800_qH_eta->SetLineStyle(5);
  hB800_qH_eta->Draw("sames");

  // **************

  TCanvas *c6 = new TCanvas("c6","6: Event Type",1000,700);
  c6->SetFillColor(10);
  c6->Divide(2,2);

  TPad *pad61 = (TCanvas*)c6->GetPad(1);
  pad61->cd(); 
  pad61->SetTopMargin(0.06);
  pad61->SetBottomMargin(0.14);
  pad61->SetLeftMargin(0.12);
  pad61->SetRightMargin(0.04);

  hT700_evtype->GetXaxis()->SetLabelFont(42);
  hT700_evtype->GetXaxis()->SetLabelSize(0.06);
  hT700_evtype->GetXaxis()->SetTitleFont(42);
  hT700_evtype->GetXaxis()->SetTitleSize(0.06);
  hT700_evtype->GetXaxis()->SetTitle("Event Type");

  hT700_evtype->GetXaxis()->SetBinLabel(1,"WbWb");
  hT700_evtype->GetXaxis()->SetBinLabel(2,"ZtZt");
  hT700_evtype->GetXaxis()->SetBinLabel(3,"HtHt");
  hT700_evtype->GetXaxis()->SetBinLabel(4,"WbZt");
  hT700_evtype->GetXaxis()->SetBinLabel(5,"WbHt");
  hT700_evtype->GetXaxis()->SetBinLabel(6,"ZtHt");

  hT700_evtype->GetYaxis()->SetLabelFont(42);
  hT700_evtype->GetYaxis()->SetLabelSize(0.04);
  hT700_evtype->GetYaxis()->SetTitleFont(42);
  hT700_evtype->GetYaxis()->SetTitleSize(0.06);
  hT700_evtype->GetYaxis()->SetTitleOffset(0.9);
  hT700_evtype->GetYaxis()->SetTitle("Fraction of Entries / Type");

  hT700_evtype->SetMaximum(0.5);
  hT700_evtype->SetMinimum(0.0);
  hT700_evtype->SetLineStyle(1);
  hT700_evtype->SetLineColor(2);
  hT700_evtype->Draw();	

  hT800_evtype->SetLineColor(2);
  hT800_evtype->SetLineStyle(2);
  hT800_evtype->Draw("sames");	

  hT900_evtype->SetLineColor(2);
  hT900_evtype->SetLineStyle(5);
  hT900_evtype->Draw("sames");

  TLegend* leg1a = new TLegend(0.30,0.60,0.60,0.90);
  leg1a->SetFillColor(10);
  leg1a->SetBorderSize(0);
  leg1a->SetTextFont(62);
  leg1a->SetTextSize(0.05);
  leg1a->AddEntry(hT700_Q_mass,"#font[42]{m_{T} = 700 GeV}","L");
  leg1a->AddEntry(hT800_Q_mass,"#font[42]{m_{T} = 800 GeV}","L");
  leg1a->AddEntry(hT900_Q_mass,"#font[42]{m_{T} = 900 GeV}","L");
  leg1a->Draw();

  TPad *pad62 = (TCanvas*)c6->GetPad(2);
  pad62->cd(); 
  pad62->SetTopMargin(0.06);
  pad62->SetBottomMargin(0.14);
  pad62->SetLeftMargin(0.12);
  pad62->SetRightMargin(0.04);

  hB700_evtype->GetXaxis()->SetLabelFont(42);
  hB700_evtype->GetXaxis()->SetLabelSize(0.06);
  hB700_evtype->GetXaxis()->SetTitleFont(42);
  hB700_evtype->GetXaxis()->SetTitleSize(0.06);
  hB700_evtype->GetXaxis()->SetTitle("Event Type");

  hB700_evtype->GetXaxis()->SetBinLabel(1,"WtWt");
  hB700_evtype->GetXaxis()->SetBinLabel(2,"ZbZb");
  hB700_evtype->GetXaxis()->SetBinLabel(3,"HbHb");
  hB700_evtype->GetXaxis()->SetBinLabel(4,"WtZb");
  hB700_evtype->GetXaxis()->SetBinLabel(5,"WtHb");
  hB700_evtype->GetXaxis()->SetBinLabel(6,"ZbHb");

  hB700_evtype->GetYaxis()->SetLabelFont(42);
  hB700_evtype->GetYaxis()->SetLabelSize(0.04);
  hB700_evtype->GetYaxis()->SetTitleFont(42);
  hB700_evtype->GetYaxis()->SetTitleSize(0.06);
  hB700_evtype->GetYaxis()->SetTitleOffset(0.9);
  hB700_evtype->GetYaxis()->SetTitle("Fraction of Entries / Type");

  hB700_evtype->SetMaximum(0.5);
  hB700_evtype->SetMinimum(0.0);
  hB700_evtype->SetLineStyle(1);
  hB700_evtype->SetLineColor(4);
  hB700_evtype->Draw();	

  hB800_evtype->SetLineColor(4);
  hB800_evtype->SetLineStyle(2);
  hB800_evtype->Draw("sames");	

  hB900_evtype->SetLineColor(4);
  hB900_evtype->SetLineStyle(5);
  hB900_evtype->Draw("sames");

  TLegend* leg2a = new TLegend(0.30,0.60,0.60,0.90);
  leg2a->SetFillColor(10);
  leg2a->SetBorderSize(0);
  leg2a->SetTextFont(62);
  leg2a->SetTextSize(0.05);
  leg2a->AddEntry(hB700_Q_mass,"#font[42]{m_{B} = 700 GeV}","L");
  leg2a->AddEntry(hB800_Q_mass,"#font[42]{m_{B} = 800 GeV}","L");
  leg2a->AddEntry(hB900_Q_mass,"#font[42]{m_{B} = 900 GeV}","L");
  leg2a->Draw();

  // **************

  TCanvas *c7 = new TCanvas("c7","7: dR boson and quark",1000,700);
  c7->SetFillColor(10);
  c7->Divide(2,2);

  TPad *pad71 = (TCanvas*)c7->GetPad(1);
  pad71->cd(); 
  pad71->SetTopMargin(0.06);
  pad71->SetBottomMargin(0.14);
  pad71->SetLeftMargin(0.12);
  pad71->SetRightMargin(0.04);

  hT700_qW_dR->GetXaxis()->SetLabelFont(42);
  hT700_qW_dR->GetXaxis()->SetLabelSize(0.05);
  hT700_qW_dR->GetXaxis()->SetTitleFont(42);
  hT700_qW_dR->GetXaxis()->SetTitleSize(0.06);
  hT700_qW_dR->GetXaxis()->SetTitle("#DeltaR(W,b)");

  hT700_qW_dR->GetYaxis()->SetLabelFont(42);
  hT700_qW_dR->GetYaxis()->SetLabelSize(0.04);
  hT700_qW_dR->GetYaxis()->SetTitleFont(42);
  hT700_qW_dR->GetYaxis()->SetTitleSize(0.06);
  hT700_qW_dR->GetYaxis()->SetTitleOffset(1.0);
  hT700_qW_dR->GetYaxis()->SetTitle("Fraction of Entries / 0.1");

  hT700_qW_dR->SetMaximum(0.08);
  hT700_qW_dR->SetLineStyle(1);
  hT700_qW_dR->SetLineColor(2);
  hT700_qW_dR->Draw();	

  hT800_qW_dR->SetLineStyle(2);
  hT800_qW_dR->SetLineColor(2);
  hT800_qW_dR->Draw("sames");	

  hT900_qW_dR->SetLineStyle(5);
  hT900_qW_dR->SetLineColor(2);
  hT900_qW_dR->Draw("sames");	

  leg1->Draw();

  TLatex* txt=new TLatex(0.65, 0.35, "T #rightarrow Wb");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.06);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  TPad *pad72 = (TCanvas*)c7->GetPad(2);
  pad72->cd(); 
  pad72->SetTopMargin(0.06);
  pad72->SetBottomMargin(0.14);
  pad72->SetLeftMargin(0.12);
  pad72->SetRightMargin(0.04);

  hB700_qW_dR->GetXaxis()->SetLabelFont(42);
  hB700_qW_dR->GetXaxis()->SetLabelSize(0.05);
  hB700_qW_dR->GetXaxis()->SetTitleFont(42);
  hB700_qW_dR->GetXaxis()->SetTitleSize(0.06);
  hB700_qW_dR->GetXaxis()->SetTitle("#DeltaR(W,t)");

  hB700_qW_dR->GetYaxis()->SetLabelFont(42);
  hB700_qW_dR->GetYaxis()->SetLabelSize(0.04);
  hB700_qW_dR->GetYaxis()->SetTitleFont(42);
  hB700_qW_dR->GetYaxis()->SetTitleSize(0.06);
  hB700_qW_dR->GetYaxis()->SetTitleOffset(1.0);
  hB700_qW_dR->GetYaxis()->SetTitle("Fraction of Entries / 0.1");

  hB700_qW_dR->SetMaximum(0.08);
  hB700_qW_dR->SetLineStyle(1);
  hB700_qW_dR->SetLineColor(4);
  hB700_qW_dR->Draw();

  hB800_qW_dR->SetLineStyle(4);
  hB800_qW_dR->SetLineColor(4);
  hB800_qW_dR->Draw("sames");

  hB900_qW_dR->SetLineStyle(5);
  hB900_qW_dR->SetLineColor(4);
  hB900_qW_dR->Draw("sames");	

  TLatex* txt=new TLatex(0.65, 0.35, "B #rightarrow Wt");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.06);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  leg2->Draw();

  TPad *pad73 = (TCanvas*)c7->GetPad(3);
  pad73->cd(); 
  pad73->SetTopMargin(0.06);
  pad73->SetBottomMargin(0.14);
  pad73->SetLeftMargin(0.12);
  pad73->SetRightMargin(0.04);

  TH1F* hcT800_qW_dR = hT800_qW_dR->Clone("hcT800_qW_dR");
  hcT800_qW_dR->GetXaxis()->SetLabelFont(42);
  hcT800_qW_dR->GetXaxis()->SetLabelSize(0.05);
  hcT800_qW_dR->GetXaxis()->SetTitleFont(42);
  hcT800_qW_dR->GetXaxis()->SetTitleSize(0.06);
  hcT800_qW_dR->GetXaxis()->SetTitle("#DeltaR(W,b;Z,t;H,t)");

  hcT800_qW_dR->GetYaxis()->SetLabelFont(42);
  hcT800_qW_dR->GetYaxis()->SetLabelSize(0.04);
  hcT800_qW_dR->GetYaxis()->SetTitleFont(42);
  hcT800_qW_dR->GetYaxis()->SetTitleSize(0.06);
  hcT800_qW_dR->GetYaxis()->SetTitleOffset(1.0);
  hcT800_qW_dR->GetYaxis()->SetTitle("Fraction of Entries / 0.1");

  hcT800_qW_dR->SetMaximum(0.08);
  hcT800_qW_dR->SetLineStyle(1);
  hcT800_qW_dR->SetLineColor(2);
  hcT800_qW_dR->Draw();	

  hT800_qZ_dR->SetLineColor(2);
  hT800_qZ_dR->SetLineStyle(2);
  hT800_qZ_dR->Draw("sames");	

  hT800_qH_dR->SetLineColor(2);
  hT800_qH_dR->SetLineStyle(5);
  hT800_qH_dR->Draw("sames");

  leg3->Draw();

  TLatex* txt=new TLatex(0.65, 0.35, "m_{T} = 800 GeV");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.06);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  TPad *pad74 = (TCanvas*)c7->GetPad(4);
  pad74->cd(); 
  pad74->SetTopMargin(0.06);
  pad74->SetBottomMargin(0.14);
  pad74->SetLeftMargin(0.12);
  pad74->SetRightMargin(0.04);

  TH1F* hcB800_qW_dR = hB800_qW_dR->Clone("hcB800_qW_dR");
  hcB800_qW_dR->GetXaxis()->SetLabelFont(42);
  hcB800_qW_dR->GetXaxis()->SetLabelSize(0.05);
  hcB800_qW_dR->GetXaxis()->SetTitleFont(42);
  hcB800_qW_dR->GetXaxis()->SetTitleSize(0.06);
  hcB800_qW_dR->GetXaxis()->SetTitle("#DeltaR(W,t;Z,b;H,b)");

  hcB800_qW_dR->GetYaxis()->SetLabelFont(42);
  hcB800_qW_dR->GetYaxis()->SetLabelSize(0.04);
  hcB800_qW_dR->GetYaxis()->SetTitleFont(42);
  hcB800_qW_dR->GetYaxis()->SetTitleSize(0.06);
  hcB800_qW_dR->GetYaxis()->SetTitleOffset(1.0);
  hcB800_qW_dR->GetYaxis()->SetTitle("Fraction of Entries / 0.1");

  hcB800_qW_dR->SetMaximum(0.08);
  hcB800_qW_dR->SetLineStyle(1);
  hcB800_qW_dR->SetLineColor(4);
  hcB800_qW_dR->Draw();	

  hB800_qZ_dR->SetLineColor(4);
  hB800_qZ_dR->SetLineStyle(2);
  hB800_qZ_dR->Draw("sames");	

  hB800_qH_dR->SetLineColor(4);
  hB800_qH_dR->SetLineStyle(5);
  hB800_qH_dR->Draw("sames");

  leg4->Draw();

  TLatex* txt=new TLatex(0.65, 0.35, "m_{T} = 800 GeV");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.06);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  // **************

  TCanvas *c8 = new TCanvas("c8","8: Lepton Multiplicity",1000,700);
  c8->SetFillColor(10);
  c8->Divide(2,2);

  TPad *pad81 = (TCanvas*)c8->GetPad(1);
  pad81->cd(); 
  pad81->SetTopMargin(0.06);
  pad81->SetBottomMargin(0.14);
  pad81->SetLeftMargin(0.12);
  pad81->SetRightMargin(0.04);
  pad81->SetLogy();

  hT700_lep_n->GetXaxis()->SetLabelFont(42);
  hT700_lep_n->GetXaxis()->SetLabelSize(0.05);
  hT700_lep_n->GetXaxis()->SetTitleFont(42);
  hT700_lep_n->GetXaxis()->SetTitleSize(0.06);
  hT700_lep_n->GetXaxis()->SetTitle("# leptons (e,#mu) from W or Z");

  hT700_lep_n->GetYaxis()->SetLabelFont(42);
  hT700_lep_n->GetYaxis()->SetLabelSize(0.04);
  hT700_lep_n->GetYaxis()->SetTitleFont(42);
  hT700_lep_n->GetYaxis()->SetTitleSize(0.06);
  hT700_lep_n->GetYaxis()->SetTitleOffset(1.0);
  hT700_lep_n->GetYaxis()->SetTitle("Fraction of Entries / N_{lep}");

  hT700_lep_n->SetMaximum(1.0);
  hT700_lep_n->SetLineStyle(1);
  hT700_lep_n->SetLineColor(2);
  hT700_lep_n->Draw();

  hT800_lep_n->SetLineStyle(4);
  hT800_lep_n->SetLineColor(2);
  hT800_lep_n->Draw("sames");

  hT900_lep_n->SetLineStyle(5);
  hT900_lep_n->SetLineColor(2);
  hT900_lep_n->Draw("sames");

  TLegend* leg1a = new TLegend(0.60,0.70,0.90,0.90);
  leg1a->SetFillColor(10);
  leg1a->SetBorderSize(0);
  leg1a->SetTextFont(62);
  leg1a->SetTextSize(0.05);
  leg1a->AddEntry(hT700_Q_mass,"#font[42]{m_{T} = 700 GeV}","L");
  leg1a->AddEntry(hT800_Q_mass,"#font[42]{m_{T} = 800 GeV}","L");
  leg1a->AddEntry(hT900_Q_mass,"#font[42]{m_{T} = 900 GeV}","L");
  leg1a->Draw();


  TPad *pad82 = (TCanvas*)c8->GetPad(2);
  pad82->cd(); 
  pad82->SetTopMargin(0.06);
  pad82->SetBottomMargin(0.14);
  pad82->SetLeftMargin(0.12);
  pad82->SetRightMargin(0.04);
  pad82->SetLogy();

  hB700_lep_n->GetXaxis()->SetLabelFont(42);
  hB700_lep_n->GetXaxis()->SetLabelSize(0.05);
  hB700_lep_n->GetXaxis()->SetTitleFont(42);
  hB700_lep_n->GetXaxis()->SetTitleSize(0.06);
  hB700_lep_n->GetXaxis()->SetTitle("# leptons (e,#mu) from W or Z");

  hB700_lep_n->GetYaxis()->SetLabelFont(42);
  hB700_lep_n->GetYaxis()->SetLabelSize(0.04);
  hB700_lep_n->GetYaxis()->SetTitleFont(42);
  hB700_lep_n->GetYaxis()->SetTitleSize(0.06);
  hB700_lep_n->GetYaxis()->SetTitleOffset(1.0);
  hB700_lep_n->GetYaxis()->SetTitle("Fraction of Entries / N_{lep}");

  hB700_lep_n->SetMaximum(1.0);
  hB700_lep_n->SetLineStyle(1);
  hB700_lep_n->SetLineColor(4);
  hB700_lep_n->Draw();

  hB800_lep_n->SetLineStyle(4);
  hB800_lep_n->SetLineColor(4);
  hB800_lep_n->Draw("sames");

  hB900_lep_n->SetLineStyle(5);
  hB900_lep_n->SetLineColor(4);
  hB900_lep_n->Draw("sames");

  TLegend* leg2a = new TLegend(0.60,0.70,0.90,0.90);
  leg2a->SetFillColor(10);
  leg2a->SetBorderSize(0);
  leg2a->SetTextFont(62);
  leg2a->SetTextSize(0.05);
  leg2a->AddEntry(hB700_Q_mass,"#font[42]{m_{B} = 700 GeV}","L");
  leg2a->AddEntry(hB800_Q_mass,"#font[42]{m_{B} = 800 GeV}","L");
  leg2a->AddEntry(hB900_Q_mass,"#font[42]{m_{B} = 900 GeV}","L");
  leg2a->Draw();

  TPad *pad83 = (TCanvas*)c8->GetPad(3);
  pad83->cd(); 
  pad83->SetTopMargin(0.06);
  pad83->SetBottomMargin(0.14);
  pad83->SetLeftMargin(0.12);
  pad83->SetRightMargin(0.04);
  pad83->SetLogy();

  hT800_lep_n_WW->GetXaxis()->SetLabelFont(42);
  hT800_lep_n_WW->GetXaxis()->SetLabelSize(0.05);
  hT800_lep_n_WW->GetXaxis()->SetTitleFont(42);
  hT800_lep_n_WW->GetXaxis()->SetTitleSize(0.06);
  hT800_lep_n_WW->GetXaxis()->SetTitle("# leptons (e,#mu) from W or Z");

  hT800_lep_n_WW->GetYaxis()->SetLabelFont(42);
  hT800_lep_n_WW->GetYaxis()->SetLabelSize(0.04);
  hT800_lep_n_WW->GetYaxis()->SetTitleFont(42);
  hT800_lep_n_WW->GetYaxis()->SetTitleSize(0.06);
  hT800_lep_n_WW->GetYaxis()->SetTitleOffset(1.0);
  hT800_lep_n_WW->GetYaxis()->SetTitle("Fraction of Entries / N_{lep}");

  hT800_lep_n_WW->SetMaximum(1.0);
  hT800_lep_n_WW->SetMinimum(0.01);
  hT800_lep_n_WW->SetLineStyle(1);
  hT800_lep_n_WW->SetLineColor(kRed);
  hT800_lep_n_WW->Draw();

  hT800_lep_n_ZZ->SetMaximum(1.0);
  hT800_lep_n_ZZ->SetLineStyle(1);
  hT800_lep_n_ZZ->SetLineColor(kRed+2);
  hT800_lep_n_ZZ->Draw("sames");

  hT800_lep_n_HH->SetMaximum(1.0);
  hT800_lep_n_HH->SetLineStyle(1);
  hT800_lep_n_HH->SetLineColor(kRed-7);
  hT800_lep_n_HH->Draw("sames");

  hT800_lep_n_WZ->SetLineStyle(2);
  hT800_lep_n_WZ->SetLineColor(kRed);
  hT800_lep_n_WZ->Draw("sames");

  hT800_lep_n_WH->SetMaximum(1.0);
  hT800_lep_n_WH->SetLineStyle(2);
  hT800_lep_n_WH->SetLineColor(kRed+2);
  hT800_lep_n_WH->Draw("sames");

  hT800_lep_n_ZH->SetMaximum(1.0);
  hT800_lep_n_ZH->SetLineStyle(2);
  hT800_lep_n_ZH->SetLineColor(kRed-7);
  hT800_lep_n_ZH->Draw("sames");

  TLegend* leg_Tfs = new TLegend(0.67,0.55,0.92,0.90);
  leg_Tfs->SetFillColor(10);
  leg_Tfs->SetBorderSize(0);
  leg_Tfs->SetTextFont(62);
  leg_Tfs->SetTextSize(0.06);
  leg_Tfs->AddEntry(hT800_lep_n_WW,"#font[42]{WbWb}","L");
  leg_Tfs->AddEntry(hT800_lep_n_ZZ,"#font[42]{ZtZt}","L");
  leg_Tfs->AddEntry(hT800_lep_n_HH,"#font[42]{HtHt}","L");
  leg_Tfs->AddEntry(hT800_lep_n_WZ,"#font[42]{WbZt}","L");
  leg_Tfs->AddEntry(hT800_lep_n_WH,"#font[42]{WbHt}","L");
  leg_Tfs->AddEntry(hT800_lep_n_ZH,"#font[42]{ZtHt}","L");
  leg_Tfs->Draw();

  TLatex* txt=new TLatex(0.70, 0.45, "m_{T} = 800 GeV");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.05);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  TPad *pad84 = (TCanvas*)c8->GetPad(4);
  pad84->cd(); 
  pad84->SetTopMargin(0.06);
  pad84->SetBottomMargin(0.14);
  pad84->SetLeftMargin(0.12);
  pad84->SetRightMargin(0.04);
  pad84->SetLogy();

  hB800_lep_n_WW->GetXaxis()->SetLabelFont(42);
  hB800_lep_n_WW->GetXaxis()->SetLabelSize(0.05);
  hB800_lep_n_WW->GetXaxis()->SetTitleFont(42);
  hB800_lep_n_WW->GetXaxis()->SetTitleSize(0.06);
  hB800_lep_n_WW->GetXaxis()->SetTitle("# leptons (e,#mu) from W or Z");

  hB800_lep_n_WW->GetYaxis()->SetLabelFont(42);
  hB800_lep_n_WW->GetYaxis()->SetLabelSize(0.04);
  hB800_lep_n_WW->GetYaxis()->SetTitleFont(42);
  hB800_lep_n_WW->GetYaxis()->SetTitleSize(0.06);
  hB800_lep_n_WW->GetYaxis()->SetTitleOffset(1.0);
  hB800_lep_n_WW->GetYaxis()->SetTitle("Fraction of Entries / N_{lep}");

  hB800_lep_n_WW->SetMaximum(1.0);
  hB800_lep_n_WW->SetMinimum(0.01);
  hB800_lep_n_WW->SetLineStyle(1);
  hB800_lep_n_WW->SetLineColor(kBlue);
  hB800_lep_n_WW->Draw();

  hB800_lep_n_ZZ->SetMaximum(1.0);
  hB800_lep_n_ZZ->SetLineStyle(1);
  hB800_lep_n_ZZ->SetLineColor(kBlue+3);
  hB800_lep_n_ZZ->Draw("sames");

  hB800_lep_n_HH->SetMaximum(1.0);
  hB800_lep_n_HH->SetLineStyle(1);
  hB800_lep_n_HH->SetLineColor(kBlue-9);
  hB800_lep_n_HH->Draw("sames");

  hB800_lep_n_WZ->SetLineStyle(2);
  hB800_lep_n_WZ->SetLineColor(kBlue);
  hB800_lep_n_WZ->Draw("sames");

  hB800_lep_n_WH->SetMaximum(1.0);
  hB800_lep_n_WH->SetLineStyle(2);
  hB800_lep_n_WH->SetLineColor(kBlue+3);
  hB800_lep_n_WH->Draw("sames");

  hB800_lep_n_ZH->SetMaximum(1.0);
  hB800_lep_n_ZH->SetLineStyle(2);
  hB800_lep_n_ZH->SetLineColor(kBlue-9);
  hB800_lep_n_ZH->Draw("sames");

  TLegend* leg_Bfs = new TLegend(0.67,0.55,0.92,0.90);
  leg_Bfs->SetFillColor(10);
  leg_Bfs->SetBorderSize(0);
  leg_Bfs->SetTextFont(62);
  leg_Bfs->SetTextSize(0.06);
  leg_Bfs->AddEntry(hB800_lep_n_WW,"#font[42]{WtWt}","L");
  leg_Bfs->AddEntry(hB800_lep_n_ZZ,"#font[42]{ZbZb}","L");
  leg_Bfs->AddEntry(hB800_lep_n_HH,"#font[42]{HbHb}","L");
  leg_Bfs->AddEntry(hB800_lep_n_WZ,"#font[42]{WtZb}","L");
  leg_Bfs->AddEntry(hB800_lep_n_WH,"#font[42]{WtHb}","L");
  leg_Bfs->AddEntry(hB800_lep_n_ZH,"#font[42]{ZbHb}","L");
  leg_Bfs->Draw();

  TLatex* txt=new TLatex(0.70, 0.45, "m_{B} = 800 GeV");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.05);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  // **************

  TCanvas *c9 = new TCanvas("c9","9: Electron pT and Eta",1000,700);
  c9->SetFillColor(10);
  c9->Divide(2,2);

  TPad *pad91 = (TCanvas*)c9->GetPad(1);
  pad91->cd(); 
  pad91->SetTopMargin(0.06);
  pad91->SetBottomMargin(0.14);
  pad91->SetLeftMargin(0.12);
  pad91->SetRightMargin(0.04);
  
  hT700_el_pt->GetXaxis()->SetLabelFont(42);
  hT700_el_pt->GetXaxis()->SetLabelSize(0.05);
  hT700_el_pt->GetXaxis()->SetTitleFont(42);
  hT700_el_pt->GetXaxis()->SetTitleSize(0.06);
  hT700_el_pt->GetXaxis()->SetTitle("p_{T}(e from W or Z) [GeV]");

  hT700_el_pt->GetYaxis()->SetLabelFont(42);
  hT700_el_pt->GetYaxis()->SetLabelSize(0.04);
  hT700_el_pt->GetYaxis()->SetTitleFont(42);
  hT700_el_pt->GetYaxis()->SetTitleSize(0.06);
  hT700_el_pt->GetYaxis()->SetTitleOffset(1.0);
  hT700_el_pt->GetYaxis()->SetTitle("Fraction of Entries / 10 GeV");

  hT700_el_pt->SetLineStyle(1);
  hT700_el_pt->SetLineColor(2);
  hT700_el_pt->Draw();

  hT800_el_pt->SetLineStyle(4);
  hT800_el_pt->SetLineColor(2);
  hT800_el_pt->Draw("sames");

  hT900_el_pt->SetLineStyle(5);
  hT900_el_pt->SetLineColor(2);
  hT900_el_pt->Draw("sames");

  leg1->Draw();

  TPad *pad92 = (TCanvas*)c9->GetPad(2);
  pad92->cd(); 
  pad92->SetTopMargin(0.06);
  pad92->SetBottomMargin(0.14);
  pad92->SetLeftMargin(0.12);
  pad92->SetRightMargin(0.04);
  
  hB700_el_pt->GetXaxis()->SetLabelFont(42);
  hB700_el_pt->GetXaxis()->SetLabelSize(0.05);
  hB700_el_pt->GetXaxis()->SetTitleFont(42);
  hB700_el_pt->GetXaxis()->SetTitleSize(0.06);
  hB700_el_pt->GetXaxis()->SetTitle("p_{T}(e from W or Z) [GeV]");

  hB700_el_pt->GetYaxis()->SetLabelFont(42);
  hB700_el_pt->GetYaxis()->SetLabelSize(0.04);
  hB700_el_pt->GetYaxis()->SetTitleFont(42);
  hB700_el_pt->GetYaxis()->SetTitleSize(0.06);
  hB700_el_pt->GetYaxis()->SetTitleOffset(1.0);
  hB700_el_pt->GetYaxis()->SetTitle("Fraction of Entries / 10 GeV");

  hB700_el_pt->SetLineStyle(1);
  hB700_el_pt->SetLineColor(4);
  hB700_el_pt->Draw();

  hB800_el_pt->SetLineStyle(4);
  hB800_el_pt->SetLineColor(4);
  hB800_el_pt->Draw("sames");

  hB900_el_pt->SetLineStyle(5);
  hB900_el_pt->SetLineColor(4);
  hB900_el_pt->Draw("sames");

  leg2->Draw();

  TPad *pad93 = (TCanvas*)c9->GetPad(3);
  pad93->cd(); 
  pad93->SetTopMargin(0.06);
  pad93->SetBottomMargin(0.14);
  pad93->SetLeftMargin(0.12);
  pad93->SetRightMargin(0.04);
  
  hT700_el_eta->GetXaxis()->SetLabelFont(42);
  hT700_el_eta->GetXaxis()->SetLabelSize(0.05);
  hT700_el_eta->GetXaxis()->SetTitleFont(42);
  hT700_el_eta->GetXaxis()->SetTitleSize(0.06);
  hT700_el_eta->GetXaxis()->SetTitle("#eta(e from W or Z)");

  hT700_el_eta->GetYaxis()->SetLabelFont(42);
  hT700_el_eta->GetYaxis()->SetLabelSize(0.04);
  hT700_el_eta->GetYaxis()->SetTitleFont(42);
  hT700_el_eta->GetYaxis()->SetTitleSize(0.06);
  hT700_el_eta->GetYaxis()->SetTitleOffset(1.0);
  hT700_el_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hT700_el_eta->SetMaximum(0.08);
  hT700_el_eta->SetLineStyle(1);
  hT700_el_eta->SetLineColor(2);
  hT700_el_eta->Draw();

  hT800_el_eta->SetLineStyle(4);
  hT800_el_eta->SetLineColor(2);
  hT800_el_eta->Draw("sames");

  hT900_el_eta->SetLineStyle(5);
  hT900_el_eta->SetLineColor(2);
  hT900_el_eta->Draw("sames");

  TPad *pad94 = (TCanvas*)c9->GetPad(4);
  pad94->cd(); 
  pad94->SetTopMargin(0.06);
  pad94->SetBottomMargin(0.14);
  pad94->SetLeftMargin(0.12);
  pad94->SetRightMargin(0.04);
  
  hB700_el_eta->GetXaxis()->SetLabelFont(42);
  hB700_el_eta->GetXaxis()->SetLabelSize(0.05);
  hB700_el_eta->GetXaxis()->SetTitleFont(42);
  hB700_el_eta->GetXaxis()->SetTitleSize(0.06);
  hB700_el_eta->GetXaxis()->SetTitle("#eta(e from W or Z)");

  hB700_el_eta->GetYaxis()->SetLabelFont(42);
  hB700_el_eta->GetYaxis()->SetLabelSize(0.04);
  hB700_el_eta->GetYaxis()->SetTitleFont(42);
  hB700_el_eta->GetYaxis()->SetTitleSize(0.06);
  hB700_el_eta->GetYaxis()->SetTitleOffset(1.0);
  hB700_el_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hB700_el_eta->SetMaximum(0.08);
  hB700_el_eta->SetLineStyle(1);
  hB700_el_eta->SetLineColor(4);
  hB700_el_eta->Draw();

  hB800_el_eta->SetLineStyle(4);
  hB800_el_eta->SetLineColor(4);
  hB800_el_eta->Draw("sames");

  hB900_el_eta->SetLineStyle(5);
  hB900_el_eta->SetLineColor(4);
  hB900_el_eta->Draw("sames");

  // **************

  TCanvas *c10 = new TCanvas("c10","10: Muon pT and Eta",1000,700);
  c10->SetFillColor(10);
  c10->Divide(2,2);

  TPad *pad101 = (TCanvas*)c10->GetPad(1);
  pad101->cd(); 
  pad101->SetTopMargin(0.06);
  pad101->SetBottomMargin(0.14);
  pad101->SetLeftMargin(0.12);
  pad101->SetRightMargin(0.04);
  
  hT700_mu_pt->GetXaxis()->SetLabelFont(42);
  hT700_mu_pt->GetXaxis()->SetLabelSize(0.05);
  hT700_mu_pt->GetXaxis()->SetTitleFont(42);
  hT700_mu_pt->GetXaxis()->SetTitleSize(0.06);
  hT700_mu_pt->GetXaxis()->SetTitle("p_{T}(#mu from W or Z) [GeV]");

  hT700_mu_pt->GetYaxis()->SetLabelFont(42);
  hT700_mu_pt->GetYaxis()->SetLabelSize(0.04);
  hT700_mu_pt->GetYaxis()->SetTitleFont(42);
  hT700_mu_pt->GetYaxis()->SetTitleSize(0.06);
  hT700_mu_pt->GetYaxis()->SetTitleOffset(1.0);
  hT700_mu_pt->GetYaxis()->SetTitle("Fraction of Entries / 10 GeV");

  hT700_mu_pt->SetLineStyle(1);
  hT700_mu_pt->SetLineColor(2);
  hT700_mu_pt->Draw();

  hT800_mu_pt->SetLineStyle(4);
  hT800_mu_pt->SetLineColor(2);
  hT800_mu_pt->Draw("sames");

  hT900_mu_pt->SetLineStyle(5);
  hT900_mu_pt->SetLineColor(2);
  hT900_mu_pt->Draw("sames");

  leg1->Draw();

  TPad *pad102 = (TCanvas*)c10->GetPad(2);
  pad102->cd(); 
  pad102->SetTopMargin(0.06);
  pad102->SetBottomMargin(0.14);
  pad102->SetLeftMargin(0.12);
  pad102->SetRightMargin(0.04);
  
  hB700_mu_pt->GetXaxis()->SetLabelFont(42);
  hB700_mu_pt->GetXaxis()->SetLabelSize(0.05);
  hB700_mu_pt->GetXaxis()->SetTitleFont(42);
  hB700_mu_pt->GetXaxis()->SetTitleSize(0.06);
  hB700_mu_pt->GetXaxis()->SetTitle("p_{T}(#mu from W or Z) [GeV]");

  hB700_mu_pt->GetYaxis()->SetLabelFont(42);
  hB700_mu_pt->GetYaxis()->SetLabelSize(0.04);
  hB700_mu_pt->GetYaxis()->SetTitleFont(42);
  hB700_mu_pt->GetYaxis()->SetTitleSize(0.06);
  hB700_mu_pt->GetYaxis()->SetTitleOffset(1.0);
  hB700_mu_pt->GetYaxis()->SetTitle("Fraction of Entries / 10 GeV");

  hB700_mu_pt->SetLineStyle(1);
  hB700_mu_pt->SetLineColor(4);
  hB700_mu_pt->Draw();

  hB800_mu_pt->SetLineStyle(4);
  hB800_mu_pt->SetLineColor(4);
  hB800_mu_pt->Draw("sames");

  hB900_mu_pt->SetLineStyle(5);
  hB900_mu_pt->SetLineColor(4);
  hB900_mu_pt->Draw("sames");

  leg2->Draw();

  TPad *pad103 = (TCanvas*)c10->GetPad(3);
  pad103->cd(); 
  pad103->SetTopMargin(0.06);
  pad103->SetBottomMargin(0.14);
  pad103->SetLeftMargin(0.12);
  pad103->SetRightMargin(0.04);
  
  hT700_mu_eta->GetXaxis()->SetLabelFont(42);
  hT700_mu_eta->GetXaxis()->SetLabelSize(0.05);
  hT700_mu_eta->GetXaxis()->SetTitleFont(42);
  hT700_mu_eta->GetXaxis()->SetTitleSize(0.06);
  hT700_mu_eta->GetXaxis()->SetTitle("#eta(#mu from W or Z)");

  hT700_mu_eta->GetYaxis()->SetLabelFont(42);
  hT700_mu_eta->GetYaxis()->SetLabelSize(0.04);
  hT700_mu_eta->GetYaxis()->SetTitleFont(42);
  hT700_mu_eta->GetYaxis()->SetTitleSize(0.06);
  hT700_mu_eta->GetYaxis()->SetTitleOffset(1.0);
  hT700_mu_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hT700_mu_eta->SetMaximum(0.08);
  hT700_mu_eta->SetLineStyle(1);
  hT700_mu_eta->SetLineColor(2);
  hT700_mu_eta->Draw();

  hT800_mu_eta->SetLineStyle(4);
  hT800_mu_eta->SetLineColor(2);
  hT800_mu_eta->Draw("sames");

  hT900_mu_eta->SetLineStyle(5);
  hT900_mu_eta->SetLineColor(2);
  hT900_mu_eta->Draw("sames");

  TPad *pad104 = (TCanvas*)c10->GetPad(4);
  pad104->cd(); 
  pad104->SetTopMargin(0.06);
  pad104->SetBottomMargin(0.14);
  pad104->SetLeftMargin(0.12);
  pad104->SetRightMargin(0.04);
  
  hB700_mu_eta->GetXaxis()->SetLabelFont(42);
  hB700_mu_eta->GetXaxis()->SetLabelSize(0.05);
  hB700_mu_eta->GetXaxis()->SetTitleFont(42);
  hB700_mu_eta->GetXaxis()->SetTitleSize(0.06);
  hB700_mu_eta->GetXaxis()->SetTitle("#eta(#mu from W or Z)");

  hB700_mu_eta->GetYaxis()->SetLabelFont(42);
  hB700_mu_eta->GetYaxis()->SetLabelSize(0.04);
  hB700_mu_eta->GetYaxis()->SetTitleFont(42);
  hB700_mu_eta->GetYaxis()->SetTitleSize(0.06);
  hB700_mu_eta->GetYaxis()->SetTitleOffset(1.0);
  hB700_mu_eta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hB700_mu_eta->SetMaximum(0.08);
  hB700_mu_eta->SetLineStyle(1);
  hB700_mu_eta->SetLineColor(4);
  hB700_mu_eta->Draw();

  hB800_mu_eta->SetLineStyle(4);
  hB800_mu_eta->SetLineColor(4);
  hB800_mu_eta->Draw("sames");

  hB900_mu_eta->SetLineStyle(5);
  hB900_mu_eta->SetLineColor(4);
  hB900_mu_eta->Draw("sames");

  // **************
  
  TCanvas *c11 = new TCanvas("c11","11: MET",1000,700);
  c11->SetFillColor(10);
  c11->Divide(2,2);

  TPad *pad111 = (TCanvas*)c11->GetPad(1);
  pad111->cd(); 
  pad111->SetTopMargin(0.06);
  pad111->SetBottomMargin(0.14);
  pad111->SetLeftMargin(0.12);
  pad111->SetRightMargin(0.04);
  
  hT700_MET->GetXaxis()->SetLabelFont(42);
  hT700_MET->GetXaxis()->SetLabelSize(0.05);
  hT700_MET->GetXaxis()->SetTitleFont(42);
  hT700_MET->GetXaxis()->SetTitleSize(0.06);
  hT700_MET->GetXaxis()->SetTitle("Missing E_{T} [GeV]");

  hT700_MET->GetYaxis()->SetLabelFont(42);
  hT700_MET->GetYaxis()->SetLabelSize(0.04);
  hT700_MET->GetYaxis()->SetTitleFont(42);
  hT700_MET->GetYaxis()->SetTitleSize(0.06);
  hT700_MET->GetYaxis()->SetTitleOffset(1.0);
  hT700_MET->GetYaxis()->SetTitle("Fraction of Entries / 10 GeV");

  hT700_MET->SetMaximum(0.20);
  hT700_MET->SetLineStyle(1);
  hT700_MET->SetLineColor(2);
  hT700_MET->Draw();

  hT800_MET->SetLineStyle(4);
  hT800_MET->SetLineColor(2);
  hT800_MET->Draw("sames");

  hT900_MET->SetLineStyle(5);
  hT900_MET->SetLineColor(2);
  hT900_MET->Draw("sames");

  leg1->Draw();

  TPad *pad112 = (TCanvas*)c11->GetPad(2);
  pad112->cd(); 
  pad112->SetTopMargin(0.06);
  pad112->SetBottomMargin(0.14);
  pad112->SetLeftMargin(0.12);
  pad112->SetRightMargin(0.04);
  
  hB700_MET->GetXaxis()->SetLabelFont(42);
  hB700_MET->GetXaxis()->SetLabelSize(0.05);
  hB700_MET->GetXaxis()->SetTitleFont(42);
  hB700_MET->GetXaxis()->SetTitleSize(0.06);
  hB700_MET->GetXaxis()->SetTitle("Missing E_{T} [GeV]");

  hB700_MET->GetYaxis()->SetLabelFont(42);
  hB700_MET->GetYaxis()->SetLabelSize(0.04);
  hB700_MET->GetYaxis()->SetTitleFont(42);
  hB700_MET->GetYaxis()->SetTitleSize(0.06);
  hB700_MET->GetYaxis()->SetTitleOffset(1.0);
  hB700_MET->GetYaxis()->SetTitle("Fraction of Entries / 10 GeV");

  hB700_MET->SetMaximum(0.20);
  hB700_MET->SetLineStyle(1);
  hB700_MET->SetLineColor(4);
  hB700_MET->Draw();

  hB800_MET->SetLineStyle(4);
  hB800_MET->SetLineColor(4);
  hB800_MET->Draw("sames");

  hB900_MET->SetLineStyle(5);
  hB900_MET->SetLineColor(4);
  hB900_MET->Draw("sames");

  leg2->Draw();


  TPad *pad113 = (TCanvas*)c11->GetPad(3);
  pad113->cd(); 
  pad113->SetTopMargin(0.06);
  pad113->SetBottomMargin(0.14);
  pad113->SetLeftMargin(0.12);
  pad113->SetRightMargin(0.04);

  hT800_MET_WW->GetXaxis()->SetLabelFont(42);
  hT800_MET_WW->GetXaxis()->SetLabelSize(0.05);
  hT800_MET_WW->GetXaxis()->SetTitleFont(42);
  hT800_MET_WW->GetXaxis()->SetTitleSize(0.06);
  hT800_MET_WW->GetXaxis()->SetTitle("Missing E_{T} [GeV]");

  hT800_MET_WW->GetYaxis()->SetLabelFont(42);
  hT800_MET_WW->GetYaxis()->SetLabelSize(0.04);
  hT800_MET_WW->GetYaxis()->SetTitleFont(42);
  hT800_MET_WW->GetYaxis()->SetTitleSize(0.06);
  hT800_MET_WW->GetYaxis()->SetTitleOffset(1.0);
  hT800_MET_WW->GetYaxis()->SetTitle("Fraction of Entries / 10 GeV");

  hT800_MET_WW->SetMaximum(0.21);
  hT800_MET_WW->SetLineStyle(1);
  hT800_MET_WW->SetLineColor(kRed);
  hT800_MET_WW->Draw();

  hT800_MET_ZZ->SetLineStyle(1);
  hT800_MET_ZZ->SetLineColor(kRed+2);
  hT800_MET_ZZ->Draw("sames");

  hT800_MET_HH->SetLineStyle(1);
  hT800_MET_HH->SetLineColor(kRed-7);
  hT800_MET_HH->Draw("sames");

  hT800_MET_WZ->SetLineStyle(2);
  hT800_MET_WZ->SetLineColor(kRed);
  hT800_MET_WZ->Draw("sames");

  hT800_MET_WH->SetLineStyle(2);
  hT800_MET_WH->SetLineColor(kRed+2);
  hT800_MET_WH->Draw("sames");

  hT800_MET_ZH->SetLineStyle(2);
  hT800_MET_ZH->SetLineColor(kRed-7);
  hT800_MET_ZH->Draw("sames");

  leg_Tfs->Draw();

  TLatex* txt=new TLatex(0.70, 0.45, "m_{T} = 800 GeV");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.05);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  TPad *pad114 = (TCanvas*)c11->GetPad(4);
  pad114->cd(); 
  pad114->SetTopMargin(0.06);
  pad114->SetBottomMargin(0.14);
  pad114->SetLeftMargin(0.12);
  pad114->SetRightMargin(0.04);

  hB800_MET_WW->GetXaxis()->SetLabelFont(42);
  hB800_MET_WW->GetXaxis()->SetLabelSize(0.05);
  hB800_MET_WW->GetXaxis()->SetTitleFont(42);
  hB800_MET_WW->GetXaxis()->SetTitleSize(0.06);
  hB800_MET_WW->GetXaxis()->SetTitle("Missing E_{T} [GeV]");

  hB800_MET_WW->GetYaxis()->SetLabelFont(42);
  hB800_MET_WW->GetYaxis()->SetLabelSize(0.04);
  hB800_MET_WW->GetYaxis()->SetTitleFont(42);
  hB800_MET_WW->GetYaxis()->SetTitleSize(0.06);
  hB800_MET_WW->GetYaxis()->SetTitleOffset(1.0);
  hB800_MET_WW->GetYaxis()->SetTitle("Fraction of Entries / 10 GeV");

  hB800_MET_WW->SetMaximum(0.21);
  hB800_MET_WW->SetLineStyle(1);
  hB800_MET_WW->SetLineColor(kBlue);
  hB800_MET_WW->Draw();

  hB800_MET_ZZ->SetLineStyle(1);
  hB800_MET_ZZ->SetLineColor(kBlue+3);
  hB800_MET_ZZ->Draw("sames");

  hB800_MET_HH->SetLineStyle(1);
  hB800_MET_HH->SetLineColor(kBlue-9);
  hB800_MET_HH->Draw("sames");

  hB800_MET_WZ->SetLineStyle(2);
  hB800_MET_WZ->SetLineColor(kBlue);
  hB800_MET_WZ->Draw("sames");

  hB800_MET_WH->SetLineStyle(2);
  hB800_MET_WH->SetLineColor(kBlue+3);
  hB800_MET_WH->Draw("sames");

  hB800_MET_ZH->SetLineStyle(2);
  hB800_MET_ZH->SetLineColor(kBlue-9);
  hB800_MET_ZH->Draw("sames");

  leg_Bfs->Draw();

  TLatex* txt=new TLatex(0.70, 0.45, "m_{B} = 800 GeV");
  txt->SetNDC(kTRUE); // <- use NDC coordinate
  txt->SetTextSize(0.05);
  txt->SetTextFont(42);
  txt->SetTextColor(1);
  txt->Draw();

  // **************
  
  TCanvas *c12 = new TCanvas("c12","12: Njet and Jet mass",1000,700);
  c12->SetFillColor(10);
  c12->Divide(2,2);

  TPad *pad121 = (TCanvas*)c12->GetPad(1);
  pad121->cd(); 
  pad121->SetTopMargin(0.06);
  pad121->SetBottomMargin(0.14);
  pad121->SetLeftMargin(0.12);
  pad121->SetRightMargin(0.04);

  hT700_njet->GetXaxis()->SetLabelFont(42);
  hT700_njet->GetXaxis()->SetLabelSize(0.05);
  hT700_njet->GetXaxis()->SetTitleFont(42);
  hT700_njet->GetXaxis()->SetTitleSize(0.06);
  hT700_njet->GetXaxis()->SetTitle("Jet Multiplicity (R=0.4)");

  hT700_njet->GetYaxis()->SetLabelFont(42);
  hT700_njet->GetYaxis()->SetLabelSize(0.04);
  hT700_njet->GetYaxis()->SetTitleFont(42);
  hT700_njet->GetYaxis()->SetTitleSize(0.06);
  hT700_njet->GetYaxis()->SetTitleOffset(1.0);
  hT700_njet->GetYaxis()->SetTitle("Fraction of Entries / N_{jet}");

  hT700_njet->SetMaximum(0.22);
  hT700_njet->SetLineStyle(1);
  hT700_njet->SetLineColor(2);
  hT700_njet->Draw();

  hT800_njet->SetLineStyle(4);
  hT800_njet->SetLineColor(2);
  hT800_njet->Draw("sames");

  hT900_njet->SetLineStyle(5);
  hT900_njet->SetLineColor(2);
  hT900_njet->Draw("sames");

  leg1->Draw();

  TPad *pad122 = (TCanvas*)c12->GetPad(2);
  pad122->cd(); 
  pad122->SetTopMargin(0.06);
  pad122->SetBottomMargin(0.14);
  pad122->SetLeftMargin(0.12);
  pad122->SetRightMargin(0.04);

  hB700_njet->SetMaximum(0.22);
  hB700_njet->GetXaxis()->SetLabelFont(42);
  hB700_njet->GetXaxis()->SetLabelSize(0.05);
  hB700_njet->GetXaxis()->SetTitleFont(42);
  hB700_njet->GetXaxis()->SetTitleSize(0.06);
  hB700_njet->GetXaxis()->SetTitle("Jet Multiplicity (R=0.4)");

  hB700_njet->GetYaxis()->SetLabelFont(42);
  hB700_njet->GetYaxis()->SetLabelSize(0.04);
  hB700_njet->GetYaxis()->SetTitleFont(42);
  hB700_njet->GetYaxis()->SetTitleSize(0.06);
  hB700_njet->GetYaxis()->SetTitleOffset(1.0);
  hB700_njet->GetYaxis()->SetTitle("Fraction of Entries / N_{jet}");

  hB700_njet->SetLineStyle(1);
  hB700_njet->SetLineColor(4);
  hB700_njet->Draw();

  hB800_njet->SetLineStyle(4);
  hB800_njet->SetLineColor(4);
  hB800_njet->Draw("sames");

  hB900_njet->SetLineStyle(5);
  hB900_njet->SetLineColor(4);
  hB900_njet->Draw("sames");

  leg2->Draw();

  TPad *pad123 = (TCanvas*)c12->GetPad(3);
  pad123->cd(); 
  pad123->SetTopMargin(0.06);
  pad123->SetBottomMargin(0.14);
  pad123->SetLeftMargin(0.12);
  pad123->SetRightMargin(0.04);
  pad123->SetLogy();

  hT700_jmass->GetXaxis()->SetLabelFont(42);
  hT700_jmass->GetXaxis()->SetLabelSize(0.05);
  hT700_jmass->GetXaxis()->SetTitleFont(42);
  hT700_jmass->GetXaxis()->SetTitleSize(0.06);
  hT700_jmass->GetXaxis()->SetTitle("Jet Mass (R=0.4) [GeV]");

  hT700_jmass->GetYaxis()->SetLabelFont(42);
  hT700_jmass->GetYaxis()->SetLabelSize(0.04);
  hT700_jmass->GetYaxis()->SetTitleFont(42);
  hT700_jmass->GetYaxis()->SetTitleSize(0.06);
  hT700_jmass->GetYaxis()->SetTitleOffset(1.0);
  hT700_jmass->GetYaxis()->SetTitle("Fraction of Entries / 5 GeV");

  hT700_jmass->SetMaximum(0.50);
  hT700_jmass->SetLineStyle(1);
  hT700_jmass->SetLineColor(2);
  hT700_jmass->Draw();

  hT800_jmass->SetLineStyle(4);
  hT800_jmass->SetLineColor(2);
  hT800_jmass->Draw("sames");

  hT900_jmass->SetLineStyle(5);
  hT900_jmass->SetLineColor(2);
  hT900_jmass->Draw("sames");
  
  TPad *pad124 = (TCanvas*)c12->GetPad(4);
  pad124->cd(); 
  pad124->SetTopMargin(0.06);
  pad124->SetBottomMargin(0.14);
  pad124->SetLeftMargin(0.12);
  pad124->SetRightMargin(0.04);
  pad124->SetLogy();

  hB700_jmass->SetMaximum(0.50);
  hB700_jmass->GetXaxis()->SetLabelFont(42);
  hB700_jmass->GetXaxis()->SetLabelSize(0.05);
  hB700_jmass->GetXaxis()->SetTitleFont(42);
  hB700_jmass->GetXaxis()->SetTitleSize(0.06);
  hB700_jmass->GetXaxis()->SetTitle("Jet Mass (R=0.4) [GeV]");

  hB700_jmass->GetYaxis()->SetLabelFont(42);
  hB700_jmass->GetYaxis()->SetLabelSize(0.04);
  hB700_jmass->GetYaxis()->SetTitleFont(42);
  hB700_jmass->GetYaxis()->SetTitleSize(0.06);
  hB700_jmass->GetYaxis()->SetTitleOffset(1.0);
  hB700_jmass->GetYaxis()->SetTitle("Fraction of Entries / 5 GeV");

  hB700_jmass->SetLineStyle(1);
  hB700_jmass->SetLineColor(4);
  hB700_jmass->Draw();

  hB800_jmass->SetLineStyle(4);
  hB800_jmass->SetLineColor(4);
  hB800_jmass->Draw("sames");

  hB900_jmass->SetLineStyle(5);
  hB900_jmass->SetLineColor(4);
  hB900_jmass->Draw("sames");

  // **************
  
  TCanvas *c13 = new TCanvas("c13","13: jet pT and eta",1000,700);
  c13->SetFillColor(10);
  c13->Divide(2,2);

  TPad *pad131 = (TCanvas*)c13->GetPad(1);
  pad131->cd(); 
  pad131->SetTopMargin(0.06);
  pad131->SetBottomMargin(0.14);
  pad131->SetLeftMargin(0.12);
  pad131->SetRightMargin(0.04);

  hT700_jpt->GetXaxis()->SetLabelFont(42);
  hT700_jpt->GetXaxis()->SetLabelSize(0.05);
  hT700_jpt->GetXaxis()->SetTitleFont(42);
  hT700_jpt->GetXaxis()->SetTitleSize(0.06);
  hT700_jpt->GetXaxis()->SetTitle("Jet p_{T} (R=0.4) [GeV]");

  hT700_jpt->GetYaxis()->SetLabelFont(42);
  hT700_jpt->GetYaxis()->SetLabelSize(0.04);
  hT700_jpt->GetYaxis()->SetTitleFont(42);
  hT700_jpt->GetYaxis()->SetTitleSize(0.06);
  hT700_jpt->GetYaxis()->SetTitleOffset(1.0);
  hT700_jpt->GetYaxis()->SetTitle("Fraction of Entries / 10 GeV");

  hT700_jpt->SetMaximum(0.09);
  hT700_jpt->SetLineStyle(1);
  hT700_jpt->SetLineColor(2);
  hT700_jpt->Draw();

  hT800_jpt->SetLineStyle(4);
  hT800_jpt->SetLineColor(2);
  hT800_jpt->Draw("sames");

  hT900_jpt->SetLineStyle(5);
  hT900_jpt->SetLineColor(2);
  hT900_jpt->Draw("sames");

  leg1->Draw();

  TPad *pad132 = (TCanvas*)c13->GetPad(2);
  pad132->cd(); 
  pad132->SetTopMargin(0.06);
  pad132->SetBottomMargin(0.14);
  pad132->SetLeftMargin(0.12);
  pad132->SetRightMargin(0.04);

  hB700_jpt->GetXaxis()->SetLabelFont(42);
  hB700_jpt->GetXaxis()->SetLabelSize(0.05);
  hB700_jpt->GetXaxis()->SetTitleFont(42);
  hB700_jpt->GetXaxis()->SetTitleSize(0.06);
  hB700_jpt->GetXaxis()->SetTitle("Jet p_{T} (R=0.4) [GeV]");

  hB700_jpt->GetYaxis()->SetLabelFont(42);
  hB700_jpt->GetYaxis()->SetLabelSize(0.04);
  hB700_jpt->GetYaxis()->SetTitleFont(42);
  hB700_jpt->GetYaxis()->SetTitleSize(0.06);
  hB700_jpt->GetYaxis()->SetTitleOffset(1.0);
  hB700_jpt->GetYaxis()->SetTitle("Fraction of Entries / 10 GeV");

  hB700_jpt->SetMaximum(0.09);
  hB700_jpt->SetLineStyle(1);
  hB700_jpt->SetLineColor(4);
  hB700_jpt->Draw();

  hB800_jpt->SetLineStyle(4);
  hB800_jpt->SetLineColor(4);
  hB800_jpt->Draw("sames");

  hB900_jpt->SetLineStyle(5);
  hB900_jpt->SetLineColor(4);
  hB900_jpt->Draw("sames");

  leg2->Draw();

  TPad *pad133 = (TCanvas*)c13->GetPad(3);
  pad133->cd(); 
  pad133->SetTopMargin(0.06);
  pad133->SetBottomMargin(0.14);
  pad133->SetLeftMargin(0.12);
  pad133->SetRightMargin(0.04);

  hT700_jeta->GetXaxis()->SetLabelFont(42);
  hT700_jeta->GetXaxis()->SetLabelSize(0.05);
  hT700_jeta->GetXaxis()->SetTitleFont(42);
  hT700_jeta->GetXaxis()->SetTitleSize(0.06);
  hT700_jeta->GetXaxis()->SetTitle("Jet #eta (R=0.4)");

  hT700_jeta->GetYaxis()->SetLabelFont(42);
  hT700_jeta->GetYaxis()->SetLabelSize(0.04);
  hT700_jeta->GetYaxis()->SetTitleFont(42);
  hT700_jeta->GetYaxis()->SetTitleSize(0.06);
  hT700_jeta->GetYaxis()->SetTitleOffset(1.0);
  hT700_jeta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hT700_jeta->SetMaximum(0.09);
  hT700_jeta->SetLineStyle(1);
  hT700_jeta->SetLineColor(2);
  hT700_jeta->Draw();

  hT800_jeta->SetLineStyle(4);
  hT800_jeta->SetLineColor(2);
  hT800_jeta->Draw("sames");

  hT900_jeta->SetLineStyle(5);
  hT900_jeta->SetLineColor(2);
  hT900_jeta->Draw("sames");

  TPad *pad134 = (TCanvas*)c13->GetPad(4);
  pad134->cd(); 
  pad134->SetTopMargin(0.06);
  pad134->SetBottomMargin(0.14);
  pad134->SetLeftMargin(0.12);
  pad134->SetRightMargin(0.04);

  hB700_jeta->GetXaxis()->SetLabelFont(42);
  hB700_jeta->GetXaxis()->SetLabelSize(0.05);
  hB700_jeta->GetXaxis()->SetTitleFont(42);
  hB700_jeta->GetXaxis()->SetTitleSize(0.06);
  hB700_jeta->GetXaxis()->SetTitle("Jet #eta (R=0.4)");

  hB700_jeta->GetYaxis()->SetLabelFont(42);
  hB700_jeta->GetYaxis()->SetLabelSize(0.04);
  hB700_jeta->GetYaxis()->SetTitleFont(42);
  hB700_jeta->GetYaxis()->SetTitleSize(0.06);
  hB700_jeta->GetYaxis()->SetTitleOffset(1.0);
  hB700_jeta->GetYaxis()->SetTitle("Fraction of Entries / 0.2");

  hB700_jeta->SetMaximum(0.09);
  hB700_jeta->SetLineStyle(1);
  hB700_jeta->SetLineColor(4);
  hB700_jeta->Draw();

  hB800_jeta->SetLineStyle(4);
  hB800_jeta->SetLineColor(4);
  hB800_jeta->Draw("sames");

  hB900_jeta->SetLineStyle(5);
  hB900_jeta->SetLineColor(4);
  hB900_jeta->Draw("sames");

  //************************************

  TCanvas* csmall = new TCanvas("csmall","Small",10,10);
  csmall->SetFillColor(10);
  csmall->cd();

}


// jet n,pT,eta,mass (4/6)

// need dphi(QQ)
// need pT(QQ)
// b-jets
// charge?