VBS/plotInterference.C at master · AndrewLevin/VBS · 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
//this should go in the directory of https://github.com/AndrewLevin/WWScattering
//change the name of the .rootlogon.C file if this macro does not work at first

#include "TROOT.h"
#include "TSystem.h"
#include "TInterpreter.h"
#include "TFile.h"
#include "TCanvas.h"
#include "TH1F.h"
#include "TStyle.h"
#include "TPad.h"
#include "Math/QuantFuncMathCore.h"
#include "TMath.h"
#include "TGraphAsymmErrors.h"
#include "StandardPlotVBS.C"

 //.x finalPlotVBS.C+(0,1,"m_{jj} (GeV)","GeV","histo_nice8TeV.root","wwss_mjj",0,19.4)


void plotInterference(){

  gStyle->SetFillStyle(0);
  gStyle->SetLegendBorderSize(0);
  gROOT->ForceStyle();


  TFile* file_ewk = new TFile("/data/blue/anlevin/VBS/21Jun2014/wwss_mjj_qed_4_qcd_0_sm.root", "read");
  TFile* file_qcd = new TFile("/data/blue/anlevin/VBS/21Jun2014/wwss_mjj_qed_2_qcd_99_sm.root", "read");
  TFile* file_ewk_plus_qcd = new TFile("/data/blue/anlevin/VBS/21Jun2014/wwss_mjj_qed_4_qcd_99_sm.root", "read");

  TH1F* hSignalEwk = (TH1F*)file_ewk->Get("histo0");
  TH1F* hSignalQCD = (TH1F*)file_qcd->Get("histo0");
  TH1F* hSignalEwkPlusQCD = (TH1F*)file_ewk_plus_qcd->Get("histo0");

  //TCanvas* c2 = new TCanvas("c2", "c2",700,50,500,500);
  //c2->cd(1);

  TH1F * hSum= (TH1F*)hSignalEwk->Clone();
  hSum->Add(hSignalQCD);

  hSum->Draw();


  TLegend* legend = new TLegend(0.5,0.74,0.8,0.9);
  legend->AddEntry(hSum,"sum of EWK and QCD");
  legend->AddEntry(hSignalEwkPlusQCD,"EWK and QCD with interference");


  legend->SetBorderSize(     0);
  legend->SetFillColor (     0);
  legend->SetTextAlign (    12);
  legend->SetTextFont  (    42);
  legend->SetTextSize(0.03);


  hSum->SetLineWidth(3);
  hSignalEwkPlusQCD->SetLineWidth(3);
  hSum->SetStats(0);


  hSum->SetTitle("");
  hSum->SetLineColor(kRed);


  hSum->Draw("e");
  hSum->GetYaxis()->SetRangeUser(0,4);
  hSum->GetXaxis()->SetTitle("m_{jj} (GeV)}");


  hSignalEwkPlusQCD->Draw("eSAME");
  legend->Draw("SAME");

}