]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG3/muon/PlotMuonQAPbPb.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Adding PWG4Base library
[u/mrichter/AliRoot.git] / PWG3 / muon / PlotMuonQAPbPb.C
CommitLineData
4916d740 1//--------------------------------------------------------------------------\r
2// Macro for QA monitoring.\r
3//\r
4// In case it is not run with full aliroot, it needs the following libraries to compile:\r
5// - libSTEERBase.so\r
6// - libESD.so\r
7// - libAOD.so\r
8// - libANALYSIS.so\r
9// - libANALYSISalice.so\r
10// - libCORRFW.so\r
11// - libPWG3muon.so\r
12//\r
13// The macro reads results of the QA task and produce monitoring plots.\r
14//\r
15// Author: Philippe Pillot - SUBATECH Nantes, Christophe Suire - IPN Orsay \r
16// and Cynthia Hadjidakis - IPN Orsay\r
17//--------------------------------------------------------------------------\r
18\r
19#if !defined(__CINT__) || defined(__MAKECINT__)\r
20\r
21#include <Riostream.h>\r
22\r
23// ROOT includes\r
24#include "TMath.h"\r
25#include "TGrid.h"\r
26#include "TFile.h"\r
27#include "TH1.h"\r
28#include "TSystem.h"\r
29#include "TStyle.h"\r
30#include "TCanvas.h"\r
31#include "TPad.h"\r
32\r
33// ALIROOT includes\r
34#include "../PWG3/base/AliCounterCollection.h"\r
35\r
36#endif\r
37\r
38//--------------------------------------------------------------------------\r
39void PlotMuonQAPbPb(const char* baseDir, const char* runList = 0x0, Bool_t selectPhysics = kFALSE)\r
40{\r
41 /// Macro for QA monitoring.\r
42 /// Examples: \r
43 /// baseDir = "alien:///alice/cern.ch/user/p/ppillot/pp7TeV/LHC10d/MuonQA/pass1/results"\r
44 /// baseDir = "results"\r
45 /// If runList != 0x0: only the given runs will be used. Otherwise use all runs found in baseDir.\r
46 /// Usage (local/grid)\r
47 /// .x PlotMuonQAPbPb.C("alien:///alice/cern.ch/user/s/suire/QAP2LHC10h/output",0,kTRUE)\r
48 \r
49\r
50 #if defined(__CINT__) && !defined(__MAKECINT__)\r
51 gSystem->Load("libTree");\r
52 gSystem->Load("libGeom");\r
53 gSystem->Load("libVMC");\r
54 gSystem->Load("libPhysics");\r
55 gSystem->Load("libSTEERBase");\r
56 gSystem->Load("libESD");\r
57 gSystem->Load("libAOD");\r
58 gSystem->Load("libANALYSIS");\r
59 gSystem->Load("libANALYSISalice");\r
60 gSystem->Load("libCORRFW");\r
61 gSystem->Load("libPWG3base");\r
62 gSystem->Load("libPWG3muon");\r
63 #endif\r
64 \r
65 // Cosmetics and configuration\r
66 gStyle->SetFillColor(10);\r
67 gStyle->SetPadGridX(kTRUE);\r
68 gStyle->SetPadGridY(kTRUE);\r
69 gStyle->SetPadRightMargin(0.01);\r
70\r
71 TString LHCPeriod = "LHC10hpass2"; \r
72\r
73 TString OutFileName = "muonQA_"; OutFileName += LHCPeriod;\r
74 TString OutFileNamePDF= OutFileName.Data(); OutFileNamePDF+= ".pdf";\r
75 TString OutFileNamePDF_open = OutFileNamePDF.Data(); OutFileNamePDF_open += "["; \r
76 TString OutFileNamePDF_close= OutFileNamePDF.Data(); OutFileNamePDF_close += "]"; \r
77 TString OutFileNameROOT= OutFileName.Data(); OutFileNameROOT+= ".root";\r
78\r
79 //TString QAFileName = "QAresults.root";\r
80 TString QAFileName = "AnalysisResults.root";\r
81 TString MergedQAFileName = "Merged"; MergedQAFileName+= QAFileName.Data();\r
82\r
83 Int_t PRINTSTAT = 1; \r
84\r
85 Int_t ALIENFASTREAD = 0 ; \r
86\r
87 if (ALIENFASTREAD){ // Equivalent to the fast read option\r
88 gEnv->SetValue("XNet.ConnectTimeout",10);\r
89 gEnv->SetValue("XNet.RequestTimeout",10);\r
90 gEnv->SetValue("XNet.MaxRedirectCount",2);\r
91 gEnv->SetValue("XNet.ReconnectTimeout",10);\r
92 gEnv->SetValue("XNet.FirstConnectMaxCnt",1);\r
93 }\r
94\r
95 TH1::AddDirectory(kFALSE);\r
96 \r
97 TString alienBaseDir = baseDir;\r
98 \r
99 if (alienBaseDir.Contains("alien:") && !TGrid::Connect("alien://")) {\r
100 Error("PlotMuonQA","cannot connect to grid");\r
101 return;\r
102 }\r
103 \r
104 Float_t LabelSize = 0.03; \r
105 \r
106 \r
107 //---------------------------------- //\r
108 // Run selection //\r
109 //---------------------------------- //\r
110 \r
111\r
112 // list runs to be analyzed\r
113 TString selectRuns = "run:";\r
114 TObjArray runs;\r
115 runs.SetOwner();\r
116 if (runList) {\r
117 \r
118 // only the ones in the runList\r
119 ifstream inFile(runList);\r
120 if (!inFile.is_open()) {\r
121 Error("PlotMuonQAPbPb",Form("unable to open file %s", runList));\r
122 return;\r
123 }\r
124 \r
125 TString currRun;\r
126 while (!inFile.eof()) {\r
127 currRun.ReadLine(inFile, kTRUE);\r
128 if (currRun.IsNull()) continue;\r
129 if (!currRun.IsDigit()) {\r
130 Error("PlotMuonQAPbPb","invalid run number: %s", currRun.Data());\r
131 return;\r
132 }\r
133 runs.AddLast(new TObjString(Form("%09d", currRun.Atoi())));\r
134 selectRuns += Form("%s,",currRun.Data());\r
135 }\r
136 selectRuns.Remove(TString::kTrailing, ',');\r
137 \r
138 inFile.close();\r
139 \r
140 } else {\r
141 \r
142 // all runs\r
143 runs.AddLast(new TObjString("*"));\r
144 \r
145 }\r
146 \r
147 // physics selection\r
148 TString select = selectPhysics ? "selected:yes" : "";\r
149 \r
150 \r
151 //---------------------------------- //\r
152 // plot global counter //\r
153 //---------------------------------- //\r
154 \r
155 TFile *globalFile = TFile::Open(Form("%s/%s", baseDir,MergedQAFileName.Data()));\r
156 if (!globalFile || ! globalFile->IsOpen()) {\r
157 Error("PlotQA", Form("failed to open file: %s/%s", baseDir,MergedQAFileName.Data()));\r
158 return;\r
159 }\r
160 globalFile->Cd("MUON_QA");\r
161 \r
162\r
163\r
164 // get counters\r
165 AliCounterCollection* eventCounters = static_cast<AliCounterCollection*>(globalFile->FindObjectAny("eventCounters"));\r
166 AliCounterCollection* trackCounters = static_cast<AliCounterCollection*>(globalFile->FindObjectAny("trackCounters"));\r
167 if (!runList) selectRuns += trackCounters->GetKeyWords("run");\r
168\r
169\r
170 //declare a default canvas c1 \r
171 TString CanvasName = "c1";\r
172 TCanvas *c1 = new TCanvas(CanvasName.Data());\r
173 c1->cd();\r
174\r
175 //TString allTriggers = "trigger:CMBAC-B,CMBS2A-B,CMBS2C-B,CMBACS2-B,C0SMH-B";\r
176 //TString minBiasTrigers = "trigger:CMBAC-B,CMBS2A-B,CMBS2C-B,CMBACS2-B";\r
177 TString allTriggers = "trigger:any"; //"trigger:CMBAC-B,CMBS2A-B,CMBS2C-B,C0SMH-B";\r
178 TString minBiasTrigers = "trigger:any"; //CMBAC-B,CMBS2A-B,CMBS2C-B";\r
179 TString bckTriggers = "trigger:CMBACS2-A,CMBACS2-C,CMBS2A-A,CMBS2A-C,CMBS2C-A,CMBS2C-C,CMBAC-A,CMBAC-A,C0SMH-A,C0SMH-C";\r
180\r
181 // Histo trigger without Phys. Sel. \r
182 TH1* hAllTriggersNoPS = eventCounters->Draw("run",Form("%s/%s",allTriggers.Data(), selectRuns.Data()));\r
183 hAllTriggersNoPS->Sumw2();\r
184 Int_t NumOfAllTriggersNoPS = hAllTriggersNoPS->Integral();\r
185 TH1* hAllTriggerslowNoPS = eventCounters->Draw("run",Form("v0mult:low/%s/%s",allTriggers.Data(), selectRuns.Data()));\r
186 hAllTriggerslowNoPS->Sumw2();\r
187 TH1* hAllTriggershighNoPS = eventCounters->Draw("run",Form("v0mult:high/%s/%s",allTriggers.Data(), selectRuns.Data()));\r
188 hAllTriggershighNoPS->Sumw2();\r
189 TH1* hMBBNoPS = eventCounters->Draw("run",Form("v0mult:low,int,high/%s/%s",minBiasTrigers.Data(), selectRuns.Data()));\r
190 hMBBNoPS->Sumw2();\r
191 Int_t NumOfMBBNoPS = hMBBNoPS->Integral();\r
192 TH1* hMBBlowNoPS = eventCounters->Draw("run",Form("v0mult:low/%s/%s",minBiasTrigers.Data(), selectRuns.Data()));\r
193 hMBBlowNoPS->Sumw2();\r
194 Int_t NumOfMBBlowNoPS = hMBBlowNoPS->Integral();\r
195 TH1* hMBBhighNoPS = eventCounters->Draw("run",Form("v0mult:high/%s/%s",minBiasTrigers.Data(), selectRuns.Data()));\r
196 hMBBhighNoPS->Sumw2();\r
197 Int_t NumOfMBBhighNoPS = hMBBhighNoPS->Integral();\r
198 //TH1* hCMUS1BNoPS = eventCounters->Draw("run",Form("trigger:CMUS1B/%s", selectRuns.Data()));\r
199 //hCMUS1BNoPS->Sumw2();\r
200 //Int_t NumOfCMUS1BNoPS = hCMUS1BNoPS->Integral();\r
201 //TH1* hCSH1BNoPS = eventCounters->Draw("run",Form("trigger:CSH1B/%s", selectRuns.Data()));\r
202 //hCSH1BNoPS->Sumw2();\r
203 //Int_t NumOfCSH1BNoPS = hCSH1BNoPS->Integral();\r
204\r
205 TH1* hBckTriggersNoPS = eventCounters->Draw("run",Form("v0mult:low,int,high/%s/%s",bckTriggers.Data(), selectRuns.Data()));\r
206 hBckTriggersNoPS->Sumw2(); \r
207 \r
208\r
209 // Histo trigger with Phys. Sel. \r
210 TH1* hAllTriggersWithPS = eventCounters->Draw("run",Form("%s/%s/selected:yes",allTriggers.Data(), selectRuns.Data()));\r
211 hAllTriggersWithPS->Sumw2();\r
212 Int_t NumOfAllTriggersWithPS = hAllTriggersWithPS->Integral();\r
213 TH1* hAllTriggersLowWithPS = eventCounters->Draw("run",Form("v0mult:low/%s/%s/selected:yes",allTriggers.Data(), selectRuns.Data()));\r
214 hAllTriggersWithPS->Sumw2();\r
215 TH1* hAllTriggersHighWithPS = eventCounters->Draw("run",Form("v0mult:high/%s/%s/selected:yes",allTriggers.Data(), selectRuns.Data()));\r
216 hAllTriggersWithPS->Sumw2();\r
217 TH1* hMBBWithPS = eventCounters->Draw("run",Form("v0mult:low,int,high/%s/%s/selected:yes",minBiasTrigers.Data(), selectRuns.Data()));\r
218 hMBBWithPS->Sumw2();\r
219 Int_t NumOfMBBWithPS = hMBBWithPS->Integral();\r
220 TH1* hMBBlowWithPS = eventCounters->Draw("run",Form("v0mult:low/%s/%s/selected:yes",minBiasTrigers.Data(), selectRuns.Data()));\r
221 hMBBlowWithPS->Sumw2();\r
222 Int_t NumOfMBBlowWithPS = hMBBlowWithPS->Integral();\r
223 TH1* hMBBhighWithPS = eventCounters->Draw("run",Form("v0mult:high/%s/%s/selected:yes",minBiasTrigers.Data(), selectRuns.Data()));\r
224 hMBBhighWithPS->Sumw2();\r
225 Int_t NumOfMBBhighWithPS = hMBBhighWithPS->Integral();\r
226 \r
227 // Histo trigger : Phys. Sel. is selected or not depending on the macro arguments\r
228 TH1* hAllTriggers = eventCounters->Draw("run",Form("%s/%s/%s",allTriggers.Data(), selectRuns.Data(), select.Data()));\r
229 hAllTriggers->Sumw2();\r
230 TH1* hMBB = eventCounters->Draw("run",Form("v0mult:low,int,high/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
231 hMBB->Sumw2();\r
232 TH1* hMBBlow = eventCounters->Draw("run",Form("v0mult:low/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
233 hMBBlow->Sumw2();\r
234 TH1* hMBBhigh = eventCounters->Draw("run",Form("v0mult:high/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
235 hMBBhigh->Sumw2();\r
236\r
237 // Histo tracking : Phys. Sel. is selected or not depending on the macro arguments\r
238 TH1* hAllTracksMBB = trackCounters->Draw("run",Form("v0mult:low,int,high/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
239 hAllTracksMBB->Sumw2();\r
240 TH1* hTriggerMBB = trackCounters->Draw("run",Form("v0mult:low,int,high/track:triggeronly/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
241 hTriggerMBB->Sumw2();\r
242 TH1* hTrackerMBB = trackCounters->Draw("run",Form("v0mult:low,int,high/track:trackeronly/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
243 hTrackerMBB->Sumw2();\r
244 TH1* hMatchedMBB = trackCounters->Draw("run",Form("v0mult:low,int,high/track:matched/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
245 hMatchedMBB->Sumw2();\r
246 TH1* hAllTracksMBBlow = trackCounters->Draw("run",Form("v0mult:low/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
247 hAllTracksMBBlow->Sumw2();\r
248 TH1* hTriggerMBBlow = trackCounters->Draw("run",Form("v0mult:low/track:triggeronly/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
249 hTriggerMBBlow->Sumw2();\r
250 TH1* hTrackerMBBlow = trackCounters->Draw("run",Form("v0mult:low/track:trackeronly/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
251 hTrackerMBBlow->Sumw2();\r
252 TH1* hMatchedMBBlow = trackCounters->Draw("run",Form("v0mult:low/track:matched/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
253 hMatchedMBBlow->Sumw2();\r
254 TH1* hAllTracksMBBhigh = trackCounters->Draw("run",Form("v0mult:high/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
255 hAllTracksMBBhigh->Sumw2();\r
256 TH1* hTriggerMBBhigh = trackCounters->Draw("run",Form("v0mult:high/track:triggeronly/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
257 hTriggerMBBhigh->Sumw2();\r
258 TH1* hTrackerMBBhigh = trackCounters->Draw("run",Form("v0mult:high/track:trackeronly/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
259 hTrackerMBBhigh->Sumw2();\r
260 TH1* hMatchedMBBhigh = trackCounters->Draw("run",Form("v0mult:high/track:matched/%s/%s/%s",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
261 hMatchedMBBhigh->Sumw2();\r
262 \r
263 TH1* hMatchedLowPtMBB = trackCounters->Draw("run",Form("v0mult:low,int,high/track:matched/acc:in/%s/%s/%s/pt:low",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
264 hMatchedLowPtMBB->Sumw2();\r
265 TH1* hMatchedHighPtMBB = trackCounters->Draw("run",Form("v0mult:low,int,high/track:matched/acc:in/%s/%s/%s/pt:high",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
266 hMatchedHighPtMBB->Sumw2();\r
267 TH1* hMatchedLowPtMBBlow = trackCounters->Draw("run",Form("v0mult:low/track:matched/acc:in/%s/%s/%s/pt:low",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
268 hMatchedLowPtMBBlow->Sumw2();\r
269 TH1* hMatchedHighPtMBBlow = trackCounters->Draw("run",Form("v0mult:low/track:matched/acc:in/%s/%s/%s/pt:high",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
270 hMatchedHighPtMBBlow->Sumw2();\r
271 TH1* hMatchedLowPtMBBhigh = trackCounters->Draw("run",Form("v0mult:high/track:matched/acc:in/%s/%s/%s/pt:low",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
272 hMatchedLowPtMBBhigh->Sumw2();\r
273 TH1* hMatchedHighPtMBBhigh = trackCounters->Draw("run",Form("v0mult:high/track:matched/acc:in/%s/%s/%s/pt:high",minBiasTrigers.Data(), selectRuns.Data(), select.Data()));\r
274 hMatchedHighPtMBBhigh->Sumw2();\r
275 \r
276 \r
277 // for MBB only, low and high mult\r
278 TH1 *hPosMatchedMBBlow = trackCounters->Draw("run",Form("v0mult:low/track:matched/charge:pos/acc:in/%s/%s/%s",minBiasTrigers.Data(),select.Data(),selectRuns.Data()));\r
279 hPosMatchedMBBlow->Sumw2();\r
280 TH1 *hNegMatchedMBBlow= trackCounters->Draw("run",Form("v0mult:low/track:matched/charge:neg/acc:in/%s/%s/%s",minBiasTrigers.Data(),select.Data(),selectRuns.Data())); \r
281 hNegMatchedMBBlow->Sumw2();\r
282 TH1 *hAllMatchedMBBlow= trackCounters->Draw("run",Form("v0mult:low/track:matched/acc:in/%s/%s/%s",minBiasTrigers.Data(),select.Data(),selectRuns.Data())); \r
283 hAllMatchedMBBlow->Sumw2();\r
284 TH1 *hPosMatchedMBBhigh = trackCounters->Draw("run",Form("v0mult:high/track:matched/charge:pos/acc:in/%s/%s/%s",minBiasTrigers.Data(),select.Data(),selectRuns.Data()));\r
285 hPosMatchedMBBhigh->Sumw2();\r
286 TH1 *hNegMatchedMBBhigh= trackCounters->Draw("run",Form("v0mult:high/track:matched/charge:neg/acc:in/%s/%s/%s",minBiasTrigers.Data(),select.Data(),selectRuns.Data())); \r
287 hNegMatchedMBBhigh->Sumw2();\r
288 TH1 *hAllMatchedMBBhigh= trackCounters->Draw("run",Form("v0mult:high/track:matched/acc:in/%s/%s/%s",minBiasTrigers.Data(),select.Data(),selectRuns.Data())); \r
289 hAllMatchedMBBhigh->Sumw2();\r
290 \r
291\r
292 TH1* hAll = eventCounters->Draw("trigger","run","run:any/selected:yes");\r
293 hAll->Draw("TEXT");\r
294 \r
295\r
296 //===================================================================================\r
297 // Put all plots in a pdf file, easier to publish (Twiki)\r
298 c1->Print(OutFileNamePDF_open.Data());\r
299\r
300 TFile *rootFileOut = TFile::Open(OutFileNameROOT.Data(),"RECREATE");\r
301 \r
302 //===================================================================================\r
303 // new canvas of the statistic wrt trigger w/ physics selection\r
304 CanvasName = LHCPeriod.Data() ; \r
305 CanvasName += "_StatByTrigger"; \r
306 TCanvas *cStatByTrigger = new TCanvas(CanvasName.Data(),"cStatByTrigger",1200,900);\r
307 cStatByTrigger->SetTopMargin(0.05);\r
308 cStatByTrigger->SetRightMargin(0.01);\r
309 cStatByTrigger->SetGridy(1);\r
310 cStatByTrigger->SetLogy(1);\r
311 cStatByTrigger->cd();\r
312\r
313 //hAllTriggersWithPS->SetLabelSize(LabelSize);\r
314 //hAllTriggersWithPS->SetTitle("Statistic per trigger w/ Phys. Sel."); \r
315 hAllTriggersWithPS->SetTitle(""); \r
316 hAllTriggersWithPS->SetLineWidth(2);\r
317 hAllTriggersWithPS->SetMinimum(10);\r
318 hAllTriggersWithPS->SetLabelSize(LabelSize);\r
319\r
320 hAllTriggersWithPS->Draw("E");\r
321 hMBBWithPS->SetLineWidth(2);\r
322 hMBBWithPS->SetLineColor(kBlue);\r
323 hMBBWithPS->Draw("ESAME");\r
324 hMBBlowWithPS->SetMarkerStyle(27);\r
325 hMBBlowWithPS->SetMarkerSize(2);\r
326 hMBBlowWithPS->SetMarkerColor(kBlue);\r
327 hMBBlowWithPS->SetLineColor(kBlue);\r
328 hMBBlowWithPS->Draw("ESAME");\r
329 hMBBhighWithPS->SetMarkerStyle(24);\r
330 hMBBhighWithPS->SetMarkerSize(2);\r
331 hMBBhighWithPS->SetMarkerColor(kBlue);\r
332 hMBBhighWithPS->SetLineColor(kBlue);\r
333 hMBBhighWithPS->Draw("ESAME");\r
334\r
335\r
336 TLegend* legcSBT = new TLegend(0.12,0.15,0.42,0.3);\r
337 legcSBT->SetHeader("Trigger Statistic w/ Phys. Sel.");\r
338 legcSBT->AddEntry(hAllTriggersWithPS,"All collisions triggers","l");\r
339 legcSBT->AddEntry(hMBBWithPS,"Min Bias [0-80]% from V0 amplitude","l");\r
340 legcSBT->AddEntry(hMBBlowWithPS,"MB low multiplicity events [60-80]%","p"); \r
341 legcSBT->AddEntry(hMBBhighWithPS,"MB high multiplicity events [0-10]%","p"); \r
342 legcSBT->Draw("same");\r
343\r
344 cStatByTrigger->Print(OutFileNamePDF.Data());\r
345 cStatByTrigger->Write();\r
346\r
347\r
348 //===================================================================================\r
349 // new canvas with the relative content of each trigger w/ and w/o physics selection\r
350 TH1 *ratioMBB = static_cast<TH1*>(hMBBWithPS->Clone("ratioMBB"));\r
351 ratioMBB->Divide(hMBBNoPS);\r
352 ratioMBB->SetLineColor(kBlack);\r
353 ratioMBB->SetLineWidth(2);\r
354 TH1 *ratioMBBlow = static_cast<TH1*>(hMBBlowWithPS->Clone("ratioMBBlow"));\r
355 ratioMBBlow->Divide(hMBBlowNoPS);\r
356 ratioMBBlow->SetLineColor(kBlue);\r
357 ratioMBBlow->SetLineWidth(2);\r
358 ratioMBBlow->SetMarkerSize(0);\r
359 ratioMBBlow->SetMarkerColor(kBlue);\r
360 TH1 *ratioMBBhigh = static_cast<TH1*>(hMBBhighWithPS->Clone("ratioMBBhigh"));\r
361 ratioMBBhigh->Divide(hMBBhighNoPS);\r
362 ratioMBBhigh->SetLineColor(kRed);\r
363 ratioMBBhigh->SetLineWidth(2);\r
364 ratioMBBhigh->SetMarkerSize(0);\r
365 ratioMBBhigh->SetMarkerColor(kRed);\r
366\r
367\r
368 TH1 *ratioMBBNoPS = static_cast<TH1*>(hMBBNoPS->Clone("ratioMBBNoPS"));\r
369 ratioMBBNoPS->Divide(hMBBNoPS);\r
370 ratioMBBNoPS->SetLineColor(kBlack);\r
371 ratioMBBNoPS->SetLineWidth(0);\r
372 ratioMBBNoPS->SetLineStyle(2);\r
373 ratioMBBNoPS->SetMarkerColor(kBlack);\r
374 ratioMBBNoPS->SetMarkerStyle(24);\r
375 ratioMBBNoPS->SetMarkerSize(1);\r
376 TH1 *ratioMBBlowNoPS = static_cast<TH1*>(hMBBlowNoPS->Clone("ratioMBBlowNoPS"));\r
377 ratioMBBlowNoPS->Divide(hMBBlowNoPS);\r
378 ratioMBBlowNoPS->SetLineColor(kBlue);\r
379 ratioMBBlowNoPS->SetLineWidth(0);\r
380 ratioMBBlowNoPS->SetLineStyle(2);\r
381 ratioMBBlowNoPS->SetMarkerColor(kBlue);\r
382 ratioMBBlowNoPS->SetMarkerStyle(24);\r
383 ratioMBBlowNoPS->SetMarkerSize(1);\r
384 TH1 *ratioMBBhighNoPS = static_cast<TH1*>(hMBBhighNoPS->Clone("ratioMBBhighNoPS"));\r
385 ratioMBBhighNoPS->Divide(hMBBhighNoPS);\r
386 ratioMBBhighNoPS->SetLineColor(kRed);\r
387 ratioMBBhighNoPS->SetLineWidth(0);\r
388 ratioMBBhighNoPS->SetLineStyle(2);\r
389 ratioMBBhighNoPS->SetMarkerColor(kRed);\r
390 ratioMBBhighNoPS->SetMarkerStyle(24);\r
391 ratioMBBhighNoPS->SetMarkerSize(1);\r
392 \r
393\r
394 TH1 *ratioBck = static_cast<TH1*>(hBckTriggersNoPS->Clone("ratioBck"));\r
395 ratioBck->Divide(hMBBNoPS);\r
396 ratioBck->SetLineColor(kBlack);\r
397 ratioBck->SetLineWidth(3);\r
398 ratioBck->SetLineStyle(3);\r
399 ratioBck->Scale(10);\r
400\r
401 CanvasName = LHCPeriod.Data() ; \r
402 CanvasName += "_RelativeTriggerContent"; \r
403 TCanvas *cRelativeTriggerContent = new TCanvas(CanvasName.Data(),"cRelativeTriggerContent",1200,900);\r
404 cRelativeTriggerContent->SetTopMargin(0.05);\r
405 cRelativeTriggerContent->SetRightMargin(0.01);\r
406 cRelativeTriggerContent->SetGridy(1);\r
407 //cRelativeTriggerContent->SetLogy(1);\r
408 cRelativeTriggerContent->cd();\r
409\r
410 ratioMBB->SetMaximum(1.0);\r
411 ratioMBB->SetMinimum(0.05);\r
412 ratioMBB->SetTitle("");\r
413 ratioMBB->SetLabelSize(LabelSize);\r
414 ratioMBB->LabelsOption("v");\r
415 ratioMBB->GetYaxis()->SetTitle("Relative trigger content w/ and w/o Phys. Sel."); \r
416\r
417 ratioMBB->Draw("E");\r
418 ratioMBBlow->Draw("ESAME");\r
419 ratioMBBhigh->Draw("ESAME");\r
420 ratioMBBNoPS->Draw("EPSAME");\r
421 ratioMBBlowNoPS->Draw("EPSAME");\r
422 ratioMBBhighNoPS->Draw("EPSAME");\r
423 ratioBck->Draw("EPSAME");\r
424\r
425 TLegend* legcRTC = new TLegend(0.4,0.25,0.70,0.45);\r
426 legcRTC->SetHeader("Physics Selection");\r
427 legcRTC->AddEntry(".","applied :","");\r
428 legcRTC->AddEntry(ratioMBB,"MBB","l");\r
429 legcRTC->AddEntry(ratioMBBlow,"MBB low mult.","l");\r
430 legcRTC->AddEntry(ratioMBBhigh,"MBB high mult.","l");\r
431 legcRTC->AddEntry(".","not applied :","");\r
432 legcRTC->AddEntry(ratioMBBNoPS,"MBB","p");\r
433 legcRTC->AddEntry(ratioMBBlowNoPS,"MBB low mult.","p");\r
434 legcRTC->AddEntry(ratioMBBhighNoPS,"MBB high mult.","p");\r
435 legcRTC->AddEntry(ratioBck,"Background (x10)","l");\r
436 legcRTC->Draw("same");\r
437 \r
438 cRelativeTriggerContent->Print(OutFileNamePDF.Data());\r
439 cRelativeTriggerContent->Write();\r
440\r
441\r
442 //===================================================================================\r
443 // new canvas with the relative ratio of multiplicity bin after physic selection\r
444 TH1 *relratioMBBlow = static_cast<TH1*>(hMBBlowWithPS->Clone("relratioMBBlow"));\r
445 relratioMBBlow->Divide(hMBBWithPS);\r
446 relratioMBBlow->SetLineColor(kBlue);\r
447 relratioMBBlow->SetLineWidth(2);\r
448 relratioMBBlow->SetMarkerSize(0);\r
449 TH1 *relratioMBBhigh = static_cast<TH1*>(hMBBhighWithPS->Clone("relratioMBBhigh"));\r
450 relratioMBBhigh->Divide(hMBBWithPS);\r
451 relratioMBBhigh->SetLineColor(kRed);\r
452 relratioMBBhigh->SetLineWidth(2);\r
453 relratioMBBhigh->SetMarkerSize(0);\r
454\r
455\r
456 TH1 *relratioMBBlowNoPS = static_cast<TH1*>(hMBBlowNoPS->Clone("relratioMBBlowNoPS"));\r
457 relratioMBBlowNoPS->Divide(hMBBNoPS);\r
458 relratioMBBlowNoPS->SetLineColor(kBlue);\r
459 relratioMBBlowNoPS->SetLineWidth(0);\r
460 relratioMBBlowNoPS->SetLineStyle(2);\r
461 relratioMBBlowNoPS->SetMarkerColor(kBlue);\r
462 relratioMBBlowNoPS->SetMarkerStyle(24);\r
463 relratioMBBlowNoPS->SetMarkerSize(1);\r
464 TH1 *relratioMBBhighNoPS = static_cast<TH1*>(hMBBhighNoPS->Clone("relratioMBBhighNoPS"));\r
465 relratioMBBhighNoPS->Divide(hMBBNoPS);\r
466 relratioMBBhighNoPS->SetLineColor(kRed);\r
467 relratioMBBhighNoPS->SetLineWidth(0);\r
468 relratioMBBhighNoPS->SetLineStyle(2);\r
469 relratioMBBhighNoPS->SetMarkerColor(kRed);\r
470 relratioMBBhighNoPS->SetMarkerStyle(24);\r
471 relratioMBBhighNoPS->SetMarkerSize(1);\r
472 \r
473 CanvasName = LHCPeriod.Data() ; \r
474 CanvasName += "_CentralityPercentileCheck"; \r
475 TCanvas *cCentralityPercentileCheck = new TCanvas(CanvasName.Data(),"CentralityPercentileCheck",1200,900);\r
476 cCentralityPercentileCheck->SetTopMargin(0.05);\r
477 cCentralityPercentileCheck->SetRightMargin(0.01);\r
478 cCentralityPercentileCheck->SetGridy(1);\r
479 //cCentralityPercentileCheck->SetLogy(1);\r
480 cCentralityPercentileCheck->cd();\r
481\r
482 relratioMBBlow->Scale(0.8);\r
483 relratioMBBhigh->Scale(0.8);\r
484 relratioMBBlowNoPS->Scale(0.8);\r
485 relratioMBBhighNoPS->Scale(0.8);\r
486\r
487 relratioMBBlow->SetMaximum(0.3);\r
488 relratioMBBlow->SetMinimum(0.01);\r
489 relratioMBBlow->SetTitle("");\r
490 relratioMBBlow->SetLabelSize(LabelSize);\r
491 relratioMBBlow->GetYaxis()->SetTitle("Centrality percentile check"); \r
492\r
493 relratioMBBlow->Draw();\r
494 relratioMBBhigh->Draw("ESAME");\r
495 relratioMBBlowNoPS->Draw("EPSAME");\r
496 relratioMBBhighNoPS->Draw("EPSAME");\r
497\r
498 TLegend* legcCPC = new TLegend(0.12,0.15,0.42,0.30);\r
499 legcCPC->SetHeader("Physics Selection");\r
500 legcCPC->AddEntry(".","applied :","");\r
501 legcCPC->AddEntry(relratioMBBlow,"MBB low mult.","l");\r
502 legcCPC->AddEntry(relratioMBBhigh,"MBB high mult.","l");\r
503 legcCPC->AddEntry(".","not applied :","");\r
504 legcCPC->AddEntry(relratioMBBlowNoPS,"MBB low mult.","p");\r
505 legcCPC->AddEntry(relratioMBBhighNoPS,"MBB high mult.","p");\r
506 legcCPC->Draw("same");\r
507 \r
508 cCentralityPercentileCheck->Print(OutFileNamePDF.Data());\r
509 cCentralityPercentileCheck->Write();\r
510\r
511\r
512 //====================================================\r
513 // Draw ratio of tracks over MBB versus run \r
514 TH1* hTrackerPerMBB= static_cast<TH1*>(hTrackerMBB->Clone("hTrackerPerMBB"));\r
515 hTrackerPerMBB->Divide(hMBB);\r
516 hTrackerPerMBB->SetLineWidth(2);\r
517 hTrackerPerMBB->SetLineColor(kRed);\r
518\r
519 TH1* hTriggerPerMBB= static_cast<TH1*>(hTriggerMBB->Clone("hTriggerPerMBB"));\r
520 hTriggerPerMBB->Divide(hMBB);\r
521 hTriggerPerMBB->SetLineWidth(2);\r
522 hTriggerPerMBB->SetLineColor(kBlue);\r
523\r
524\r
525 TH1* hMatchedPerMBB= static_cast<TH1*>(hMatchedMBB->Clone("hMatchedPerMBB"));\r
526 hMatchedPerMBB->Divide(hMBB);\r
527 hMatchedPerMBB->SetLineWidth(2);\r
528 hMatchedPerMBB->SetLineColor(kViolet);\r
529\r
530\r
531 TH1* hAllTracksPerMBB= static_cast<TH1*>(hAllTracksMBB->Clone("hAllTracksPerMBB"));\r
532 hAllTracksPerMBB->Divide(hMBB);\r
533 hAllTracksPerMBB->SetLineWidth(3);\r
534 hAllTracksPerMBB->SetLineColor(kBlack);\r
535\r
536\r
537 CanvasName = LHCPeriod.Data() ; \r
538 CanvasName += "_RatioTrackTypesMBB"; \r
539 TCanvas *cRatioTrackTypesMBB = new TCanvas(CanvasName.Data(),"cRatioTrackTypesMBB",1200,900);\r
540 cRatioTrackTypesMBB->SetRightMargin(0.01);\r
541 cRatioTrackTypesMBB->SetGridy(1);\r
542 //cRatioTrackTypesMBB->SetLogy(1);\r
543 cRatioTrackTypesMBB->cd();\r
544\r
545 hAllTracksPerMBB->SetTitle("Ratio (Number of Tracks)/MBB for [0-80]% centrality");\r
546 hAllTracksPerMBB->SetMaximum(4);\r
547 hAllTracksPerMBB->SetMinimum(0.1);\r
548 hAllTracksPerMBB->SetLabelSize(LabelSize);\r
549 hAllTracksPerMBB->Draw("E");\r
550\r
551 hTrackerPerMBB->Draw("Esame");\r
552 hMatchedPerMBB->Draw("Esame");\r
553 hTriggerPerMBB->Draw("Esame");\r
554\r
555 TLegend* legcTTMBB = new TLegend(0.15,0.50,0.35,0.70);\r
556 legcTTMBB->AddEntry(hAllTracksPerMBB,"All tracks","l");\r
557 legcTTMBB->AddEntry(hTrackerPerMBB,"Tracking (only) tracks","l");\r
558 legcTTMBB->AddEntry(hMatchedPerMBB,"Matched tracks","l");\r
559 legcTTMBB->AddEntry(hTriggerPerMBB,"Trigger (only) tracks","l");\r
560 legcTTMBB->Draw("same");\r
561\r
562 cRatioTrackTypesMBB->Print(OutFileNamePDF.Data());\r
563 cRatioTrackTypesMBB->Write();\r
564\r
565 //==========================================================================\r
566 // Draw ratio of tracks over MBB versus runs and low and high multiplicities \r
567\r
568 TH1* hTrackerPerMBBlow= static_cast<TH1*>(hTrackerMBBlow->Clone("hTrackerPerMBBlow"));\r
569 hTrackerPerMBBlow->Divide(hMBBlow);\r
570 hTrackerPerMBBlow->SetLineWidth(2);\r
571 hTrackerPerMBBlow->SetLineStyle(2);\r
572 hTrackerPerMBBlow->SetLineColor(kRed);\r
573 hTrackerPerMBBlow->SetMarkerColor(kRed);\r
574 hTrackerPerMBBlow->SetMarkerSize(2);\r
575 hTrackerPerMBBlow->SetMarkerStyle(27);\r
576 TH1* hTrackerPerMBBhigh= static_cast<TH1*>(hTrackerMBBhigh->Clone("hTrackerPerMBBhigh"));\r
577 hTrackerPerMBBhigh->Divide(hMBBhigh);\r
578 hTrackerPerMBBhigh->SetLineWidth(2);\r
579 hTrackerPerMBBhigh->SetLineStyle(2);\r
580 hTrackerPerMBBhigh->SetLineColor(kRed);\r
581 hTrackerPerMBBhigh->SetMarkerColor(kRed);\r
582 hTrackerPerMBBhigh->SetMarkerSize(2);\r
583 hTrackerPerMBBhigh->SetMarkerStyle(24);\r
584\r
585\r
586 TH1* hTriggerPerMBBlow= static_cast<TH1*>(hTriggerMBBlow->Clone("hTriggerPerMBBlow"));\r
587 hTriggerPerMBBlow->Divide(hMBBlow);\r
588 hTriggerPerMBBlow->SetLineWidth(2);\r
589 hTriggerPerMBBlow->SetLineStyle(2);\r
590 hTriggerPerMBBlow->SetLineColor(kBlue);\r
591 hTriggerPerMBBlow->SetMarkerColor(kBlue);\r
592 hTriggerPerMBBlow->SetMarkerSize(2);\r
593 hTriggerPerMBBlow->SetMarkerStyle(27);\r
594 TH1* hTriggerPerMBBhigh= static_cast<TH1*>(hTriggerMBBhigh->Clone("hTriggerPerMBBhigh"));\r
595 hTriggerPerMBBhigh->Divide(hMBBhigh);\r
596 hTriggerPerMBBhigh->SetLineWidth(2);\r
597 hTriggerPerMBBhigh->SetLineStyle(2);\r
598 hTriggerPerMBBhigh->SetLineColor(kBlue);\r
599 hTriggerPerMBBhigh->SetMarkerColor(kBlue);\r
600 hTriggerPerMBBhigh->SetMarkerSize(2);\r
601 hTriggerPerMBBhigh->SetMarkerStyle(24);\r
602\r
603\r
604 TH1* hMatchedPerMBBlow= static_cast<TH1*>(hMatchedMBBlow->Clone("hMatchedPerMBBlow"));\r
605 hMatchedPerMBBlow->Divide(hMBBlow);\r
606 hMatchedPerMBBlow->SetLineWidth(2);\r
607 hMatchedPerMBBlow->SetLineStyle(2);\r
608 hMatchedPerMBBlow->SetLineColor(kViolet);\r
609 hMatchedPerMBBlow->SetMarkerColor(kViolet);\r
610 hMatchedPerMBBlow->SetMarkerSize(2);\r
611 hMatchedPerMBBlow->SetMarkerStyle(27);\r
612 TH1* hMatchedPerMBBhigh= static_cast<TH1*>(hMatchedMBBhigh->Clone("hMatchedPerMBBhigh"));\r
613 hMatchedPerMBBhigh->Divide(hMBBhigh);\r
614 hMatchedPerMBBhigh->SetLineWidth(2);\r
615 hMatchedPerMBBhigh->SetLineStyle(2);\r
616 hMatchedPerMBBhigh->SetLineColor(kViolet);\r
617 hMatchedPerMBBhigh->SetMarkerColor(kViolet);\r
618 hMatchedPerMBBhigh->SetMarkerSize(2);\r
619 hMatchedPerMBBhigh->SetMarkerStyle(24);\r
620\r
621\r
622 TH1* hAllTracksPerMBBlow= static_cast<TH1*>(hAllTracksMBBlow->Clone("hAllTracksPerMBBlow"));\r
623 hAllTracksPerMBBlow->Divide(hMBBlow);\r
624 hAllTracksPerMBBlow->SetLineWidth(2);\r
625 hAllTracksPerMBBlow->SetLineStyle(2);\r
626 hAllTracksPerMBBlow->SetLineColor(kBlack);\r
627 hAllTracksPerMBBlow->SetMarkerColor(kBlack);\r
628 hAllTracksPerMBBlow->SetMarkerSize(2);\r
629 hAllTracksPerMBBlow->SetMarkerStyle(27);\r
630 TH1* hAllTracksPerMBBhigh= static_cast<TH1*>(hAllTracksMBBhigh->Clone("hAllTracksPerMBBhigh"));\r
631 hAllTracksPerMBBhigh->Divide(hMBBhigh);\r
632 hAllTracksPerMBBhigh->SetLineWidth(2);\r
633 hAllTracksPerMBBhigh->SetLineStyle(2);\r
634 hAllTracksPerMBBhigh->SetLineColor(kBlack);\r
635 hAllTracksPerMBBhigh->SetMarkerColor(kBlack);\r
636 hAllTracksPerMBBhigh->SetMarkerSize(2);\r
637 hAllTracksPerMBBhigh->SetMarkerStyle(24);\r
638\r
639\r
640 CanvasName = LHCPeriod.Data() ; \r
641 CanvasName += "_RatioTrackTypesMBBVsMult"; \r
642 TCanvas *cRatioTrackTypesMBBVsMult = new TCanvas(CanvasName.Data(),"cRatioTrackTypesMBBVsMult",1200,900);\r
643 cRatioTrackTypesMBBVsMult->SetRightMargin(0.01);\r
644 cRatioTrackTypesMBBVsMult->SetGridy(1);\r
645 //cRatioTrackTypesMBBVsMult->SetLogy(1);\r
646\r
647\r
648 cRatioTrackTypesMBBVsMult->Divide(1,2);\r
649\r
650\r
651 cRatioTrackTypesMBBVsMult->cd(1);\r
652\r
653 hAllTracksPerMBBlow->SetTitle("Ratio (Number of Tracks)/MBB low mult");\r
654 hAllTracksPerMBBlow->SetMaximum(0.30);\r
655 hAllTracksPerMBBlow->SetMinimum(0.001);\r
656 hAllTracksPerMBBlow->SetLabelSize(LabelSize);\r
657 hAllTracksPerMBBlow->Draw("E");\r
658 hTrackerPerMBBlow->Draw("Esame");\r
659 hTriggerPerMBBlow->Draw("Esame");\r
660 hMatchedPerMBBlow->Draw("Esame");\r
661\r
662 TLegend* legcTTMBBlow = new TLegend(0.40,0.8,0.60,1.);\r
663 legcTTMBBlow->AddEntry(hAllTracksPerMBBlow,"All tracks","p");\r
664 legcTTMBBlow->AddEntry(hTrackerPerMBBlow,"Tracking (only) tracks","p");\r
665 legcTTMBBlow->AddEntry(hMatchedPerMBBlow,"Matched tracks","p");\r
666 legcTTMBBlow->AddEntry(hTriggerPerMBBlow,"Trigger (only) tracks","p");\r
667 legcTTMBBlow->Draw("same");\r
668 \r
669\r
670 cRatioTrackTypesMBBVsMult->cd(2);\r
671 hAllTracksPerMBBhigh->SetTitle("Ratio (Number of Tracks)/MBB high mult");\r
672 hAllTracksPerMBBhigh->SetMaximum(15);\r
673 hAllTracksPerMBBhigh->SetMinimum(0.01);\r
674 hAllTracksPerMBBhigh->SetLabelSize(LabelSize);\r
675 hAllTracksPerMBBhigh->Draw("E");\r
676 hTrackerPerMBBhigh->Draw("Esame");\r
677 hTriggerPerMBBhigh->Draw("Esame");\r
678 hMatchedPerMBBhigh->Draw("Esame");\r
679\r
680\r
681\r
682 TLegend* legcTTMBBhigh = new TLegend(0.40,0.8,0.60,1.);\r
683 legcTTMBBhigh->AddEntry(hAllTracksPerMBBhigh,"All tracks","p");\r
684 legcTTMBBhigh->AddEntry(hTrackerPerMBBhigh,"Tracking (only) tracks","p");\r
685 legcTTMBBhigh->AddEntry(hMatchedPerMBBhigh,"Matched tracks","p");\r
686 legcTTMBBhigh->AddEntry(hTriggerPerMBBhigh,"Trigger (only) tracks","p");\r
687 legcTTMBBhigh->Draw("same");\r
688 \r
689\r
690 // cRatioTrackTypesMBBVsMult->cd();\r
691 // TPaveText pt1(0.2,0.2,0.7,0.9);\r
692 // TText* t1=pt1.AddText("All tracks (tracking-only + trigger-only + matched)");\r
693 // t1->SetTextColor(kBlack);\r
694 // pt1.Draw();\r
695 // break;\r
696\r
697 cRatioTrackTypesMBBVsMult->Print(OutFileNamePDF.Data());\r
698 cRatioTrackTypesMBBVsMult->Write();\r
699\r
700\r
701 //===========================================\r
702 // draw ratio of track (summed, i.e all trigger=(matched+trigger-only)) for MB versus run for low mult\r
703 CanvasName = LHCPeriod.Data() ; \r
704 CanvasName += "_TrackMultBBlow"; \r
705 TCanvas* cTrackMultBBlow = new TCanvas("cTrackMultBBlow","cTrackMultBBlow",1200,900);\r
706 // must be applied on the pads created by divide\r
707 cTrackMultBBlow->Divide(1,2);\r
708 cTrackMultBBlow->cd(1);\r
709\r
710 TH1* hSumTriggerOverMBBlow = static_cast<TH1*>(hTriggerMBBlow->Clone("hSumTriggerOverMBBlow"));\r
711 hSumTriggerOverMBBlow->Add(hMatchedMBBlow);\r
712 hSumTriggerOverMBBlow->Divide(hMBBlow);\r
713 hSumTriggerOverMBBlow->SetTitle("Sum of trigger tracks (matched+trigger-only) in MBB events / # MBB events in low mult.");\r
714 //hSumTriggerOverMBBlow->LabelsOption("u");\r
715 hSumTriggerOverMBBlow->SetLabelSize(LabelSize);\r
716 hSumTriggerOverMBBlow->SetLineWidth(2);\r
717 hSumTriggerOverMBBlow->SetLineColor(kBlue);\r
718 hSumTriggerOverMBBlow->Draw("e");\r
719 cTrackMultBBlow->cd(2);\r
720 TH1* hSumTrackerOverMBBlow = static_cast<TH1*>(hTrackerMBBlow->Clone("hSumTrackerOverMBBlow"));\r
721 hSumTrackerOverMBBlow->Add(hMatchedMBBlow);\r
722 hSumTrackerOverMBBlow->Divide(hMBBlow);\r
723 hSumTrackerOverMBBlow->SetTitle("Sum tracker tracks (matched+tracker-only) in MBB events / # MBB events in low mult.");\r
724 //hSumTrackerOverMBBlow->LabelsOption("u");\r
725 hSumTrackerOverMBBlow->SetLabelSize(LabelSize);\r
726 hSumTrackerOverMBBlow->SetLineWidth(2);\r
727 hSumTrackerOverMBBlow->SetLineColor(kBlue);\r
728 hSumTrackerOverMBBlow->Draw("e");\r
729\r
730\r
731 cTrackMultBBlow->Print(OutFileNamePDF.Data());\r
732 cTrackMultBBlow->Write();\r
733\r
734\r
735\r
736 //===========================================\r
737 // draw ratio of track (summed, i.e all trigger=(matched+trigger-only)) for MB versus run for high mult\r
738 CanvasName = LHCPeriod.Data() ; \r
739 CanvasName += "_TrackMultBBhigh"; \r
740 TCanvas* cTrackMultBBhigh = new TCanvas("cTrackMultBBhigh","cTrackMultBBhigh",1200,900);\r
741 // must be applied on the pads created by divide\r
742 cTrackMultBBhigh->Divide(1,2);\r
743 cTrackMultBBhigh->cd(1);\r
744\r
745 TH1* hSumTriggerOverMBBhigh = static_cast<TH1*>(hTriggerMBBhigh->Clone("hSumTriggerOverMBBhigh"));\r
746 hSumTriggerOverMBBhigh->Add(hMatchedMBBhigh);\r
747 hSumTriggerOverMBBhigh->Divide(hMBBhigh);\r
748 hSumTriggerOverMBBhigh->SetTitle("Sum of trigger tracks (matched+trigger-only) in MBB events / # MBB events in high mult.");\r
749 //hSumTriggerOverMBBhigh->LabelsOption("u");\r
750 hSumTriggerOverMBBhigh->SetLabelSize(LabelSize);\r
751 hSumTriggerOverMBBhigh->SetLineWidth(2);\r
752 hSumTriggerOverMBBhigh->SetLineColor(kBlue);\r
753 hSumTriggerOverMBBhigh->Draw("e");\r
754 cTrackMultBBhigh->cd(2);\r
755 TH1* hSumTrackerOverMBBhigh = static_cast<TH1*>(hTrackerMBBhigh->Clone("hSumTrackerOverMBBhigh"));\r
756 hSumTrackerOverMBBhigh->Add(hMatchedMBBhigh);\r
757 hSumTrackerOverMBBhigh->Divide(hMBBhigh);\r
758 hSumTrackerOverMBBhigh->SetTitle("Sum tracker tracks (matched+tracker-only) in MBB events / # MBB events in high mult.");\r
759 //hSumTrackerOverMBBhigh->LabelsOption("u");\r
760 hSumTrackerOverMBBhigh->SetLabelSize(LabelSize);\r
761 hSumTrackerOverMBBhigh->SetLineWidth(2);\r
762 hSumTrackerOverMBBhigh->SetLineColor(kBlue);\r
763 hSumTrackerOverMBBhigh->Draw("e");\r
764\r
765\r
766 cTrackMultBBhigh->Print(OutFileNamePDF.Data());\r
767 cTrackMultBBhigh->Write();\r
768\r
769 \r
770 \r
771 //===========================================\r
772 // draw mixed ratio of track over track versus run for MBB in low mult \r
773 CanvasName = LHCPeriod.Data() ; \r
774 CanvasName += "_RatioTrackMBBlow"; \r
775 TCanvas* cRatioTrackMBBlow = new TCanvas("cRatioTrackMBBlow","cRatioTrackMBBlow",1200,900);\r
776 cRatioTrackMBBlow->Divide(1,3);\r
777\r
778 cRatioTrackMBBlow->cd(1);\r
779 TH1* hTrackerOverTriggerMBBlow = static_cast<TH1*>(hTrackerMBBlow->Clone("hTrackerOverTriggerMBBlow"));\r
780 hTrackerOverTriggerMBBlow->Divide(hTriggerMBBlow);\r
781 hTrackerOverTriggerMBBlow->SetTitle("# tracker tracks / # trigger tracks in MBB low mult.");\r
782 //hTrackerOverTriggerMBBlow->LabelsOption("u");\r
783 hTrackerOverTriggerMBBlow->SetLabelSize(LabelSize);\r
784 hTrackerOverTriggerMBBlow->SetLineWidth(2);\r
785 hTrackerOverTriggerMBBlow->SetLineColor(kBlue);\r
786 hTrackerOverTriggerMBBlow->Draw("e");\r
787\r
788 cRatioTrackMBBlow->cd(2);\r
789 TH1* hMatchedOverTriggerMBBlow = static_cast<TH1*>(hMatchedMBBlow->Clone("hMatchedOverTriggerMBBlow"));\r
790 hMatchedOverTriggerMBBlow->Divide(hTriggerMBBlow);\r
791 hMatchedOverTriggerMBBlow->SetTitle("# matched tracks / # trigger tracks in MBB low mult.");\r
792 //hMatchedOverTriggerMBBlow->LabelsOption("u");\r
793 hMatchedOverTriggerMBBlow->SetLabelSize(LabelSize);\r
794 hMatchedOverTriggerMBBlow->SetLineWidth(2);\r
795 hMatchedOverTriggerMBBlow->SetLineColor(kBlue);\r
796 hMatchedOverTriggerMBBlow->Draw("e");\r
797\r
798 cRatioTrackMBBlow->cd(3);\r
799 TH1* hMatchedOverTrackerMBBlow = static_cast<TH1*>(hMatchedMBBlow->Clone("hMatchedOverTrackerMBBlow"));\r
800 hMatchedOverTrackerMBBlow->Divide(hTrackerMBBlow);\r
801 hMatchedOverTrackerMBBlow->SetTitle("# matched tracks / # tracker tracks in MBB low mult.");\r
802 //hMatchedOverTrackerMBBlow->LabelsOption("u");\r
803 hMatchedOverTrackerMBBlow->SetLabelSize(LabelSize);\r
804 hMatchedOverTrackerMBBlow->SetLineWidth(2);\r
805 hMatchedOverTrackerMBBlow->SetLineColor(kBlue);\r
806 hMatchedOverTrackerMBBlow->Draw("e");\r
807 \r
808 cRatioTrackMBBlow->Print(OutFileNamePDF.Data());\r
809 cRatioTrackMBBlow->Write();\r
810\r
811 \r
812 //===========================================\r
813 // draw mixed ratio of track over track versus run for MBB in high mult \r
814 CanvasName = LHCPeriod.Data() ; \r
815 CanvasName += "_RatioTrackMBBhigh"; \r
816 TCanvas* cRatioTrackMBBhigh = new TCanvas("cRatioTrackMBBhigh","cRatioTrackMBBhigh",1200,900);\r
817 cRatioTrackMBBhigh->Divide(1,3);\r
818\r
819 cRatioTrackMBBhigh->cd(1);\r
820 TH1* hTrackerOverTriggerMBBhigh = static_cast<TH1*>(hTrackerMBBhigh->Clone("hTrackerOverTriggerMBBhigh"));\r
821 hTrackerOverTriggerMBBhigh->Divide(hTriggerMBBhigh);\r
822 hTrackerOverTriggerMBBhigh->SetTitle("# tracker tracks / # trigger tracks in MBB high mult.");\r
823 //hTrackerOverTriggerMBBhigh->LabelsOption("u");\r
824 hTrackerOverTriggerMBBhigh->SetLabelSize(LabelSize);\r
825 hTrackerOverTriggerMBBhigh->SetLineWidth(2);\r
826 hTrackerOverTriggerMBBhigh->SetLineColor(kBlue);\r
827 hTrackerOverTriggerMBBhigh->Draw("e");\r
828\r
829 cRatioTrackMBBhigh->cd(2);\r
830 TH1* hMatchedOverTriggerMBBhigh = static_cast<TH1*>(hMatchedMBBhigh->Clone("hMatchedOverTriggerMBBhigh"));\r
831 hMatchedOverTriggerMBBhigh->Divide(hTriggerMBBhigh);\r
832 hMatchedOverTriggerMBBhigh->SetTitle("# matched tracks / # trigger tracks in MBB high mult.");\r
833 //hMatchedOverTriggerMBBhigh->LabelsOption("u");\r
834 hMatchedOverTriggerMBBhigh->SetLabelSize(LabelSize);\r
835 hMatchedOverTriggerMBBhigh->SetLineWidth(2);\r
836 hMatchedOverTriggerMBBhigh->SetLineColor(kBlue);\r
837 hMatchedOverTriggerMBBhigh->Draw("e");\r
838\r
839 cRatioTrackMBBhigh->cd(3);\r
840 TH1* hMatchedOverTrackerMBBhigh = static_cast<TH1*>(hMatchedMBBhigh->Clone("hMatchedOverTrackerMBBhigh"));\r
841 hMatchedOverTrackerMBBhigh->Divide(hTrackerMBBhigh);\r
842 hMatchedOverTrackerMBBhigh->SetTitle("# matched tracks / # tracker tracks in MBB high mult.");\r
843 //hMatchedOverTrackerMBBhigh->LabelsOption("u");\r
844 hMatchedOverTrackerMBBhigh->SetLabelSize(LabelSize);\r
845 hMatchedOverTrackerMBBhigh->SetLineWidth(2);\r
846 hMatchedOverTrackerMBBhigh->SetLineColor(kBlue);\r
847 hMatchedOverTrackerMBBhigh->Draw("e");\r
848 \r
849 cRatioTrackMBBhigh->Print(OutFileNamePDF.Data());\r
850 cRatioTrackMBBhigh->Write();\r
851\r
852\r
853\r
854\r
855\r
856// //===========================================\r
857// // Draw ratio of tracks over CMUS1B versus run\r
858// TH1* hTrackerPerCMUS1B= static_cast<TH1*>(hTrackerCMUS1B->Clone("hTrackerPerCMUS1B"));\r
859// hTrackerPerCMUS1B->Divide(hCMUS1B);\r
860// hTrackerPerCMUS1B->SetLineWidth(2);\r
861// hTrackerPerCMUS1B->SetLineColor(kRed);\r
862\r
863// TH1* hTriggerPerCMUS1B= static_cast<TH1*>(hTriggerCMUS1B->Clone("hTriggerPerCMUS1B"));\r
864// hTriggerPerCMUS1B->Divide(hCMUS1B);\r
865// hTriggerPerCMUS1B->SetLineWidth(2);\r
866// hTriggerPerCMUS1B->SetLineColor(kBlue);\r
867\r
868// TH1* hMatchedPerCMUS1B= static_cast<TH1*>(hMatchedCMUS1B->Clone("hMatchedPerCMUS1B"));\r
869// hMatchedPerCMUS1B->Divide(hCMUS1B);\r
870// hMatchedPerCMUS1B->SetLineWidth(2);\r
871// hMatchedPerCMUS1B->SetLineColor(kViolet);\r
872\r
873\r
874// TH1* hAllTracksPerCMUS1B= static_cast<TH1*>(hAllTracksCMUS1B->Clone("hAllTracksPerCMUS1B"));\r
875// hAllTracksPerCMUS1B->Divide(hCMUS1B);\r
876// hAllTracksPerCMUS1B->SetLineWidth(3);\r
877// hAllTracksPerCMUS1B->SetLineColor(kBlack);\r
878\r
879\r
880// CanvasName = LHCPeriod.Data() ; \r
881// CanvasName += "_RatioTrackTypesCMUS1B"; \r
882// TCanvas *cRatioTrackTypesCMUS1B = new TCanvas(CanvasName.Data(),"cRatioTrackTypesCMUS1B",1200,900);\r
883// cRatioTrackTypesCMUS1B->SetRightMargin(0.01);\r
884// cRatioTrackTypesCMUS1B->SetGridy(1);\r
885// cRatioTrackTypesCMUS1B->cd();\r
886\r
887// hAllTracksPerCMUS1B->SetTitle("Ratio (Number of Tracks)/CMUS1B");\r
888// hAllTracksPerCMUS1B->SetMinimum(0.01);\r
889// hAllTracksPerCMUS1B->SetLabelSize(LabelSize);\r
890// hAllTracksPerCMUS1B->Draw("E");\r
891\r
892// hTrackerPerCMUS1B->Draw("Esame");\r
893// hMatchedPerCMUS1B->Draw("Esame");\r
894// hTriggerPerCMUS1B->Draw("Esame");\r
895\r
896// TLegend* legcTTCMUS1B = new TLegend(0.75,0.55,0.90,0.75);\r
897// legcTTCMUS1B->AddEntry(hAllTracksPerCMUS1B,"All tracks","l");\r
898// legcTTCMUS1B->AddEntry(hTrackerPerCMUS1B,"Tracking (only) tracks","l");\r
899// legcTTCMUS1B->AddEntry(hMatchedPerCMUS1B,"Matched tracks","l");\r
900// legcTTCMUS1B->AddEntry(hTriggerPerCMUS1B,"Trigger (only) tracks","l");\r
901// legcTTCMUS1B->Draw("same");\r
902\r
903// cRatioTrackTypesCMUS1B->Print(OutFileNamePDF.Data());\r
904// cRatioTrackTypesCMUS1B->Write();\r
905\r
906\r
907// //===========================================\r
908// // draw ratio of track (summed, i.e all trigger=(matched+trigger-only)) over CMUS1B versus run\r
909// TCanvas* cTrackMultCMUS1B = new TCanvas("cTrackMultCMUS1B","cTrackMultCMUS1B",1200,900);\r
910// cTrackMultCMUS1B->Divide(1,2);\r
911// cTrackMultCMUS1B->cd(1);\r
912// TH1* hSumTriggerOverCMUS1B = static_cast<TH1*>(hTriggerCMUS1B->Clone("hSumTriggerOverCMUS1B"));\r
913// hSumTriggerOverCMUS1B->Add(hMatchedCMUS1B);\r
914// hSumTriggerOverCMUS1B->Divide(hCMUS1B);\r
915// hSumTriggerOverCMUS1B->SetTitle("Sum of trigger tracks (matched+trigger-only) in CMUS1B events / # CMUS1B events");\r
916// //hSumTriggerOverCMUS1B->LabelsOption("u");\r
917// hSumTriggerOverCMUS1B->SetLabelSize(LabelSize);\r
918// hSumTriggerOverCMUS1B->SetLineWidth(2);\r
919// hSumTriggerOverCMUS1B->SetLineColor(kRed);\r
920// hSumTriggerOverCMUS1B->Draw("e");\r
921\r
922// cTrackMultCMUS1B->cd(2);\r
923// TH1* hSumTrackerOverCMUS1B = static_cast<TH1*>(hTrackerCMUS1B->Clone("hSumTrackerOverCMUS1B"));\r
924// hSumTrackerOverCMUS1B->Add(hMatchedCMUS1B);\r
925// hSumTrackerOverCMUS1B->Divide(hCMUS1B);\r
926// hSumTrackerOverCMUS1B->SetTitle("Sum of tracker tracks (matched+tracker-only) in CMUS1B events / # CMUS1B events");\r
927// //hSumTrackerOverCMUS1B->LabelsOption("u");\r
928// hSumTrackerOverCMUS1B->SetLabelSize(LabelSize);\r
929// hSumTrackerOverCMUS1B->SetLineWidth(2);\r
930// hSumTrackerOverCMUS1B->SetLineColor(kRed);\r
931// hSumTrackerOverCMUS1B->Draw("e");\r
932\r
933// cTrackMultCMUS1B->Print(OutFileNamePDF.Data());\r
934// cTrackMultCMUS1B->Write();\r
935 \r
936// //===========================================\r
937// // draw mixed ratio of track over track versus run for CMUS1B\r
938// TCanvas* cRatioTrackCMUS1B = new TCanvas("cRatioTrackCMUS1B","cRatioTrackCMUS1B",1200,900);\r
939// cRatioTrackCMUS1B->Divide(1,3);\r
940\r
941// cRatioTrackCMUS1B->cd(1);\r
942// TH1* hTrackerOverTriggerCMUS1B = static_cast<TH1*>(hTrackerCMUS1B->Clone("hTrackerOverTriggerCMUS1B"));\r
943// hTrackerOverTriggerCMUS1B->Divide(hTriggerCMUS1B);\r
944// hTrackerOverTriggerCMUS1B->SetTitle("# tracker tracks / # trigger tracks in CMUS1B");\r
945// //hTrackerOverTriggerCMUS1B->LabelsOption("u");\r
946// hTrackerOverTriggerCMUS1B->SetLabelSize(LabelSize);\r
947// hTrackerOverTriggerCMUS1B->SetLineWidth(2);\r
948// hTrackerOverTriggerCMUS1B->SetLineColor(kRed);\r
949// hTrackerOverTriggerCMUS1B->Draw("e");\r
950\r
951// cRatioTrackCMUS1B->cd(2);\r
952// TH1* hMatchedOverTriggerCMUS1B = static_cast<TH1*>(hMatchedCMUS1B->Clone("hMatchedOverTriggerCMUS1B"));\r
953// hMatchedOverTriggerCMUS1B->Divide(hTriggerCMUS1B);\r
954// hMatchedOverTriggerCMUS1B->SetTitle("# matched tracks / # trigger tracks in CMUS1B");\r
955// //hMatchedOverTriggerCMUS1B->LabelsOption("u");\r
956// hMatchedOverTriggerCMUS1B->SetLabelSize(LabelSize);\r
957// hMatchedOverTriggerCMUS1B->SetLineWidth(2);\r
958// hMatchedOverTriggerCMUS1B->SetLineColor(kRed);\r
959// hMatchedOverTriggerCMUS1B->Draw("e");\r
960\r
961// cRatioTrackCMUS1B->cd(3);\r
962// TH1* hMatchedOverTrackerCMUS1B = static_cast<TH1*>(hMatchedCMUS1B->Clone("hMatchedOverTrackerCMUS1B"));\r
963// hMatchedOverTrackerCMUS1B->Divide(hTrackerCMUS1B);\r
964// hMatchedOverTrackerCMUS1B->SetTitle("# matched tracks / # tracker tracks in CMUS1B");\r
965// //hMatchedOverTrackerCMUS1B->LabelsOption("u");\r
966// hMatchedOverTrackerCMUS1B->SetLabelSize(LabelSize);\r
967// hMatchedOverTrackerCMUS1B->SetLineWidth(2);\r
968// hMatchedOverTrackerCMUS1B->SetLineColor(kRed);\r
969// hMatchedOverTrackerCMUS1B->Draw("e");\r
970 \r
971// cRatioTrackCMUS1B->Print(OutFileNamePDF.Data());\r
972// cRatioTrackCMUS1B->Write();\r
973\r
974 //==================================================\r
975 // Draw matched tracks charge asymetry in low mult for MBB triggers \r
976 TH1 *hDiffMatchedMBBlow= static_cast<TH1*>(hPosMatchedMBBlow->Clone("hDiffMatchedMBBlow"));\r
977 hDiffMatchedMBBlow->Add(hNegMatchedMBBlow,-1);\r
978 hDiffMatchedMBBlow->Sumw2();\r
979 \r
980 TH1 *hAsymMatchedMBBlow= static_cast<TH1*>(hDiffMatchedMBBlow->Clone("hAsymMatchedMBBlow"));\r
981 hAsymMatchedMBBlow->Divide(hAllMatchedMBBlow);\r
982 hAsymMatchedMBBlow->SetLineColor(kBlue);\r
983 hAsymMatchedMBBlow->SetLineWidth(2);\r
984 hAsymMatchedMBBlow->SetMinimum(-0.2);\r
985 hAsymMatchedMBBlow->SetMaximum(0.2);\r
986 hAsymMatchedMBBlow->SetLabelSize(LabelSize);\r
987 hAsymMatchedMBBlow->SetTitle("Matched tracks charge asymmetry in Physics Selected events for MBB low (blue) and high (red) mult.");\r
988\r
989 TH1 *hDiffMatchedMBBhigh= static_cast<TH1*>(hPosMatchedMBBhigh->Clone("hDiffMatchedMBBhigh"));\r
990 hDiffMatchedMBBhigh->Add(hNegMatchedMBBhigh,-1);\r
991 hDiffMatchedMBBhigh->Sumw2();\r
992 \r
993 TH1 *hAsymMatchedMBBhigh= static_cast<TH1*>(hDiffMatchedMBBhigh->Clone("hAsymMatchedMBBhigh"));\r
994 hAsymMatchedMBBhigh->Divide(hAllMatchedMBBhigh);\r
995 hAsymMatchedMBBhigh->SetLineColor(kRed);\r
996 hAsymMatchedMBBhigh->SetLineWidth(2);\r
997 hAsymMatchedMBBhigh->SetMinimum(-0.2);\r
998 hAsymMatchedMBBhigh->SetMaximum(0.2);\r
999 hAsymMatchedMBBhigh->SetLabelSize(LabelSize);\r
1000 hAsymMatchedMBBhigh->SetTitle("Matched tracks asymetry in Physics Selected events");\r
1001 \r
1002\r
1003 CanvasName = LHCPeriod.Data() ; \r
1004 CanvasName += "_AsymMatched"; \r
1005 TCanvas *cAsymMatched = new TCanvas(CanvasName.Data(),"cAsymMatched",1200,900);\r
1006 cAsymMatched->SetRightMargin(0.01);\r
1007 cAsymMatched->SetGridy(1);\r
1008 cAsymMatched->cd();\r
1009 hAsymMatchedMBBlow->GetYaxis()->SetTitle("Charge asymmetry"); \r
1010 hAsymMatchedMBBlow->Draw("EH");\r
1011 hAsymMatchedMBBhigh->Draw("EHsame");\r
1012 \r
1013 TLegend* legcAMT = new TLegend(0.60,0.25,0.98,0.45);\r
1014 legcAMT->SetHeader("Charge asymmetry of matched track per MBB (include Vtx, #eta and R_{Abs} cuts) ");\r
1015 legcAMT->AddEntry(hAsymMatchedMBBlow," Low mult. events","l");\r
1016 legcAMT->AddEntry(hAsymMatchedMBBhigh," High mult. events ","l");\r
1017 legcAMT->Draw("same");\r
1018\r
1019 cAsymMatched->Print(OutFileNamePDF.Data());\r
1020 cAsymMatched->Write();\r
1021\r
1022\r
1023 //=========================================================\r
1024 // Draw low/high pt tracks in acceptance per MBB all mult \r
1025 TH1* hMatchedLowPtPerMBB = static_cast<TH1*> (hMatchedLowPtMBB->Clone("hMatchedLowPtPerMBB"));\r
1026 hMatchedLowPtPerMBB->Sumw2();\r
1027 hMatchedLowPtPerMBB->Divide(hMBB);\r
1028 hMatchedLowPtPerMBB->SetLineWidth(2);\r
1029 hMatchedLowPtPerMBB->SetLineColor(kBlue);\r
1030\r
1031 TH1* hMatchedHighPtPerMBB = static_cast<TH1*> (hMatchedHighPtMBB->Clone("hMatchedHighPtPerMBB"));\r
1032 hMatchedHighPtPerMBB->Sumw2();\r
1033 hMatchedHighPtPerMBB->Divide(hMBB);\r
1034 hMatchedHighPtPerMBB->SetLineWidth(2);\r
1035 hMatchedHighPtPerMBB->SetLineColor(kRed);\r
1036\r
1037 TH1* hMatchedLowPtPerMBBlow = static_cast<TH1*> (hMatchedLowPtMBBlow->Clone("hMatchedLowPtPerMBBlow"));\r
1038 hMatchedLowPtPerMBBlow->Sumw2();\r
1039 hMatchedLowPtPerMBBlow->Divide(hMBBlow);\r
1040 hMatchedLowPtPerMBBlow->SetLineWidth(2);\r
1041 hMatchedLowPtPerMBBlow->SetMarkerSize(2);\r
1042 hMatchedLowPtPerMBBlow->SetMarkerStyle(27);\r
1043 hMatchedLowPtPerMBBlow->SetMarkerColor(kBlue);\r
1044 hMatchedLowPtPerMBBlow->SetLineWidth(2);\r
1045 hMatchedLowPtPerMBBlow->SetLineColor(kBlue);\r
1046\r
1047 TH1* hMatchedHighPtPerMBBlow = static_cast<TH1*> (hMatchedHighPtMBBlow->Clone("hMatchedHighPtPerMBBlow"));\r
1048 hMatchedHighPtPerMBBlow->Sumw2();\r
1049 hMatchedHighPtPerMBBlow->Divide(hMBB);\r
1050 hMatchedHighPtPerMBBlow->SetLineWidth(2);\r
1051 hMatchedHighPtPerMBBlow->SetLineColor(kRed);\r
1052 hMatchedHighPtPerMBBlow->SetMarkerSize(2);\r
1053 hMatchedHighPtPerMBBlow->SetMarkerStyle(27);\r
1054 hMatchedHighPtPerMBBlow->SetMarkerColor(kRed);\r
1055\r
1056 TH1* hMatchedLowPtPerMBBhigh = static_cast<TH1*> (hMatchedLowPtMBBhigh->Clone("hMatchedLowPtPerMBBhigh"));\r
1057 hMatchedLowPtPerMBBhigh->Sumw2();\r
1058 hMatchedLowPtPerMBBhigh->Divide(hMBBhigh);\r
1059 hMatchedLowPtPerMBBhigh->SetLineWidth(2);\r
1060 hMatchedLowPtPerMBBhigh->SetLineColor(kBlue);\r
1061 hMatchedLowPtPerMBBhigh->SetMarkerSize(2);\r
1062 hMatchedLowPtPerMBBhigh->SetMarkerStyle(24);\r
1063 hMatchedLowPtPerMBBhigh->SetMarkerColor(kBlue);\r
1064\r
1065\r
1066 TH1* hMatchedHighPtPerMBBhigh = static_cast<TH1*> (hMatchedHighPtMBBhigh->Clone("hMatchedHighPtPerMBBhigh"));\r
1067 hMatchedHighPtPerMBBhigh->Sumw2();\r
1068 hMatchedHighPtPerMBBhigh->Divide(hMBB);\r
1069 hMatchedHighPtPerMBBhigh->SetLineWidth(2);\r
1070 hMatchedHighPtPerMBBhigh->SetLineColor(kRed);\r
1071 hMatchedHighPtPerMBBhigh->SetMarkerSize(2);\r
1072 hMatchedHighPtPerMBBhigh->SetMarkerStyle(24);\r
1073 hMatchedHighPtPerMBBhigh->SetMarkerColor(kRed);\r
1074\r
1075\r
1076 CanvasName = LHCPeriod.Data() ; \r
1077 CanvasName += "_HighPtMuons"; \r
1078 TCanvas *cHighPtMuons = new TCanvas(CanvasName.Data(),"cHighPtMuons",1200,900);\r
1079 cHighPtMuons->SetTopMargin(0.05);\r
1080 cHighPtMuons->SetRightMargin(0.01);\r
1081 //cHighPtMuons->SetLogy(1);\r
1082 cHighPtMuons->SetGridy(1);\r
1083 cHighPtMuons->cd();\r
1084\r
1085 hMatchedLowPtPerMBB->SetTitle("");\r
1086 hMatchedLowPtPerMBB->GetYaxis()->SetTitle("Ratio per MBB");\r
1087 hMatchedLowPtPerMBB->SetMaximum(0.13);\r
1088 hMatchedLowPtPerMBB->SetMinimum(0.02);\r
1089 hMatchedLowPtPerMBB->SetLabelSize(LabelSize);\r
1090\r
1091 hMatchedLowPtPerMBB->Draw("E");\r
1092 hMatchedHighPtPerMBB->Scale(5);\r
1093 hMatchedHighPtPerMBB->Draw("Esame");\r
1094 \r
1095 TLegend* legcHPM = new TLegend(0.60,0.50,0.98,0.65);\r
1096 legcHPM->SetHeader("Number of matched track per MBB (include Vtx, #eta and R_{Abs} cuts) ");\r
1097 legcHPM->AddEntry(hMatchedLowPtPerMBB," p_{T} > 1 GeV/c ","l");\r
1098 legcHPM->AddEntry(hMatchedHighPtPerMBB," (x5) p_{T} > 2 GeV/c ","l");\r
1099 legcHPM->Draw("same");\r
1100\r
1101 cHighPtMuons->Print(OutFileNamePDF.Data());\r
1102 cHighPtMuons->Write(); \r
1103\r
1104\r
1105 //=====================================================================\r
1106 // Draw low/high pt tracks in acceptance per MBB in low and high mult \r
1107 CanvasName = LHCPeriod.Data() ; \r
1108 CanvasName += "_HighPtMuonsVsMult"; \r
1109 TCanvas *cHighPtMuonsVsMult = new TCanvas(CanvasName.Data(),"cHighPtMuonsVsMult",1200,900);\r
1110 cHighPtMuonsVsMult->SetTopMargin(0.05);\r
1111 cHighPtMuonsVsMult->SetRightMargin(0.01);\r
1112 cHighPtMuonsVsMult->SetLogy(1);\r
1113 cHighPtMuonsVsMult->SetGridy(1);\r
1114 cHighPtMuonsVsMult->Divide(1,2);\r
1115\r
1116\r
1117 cHighPtMuonsVsMult->cd(1);\r
1118\r
1119 hMatchedLowPtPerMBBlow->SetTitle("");\r
1120 hMatchedLowPtPerMBBlow->GetYaxis()->SetTitle("Ratio per MBB in low mult");\r
1121 hMatchedLowPtPerMBBlow->SetMaximum(0.015);\r
1122 hMatchedLowPtPerMBBlow->SetMinimum(0.000001);\r
1123 hMatchedLowPtPerMBBlow->SetLabelSize(LabelSize);\r
1124\r
1125 hMatchedLowPtPerMBBlow->Draw("E");\r
1126 hMatchedHighPtPerMBBlow->Scale(10);\r
1127 hMatchedHighPtPerMBBlow->Draw("Esame");\r
1128 TLegend* legcHPMlow = new TLegend(0.60,0.35,0.98,0.5);\r
1129 legcHPMlow->SetHeader("Number of matched track per MBB (include Vtx, #eta and R_{Abs} cuts) ");\r
1130 legcHPMlow->AddEntry(hMatchedLowPtPerMBBlow," p_{T} > 1 GeV/c in low mult","p");\r
1131 legcHPMlow->AddEntry(hMatchedHighPtPerMBBlow,"(x10) p_{T} > 2 GeV/c in low mult","p");\r
1132 legcHPMlow->Draw("same");\r
1133\r
1134\r
1135 cHighPtMuonsVsMult->cd(2);\r
1136 hMatchedLowPtPerMBBhigh->SetTitle("");\r
1137 hMatchedLowPtPerMBBhigh->GetYaxis()->SetTitle("Ratio per MBB in high mult");\r
1138 hMatchedLowPtPerMBBhigh->SetMinimum(0.005);\r
1139 hMatchedLowPtPerMBBhigh->Draw("E");\r
1140 hMatchedHighPtPerMBBhigh->Scale(10);\r
1141 hMatchedHighPtPerMBBhigh->Draw("Esame");\r
1142 \r
1143 TLegend* legcHPMhigh = new TLegend(0.60,0.4,0.98,0.55);\r
1144 legcHPMhigh->SetHeader("Number of matched track per MBB (include Vtx, #eta and R_{Abs} cuts) ");\r
1145 legcHPMhigh->AddEntry(hMatchedLowPtPerMBBhigh," p_{T} > 1 GeV/c in high mult","p");\r
1146 legcHPMhigh->AddEntry(hMatchedHighPtPerMBBhigh,"(x10) p_{T} > 2 GeV/c in high mult","p");\r
1147 legcHPMhigh->Draw("same");\r
1148 \r
1149 cHighPtMuonsVsMult->Print(OutFileNamePDF.Data());\r
1150 cHighPtMuonsVsMult->Write(); \r
1151\r
1152 // close merged file\r
1153 globalFile->Close();\r
1154\r
1155\r
1156 //--------------------------------------------- //\r
1157 // monitor quantities run per run //\r
1158 //--------------------------------------------- //\r
1159 \r
1160\r
1161 // Are the runs stored locally or in alien?\r
1162 Int_t isAlienFile = 0;\r
1163 if(alienBaseDir.Contains("alien:")){\r
1164 isAlienFile = 1;\r
1165 alienBaseDir.ReplaceAll("alien://","");\r
1166 }\r
1167\r
1168 /* if ( ! alienBaseDir.Contains("alien:")){\r
1169 Info("PlotMuonQA","Working locally, stopping here");\r
1170 rootFileOut->Close();\r
1171 return;\r
1172 }\r
1173 */ \r
1174 // else (!TGrid::Connect("alien://")) {\r
1175 // Error("PlotMuonQA","cannot connect to grid. It's really needed here");\r
1176 // c1->Print(OutFileNamePDF_close.Data());\r
1177 // rootFileOut->Close();\r
1178 // return;\r
1179 // }\r
1180\r
1181 TH1F* hNClustersPerTrackVsRun_Mean = new TH1F("hNClustersPerTrackVsRun_Mean", "averaged number of associated clusters per track;run;<n_{clusters}>",10000,1,10000);\r
1182 TH1F* hNClustersPerTrackVsRun_Sigma = new TH1F("hNClustersPerTrackVsRun_Sigma", "dispersion of the number of associated clusters per track;run;#sigma_{n_{clusters}}",10000,1,10000);\r
1183 TH1F* hNChamberHitPerTrack_Mean = new TH1F("hNChamberHitPerTrack_Mean", "averaged number of chambers hit per track;run;<n_{chamber hit}>",10000,1,10000);\r
1184 TH1F* hNChamberHitPerTrack_Sigma = new TH1F("hNChamberHitPerTrack_Sigma", "dispersion of the number of chambers hit per track;run;#sigma_{n_{chamber hit}}",10000,1,10000);\r
1185 TH1F* hChi2_Mean = new TH1F("hChi2_Mean", "averaged normalized #chi^{2} distribution;run;<#chi^{2} / ndf>",10000,1,10000);\r
1186 TH1F* hChi2_Sigma = new TH1F("hChi2_Sigma", "dispersion of normalized #chi^{2} distribution;run;#sigma_{#chi^{2} / ndf}",10000,1,10000);\r
1187 TH1F* hNClustersInCh[10];\r
1188 for (Int_t ich=0; ich<10; ich++){\r
1189 hNClustersInCh[ich] = new TH1F(Form("hNClustersInCh%d",ich+1), Form("averaged number of clusters in chamber %d per track;run;<n_{clusters}>",ich+1),10000,1,10000);\r
1190 }\r
1191 TH1F* hClusterHitMapXInCh[10];\r
1192 for (Int_t ich=0; ich<10; ich++){\r
1193 hClusterHitMapXInCh[ich] = new TH1F(Form("hClusterHitMapXInCh%d",ich+1), Form("averaged cluster position distribution in chamber %d;X (cm)",ich+1),10000,1,10000);\r
1194 }\r
1195 TH1F* hClusterHitMapYInCh[10];\r
1196 for (Int_t ich=0; ich<10; ich++){\r
1197 hClusterHitMapYInCh[ich] = new TH1F(Form("hClusterHitMapYInCh%d",ich+1), Form("averaged cluster position distribution in chamber %d;Y (cm)",ich+1),10000,1,10000);\r
1198 }\r
1199\r
1200 Int_t ibin = 1;\r
1201 \r
1202 // loop over runs\r
1203 for ( Int_t irun=0; irun<runs.GetEntriesFast(); irun++ ) { //loop over runs \r
1204 \r
1205 TString run = ((TObjString*)runs.UncheckedAt(irun))->GetString();\r
1206 // get the file (or list of files) to be analyzed\r
1207 TString command;\r
1208 TGridResult *res = 0;\r
1209 TObjString *objs = 0;\r
1210 \r
1211 if(isAlienFile){\r
1212 // get the file (or list of files) to be analyzed\r
1213 command = Form("find %s/ %s/%s", alienBaseDir.Data(), run.Data(), QAFileName.Data());\r
1214 res = gGrid->Command(command);\r
1215 if (!res) {\r
1216 Error("PlotMuonQAPbPb",Form("no result for the command: %s",command.Data()));\r
1217 return;\r
1218 }\r
1219 }\r
1220 else{\r
1221 res = new TGridResult();\r
1222\r
1223 if(runList){\r
1224 objs = new TObjString(Form("%s/%s/AnalysisResults.root", alienBaseDir.Data(), run.Data()));\r
1225 res->Add(objs);\r
1226 }\r
1227 else {\r
1228 //loop over the directory to find the root files \r
1229 void *dir = gSystem->OpenDirectory(alienBaseDir.Data());\r
1230 TString sDirFilename;\r
1231 Int_t iEntry=0, iFile=0;\r
1232\r
1233 while(kTRUE){\r
1234 iEntry++;\r
1235 const char* dirFilename = gSystem->GetDirEntry(dir);\r
1236 if(!dirFilename) break;\r
1237 sDirFilename = dirFilename;\r
1238 if(!sDirFilename.IsDigit()) continue;\r
1239 iFile++;\r
1240 objs = new TObjString(Form("%s/%s/AnalysisResults.root", alienBaseDir.Data(), sDirFilename.Data()));\r
1241 res->Add(objs);\r
1242 }\r
1243 }\r
1244 }\r
1245\r
1246 // Loop over the 'find' results and get next LFN\r
1247 TIter nextmap(res);\r
1248 TMap *map = 0;\r
1249 TObjString *objs;\r
1250\r
1251 //some checks\r
1252 Int_t iLoop=0, iLoopMax=1000;\r
1253\r
1254 while (kTRUE){\r
1255 // get the current file url\r
1256 if(isAlienFile){\r
1257 map=(TMap*)nextmap();\r
1258 if(!map) break;\r
1259 objs = dynamic_cast<TObjString*>(map->GetValue("turl"));\r
1260 }\r
1261 else{\r
1262 objs=(TObjString*)nextmap();\r
1263 if(!objs) break;\r
1264 }\r
1265\r
1266 //in case of infinite loop\r
1267 iLoop++;\r
1268 if(iLoop>iLoopMax) break;\r
1269\r
1270\r
1271 if (!objs || !objs->GetString().Length()) {\r
1272 Error("PlotMuonQA","turl not found for the run %s... SKIPPING", run.Data());\r
1273 continue;\r
1274 }\r
1275 \r
1276 // open the outfile for this run\r
1277 TFile *runFile = TFile::Open(objs->GetString());\r
1278 if (!runFile || ! runFile->IsOpen()) {\r
1279 Error("PlotMuonQA", Form("failed to open file: %s", objs->GetName()));\r
1280 return;\r
1281 }\r
1282\r
1283 // change here\r
1284 runFile->Cd("MUON_QA");\r
1285 \r
1286 // get interesting histos\r
1287 TObjArray* general1 = static_cast<TObjArray*>(runFile->FindObjectAny("general1"));\r
1288 TObjArray* general2 = static_cast<TObjArray*>(runFile->FindObjectAny("general2"));\r
1289 TObjArray* expert = static_cast<TObjArray*>(runFile->FindObjectAny("expert"));\r
1290\r
1291 if (!general1 || !general2 || !expert){\r
1292 Error("PlotMUONQAPbPb", Form("All objects not here !!! ===> Skipping...for %s",objs->GetName()));\r
1293 continue;\r
1294 }\r
1295\r
1296 TH1* hNClustersPerTrack = static_cast<TH1*>(general1->FindObject("hNClustersPerTrack"));\r
1297 TH1* hNChamberHitPerTrack = static_cast<TH1*>(general1->FindObject("hNChamberHitPerTrack"));\r
1298 TH1* hChi2 = static_cast<TH1*>(general1->FindObject("hChi2"));\r
1299 TH1* hNClustersPerCh = static_cast<TH1*>(general2->FindObject("hNClustersPerCh"));\r
1300\r
1301 TH2* hClusterHitMapInCh[10];\r
1302 for(Int_t ich=0; ich<10; ich++) hClusterHitMapInCh[ich] = static_cast<TH2*>(expert->FindObject(Form("hClusterHitMapInCh%d",ich+1)));\r
1303\r
1304 // do not skip empty runs butrather set values to 0\r
1305 if (!hNClustersPerCh) {\r
1306 Warning("PlotQA", Form("File: %s has empty histograms !", objs->GetName()));\r
1307 hNClustersPerTrackVsRun_Mean->SetBinContent(ibin, 0.);\r
1308 hNClustersPerTrackVsRun_Mean->SetBinError(ibin, 1.);\r
1309 hNClustersPerTrackVsRun_Sigma->SetBinContent(ibin, 0.);\r
1310 hNClustersPerTrackVsRun_Sigma->SetBinError(ibin, 1.);\r
1311 hNChamberHitPerTrack_Mean->SetBinContent(ibin, 0.);\r
1312 hNChamberHitPerTrack_Mean->SetBinError(ibin, 1.);\r
1313 hNChamberHitPerTrack_Sigma->SetBinContent(ibin, 0.);\r
1314 hNChamberHitPerTrack_Sigma->SetBinError(ibin, 1.);\r
1315 hChi2_Mean->SetBinContent(ibin, 0.);\r
1316 hChi2_Mean->SetBinError(ibin, 1.);\r
1317 hChi2_Sigma->SetBinContent(ibin, 0.);\r
1318 hChi2_Sigma->SetBinError(ibin, 1.);\r
1319 for (Int_t ich=0; ich<10; ich++) {\r
1320 hNClustersInCh[ich]->SetBinContent(ibin,0.);\r
1321 hNClustersInCh[ich]->SetBinError(ibin,1.);\r
1322 hClusterHitMapXInCh[ich]->SetBinContent(ibin,0.);\r
1323 hClusterHitMapXInCh[ich]->SetBinError(ibin,1.);\r
1324 hClusterHitMapYInCh[ich]->SetBinContent(ibin,0.);\r
1325 hClusterHitMapYInCh[ich]->SetBinError(ibin,1.); \r
1326 }\r
1327 //runFile->Close();\r
1328 //continue;\r
1329 }\r
1330 else {\r
1331 // fill monitoring plots\r
1332 hNClustersPerTrackVsRun_Mean->SetBinContent(ibin, hNClustersPerTrack->GetMean());\r
1333 hNClustersPerTrackVsRun_Mean->SetBinError(ibin, hNClustersPerTrack->GetMeanError());\r
1334 hNClustersPerTrackVsRun_Sigma->SetBinContent(ibin, hNClustersPerTrack->GetRMS());\r
1335 hNClustersPerTrackVsRun_Sigma->SetBinError(ibin, hNClustersPerTrack->GetRMSError());\r
1336 hNChamberHitPerTrack_Mean->SetBinContent(ibin, hNChamberHitPerTrack->GetMean());\r
1337 hNChamberHitPerTrack_Mean->SetBinError(ibin, hNChamberHitPerTrack->GetMeanError());\r
1338 hNChamberHitPerTrack_Sigma->SetBinContent(ibin, hNChamberHitPerTrack->GetRMS());\r
1339 hNChamberHitPerTrack_Sigma->SetBinError(ibin, hNChamberHitPerTrack->GetRMSError());\r
1340 hChi2_Mean->SetBinContent(ibin, hChi2->GetMean());\r
1341 hChi2_Mean->SetBinError(ibin, hChi2->GetMeanError());\r
1342 hChi2_Sigma->SetBinContent(ibin, hChi2->GetRMS());\r
1343 hChi2_Sigma->SetBinError(ibin, hChi2->GetRMSError());\r
1344 for (Int_t ich=0; ich<10; ich++) {\r
1345 hNClustersInCh[ich]->SetBinContent(ibin,hNClustersPerCh->GetBinContent(ich+1));\r
1346 hNClustersInCh[ich]->SetBinError(ibin,hNClustersPerCh->GetBinError(ich+1));\r
1347 hClusterHitMapXInCh[ich]->SetBinContent(ibin,hClusterHitMapInCh[ich]->GetMean(1));\r
1348 hClusterHitMapXInCh[ich]->SetBinError(ibin,hClusterHitMapInCh[ich]->GetMeanError(1));\r
1349 hClusterHitMapYInCh[ich]->SetBinContent(ibin,hClusterHitMapInCh[ich]->GetMean(2));\r
1350 hClusterHitMapYInCh[ich]->SetBinError(ibin,hClusterHitMapInCh[ich]->GetMeanError(2));\r
1351 }\r
1352 }\r
1353\r
1354 // set labels\r
1355 run = objs->GetString();\r
1356 run.ReplaceAll(baseDir, "");\r
1357 run.Remove(TString::kLeading, '/');\r
1358 run.Remove(TString::kLeading, '0');\r
1359 run.ReplaceAll(QAFileName.Data(),"");\r
1360 run.Remove(TString::kTrailing, '/');\r
1361 hNClustersPerTrackVsRun_Mean->GetXaxis()->SetBinLabel(ibin, run.Data());\r
1362 hNClustersPerTrackVsRun_Sigma->GetXaxis()->SetBinLabel(ibin, run.Data());\r
1363 hNChamberHitPerTrack_Mean->GetXaxis()->SetBinLabel(ibin, run.Data());\r
1364 hNChamberHitPerTrack_Sigma->GetXaxis()->SetBinLabel(ibin, run.Data());\r
1365 hChi2_Mean->GetXaxis()->SetBinLabel(ibin, run.Data());\r
1366 hChi2_Sigma->GetXaxis()->SetBinLabel(ibin, run.Data());\r
1367 for (Int_t ich=0; ich<10; ich++) {\r
1368 hNClustersInCh[ich]->GetXaxis()->SetBinLabel(ibin, run.Data());\r
1369 hClusterHitMapXInCh[ich]->GetXaxis()->SetBinLabel(ibin, run.Data());\r
1370 hClusterHitMapYInCh[ich]->GetXaxis()->SetBinLabel(ibin, run.Data());\r
1371 }\r
1372 // close outfile for this run\r
1373 runFile->Close();\r
1374 ibin++; \r
1375 }\r
1376 \r
1377 delete res;\r
1378 }\r
1379\r
1380 //sort label\r
1381 hNClustersPerTrackVsRun_Mean->LabelsOption("a");\r
1382 hNClustersPerTrackVsRun_Sigma->LabelsOption("a");\r
1383 hNChamberHitPerTrack_Mean->LabelsOption("a");\r
1384 hNChamberHitPerTrack_Sigma->LabelsOption("a");\r
1385 hChi2_Mean->LabelsOption("a");\r
1386 hChi2_Sigma->LabelsOption("a");\r
1387\r
1388 for(Int_t ich=0; ich<10; ich++){\r
1389 hNClustersInCh[ich]->LabelsOption("a");\r
1390 hClusterHitMapXInCh[ich]->LabelsOption("a");\r
1391 hClusterHitMapYInCh[ich]->LabelsOption("a");\r
1392 }\r
1393 \r
1394 TString dirToGo = OutFileNameROOT.Data(); dirToGo+=":/";\r
1395 gDirectory->Cd(dirToGo.Data());\r
1396 //==================================================\r
1397 //Display Mean and Sigma of the number of associated clusters to a track \r
1398 TLegend *lNClusters = new TLegend(0.75,0.85,0.99,0.99);\r
1399 lNClusters->AddEntry(hNClustersPerTrackVsRun_Mean,"clusters","PL");\r
1400 lNClusters->AddEntry(hNChamberHitPerTrack_Mean,"chamber hit","PL");\r
1401 \r
1402 TCanvas* cNClusters = new TCanvas("cNClusters","cNClusters",1200,900);\r
1403 cNClusters->SetLogy(0);\r
1404 cNClusters->Divide(1,2);\r
1405 cNClusters->cd(1);\r
1406 //hNClustersPerTrackVsRun_Mean->SetMaximum(11);\r
1407 hNClustersPerTrackVsRun_Mean->SetMinimum(7);\r
1408 hNClustersPerTrackVsRun_Mean->SetStats(kFALSE);\r
1409 hNClustersPerTrackVsRun_Mean->GetXaxis()->SetRange(1,ibin-1);\r
1410 hNClustersPerTrackVsRun_Mean->GetXaxis()->SetNdivisions(1,kFALSE);\r
1411 //hNClustersPerTrackVsRun_Mean->LabelsOption("u");\r
1412 hNClustersPerTrackVsRun_Mean->SetLabelSize(LabelSize);\r
1413 hNClustersPerTrackVsRun_Mean->SetTitle("averaged number of associated clusters or of the number of chamber hit per track");\r
1414 hNClustersPerTrackVsRun_Mean->SetLineWidth(2);\r
1415 hNClustersPerTrackVsRun_Mean->Draw("e");\r
1416 hNChamberHitPerTrack_Mean->SetLineColor(kRed);\r
1417 hNChamberHitPerTrack_Mean->SetLineWidth(2);\r
1418 hNChamberHitPerTrack_Mean->Draw("esame");\r
1419 lNClusters->Draw("same");\r
1420\r
1421 cNClusters->cd(2);\r
1422 //hNClustersPerTrackVsRun_Sigma->SetMaximum(1.1);\r
1423 hNClustersPerTrackVsRun_Sigma->SetMinimum(0.4);\r
1424 hNClustersPerTrackVsRun_Sigma->SetStats(kFALSE);\r
1425 hNClustersPerTrackVsRun_Sigma->GetXaxis()->SetRange(1,ibin-1);\r
1426 hNClustersPerTrackVsRun_Sigma->GetXaxis()->SetNdivisions(1,kFALSE);\r
1427 //hNClustersPerTrackVsRun_Sigma->LabelsOption("u");\r
1428 hNClustersPerTrackVsRun_Sigma->SetLabelSize(LabelSize);\r
1429 hNClustersPerTrackVsRun_Sigma->SetTitle("dispersion of the number of associated clusters or of the number of chamber hit per track");\r
1430 hNClustersPerTrackVsRun_Sigma->SetLineWidth(2);\r
1431 hNClustersPerTrackVsRun_Sigma->Draw("e");\r
1432 hNChamberHitPerTrack_Sigma->SetLineWidth(2);\r
1433 hNChamberHitPerTrack_Sigma->SetLineColor(kRed);\r
1434 hNChamberHitPerTrack_Sigma->Draw("esame");\r
1435 lNClusters->Draw("same");\r
1436\r
1437 cNClusters->Print(OutFileNamePDF.Data());\r
1438 cNClusters->Write();\r
1439\r
1440\r
1441 //==================================================\r
1442 // Display average number of cluster per chamber\r
1443 TLegend *lNClustersPerCh = new TLegend(0.92,0.45,0.99,0.99);\r
1444 TCanvas* cNClustersPerCh = new TCanvas("cNClustersPerCh","cNClustersPerCh",1200,900);\r
1445 cNClustersPerCh->cd();\r
1446 cNClustersPerCh->SetRightMargin(0.1);\r
1447 hNClustersInCh[0]->SetStats(kFALSE);\r
1448 hNClustersInCh[0]->GetXaxis()->SetRange(1,ibin-1);\r
1449 hNClustersInCh[0]->GetXaxis()->SetNdivisions(1,kFALSE);\r
1450 //hNClustersInCh[0]->LabelsOption("u");\r
1451 hNClustersInCh[0]->SetLabelSize(LabelSize);\r
1452 hNClustersInCh[0]->SetTitle("averaged number of clusters in chamber i per track");\r
1453 hNClustersInCh[0]->SetMaximum(1.2);\r
1454 hNClustersInCh[0]->SetMinimum(0.01);\r
1455 for (Int_t ich=0; ich<10; ich++) {\r
1456 hNClustersInCh[ich]->SetLineColor(ich+1+ich/9);\r
1457 hNClustersInCh[ich]->SetLineWidth(2);\r
1458 if (ich == 0) hNClustersInCh[ich]->Draw("e");\r
1459 else hNClustersInCh[ich]->Draw("esame");\r
1460 lNClustersPerCh->AddEntry(hNClustersInCh[ich],Form("ch%d",ich+1),"PL");\r
1461 }\r
1462 lNClustersPerCh->Draw("same");\r
1463 \r
1464 cNClustersPerCh->Print(OutFileNamePDF.Data());\r
1465 cNClustersPerCh->Write();\r
1466\r
1467 //==================================================\r
1468 // Display average X and Y position of clusters per chamber\r
1469 TLegend *lClusterHitMapPerCh = new TLegend(0.92,0.45,0.99,0.99);\r
1470 TCanvas* cClusterHitMapPerCh = new TCanvas("cClusterHitMapPerCh","cClusterHitMapPerCh",1200,900);\r
1471 cClusterHitMapPerCh->Divide(1,2);\r
1472 cClusterHitMapPerCh->GetPad(1)->SetRightMargin(0.1);\r
1473 cClusterHitMapPerCh->GetPad(2)->SetRightMargin(0.1);\r
1474 \r
1475 cClusterHitMapPerCh->cd(1);\r
1476 hClusterHitMapXInCh[0]->SetStats(kFALSE);\r
1477 hClusterHitMapXInCh[0]->GetXaxis()->SetRange(1,ibin-1);\r
1478 hClusterHitMapXInCh[0]->GetXaxis()->SetNdivisions(1,kFALSE);\r
1479 //hNClustersInCh[0]->LabelsOption("u");\r
1480 hClusterHitMapXInCh[0]->SetLabelSize(LabelSize);\r
1481 hClusterHitMapXInCh[0]->SetTitle("<X> of clusters - associated to a track - in chamber i");\r
1482 hClusterHitMapXInCh[0]->SetMaximum(30);\r
1483 hClusterHitMapXInCh[0]->SetMinimum(-30);\r
1484 for (Int_t ich=0; ich<10; ich++) {\r
1485 hClusterHitMapXInCh[ich]->SetLineColor(ich+1+ich/9);\r
1486 hClusterHitMapXInCh[ich]->SetLineWidth(2);\r
1487 if (ich == 0) hClusterHitMapXInCh[ich]->Draw("e");\r
1488 else hClusterHitMapXInCh[ich]->Draw("esame");\r
1489 \r
1490 lClusterHitMapPerCh->AddEntry(hClusterHitMapXInCh[ich],Form("ch%d",ich+1),"PL");\r
1491 }\r
1492 lClusterHitMapPerCh->Draw("same");\r
1493 \r
1494 cClusterHitMapPerCh->cd(2);\r
1495 hClusterHitMapYInCh[0]->SetStats(kFALSE);\r
1496 hClusterHitMapYInCh[0]->GetXaxis()->SetRange(1,ibin-1);\r
1497 hClusterHitMapYInCh[0]->GetXaxis()->SetNdivisions(1,kFALSE);\r
1498 //hNClustersInCh[0]->LabelsOption("u");\r
1499 hClusterHitMapYInCh[0]->SetLabelSize(LabelSize);\r
1500 hClusterHitMapYInCh[0]->SetTitle("<Y> of clusters - associated to a track - in chamber i");\r
1501 hClusterHitMapYInCh[0]->SetMaximum(30);\r
1502 hClusterHitMapYInCh[0]->SetMinimum(-30);\r
1503 for (Int_t ich=0; ich<10; ich++) {\r
1504 hClusterHitMapYInCh[ich]->SetLineColor(ich+1+ich/9);\r
1505 hClusterHitMapYInCh[ich]->SetLineWidth(2);\r
1506 if (ich == 0) hClusterHitMapYInCh[ich]->Draw("e");\r
1507 else hClusterHitMapYInCh[ich]->Draw("esame");\r
1508 }\r
1509 lClusterHitMapPerCh->Draw("same");\r
1510 \r
1511 cClusterHitMapPerCh->Print(OutFileNamePDF.Data());\r
1512 cClusterHitMapPerCh->Write();\r
1513\r
1514\r
1515 //==================================================\r
1516 // Display tracks ChiSquare \r
1517 TCanvas* cChi2 = new TCanvas("cChi2","cChi2",1200,900);\r
1518 cChi2->Divide(1,2);\r
1519 cChi2->cd(1);\r
1520 hChi2_Mean->SetStats(kFALSE);\r
1521 hChi2_Mean->GetXaxis()->SetRange(1,ibin-1);\r
1522 hChi2_Mean->GetXaxis()->SetNdivisions(1,kFALSE);\r
1523 //hChi2_Mean->LabelsOption("u");\r
1524 hChi2_Mean->SetLabelSize(LabelSize);\r
1525 hChi2_Mean->SetLineWidth(2);\r
1526 hChi2_Mean->Draw("e");\r
1527\r
1528 cChi2->cd(2);\r
1529 hChi2_Sigma->SetStats(kFALSE);\r
1530 hChi2_Sigma->GetXaxis()->SetRange(1,ibin-1);\r
1531 hChi2_Sigma->GetXaxis()->SetNdivisions(1,kFALSE);\r
1532 //hChi2_Sigma->LabelsOption("u");\r
1533 hChi2_Sigma->SetLabelSize(LabelSize);\r
1534 hChi2_Sigma->SetLineWidth(2);\r
1535 hChi2_Sigma->Draw("e");\r
1536\r
1537 cChi2->Print(OutFileNamePDF.Data());\r
1538 cChi2->Write();\r
1539\r
1540 // close the PDF file\r
1541 c1->Print(OutFileNamePDF_close.Data());\r
1542 rootFileOut->Close();\r
1543 \r
1544 //====================================================\r
1545 if (PRINTSTAT){\r
1546 // set the format to print labels\r
1547 THashList* labels = hMBB->GetXaxis()->GetLabels();\r
1548 TString format(Form("\n%%%ds %%9d",0));\r
1549 \r
1550 // print value for each label\r
1551 TObjString* label = 0x0;\r
1552 TIter nextLabel(labels);\r
1553 cout << "-------------------------------------------------" << endl;\r
1554 cout << "Run Number" << "\t Number of CMBB after Phys. Sel. " << endl ; \r
1555 while ((label = static_cast<TObjString*>(nextLabel()))) {\r
1556 Int_t bin = (Int_t) label->GetUniqueID();\r
1557 printf(format.Data(), label->String().Data(), (Int_t) hMBB->GetBinContent(bin));\r
1558 }\r
1559 printf("\n\n");\r
1560\r
1561\r
1562 cout << "-------------------------------------------------" << endl;\r
1563 cout << "Total statistic" << endl; \r
1564 cout << " " << endl ; \r
1565 cout << "Number of NumOfAllTriggers " << endl ;\r
1566 cout << "\t before selection " << NumOfAllTriggersNoPS << "\t after selection " << NumOfAllTriggersWithPS << " --> rejection = " << (Double_t) (NumOfAllTriggersNoPS-NumOfAllTriggersWithPS)/(NumOfAllTriggersNoPS)*100. << "%" << endl ; \r
1567 cout << "Number of NumOfCMBB " << endl ;\r
1568 cout << "\t before selection " << NumOfMBBNoPS << "\t after selection " << NumOfMBBWithPS << " --> rejection = " << (Double_t) (NumOfMBBNoPS-NumOfMBBWithPS)/(NumOfMBBNoPS)*100. << "%" << endl ; \r
1569 cout << " " << endl ; \r
1570 \r
1571 }\r
1572}\r