source: huonglan/MYSIMULATION/mttggqq.C @ 15

#include <iostream>
#include <cstdio>
#include <math.h>
#include <TROOT.h>
#include <TFile.h>
#include <TTree.h>
#include <TCanvas.h>

void mtt()
{
  gROOT->SetStyle("Plain");
  TCanvas *c1 = new TCanvas("mycanvas1","",900,450);
  c1 -> Divide(2,1);

  TFile *file1 = new TFile("./pyres400_15_6.root");
  TString mttgg = "pymtt_gg";
  TString mttqq = "pymtt_qq";
  TH1F  *hmttgg = (TH1F*)file1-> Get(mttgg);
  TH1F  *hmttqq = (TH1F*)file1-> Get(mttqq);
  hmttgg->Scale(1./(hmttgg->Integral()));
  hmttqq->Scale(1./(hmttqq->Integral()));

  TFile *file2 = new TFile("./res400done_15_6.root");
  TString mttggtr = "truth_sqrtshatGG";
  TString mttqqtr = "truth_sqrtshatQQ";
  TString mttggan = "ana_mttGG";
  TString mttqqan = "ana_mttQQ";
  TH1F  *hmttggtr = (TH1F*)file2-> Get(mttggtr);
  TH1F  *hmttqqtr = (TH1F*)file2-> Get(mttqqtr);
  TH1F  *hmttggan = (TH1F*)file2-> Get(mttggan);
  TH1F  *hmttqqan = (TH1F*)file2-> Get(mttqqan);
  hmttggtr->Scale(1./(hmttggtr->Integral()));
  hmttqqtr->Scale(1./(hmttqqtr->Integral()));
  hmttggan->Scale(1./(hmttggan->Integral()));
  hmttqqan->Scale(1./(hmttqqan->Integral()));

  TLegend *legend1 = new TLegend(0.85,0.6,0.5,0.75);
  legend1 -> AddEntry(hmttgg,"Pythia");
  legend1 -> AddEntry(hmttggtr,"My analysis - truth");
  legend1 -> AddEntry(hmttggan,"My analysis - reco");

  TLegend *legend2 = new TLegend(0.85,0.6,0.5,0.75);
  legend2 -> AddEntry(hmttqq,"Pythia");
  legend2 -> AddEntry(hmttqqtr,"My analysis - truth");
  legend2 -> AddEntry(hmttqqan,"My analysis - reco");

  c1 -> cd(1);
  hmttgg->Draw();
  hmttggtr->SetLineColor(kRed);
  hmttggtr->Draw("same");
  hmttggan->SetLineColor(kBlue);
  hmttggan->Draw("same");
  legend1 -> Draw("same");

  c1 -> cd(2);
  hmttqq->Draw();
  hmttqqtr->SetLineColor(kRed);
  hmttqqtr->Draw("same");
  hmttqqan->SetLineColor(kBlue);
  hmttqqan->Draw("same");
  legend2 -> Draw("same");

  c1 -> Print("mtt_ggqq.ps");
}