ROOTOutputAnalyzers/makeQCDSuperimposePlot.C at main · emcannaert/ROOTOutputAnalyzers · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
#include <iostream>
#include <string>
#include "TLorentzVector.h"
#include "TRatioPlot.h"
#include "THStack.h"
using namespace std;

void makeQCDSuperimposePlot()
{

  TFile *f1 = new TFile("/home/ethan/Documents/clustAlg_QCD2000toInf_processed_cos70.root");
  TFile *f2 = new TFile("/home/ethan/Documents/clustAlg_QCD2000toInf_processed_cos75.root");
  TFile *f3 = new TFile("/home/ethan/Documents/clustAlg_QCD2000toInf_processed_cos80.root");
  TFile *f4 = new TFile("/home/ethan/Documents/clustAlg_QCD2000toInf_processed_cos85.root");

  TCanvas *c1 = new TCanvas("c1","",400,20, 1500,1500);
  gPad->SetLogy();
  gStyle->SetOptStat(0);

  TH1F *h_mass_cos70 = (TH1F*)f1->Get("h_SJ_mass_100");
  TH1F *h_mass_cos75 = (TH1F*)f2->Get("h_SJ_mass_100");
  TH1F *h_mass_cos80 = (TH1F*)f3->Get("h_SJ_mass_100");
  TH1F *h_mass_cos85 = (TH1F*)f4->Get("h_SJ_mass_100");


  h_mass_cos85->SetTitle("Mass of QCD SuperJets (Different Sorting cos thresholds) [59.83 fb^{-1}]");
  h_mass_cos85->SetLineColor(kBlue);
  h_mass_cos85->SetMarkerColor(kBlue);
  h_mass_cos85->SetLineWidth(3);
  h_mass_cos85->Draw("HIST");

  h_mass_cos80->SetLineColor(kOrange);
  h_mass_cos80->SetMarkerColor(kOrange);
  h_mass_cos80->SetLineWidth(3);
  h_mass_cos80->Draw("HIST same");

  h_mass_cos75->SetLineColor(kRed);
  h_mass_cos75->SetMarkerColor(kRed);
  h_mass_cos75->SetLineWidth(3);
  h_mass_cos75->Draw("HIST same");

  h_mass_cos70->SetLineColor(kGreen);
  h_mass_cos70->SetMarkerColor(kGreen);
  h_mass_cos70->SetLineWidth(3);
  h_mass_cos70->Draw("HIST same");

  TLegend *legend = new TLegend(1.0,0.8,0.48,0.9);
  legend->AddEntry(h_mass_cos70,"QCD MC cos 70","l");
  legend->AddEntry(h_mass_cos75,"QCD MC cos 75","l");
  legend->AddEntry(h_mass_cos80,"QCD MC cos 80","l");
  legend->AddEntry(h_mass_cos85,"QCD MC cos 85","l");
  legend->Draw();

  c1->SaveAs("QCD_SJ_mass_100_superimpose.png");


}