ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/dasu/UltraFastSim/Plot.C
(Generate patch)

Comparing UserCode/dasu/UltraFastSim/Plot.C (file contents):
Revision 1.2 by aarond, Thu Feb 24 16:36:34 2011 UTC vs.
Revision 1.3 by dasu, Fri Mar 11 19:44:55 2011 UTC

# Line 1 | Line 1
1 + #include <string>
2 +
3 + void Plot()
4   {
5 < // Make plots
6 <  TFile zh("ZHmmbbOut.root");
7 <  TFile zj("ZmmOut.root");
8 <  TFile zz("ZZOut.root");
5 >
6 >  char fName[1024];
7 >
8 >  std::string zhName("ZH-plots-v5.root");
9 >  std::string zzName("ZZ-plots-v5.root");
10 >  std::string zjName("Zmm-v5-Plots.root");
11 >
12 >  // Make plots
13 >
14 >  TFile zh(zhName.c_str());
15 >  TFile zj(zjName.c_str());
16 >  TFile zz(zzName.c_str());
17  
18    TH1F *zhmuonpt=zh.Get("Muons Pt");
19    TH1F *zjmuonpt=zj.Get("Muons Pt");
# Line 18 | Line 29
29  
30    TCanvas c1;
31  
32 <  zjmuonpt->Draw();
32 >  zhmuonpt->Draw();
33    zzmuonpt->Draw("same");
34 <  zhmuonpt->Draw("same");
34 >  zjmuonpt->Draw("same");
35  
36    TLegend *lg=new TLegend(0.5,0.5,0.7,0.7);
37    lg->SetFillColor(kWhite);
# Line 29 | Line 40
40    lg->AddEntry(zzmuonpt,"ZZ","l");
41    lg->Draw();
42  
43 <  /////////////////////////////
33 <
34 <  TH1F *zhbjetpt=zh.Get("B-jets Et");
35 <  TH1F *zjbjetpt=zj.Get("B-jets Et");
36 <  TH1F *zzbjetpt=zz.Get("B-jets Et");
37 <
38 <  zjbjetpt->SetLineColor(kGreen);
39 <  zzbjetpt->SetLineColor(kOrange+3);
40 <  zhbjetpt->SetLineColor(kRed);
41 <
42 <  zjbjetpt->Scale(1./zjbjetpt->Integral());
43 <  zzbjetpt->Scale(1./zzbjetpt->Integral());
44 <  zhbjetpt->Scale(1./zhbjetpt->Integral());
45 <
46 <  TCanvas c2;
47 <
48 <  zjbjetpt->Draw();
49 <  zzbjetpt->Draw("same");
50 <  zhbjetpt->Draw("same");
51 <
52 <  TLegend *lg2=new TLegend(0.5,0.5,0.7,0.7);
53 <  lg2->SetFillColor(kWhite);
54 <  lg2->AddEntry(zhbjetpt,"ZH, m(H)=120 GeV","l");
55 <  lg2->AddEntry(zjbjetpt,"Z+jets","l");
56 <  lg2->AddEntry(zzbjetpt,"ZZ","l");
57 <  lg2->Draw();
58 <
59 <  /////////////////////////////
60 <
61 <  TH1F *zhzpt=zh.Get("Z Pt");
62 <  TH1F *zjzpt=zj.Get("Z Pt");
63 <  TH1F *zzzpt=zz.Get("Z Pt");
64 <
65 <  zjzpt->SetLineColor(kGreen);
66 <  zzzpt->SetLineColor(kOrange+3);
67 <  zhzpt->SetLineColor(kRed);
68 <
69 <  zjzpt->Scale(1./zjzpt->Integral());
70 <  zzzpt->Scale(1./zzzpt->Integral());
71 <  zhzpt->Scale(1./zhzpt->Integral());
72 <
73 <  TCanvas c3;
74 <
75 <  zjzpt->Draw();
76 <  zzzpt->Draw("same");
77 <  zhzpt->Draw("same");
78 <
79 <  TLegend *lg3=new TLegend(0.5,0.5,0.7,0.7);
80 <  lg3->SetFillColor(kWhite);
81 <  lg3->AddEntry(zhzpt,"ZH, m(H)=120 GeV","l");
82 <  lg3->AddEntry(zjzpt,"Z+jets","l");
83 <  lg3->AddEntry(zzzpt,"ZZ","l");
84 <  lg3->Draw();
85 <
86 <
87 <  /////////////////////////////
88 <
89 <  TH1F *zhhpt=zh.Get("H Pt");
90 <  TH1F *zjhpt=zj.Get("H Pt");
91 <  TH1F *zzhpt=zz.Get("H Pt");
92 <
93 <  zjhpt->SetLineColor(kGreen);
94 <  zzhpt->SetLineColor(kOrange+3);
95 <  zhhpt->SetLineColor(kRed);
96 <
97 <  zjhpt->Rebin(5);
98 <  zzhpt->Rebin(5);
99 <  zhhpt->Rebin(5);
100 <
101 <  zjhpt->Scale(1./zjhpt->Integral());
102 <  zzhpt->Scale(1./zzhpt->Integral());
103 <  zhhpt->Scale(1./zhhpt->Integral());
104 <
105 <  TCanvas c4;
106 <
107 <  zjhpt->Draw();
108 <  zzhpt->Draw("same");
109 <  zhhpt->Draw("same");
110 <
111 <  TLegend *lg4=new TLegend(0.5,0.5,0.7,0.7);
112 <  lg4->SetFillColor(kWhite);
113 <  lg4->AddEntry(zhhpt,"ZH, m(H)=120 GeV","l");
114 <  lg4->AddEntry(zjhpt,"Z+jets","l");
115 <  lg4->AddEntry(zzhpt,"ZZ","l");
116 <  lg4->Draw();
117 <
43 >  c1.SaveAs("muonPt.png");;
44  
45    /////////////////////////////
46  
# Line 132 | Line 58
58  
59    TCanvas c5;
60  
61 <  zjzinvmass->Draw();
61 >  zhzinvmass->Draw();
62    zzzinvmass->Draw("same");
63 <  zhzinvmass->Draw("same");
63 >  zjzinvmass->Draw("same");
64  
65    TLegend *lg5=new TLegend(0.5,0.5,0.7,0.7);
66    lg5->SetFillColor(kWhite);
# Line 143 | Line 69
69    lg5->AddEntry(zzzinvmass,"ZZ","l");
70    lg5->Draw();
71  
72 +  c5.SaveAs("ZMass.png");;
73  
74    /////////////////////////////
75  
76 <  TH1F *zhhinvmass=zh.Get("H invmass preselection");
77 <  TH1F *zjhinvmass=zj.Get("H invmass preselection");
78 <  TH1F *zzhinvmass=zz.Get("H invmass preselection");
79 <
80 <  zjhinvmass->SetLineColor(kGreen);
81 <  zzhinvmass->SetLineColor(kOrange+3);
82 <  zhhinvmass->SetLineColor(kRed);
83 <
84 <  zjhinvmass->Rebin(5);
85 <  zzhinvmass->Rebin(5);
86 <  zhhinvmass->Rebin(5);
87 <
88 <  zjhinvmass->Scale(1./zjhinvmass->Integral());
89 <  zzhinvmass->Scale(1./zzhinvmass->Integral());
90 <  zhhinvmass->Scale(1./zhhinvmass->Integral());
91 <
92 <  TCanvas c6;
93 <
94 <  zjhinvmass->Draw();
95 <  zzhinvmass->Draw("same");
96 <  zhhinvmass->Draw("same");
97 <
98 <  TLegend *lg6=new TLegend(0.5,0.5,0.7,0.7);
99 <  lg6->SetFillColor(kWhite);
100 <  lg6->AddEntry(zhhinvmass,"ZH, m(H)=120 GeV","l");
101 <  lg6->AddEntry(zjhinvmass,"Z+jets","l");
102 <  lg6->AddEntry(zzhinvmass,"ZZ","l");
103 <  lg6->Draw();
104 <
105 <  /////////////////////////////
106 <
107 <  TH1F *zhhinvmass2=zh.Get("H invmass");
108 <  TH1F *zjhinvmass2=zj.Get("H invmass");
109 <  TH1F *zzhinvmass2=zz.Get("H invmass");
110 <
111 <  zjhinvmass2->SetLineColor(kGreen);
112 <  zzhinvmass2->SetLineColor(kOrange+3);
113 <  zhhinvmass2->SetLineColor(kRed);
114 <
115 <  zjhinvmass2->Rebin(5);
116 <  zzhinvmass2->Rebin(5);
117 <  zhhinvmass2->Rebin(5);
118 <
119 <  //zjhinvmass2->Scale(1./zjhinvmass2->Integral());
120 <  zzhinvmass2->Scale(1./zzhinvmass2->Integral());
121 <  zhhinvmass2->Scale(1./zhhinvmass2->Integral());
122 <
123 <  TCanvas c7;
124 <
125 <  zjhinvmass2->Draw();
126 <  zzhinvmass2->Draw("same");
127 <  zhhinvmass2->Draw("same");
128 <
129 <  TLegend *lg7=new TLegend(0.5,0.5,0.7,0.7);
130 <  lg7->SetFillColor(kWhite);
131 <  lg7->AddEntry(zhhinvmass2,"ZH, m(H)=120 GeV","l");
132 <  lg7->AddEntry(zjhinvmass2,"Z+jets","l");
133 <  lg7->AddEntry(zzhinvmass2,"ZZ","l");
134 <  lg7->Draw();
135 <
136 <  /////////////////////////////
137 <
138 <  TH1F *zhdphi=zh.Get("Dphi(Z,H)");
139 <  TH1F *zjdphi=zj.Get("Dphi(Z,H)");
140 <  TH1F *zzdphi=zz.Get("Dphi(Z,H)");
141 <
142 <  zjdphi->SetLineColor(kGreen);
143 <  zzdphi->SetLineColor(kOrange+3);
144 <  zhdphi->SetLineColor(kRed);
145 <
146 <  zjdphi->Scale(1./zjdphi->Integral());
147 <  zzdphi->Scale(1./zzdphi->Integral());
148 <  zhdphi->Scale(1./zhdphi->Integral());
149 <
150 <  TCanvas c8;
151 <
152 <  zjdphi->Draw();
153 <  zzdphi->Draw("same");
227 <  zhdphi->Draw("same");
228 <
229 <  TLegend *lg8=new TLegend(0.5,0.5,0.7,0.7);
230 <  lg8->SetFillColor(kWhite);
231 <  lg8->AddEntry(zhdphi,"ZH, m(H)=120 GeV","l");
232 <  lg8->AddEntry(zjdphi,"Z+jets","l");
233 <  lg8->AddEntry(zzdphi,"ZZ","l");
234 <  lg8->Draw();
235 <
236 <
237 <
238 <
76 >  vector<string> geom;
77 >  geom.push_back("Phase 1");
78 >  geom.push_back("StdGeom");
79 >
80 >  vector<string> bJetOP;
81 >  bJetOP.push_back("Loose");
82 >  bJetOP.push_back("Medium");
83 >  bJetOP.push_back("Tight");
84 >
85 >  vector<string> hNames;
86 >  hNames.push_back("B-jets Et ");
87 >  hNames.push_back("B-jets Eta ");
88 >  hNames.push_back("B-jets Mult. ");
89 >  hNames.push_back("B-jets Phi ");
90 >  hNames.push_back("Z Pt ");
91 >  hNames.push_back("H Phi ");
92 >  hNames.push_back("H Pt ");
93 >  hNames.push_back("H invmass ");
94 >  hNames.push_back("H invmass preselection ");
95 >  hNames.push_back("Dphi(Z,H) ");
96 >  
97 >  for(int i = 0; i < geom.size(); i++) {
98 >    for(int j = 0; j < bJetOP.size(); j++) {
99 >      for(int k = 0; k < hNames.size(); k++) {
100 >        
101 >        cout << hNames[k] << endl;
102 >
103 >        string hName(hNames[k].c_str());
104 >        hName += bJetOP[j];
105 >        hName += " ";
106 >        hName += geom[i];
107 >
108 >        cout << hName << endl;
109 >        
110 >        TH1F *zhHist=zh.Get(hName.c_str());
111 >        TH1F *zjHist=zj.Get(hName.c_str());
112 >        TH1F *zzHist=zz.Get(hName.c_str());
113 >
114 >        if(zhHist == 0 || zjHist == 0 || zzHist == 0) {
115 >          cout << "Failed to find " << hName << endl;
116 >          break;
117 >        }
118 >        
119 >        zjHist->SetLineColor(kGreen);
120 >        zzHist->SetLineColor(kOrange+3);
121 >        zhHist->SetLineColor(kRed);
122 >        
123 >        zjHist->Scale(1./zjHist->Integral());
124 >        zzHist->Scale(1./zzHist->Integral());
125 >        zhHist->Scale(1./zhHist->Integral());
126 >        
127 >        TCanvas *c = new TCanvas();
128 >        
129 >        zhHist->Draw();
130 >        zzHist->Draw("same");
131 >        zjHist->Draw("same");
132 >        
133 >        TLegend *lg2=new TLegend(0.5,0.5,0.7,0.7);
134 >        lg2->SetFillColor(kWhite);
135 >        lg2->AddEntry(zhHist,"ZH, m(H)=120 GeV","l");
136 >        lg2->AddEntry(zjHist,"Z+jets","l");
137 >        lg2->AddEntry(zzHist,"ZZ","l");
138 >        lg2->Draw();
139 >        
140 >        strcpy(fName, hName.c_str());
141 >        for(int l = 0; l < strlen(fName); l++) {
142 >          if(fName[l] == ' ') fName[l] = '-';
143 >        }
144 >        strcat(fName, ".png");
145 >        c->SaveAs(fName);
146 >
147 >        //      delete lg2;
148 >        //      delete c;
149 >        
150 >      }
151 >    }
152 >  }
153 >  
154   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines