ROOTOutputAnalyzers/divideTH2.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
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
#include <iostream>
#include <string>
#include "TLorentzVector.h"
#include "TRatioPlot.h"
#include "THStack.h"
using namespace std;

void divideTH2()
{

  double QCD1000to1500_SF = 4.289571744;
  double QCD1500to2000_SF = 0.6042726685;
  double QCD2000toInf_SF = 0.2132134533;
  double TTJets_inclusive_SF = 4.742226033;
  double h_WJetsToLNu_SF = 0.02459688845;
  double h_WJet_QQ_SF = 0.1266526072;
  double h_TTToHadronic_SF = 0.06588049107;
  double h_TTTo2l2nu_SF =0.03616639075 ;
  double h_TTToSemiLeptonic_SF = 0.04563489275;
  //double h_TTJets2500toInf_SF  = 0.00008408965681;


  TFile *f1 = new TFile("/home/ethan/Documents/QCD_HT1000to1500_combined_processed.root");
  TFile *f2 = new TFile("/home/ethan/Documents/QCD_HT1500to2000_combined_processed.root");
  TFile *f3 = new TFile("/home/ethan/Documents/QCD_HT2000toInf_combined_processed.root");   //need to add these two together
  TFile *f4 = new TFile("/home/ethan/Documents/TTJets_inclusive_combined_processed.root");
  TFile *f5 = new TFile("/home/ethan/Documents/dataMistagRates_JetHT2018A.root");
  TFile *f6 = new TFile("/home/ethan/Documents/dataMistagRates_JetHT2018B.root");
  TFile *f7 = new TFile("/home/ethan/Documents/dataMistagRates_JetHT2018C.root");
  TFile *f8 = new TFile("/home/ethan/Documents/dataMistagRates_JetHT2018D.root");

  TFile *f9 = new TFile("/home/ethan/Documents/WJetsToLNu_combined_processed.root");
  TFile *f10 = new TFile("/home/ethan/Documents/WJet_QQ_combined_processed.root");
  TFile *f11 = new TFile("/home/ethan/Documents/TTToHadronic_combined_processed.root");
  TFile *f12 = new TFile("/home/ethan/Documents/TTTo2l2nu_combine_processed.root");
  TFile *f13 = new TFile("/home/ethan/Documents/TTtoSemiLeptonic_combined_processed.root");


  TH2F *h_HT1000to1500_HT = (TH2F*)f1->Get("h_MSJ_mass_vs_MdSJ");
  TH2F *h_HT1500to2000_HT = (TH2F*)f2->Get("h_MSJ_mass_vs_MdSJ");
  TH2F *h_HT2000toInf_HT = (TH2F*)f3->Get("h_MSJ_mass_vs_MdSJ");

  h_HT1000to1500_HT->Scale(QCD1000to1500_SF);
  h_HT1500to2000_HT->Scale(QCD1500to2000_SF);
  h_HT2000toInf_HT->Scale(QCD2000toInf_SF);


  TH2F *h_WJetsToLNu = (TH2F*)f9->Get("h_MSJ_mass_vs_MdSJ");
  TH2F *h_WJet_QQ = (TH2F*)f10->Get("h_MSJ_mass_vs_MdSJ");
  TH2F *h_TTToHadronic = (TH2F*)f11->Get("h_MSJ_mass_vs_MdSJ");
  TH2F *h_TTTo2l2nu = (TH2F*)f12->Get("h_MSJ_mass_vs_MdSJ");
  TH2F *h_TTToSemiLeptonic = (TH2F*)f13->Get("h_MSJ_mass_vs_MdSJ");

  h_WJetsToLNu->Scale(h_WJetsToLNu_SF);
  h_WJet_QQ->Scale(h_WJet_QQ_SF);
  h_TTToHadronic->Scale(h_TTToHadronic_SF);
  h_TTTo2l2nu->Scale(h_TTTo2l2nu_SF);
  h_TTToSemiLeptonic->Scale(h_TTToSemiLeptonic_SF);

  TH2F *h_TTbar_combined = new TH2F(*h_TTToHadronic);
  h_TTbar_combined->Add(h_TTTo2l2nu);
  h_TTbar_combined->Add(h_TTToSemiLeptonic);

  TH2F *h_WJets_combined = new TH2F(*h_WJetsToLNu);
  h_WJets_combined->Add(h_WJet_QQ);


  TH2F *h_JetHTA_HT = (TH2F*)f5->Get("h_MSJ_mass_vs_MdSJ");
  TH2F *h_JetHTB_HT = (TH2F*)f6->Get("h_MSJ_mass_vs_MdSJ");
  TH2F *h_JetHTC_HT = (TH2F*)f7->Get("h_MSJ_mass_vs_MdSJ");
  TH2F *h_JetHTD_HT = (TH2F*)f8->Get("h_MSJ_mass_vs_MdSJ");
  int totalNDataEvents = h_JetHTA_HT->GetEntries() + h_JetHTB_HT->GetEntries() + h_JetHTC_HT->GetEntries() + h_JetHTD_HT->GetEntries();
  std::cout << "Total number of events is " << totalNDataEvents << std::endl;


  TCanvas *c1 = new TCanvas("c1","",400,20, 2000,2000);
  gPad->SetLogz();

  gStyle->SetOptStat(0);

  TH2F *h_JetHT_combined = new TH2F(*h_JetHTA_HT);
  h_JetHT_combined->GetYaxis()->SetTitleOffset(1.35);
  h_JetHT_combined->GetYaxis()->SetLabelSize(0.015);
  h_JetHT_combined->Add(h_JetHTB_HT);
  h_JetHT_combined->Add(h_JetHTC_HT);
  h_JetHT_combined->Add(h_JetHTD_HT);
  h_JetHT_combined->SetTitle("M_{SJ} vs M_{diSJ} for Data (JetHT 2018);M_{diSJ};M_{SJ}");

  h_JetHT_combined->Draw("colz");
  c1->SaveAs("mSJ_vs_mdiSJ_DATA.png");


  TH2F *h_QCD_combined = new TH2F(*h_HT1000to1500_HT);
  h_QCD_combined->Add(h_HT1500to2000_HT);
  h_QCD_combined->Add(h_HT2000toInf_HT);


  h_QCD_combined->SetFillColor(kRed);
  h_TTbar_combined->SetFillColor(kBlue);
  h_WJets_combined->SetFillColor(kGreen);


  TH2F *h_BR_combined = new TH2F(*h_QCD_combined);
  h_BR_combined->Add(h_TTbar_combined);
  h_BR_combined->Add(h_WJets_combined);
  h_BR_combined->SetTitle("M_{SJ} vs M_{diSJ} for Combined MC BR;M_{diSJ};M_{SJ}");
  h_BR_combined->GetYaxis()->SetTitleOffset(1.35);
  h_BR_combined->GetYaxis()->SetLabelSize(0.015);

  h_BR_combined->Draw("colz");
  c1->SaveAs("mSJ_vs_mdiSJ_BR.png");


  h_JetHT_combined->Divide(h_BR_combined);
  h_JetHT_combined->SetTitle("data/MC ratio of M_{SJ} vs M_{diSJ};M_{diSJ};M_{SJ}");
  h_JetHT_combined->GetYaxis()->SetTitleOffset(1.35);
  h_JetHT_combined->GetYaxis()->SetLabelSize(0.015);

  h_JetHT_combined->Draw("colz");
  c1->SaveAs("mSJ_vs_mdiSJ_ratio.png");


}


//There are 3338 total events, with 97 anti-tagged, and of these events there are 82 with second superjet tagged.