#include <TFile.h>
#include <TTree.h>
#include <TROOT.h>
#include <TH1.h>
#include <TH2.h>
#include <TSystem.h>
#include <vector>
#include <TLorentzVector.h>
#include <math.h>
#include "TMath.h"
#include <iostream>
#include <fstream>
#include <algorithm>
#include <map>
#include <TFileMerger.h>

void nt2()

{

   TFile *f1 = new TFile("/afs/cern.ch/work/n/nmajeed/private/CMSSW_7_3_0/my_work/ttbar/tt_final_new.root");
   TH2F *h2 = (TH2F*)gDirectory->Get("third_loose_mu_pt_eta");
   TH2F *h1 = (TH2F*)gDirectory->Get("third_tight_mu_pt_eta");
   TH2F *nt2 = (TH2F*)gDirectory->Get("nt2_pt_eta");

   //TH2F *f = new TH2F("f","ratio",100,0,300,100,-3,3);

  double bins_x[] = {10, 20, 40, 70, 100};
  double bins_y[] = {0, 1, 1.5, 2.5};

   TH2F *h1 = new TH2F("h1", "tight muons",(sizeof(bins_x) / sizeof(double) - 1), bins_x, (sizeof(bins_y) / sizeof(double) - 1), bins_y);
   h1->Sumw2(kTRUE);

   TH2F *h2 = new TH2F("h2", "loose muons",(sizeof(bins_x) / sizeof(double) - 1), bins_x, (sizeof(bins_y) / sizeof(double) - 1), bins_y);
   h2->Sumw2(kTRUE);

   TH2F *nt2 = new TH2F("nt2", "b-jets",(sizeof(bins_x) / sizeof(double) - 1), bins_x, (sizeof(bins_y) / sizeof(double) - 1), bins_y);
   h2->Sumw2(kTRUE);

   TH2F *f = new TH2F("f", "tight/loose",(sizeof(bins_x) / sizeof(double) - 1), bins_x, (sizeof(bins_y) / sizeof(double) - 1), bins_y);
   f->Sumw2(kTRUE);

   TH2F *f = (TH2F*)h1->Clone();
  
   f->Divide(h2);

   f->Draw("colz text e");
   h1->Draw("colz text e");

   for(int i; i<=bins_x; i++) {
      for(int j; j<=bins_y; ij++) {

     double Num  = h1->GetBinContent(i,j);
     double Den    = h2->GetBinContent(i,j);
     if (Num==0 &&  Den==0 )  continue ;
     double ratio = Num / Den ;
     cout << " ratio = " <<  ratio << endl;
     double mins = (1-ratio) ;
     double mul = ratio * mins;

      }
      }
   TH2D *Nt3;
   Nt3->Divide(mul,nt2);

   TCanvas* c1 = new TCanvas("c1","c1",800,600) ;
   c1->Divide(2,2);
   c1->cd(1);
   f->Draw("colz text e");
   c1->cd(2);
   h1->Draw("colz text e");
   c1->cd(3);
   nt2->Draw("colz text e");
   c1->cd(4);
   Nt3->Draw("colz text e");

}