]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG4/JetTasks/AliHistogramsUE.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
new UE Task from Sara
[u/mrichter/AliRoot.git] / PWG4 / JetTasks / AliHistogramsUE.cxx
CommitLineData
1f329128 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: A.Abrahantes, E.Lopez, S.Vallero *
5 * Version 1.0 *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16/* $Id:$ */
17////////////////////////////////////////////////
18//---------------------------------------------
19// Class to handle histograms for UE analysis
20//---------------------------------------------
21////////////////////////////////////////////////
22#include <TROOT.h>
23#include <TBranch.h>
24#include <TCanvas.h>
25#include <TChain.h>
26#include <TFile.h>
27#include <TH1F.h>
28#include <TH1I.h>
29#include <TH2F.h>
30#include <TLegend.h>
31#include <TList.h>
32#include <TLorentzVector.h>
33#include <TMath.h>
34#include <TProfile.h>
35#include <TRandom.h>
36#include <TStyle.h>
37#include <TSystem.h>
38#include <TTree.h>
39#include <TVector3.h>
40
41#include "AliHistogramsUE.h"
42#include "AliAnalyseUE.h"
43#include "AliLog.h"
44
45ClassImp( AliHistogramsUE)
46
47//____________________________________________________________________
48AliHistogramsUE:: AliHistogramsUE():
49TObject(),
50fBinsPtInHist(0),
51fMinJetPtInHist(0.),
52fMaxJetPtInHist(0.),
53fTrackEtaCut(0.),
54fListOfHistos(0x0),
55fhNJets(0x0),
56fhEleadingPt(0x0),
57fhMinRegPtDist(0x0),
58fhRegionMultMin(0x0),
59fhMinRegAvePt(0x0),
60fhMinRegSumPt(0x0),
61fhMinRegMaxPtPart(0x0),
62fhMinRegSumPtvsMult(0x0),
63fhdNdEtaPhiDist(0x0),
64fhFullRegPartPtDistVsEt(0x0),
65fhTransRegPartPtDistVsEt(0x0),
66fhRegionSumPtMaxVsEt(0x0),
67fhRegionMultMax(0x0),
68fhRegionMultMaxVsEt(0x0),
69fhRegionSumPtMinVsEt(0x0), //fhRegionMultMin(0x0),
70fhRegionMultMinVsEt(0x0),
71fhRegionAveSumPtVsEt(0x0),
72fhRegionDiffSumPtVsEt(0x0),
73fhRegionAvePartPtMaxVsEt(0x0),
74fhRegionAvePartPtMinVsEt(0x0),
75fhRegionMaxPartPtMaxVsEt(0x0),
76fhRegForwardMult(0x0),
77fhRegForwardSumPtvsMult(0x0),
78fhRegBackwardMult(0x0),
79fhRegBackwardSumPtvsMult(0x0),
80fhRegForwardPartPtDistVsEt(0x0),
81fhRegBackwardPartPtDistVsEt(0x0),
82fhRegTransMult(0x0),
83fhRegTransSumPtVsMult(0x0),
84fhMinRegSumPtJetPtBin(0x0),
85fhMaxRegSumPtJetPtBin(0x0),
86fhVertexMult(0x0),
87fh1Xsec(0x0),
88fh1Trials(0x0)
89{
90 // Default constructor
91
92}
93
94//____________________________________________________________________
95AliHistogramsUE:: AliHistogramsUE(TList *list):
96TObject(),
97fBinsPtInHist(0),
98fMinJetPtInHist(0.),
99fMaxJetPtInHist(0.),
100fTrackEtaCut(0.),
101fListOfHistos(0x0),
102fhNJets(0x0),
103fhEleadingPt(0x0),
104fhMinRegPtDist(0x0),
105fhRegionMultMin(0x0),
106fhMinRegAvePt(0x0),
107fhMinRegSumPt(0x0),
108fhMinRegMaxPtPart(0x0),
109fhMinRegSumPtvsMult(0x0),
110fhdNdEtaPhiDist(0x0),
111fhFullRegPartPtDistVsEt(0x0),
112fhTransRegPartPtDistVsEt(0x0),
113fhRegionSumPtMaxVsEt(0x0),
114fhRegionMultMax(0x0),
115fhRegionMultMaxVsEt(0x0),
116fhRegionSumPtMinVsEt(0x0), //fhRegionMultMin(0x0),
117fhRegionMultMinVsEt(0x0),
118fhRegionAveSumPtVsEt(0x0),
119fhRegionDiffSumPtVsEt(0x0),
120fhRegionAvePartPtMaxVsEt(0x0),
121fhRegionAvePartPtMinVsEt(0x0),
122fhRegionMaxPartPtMaxVsEt(0x0),
123fhRegForwardMult(0x0),
124fhRegForwardSumPtvsMult(0x0),
125fhRegBackwardMult(0x0),
126fhRegBackwardSumPtvsMult(0x0),
127fhRegForwardPartPtDistVsEt(0x0),
128fhRegBackwardPartPtDistVsEt(0x0),
129fhRegTransMult(0x0),
130fhRegTransSumPtVsMult(0x0),
131fhMinRegSumPtJetPtBin(0x0),
132fhMaxRegSumPtJetPtBin(0x0),
133fhVertexMult(0x0),
134fh1Xsec(0x0),
135fh1Trials(0x0)
136{
137 // Constructor, initialize members from list
138 fhNJets = (TH1F*)list->FindObject("hNJets");
139 fhEleadingPt = (TH1F*)list->FindObject("hEleadingPt");
140 fhMinRegPtDist = (TH1F*)list->FindObject("hMinRegPtDist");
141 fhRegionMultMin = (TH1F*)list->FindObject("hRegionMultMin");
142 fhMinRegAvePt = (TH1F*)list->FindObject("hMinRegAvePt");
143 fhMinRegSumPt = (TH1F*)list->FindObject("hMinRegSumPt");
144 fhMinRegMaxPtPart = (TH1F*)list->FindObject("hMinRegMaxPtPart");
145 fhMinRegSumPtvsMult = (TH1F*)list->FindObject("hMinRegSumPtvsMult");
146 fhdNdEtaPhiDist = (TH2F*)list->FindObject("hdNdEtaPhiDist");
147 fhFullRegPartPtDistVsEt = (TH2F*)list->FindObject("hFullRegPartPtDistVsEt");
148 fhTransRegPartPtDistVsEt = (TH2F*)list->FindObject("hTransRegPartPtDistVsEt");
149 fhRegionSumPtMaxVsEt = (TH1F*)list->FindObject("hRegionSumPtMaxVsEt");
150 fhRegionMultMax = (TH1I*)list->FindObject("hRegionMultMax");
151 fhRegionMultMaxVsEt = (TH1F*)list->FindObject("hRegionMultMaxVsEt");
152 fhRegionSumPtMinVsEt = (TH1F*)list->FindObject("hRegionSumPtMinVsEt");
153 fhRegionMultMinVsEt = (TH1F*)list->FindObject("hRegionMultMinVsEt");
154 fhRegionAveSumPtVsEt = (TH1F*)list->FindObject("hRegionAveSumPtVsEt");
155 fhRegionDiffSumPtVsEt = (TH1F*)list->FindObject("hRegionDiffSumPtVsEt");
156 fhRegionAvePartPtMaxVsEt = (TH1F*)list->FindObject("hRegionAvePartPtMaxVsEt");
157 fhRegionAvePartPtMinVsEt = (TH1F*)list->FindObject("hRegionAvePartPtMinVsEt");
158 fhRegionMaxPartPtMaxVsEt = (TH1F*)list->FindObject("hRegionMaxPartPtMaxVsEt");
159 fhRegForwardMult = (TH2F*)list->FindObject("hRegForwardMult");
160 fhRegForwardSumPtvsMult = (TH2F*)list->FindObject("hRegForwardSumPtvsMult");
161 fhRegBackwardMult = (TH2F*)list->FindObject("hRegBackwardMult");
162 fhRegBackwardSumPtvsMult = (TH2F*)list->FindObject("hRegBackwardSumPtvsMult");
163 fhRegForwardPartPtDistVsEt = (TH2F*)list->FindObject("hRegForwardPartPtDistVsEt");
164 fhRegBackwardPartPtDistVsEt = (TH2F*)list->FindObject("hRegBackwardPartPtDistVsEt");
165 fhRegTransMult = (TH2F*)list->FindObject("hRegTransMult");
166 fhRegTransSumPtVsMult = (TH2F*)list->FindObject("hRegTransSumPtVsMult");
167 fhMinRegSumPtJetPtBin = (TH2F*)list->FindObject("hMinRegSumPtJetPtBin");
168 fhMaxRegSumPtJetPtBin = (TH2F*)list->FindObject("hMaxRegSumPtJetPtBin");
169 fhVertexMult = (TH1F*)list->FindObject("hVertexMult");
170 fh1Xsec = (TProfile*)list->FindObject("h1Xsec");
171 fh1Trials = (TH1F*)list->FindObject("h1Trials");
172
173}
174//____________________________________________________________________
175AliHistogramsUE:: AliHistogramsUE(const AliHistogramsUE & original):
176TObject(),
177fBinsPtInHist(original.fBinsPtInHist),
178fMinJetPtInHist(original.fMinJetPtInHist),
179fMaxJetPtInHist(original.fMaxJetPtInHist),
180fTrackEtaCut(original.fTrackEtaCut),
181fListOfHistos(original.fListOfHistos),
182fhNJets(original.fhNJets),
183fhEleadingPt(original.fhEleadingPt),
184fhMinRegPtDist(original.fhMinRegPtDist),
185fhRegionMultMin(original.fhRegionMultMin),
186fhMinRegAvePt(original.fhMinRegAvePt),
187fhMinRegSumPt(original.fhMinRegSumPt),
188fhMinRegMaxPtPart(original.fhMinRegMaxPtPart),
189fhMinRegSumPtvsMult(original.fhMinRegSumPtvsMult),
190fhdNdEtaPhiDist(original.fhdNdEtaPhiDist),
191fhFullRegPartPtDistVsEt(original.fhFullRegPartPtDistVsEt),
192fhTransRegPartPtDistVsEt(original.fhTransRegPartPtDistVsEt),
193fhRegionSumPtMaxVsEt(original.fhRegionSumPtMaxVsEt),
194fhRegionMultMax(original.fhRegionMultMax),
195fhRegionMultMaxVsEt(original.fhRegionMultMaxVsEt),
196fhRegionSumPtMinVsEt(original.fhRegionSumPtMinVsEt),
197fhRegionMultMinVsEt(original.fhRegionMultMinVsEt),
198fhRegionAveSumPtVsEt(original.fhRegionAveSumPtVsEt),
199fhRegionDiffSumPtVsEt(original.fhRegionDiffSumPtVsEt),
200fhRegionAvePartPtMaxVsEt(original.fhRegionAvePartPtMaxVsEt),
201fhRegionAvePartPtMinVsEt(original.fhRegionAvePartPtMinVsEt),
202fhRegionMaxPartPtMaxVsEt(original.fhRegionMaxPartPtMaxVsEt),
203fhRegForwardMult(original.fhRegForwardMult),
204fhRegForwardSumPtvsMult(original.fhRegForwardSumPtvsMult),
205fhRegBackwardMult(original.fhRegBackwardMult),
206fhRegBackwardSumPtvsMult(original.fhRegBackwardSumPtvsMult),
207fhRegForwardPartPtDistVsEt(original.fhRegForwardPartPtDistVsEt),
208fhRegBackwardPartPtDistVsEt(original.fhRegBackwardPartPtDistVsEt),
209fhRegTransMult(original.fhRegTransMult),
210fhRegTransSumPtVsMult(original.fhRegTransSumPtVsMult),
211fhMinRegSumPtJetPtBin(original.fhMinRegSumPtJetPtBin),
212fhMaxRegSumPtJetPtBin(original.fhMaxRegSumPtJetPtBin),
213fhVertexMult(original.fhVertexMult),
214fh1Xsec(original.fh1Xsec),
215fh1Trials(original.fh1Trials)
216{
217
218 // Copy constructor
219
220}
221
222
223//______________________________________________________________
224AliHistogramsUE & AliHistogramsUE::operator = (const AliHistogramsUE & /*source*/)
225{
226
227 // assignment operator
228 return *this;
229
230}
231
232//______________________________________________________________
233TObjArray* AliHistogramsUE::CreateCanvas(const Int_t ncanv){
234
235 // Create canvas for plotting
236 printf("Creating %d canvas ... \n",ncanv);
237 TObjArray *arr=new TObjArray;
238 TString name;
239 for(Int_t i=0;i<ncanv;i++ ){
240 name=Form("Canvas %d",i);
241 TCanvas *c=new TCanvas(name,name);
242 c->SetFillColor(0);
243 gStyle->SetOptStat(0);
244 gStyle->SetOptTitle(0);
245 arr->Add(c);
246 }
247return arr;
248}
249
250//____________________________________________________________________
251TList* AliHistogramsUE::CreateHistos(Int_t bins, Double_t min, Double_t max, Double_t etacut)
252{
253
254 // Create all histograms necessary for UE analysis
255 fBinsPtInHist = bins;
256 fMinJetPtInHist = min;
257 fMaxJetPtInHist = max;
258 fTrackEtaCut= etacut;
259
260 fListOfHistos = new TList();
261
262 //Number of reconstructed clusters
263 fhNJets = new TH1F("hNJets", "Number of clusters", 20, 0, 20);
264 fhNJets->SetXTitle("Number of reconstructed clusters");
265 fhNJets->SetYTitle("#");
266 fhNJets->Sumw2();
267 fListOfHistos->Add( fhNJets ); // At(0)
268
269 //pT distribution of leading clusters
270 fhEleadingPt = new TH1F("hEleadingPt", "Leading cluster p_{T}", bins, min, max);
271 fhEleadingPt->SetXTitle("p_{T} of cluster (GeV/c)");
272 fhEleadingPt->SetYTitle("1/N_{ev} dN/dp_{T} (|#eta|<0.5)");
273 fhEleadingPt->Sumw2();
274 fListOfHistos->Add( fhEleadingPt ); // At(1)
275
276 //Track pT distribution in MIN zone
277 fhMinRegPtDist = new TH1F("hMinRegPtDist", "p_{T} distribution in MIN zone", 50,0.,20.);
278 fhMinRegPtDist->SetXTitle("Track p_{T} (GeV/c)");
279 fhMinRegPtDist->SetYTitle("dN/dp_{T}");
280 fhMinRegPtDist->Sumw2();
281 fListOfHistos->Add( fhMinRegPtDist ); // At(2)
282
283 //Multiplicity in MIN zone
284 fhRegionMultMin = new TH1F("hRegionMultMin", "N_{ch}^{90, min}", 21, -0.5, 20.5);
285 fhRegionMultMin->SetXTitle("N_{ch tracks}");
286 fhRegionMultMin->Sumw2();
287 fListOfHistos->Add( fhRegionMultMin ); // At(3)
288
289 //Average pT in MIN region
290 fhMinRegAvePt = new TH1F("hMinRegAvePt", "#LTp_{T}#GT", 50, 0., 20.);
291 fhMinRegAvePt->SetXTitle("p_{T} (GeV/c)");
292 fhMinRegAvePt->Sumw2();
293 fListOfHistos->Add( fhMinRegAvePt ); // At(4)
294
295 //Sum pT in MIN region
296 fhMinRegSumPt = new TH1F("hMinRegSumPt", "#Sigma p_{T} ", 50, 0., 20.);
297 fhMinRegSumPt->SetYTitle("Ed^{3}N_{tracks}/dp^{3} (c^{3}/GeV^{2})");
298 fhMinRegSumPt->SetXTitle("#Sigma p_{T} (GeV/c)");
299 fhMinRegSumPt->Sumw2();
300 fListOfHistos->Add( fhMinRegSumPt ); // At(5)
301
302 //Track with maximum pT in MIN region
303 fhMinRegMaxPtPart = new TH1F("hMinRegMaxPtPart", "max(p_{T})|_{event} ", 50, 0., 20.);
304 fhMinRegMaxPtPart->SetYTitle("Ed^{3}N_{tracks}/dp^{3} (c^{3}/GeV^{2})");
305 fhMinRegMaxPtPart->SetXTitle("p_{T} (GeV/c)");
306 fhMinRegMaxPtPart->Sumw2();
307 fListOfHistos->Add( fhMinRegMaxPtPart ); // At(6)
308
309 //Sum pT vs. multiplicity in MIN region
310 fhMinRegSumPtvsMult = new TH1F("hMinRegSumPtvsMult", "#Sigma p_{T} vs. Multiplicity ", 21, -0.5, 20.5);
311 fhMinRegSumPtvsMult->SetYTitle("#Sigma p_{T} (GeV/c)");
312 fhMinRegSumPtvsMult->SetXTitle("N_{charge}");
313 fhMinRegSumPtvsMult->Sumw2();
314 fListOfHistos->Add( fhMinRegSumPtvsMult ); // At(7);
315
316 //Phi correlation track-cluster vs. leading cluster pT
317 fhdNdEtaPhiDist = new TH2F("hdNdEtaPhiDist", Form("Charge particle density |#eta|<%3.1f vs #Delta#phi", fTrackEtaCut),62, 0., 2.*TMath::Pi(), bins, min, max);
318 fhdNdEtaPhiDist->SetXTitle("#Delta#phi");
319 fhdNdEtaPhiDist->SetYTitle("Leading cluster p_{T}");
320 fhdNdEtaPhiDist->Sumw2();
321 fListOfHistos->Add( fhdNdEtaPhiDist ); // At(8)
322
323 //Can be used to get track pT distribution for different cluster pT bins (full region)
324 fhFullRegPartPtDistVsEt = new TH2F("hFullRegPartPtDistVsEt", Form( "dN/dp_{T} |#eta|<%3.1f vs Leading cluster p_{T}", fTrackEtaCut),100,0.,50., bins, min, max);
325 fhFullRegPartPtDistVsEt->SetYTitle("Leading cluster p_{T}");
326 fhFullRegPartPtDistVsEt->SetXTitle("p_{T}");
327 fhFullRegPartPtDistVsEt->Sumw2();
328 fListOfHistos->Add( fhFullRegPartPtDistVsEt ); // At(9)
329
330 //Can be used to get part pT distribution for different cluster pT bins (transverse region)
331 fhTransRegPartPtDistVsEt = new TH2F("hTransRegPartPtDistVsEt", Form( "dN/dp_{T} in tranvese regions |#eta|<%3.1f vs Leading cluster p_{T}", fTrackEtaCut),100,0.,50., bins, min, max);
332 fhTransRegPartPtDistVsEt->SetYTitle("Leading cluster p_{T}");
333 fhTransRegPartPtDistVsEt->SetXTitle("p_{T}");
334 fhTransRegPartPtDistVsEt->Sumw2();
335 fListOfHistos->Add( fhTransRegPartPtDistVsEt ); // At(10)
336
337 //Sum pT in MAX region vs. leading-cluster pT
338 fhRegionSumPtMaxVsEt = new TH1F("hRegionSumPtMaxVsEt", "P_{T}^{90, max} vs Leading cluster p_{T}", bins, min, max);
339 fhRegionSumPtMaxVsEt->SetXTitle("p_{T} (GeV/c)");
340 fhRegionSumPtMaxVsEt->Sumw2();
341 fListOfHistos->Add( fhRegionSumPtMaxVsEt ); // At(11)
342
343
344 //Sum pT in MIN region vs. leading-cluster pT
345 fhRegionSumPtMinVsEt = new TH1F("hRegionSumPtMinVsEt", "P_{T}^{90, min} vs Leading cluster p_{T}", bins, min, max);
346 fhRegionSumPtMinVsEt->SetXTitle("p_{T} (GeV/c)");
347 fhRegionSumPtMinVsEt->Sumw2();
348 fListOfHistos->Add( fhRegionSumPtMinVsEt ); // At(12)
349
350 //Multiplicity in MAX region
351 fhRegionMultMax = new TH1I("hRegionMultMax", "N_{ch}^{90, max}", 21, -0.5, 20.5);
352 fhRegionMultMax->SetXTitle("N_{ch tracks}");
353 fhRegionMultMax->Sumw2();
354 fListOfHistos->Add( fhRegionMultMax ); // At(13)
355
356 //Multiplicity in MAX region vs. leading-cluster pT
357 fhRegionMultMaxVsEt = new TH1F("hRegionMultMaxVsEt", "N_{ch}^{90, max} vs Leading cluster p_{T}", bins, min, max);
358 fhRegionMultMaxVsEt->SetXTitle("p_{T} (GeV/c)");
359 fhRegionMultMaxVsEt->Sumw2();
360 fListOfHistos->Add( fhRegionMultMaxVsEt ); // At(14)
361
362 //Multiplicity in MIN region vs. leading-cluster pT
363 fhRegionMultMinVsEt = new TH1F("hRegionMultMinVsEt", "N_{ch}^{90, min} vs Leading cluster p_{T}", bins, min, max);
364 fhRegionMultMinVsEt->SetXTitle("p_{T} (GeV/c)");
365 fhRegionMultMinVsEt->Sumw2();
366 fListOfHistos->Add( fhRegionMultMinVsEt ); // At(15)
367
368 //Average sum pT in TRANSVERSE(MIN+MAX) region vs. leading-cluster pT
369 fhRegionAveSumPtVsEt = new TH1F("hRegionAveSumPtVsEt", "(P_{T}^{90, max} + P_{T}^{90, min})/2 vs Leading cluster p_{T}", bins, min, max);
370 fhRegionAveSumPtVsEt->SetXTitle("p_{T} (GeV/c)");
371 fhRegionAveSumPtVsEt->Sumw2();
372 fListOfHistos->Add( fhRegionAveSumPtVsEt ); // At(16)
373
374 //Difference sum pT (MAX-MIN) vs. leading-cluster pT
375 fhRegionDiffSumPtVsEt= new TH1F("hRegionDiffSumPtVsEt", "(P_{T}^{90, max} - P_{T}^{90, min}) vs Leading cluster p_{T}", bins, min, max);
376 fhRegionDiffSumPtVsEt->SetXTitle("P_{T} (GeV/c)");
377 fhRegionDiffSumPtVsEt->Sumw2();
378 fListOfHistos->Add( fhRegionDiffSumPtVsEt ); // At(17)
379
380 //Average track pT in MAX region vs. leading-cluster pT
381 fhRegionAvePartPtMaxVsEt = new TH1F("hRegionAvePartPtMaxVsEt", "#LTp_{T}#GT^{90, max} vs Leading cluster p_{T}", bins, min, max);
382 fhRegionAvePartPtMaxVsEt->SetXTitle("p_{T} (GeV/c)");
383 fhRegionAvePartPtMaxVsEt->Sumw2();
384 fListOfHistos->Add( fhRegionAvePartPtMaxVsEt ); // At(18)
385
386 //Average track pT in MIN region vs. leading-cluster pT
387 fhRegionAvePartPtMinVsEt = new TH1F("hRegionAvePartPtMinVsEt", "#LTp_{T}#GT^{90, min} vs Leading cluster p_{T}", bins, min, max);
388 fhRegionAvePartPtMinVsEt->SetXTitle("p_{T} (GeV/c)");
389 fhRegionAvePartPtMinVsEt->Sumw2();
390 fListOfHistos->Add( fhRegionAvePartPtMinVsEt ); // At(19)
391
392 //Maximum track pT in MAX region vs. leading-cluster pT
393 fhRegionMaxPartPtMaxVsEt = new TH1F("hRegionMaxPartPtMaxVsEt", "max(p_{T})^{90} vs Leading cluster p_{T}", bins, min, max);
394 fhRegionMaxPartPtMaxVsEt->SetXTitle("p_{T} (GeV/c)");
395 fhRegionMaxPartPtMaxVsEt->Sumw2();
396 fListOfHistos->Add( fhRegionMaxPartPtMaxVsEt ); // At(20)
397
398 //Multiplicity in FORWARD region
399 fhRegForwardMult = new TH2F("hRegForwardMult", "N_{ch}^{forward}", bins, min, max, 21, -0.5, 20.5);
400 fhRegForwardMult->SetXTitle("N_{ch tracks}");
401 fhRegForwardMult->Sumw2();
402 fListOfHistos->Add( fhRegForwardMult ); // At(25)
403
404 //Sum pT in FORWARD region vs. multiplicity
405 fhRegForwardSumPtvsMult = new TH2F("hRegForwardSumPtvsMult", "Forward #Sigma p_{T} vs. Multiplicity ", bins, min, max, 21, -0.5, 20.5);
406 fhRegForwardSumPtvsMult->SetYTitle("#Sigma p_{T} (GeV/c)");
407 fhRegForwardSumPtvsMult->SetXTitle("N_{charge}");
408 fhRegForwardSumPtvsMult->Sumw2();
409 fListOfHistos->Add( fhRegForwardSumPtvsMult ); // At(26);
410
411 //Multiplicity in BACKWARD region
412 fhRegBackwardMult = new TH2F("hRegBackwardMult", "N_{ch}^{backward}", bins, min, max, 21, -0.5, 20.5);
413 fhRegBackwardMult->SetXTitle("N_{ch tracks}");
414 fhRegBackwardMult->Sumw2();
415 fListOfHistos->Add( fhRegBackwardMult ); // At(27)
416
417 //Sum pT in BACKWARD region vs. multiplicity
418 fhRegBackwardSumPtvsMult = new TH2F("hRegBackwardSumPtvsMult", "Backward #Sigma p_{T} vs. Multiplicity ", bins, min, max, 21, -0.5, 20.5);
419 fhRegBackwardSumPtvsMult->SetYTitle("#Sigma p_{T} (GeV/c)");
420 fhRegBackwardSumPtvsMult->SetXTitle("N_{charge}");
421 fhRegBackwardSumPtvsMult->Sumw2();
422 fListOfHistos->Add( fhRegBackwardSumPtvsMult ); // At(28);
423
424 //Track pT distribution in FORWARD region vs. leading-cluster pT
425 fhRegForwardPartPtDistVsEt = new TH2F("hRegForwardPartPtDistVsEt", Form( "dN/dP_{T} |#eta|<%3.1f vs Leading cluster p_{T}", fTrackEtaCut), 100,0.,50., bins, min, max);
426 fhRegForwardPartPtDistVsEt->SetYTitle("Leading cluster p_{T}");
427 fhRegForwardPartPtDistVsEt->SetXTitle("p_{T} (GeV/c)");
428 fhRegForwardPartPtDistVsEt->Sumw2();
429 fListOfHistos->Add( fhRegForwardPartPtDistVsEt ); // At(29)
430
431 //Track pT distribution in BACKWARD region vs. leading-cluster pT
432 fhRegBackwardPartPtDistVsEt = new TH2F("hRegBackwardPartPtDistVsEt", Form( "dN/dP_{T} |#eta|<%3.1f vs Leading cluster p_{T}", fTrackEtaCut), 100,0.,50., bins, min, max);
433 fhRegBackwardPartPtDistVsEt->SetYTitle("Leading cluster p_{T}");
434 fhRegBackwardPartPtDistVsEt->SetXTitle("p_{T}");
435 fhRegBackwardPartPtDistVsEt->Sumw2();
436 fListOfHistos->Add( fhRegBackwardPartPtDistVsEt ); // At(30)
437
438 //Multiplicity in TRANSVERSE (MIN+MAX) region
439 fhRegTransMult = new TH2F("hRegTransMult", "N_{ch}^{transv}", bins, min, max, 21, -0.5, 20.5);
440 fhRegTransMult->SetXTitle("N_{ch tracks}");
441 fhRegTransMult->Sumw2();
442 fListOfHistos->Add( fhRegTransMult ); // At(31)
443
444 //Sum pT in TRANSVERSE (MIN+MAX) region vs. multiplicity
445 fhRegTransSumPtVsMult = new TH2F("hRegTransSumPtVsMult", "Transverse #Sigma p_{T} vs. Multiplicity ",bins, min, max, 21, -0.5, 20.5);
446 fhRegTransSumPtVsMult->SetYTitle("#Sigma p_{T} (GeV/c)");
447 fhRegTransSumPtVsMult->SetXTitle("N_{charge}");
448 fhRegTransSumPtVsMult->Sumw2();
449 fListOfHistos->Add( fhRegTransSumPtVsMult ); // At(32);
450
451 //Sum pT in MIN region per cluster pT bin
452 fhMinRegSumPtJetPtBin = new TH2F("hMinRegSumPtJetPtBin", "Transverse Min Reg #Sigma p_{T} per cluster pT bin",bins, min, max, 50, 0., 20.);
453 fhMinRegSumPtJetPtBin->SetXTitle("Leading cluster p_{T}");
454 fhMinRegSumPtJetPtBin->Sumw2();
455 fListOfHistos->Add( fhMinRegSumPtJetPtBin ); // At(33)
456
457 //Sum pT in MAX region per cluster pT bin
458 fhMaxRegSumPtJetPtBin = new TH2F("hMaxRegSumPtJetPtBin", "Transverse Max Reg #Sigma p_{T} per cluster pT bin", bins, min, max, 50, 0., 20.);
459 fhMaxRegSumPtJetPtBin->SetXTitle("Leading cluster p_{T}");
460 fhMaxRegSumPtJetPtBin->Sumw2();
461 fListOfHistos->Add( fhMaxRegSumPtJetPtBin ); // At(34)
462
463 //Multiplicity in main vertex
464 fhVertexMult = new TH1F("hVertexMult", "Multiplicity in Main Vertex", 81, -0.5 , 80.5);
465 fhVertexMult->SetXTitle("Main Vertex Multiplicity");
466 fhVertexMult->Sumw2();
467 fListOfHistos->Add( fhVertexMult ); //At(35)
468
469 fh1Xsec = new TProfile("h1Xsec","xsec from pyxsec.root",1,0,1);
470 fh1Xsec->GetXaxis()->SetBinLabel(1,"<#sigma>");
471 fh1Xsec->Sumw2();
472 fListOfHistos->Add( fh1Xsec ); //At(36)
473
474 fh1Trials = new TH1F("h1Trials","trials from pyxsec.root",1,0,1);
475 fh1Trials->GetXaxis()->SetBinLabel(1,"#sum{ntrials}");
476 fh1Trials->Sumw2();
477 fListOfHistos->Add( fh1Trials ); //At(37)
478
479 return fListOfHistos;
480}
481
482
483//____________________________________________________________________
484void AliHistogramsUE::DrawUE(Int_t debug){
485
486 // To draw histograms at the end of task running
487 // Normalize histos to region area TODO:
488 // Normalization done at Analysis, taking into account
489 // area variations on per-event basis (cone case)
490
491 //HIGH WARNING!!!!!: DO NOT SCALE ANY OF THE ORIGINAL HISTOGRAMS
492 //MAKE A COPY, DRAW IT, And later sacale that copy. CAF Issue!!!!!
493
494 Int_t binsPtInHist = fhEleadingPt->GetNbinsX();
495 Double_t minJetPtInHist = fhEleadingPt->GetXaxis()->GetBinLowEdge(1);
496 Double_t maxJetPtInHist = fhEleadingPt->GetXaxis()->GetBinUpEdge(binsPtInHist);
497
498 //Sum pT
499 TCanvas* c1 = new TCanvas("c1",Form("sumPt dist (%s)", GetTitle()),60,60,1100,700);
500 c1->Divide(2,2);
501 c1->cd(1);
502 TH1F *h1r = new TH1F("hRegionEtvsSumPtMax" , "", binsPtInHist, minJetPtInHist, maxJetPtInHist);
503 //TH1F *h1r = new TH1F();
504 h1r->Divide(fhRegionSumPtMaxVsEt,fhEleadingPt,1,1);
505 //h1r->Scale( areafactor );
506 h1r->SetMarkerStyle(20);
507 h1r->SetXTitle("P_{T} of Leading Jet (GeV/c)");
508 h1r->SetYTitle("P_{T}^{90, max}");
509 h1r->DrawCopy("p");
510
511 c1->cd(2);
512 TH1F *h2r = new TH1F("hRegionEtvsSumPtMin" , "", binsPtInHist, minJetPtInHist, maxJetPtInHist);
513 h2r->Divide(fhRegionSumPtMinVsEt,fhEleadingPt,1,1);
514 //h2r->Scale( areafactor );
515 h2r->SetMarkerStyle(20);
516 h2r->SetXTitle("P_{T} of Leading Jet (GeV/c)");
517 h2r->SetYTitle("P_{T}^{90, min}");
518 h2r->DrawCopy("p");
519
520 c1->cd(3);
521 TH1F *h4r = new TH1F("hRegionEtvsDiffPt" , "", binsPtInHist, minJetPtInHist, maxJetPtInHist);
522 //TH1F *h41r = new TH1F("hRegForwvsDiffPt" , "", fbinsPtInHist, fMinJetPtInHist, fMaxJetPtInHist);
523 //TH1F *h42r = new TH1F("hRegBackvsDiffPt" , "", fbinsPtInHist, fMinJetPtInHist, fMaxJetPtInHist);
524 //h41r->Divide(fhRegForwardSumPtVsEt,fhEleadingPt,1,1);
525 //h42r->Divide(fhRegBackwardSumPtVsEt,fhEleadingPt,1,1);
526 h4r->Divide(fhRegionAveSumPtVsEt,fhEleadingPt,1,1);
527 //h4r->Scale(2.); // make average
528 //h4r->Scale( areafactor );
529 h4r->SetYTitle("#DeltaP_{T}^{90}");
530 h4r->SetXTitle("P_{T} of Leading Jet (GeV/c)");
531 h4r->SetMarkerStyle(20);
532 h4r->DrawCopy("p");
533
534 c1->cd(4);
535 TH1F *h5r = new TH1F("hRegionMultMaxVsEtleading", "", binsPtInHist, minJetPtInHist, maxJetPtInHist);
536 TH1F *h6r = new TH1F("hRegionMultMinVsEtleading", "", binsPtInHist, minJetPtInHist, maxJetPtInHist);
537 h5r->Divide(fhRegionMultMaxVsEt,fhEleadingPt,1,1);
538 h6r->Divide(fhRegionMultMinVsEt,fhEleadingPt,1,1);
539 //h5r->Scale( areafactor );
540 h5r->SetYTitle("N_{Tracks}^{90}");
541 h5r->SetXTitle("P_{T} of Leading Jet (GeV/c)");
542 h5r->SetMarkerStyle(20);
543 h5r->DrawCopy("p");
544 h6r->SetMarkerStyle(21);
545 h6r->SetMarkerColor(2);
546 h6r->SetYTitle("N_{Tracks}^{90}");
547 h6r->SetXTitle("P_{T} of Leading Jet (GeV/c)");
548 h6r->DrawCopy("p same");
549 c1->Update();
550
551 //Get Normalization
552 //Double_t xsec = fh1Xsec->GetBinContent(1);
553 Double_t xsec = fh1Xsec->GetBinContent(1);
554 Double_t ntrials = fh1Trials->GetBinContent(1);
555 Double_t normFactor = xsec/ntrials;
556 if(debug > 1)Printf("xSec %f nTrials %f Norm %f \n",xsec,ntrials,normFactor);
557
558
559 //Jet pT distribution
560 TCanvas* c2 = new TCanvas("c2","Jet Pt dist",160,160,1200,800);
561 TH1 * copy = 0;
562 c2->Divide(2,2);
563 c2->cd(1);
564 fhEleadingPt->SetMarkerStyle(20);
565 fhEleadingPt->SetMarkerColor(2);
566 //if( normFactor > 0.) fhEleadingPt->Scale(normFactor);
567 //fhEleadingPt->Draw("p");
568 copy = fhEleadingPt->DrawCopy("p");
569 if( normFactor > 0.) copy->Scale(normFactor);
570 gPad->SetLogy();
571
572 c2->cd(2);
573 Int_t xbin1 = fhdNdEtaPhiDist->GetYaxis()->FindFixBin(minJetPtInHist);
574 Int_t xbin2 = fhdNdEtaPhiDist->GetYaxis()->FindFixBin(maxJetPtInHist);
575 TH1D* dNdEtaPhiDistAllJets = fhdNdEtaPhiDist->ProjectionX("dNdEtaPhiDistAllJets",xbin1,xbin2);
576 dNdEtaPhiDistAllJets->SetMarkerStyle(20);
577 dNdEtaPhiDistAllJets->SetMarkerColor(2);
578 dNdEtaPhiDistAllJets->DrawCopy("p");
579 gPad->SetLogy();
580
581 c2->cd(3);
582 fhNJets->DrawCopy();
583 //c2->cd(4);
584 //fhValidRegion->DrawCopy("p");
585 //fhTransRegPartPtDist = (TH1F*)fListOfHistos->At(2);
586 //fhRegionMultMin = (TH1F*)fListOfHistos->At(3);
587 //fhMinRegAvePt = (TH1F*)fListOfHistos->At(4);
588 //fhMinRegSumPt = (TH1F*)fListOfHistos->At(5);
589 //fhMinRegMaxPtPart = (TH1F*)fListOfHistos->At(6);
590 //fhMinRegSumPtvsMult = (TH1F*)fListOfHistos->At(7);
591 c2->Update();
592
593 //pT distributions
594 TCanvas* c3 = new TCanvas("c3"," pT dist",160,160,1200,800);
595 c3->Divide(2,2);
596 c3->cd(1);
597 //fhTransRegPartPtDist->SetMarkerStyle(20);
598 //fhTransRegPartPtDist->SetMarkerColor(2);
599 //fhTransRegPartPtDist->Scale(areafactor/fhTransRegPartPtDist->GetEntries());
600 //fhTransRegPartPtDist->DrawCopy("p");
601 //gPad->SetLogy();
602
603
604 c3->cd(2);
605 fhMinRegSumPt->SetMarkerStyle(20);
606 fhMinRegSumPt->SetMarkerColor(2);
607 //fhMinRegSumPt->Scale(areafactor);
608 fhMinRegSumPt->DrawCopy("p");
609 gPad->SetLogy();
610
611 c3->cd(3);
612 fhMinRegAvePt->SetMarkerStyle(20);
613 fhMinRegAvePt->SetMarkerColor(2);
614 //fhMinRegAvePt->Scale(areafactor);
615 fhMinRegAvePt->DrawCopy("p");
616 gPad->SetLogy();
617
618 c3->cd(4);
619 TH1F *h7r = new TH1F("hRegionMultMinVsMult", "", 21, -0.5, 20.5);
620 h7r->Divide(fhMinRegSumPtvsMult,fhRegionMultMin,1,1);
621 h7r->SetMarkerStyle(20);
622 h7r->SetMarkerColor(2);
623 h7r->DrawCopy("p");
624 c3->Update();
625
626
627
628 //Save canvas
629 c1->SaveAs("c1.pdf");
630 AliInfo("Canvas 1 saved");
631 c2->SaveAs("c2.pdf");
632 AliInfo("Canvas 2 saved");
633 c3->SaveAs("c3.pdf");
634 AliInfo("Canvas 3 saved");
635
636}
637
638//____________________________________________________________________
639void AliHistogramsUE::FillHistogram(const char* name, Double_t fillX){
640
641 // Fill 1D histogram with double
642 ((TH1F*)fListOfHistos->FindObject(name))->Fill(fillX);
643
644}
645
646//____________________________________________________________________
647void AliHistogramsUE::FillHistogram(const char* name, Int_t fillX){
648
649 // Fill 1D histogram with integer
650 ((TH1F*)fListOfHistos->FindObject(name))->Fill(fillX);
651
652}
653
654//____________________________________________________________________
655void AliHistogramsUE::FillHistogram(const char* name, Double_t fillX, Double_t fillY){
656
657 // Case of TH1F with weight or TH2F w/o weight
658 TObject *obj = fListOfHistos->FindObject(name);
659 if (obj->InheritsFrom("TH1F")){
660 ((TH1F*)fListOfHistos->FindObject(name))->Fill(fillX, fillY);
661 } else {
662 ((TH2F*)fListOfHistos->FindObject(name))->Fill(fillX, fillY);
663 }
664
665
666}
667
668//____________________________________________________________________
669void AliHistogramsUE::FillHistogram(const char* name, Double_t fillX, Double_t fillY, Double_t weight){
670
671 // Fill 2D histogram with double and weight
672 ((TH2F*)fListOfHistos->FindObject(name))->Fill(fillX, fillY, weight);
673
674}
675
676//____________________________________________________________________
677void AliHistogramsUE::FillHistogram(const char* name, Double_t fillX, Int_t fillY, Double_t weight){
678
679 // Fill 2D histogram with integer and weight
680 ((TH2F*)fListOfHistos->FindObject(name))->Fill(fillX, fillY, weight);
681
682}
683
684//____________________________________________________________________
685TObjArray* AliHistogramsUE::GetHistosForPlotting(TString data, TString branches){
686
687 // Instance filled histos for plotting purpose
688 printf("Creating histograms ... \n");
689
690 printf("Reading file: %s\n",data.Data());
691
692 // Read input files -----------------------------------------
693 TFile *fdata = new TFile(data.Data());
694 TDirectoryFile *ddata[20];
695 TList *ldata[20];
696
697 TObjArray *arrb=branches.Tokenize(";");
698 TIter next(arrb);
699 TObject *o=0;
700 Int_t br=0;
701 while ( (o=next()) ){
702 ddata[br] = (TDirectoryFile*)fdata->Get(Form("PWG4_UE%s",o->GetName()));
703 if(!ddata[br]) printf("ERROR: No histo dir found! \n");
704 ldata[br] = (TList*)ddata[br]->Get(Form("histosUE%s",o->GetName()));
705 printf("Reading branch: %s\n",o->GetName());
706 if(!ldata[br]) printf("ERROR: No histo list found! \n");
707 br++;
708 }
709
710 TObjArray *arr=new TObjArray;
711
712 TH1F *hjets[20]; // accepted leading jets
713 TH1F *hnjets[20]; // number of accepted jets
714 TH2F *hetaphi[20]; // delta-phi particle-jet correlation
715 TH2F *hptfull[20]; // particle pT all regions vs. jet pT
716 TH2F *hpttransv[20]; // particle pT transv. regions vs. jet pT
717 TH1F *hmax[20]; // sum pT in MAX region
718 TH1F *hmin[20]; // sum pT in MIN region
719 TH1F *hmultmax[20]; // multiplicity in MAX region
720 TH1F *hmultmin[20]; // multiplicity in MIN region
721
722 for (Int_t i =0; i<br; i++){
723
724 //Number of jets --------------------------------------------
725 hnjets[i] = (TH1F*) ldata[i]->FindObject("fhNJets");
726 hnjets[i]->GetXaxis()->SetTitle("Number of jets");
727 hnjets[i]->GetYaxis()->SetTitle("1/n_{ev} dN/dN_{jets}");
728 hnjets[i]->SetMarkerStyle(20);
729 hnjets[i]->SetMarkerColor(1+i);
730
731
732 //Leading jets ----------------------------------------------
733 hjets[i] = (TH1F*) ldata[i]->FindObject("hEleadingPt");
734 hjets[i]->GetXaxis()->SetTitle("P_{T} (GeV/c)");
735 hjets[i]->GetYaxis()->SetTitle("1/n_{ev} dN/dp_{T} (|#eta<0.5|)");
736 hjets[i]->SetMarkerStyle(20);
737 hjets[i]->SetMarkerColor(1+i);
738 hjets[i]->SetMinimum(0.1);
739 hjets[i]->SetMaximum(1000.);
740
741
742 //Transverse Region MAX -------------------------------------
743 hmax[i] = (TH1F*) ldata[i]->FindObject("hRegionSumPtMaxVsEt");
744 if (!hmax[i])AliInfo("Histo not found!!!");
745 hmax[i]->GetXaxis()->SetTitle("Leading jet P_{T} (GeV/c)");
746 hmax[i]->GetYaxis()->SetTitle("P_{T}^{90,max} (GeV/c)");
747 hmax[i]->SetMarkerStyle(20);
748 hmax[i]->SetMarkerColor(1+i);
749 hmax[i]->Divide(hjets[i]); // normalize for jet spectrum
750 hmax[i]->SetMaximum(5.);
751
752
753 //Transverse Region MIN -------------------------------------
754 hmin[i] = (TH1F*) ldata[i]->FindObject("hRegionSumPtMinVsEt");
755 hmin[i]->GetXaxis()->SetTitle("Leading jet P_{T} (GeV/c)");
756 hmin[i]->GetYaxis()->SetTitle("P_{T}^{90,min} (GeV/c)");
757 hmin[i]->SetMarkerStyle(20);
758 hmin[i]->SetMarkerColor(1+i);
759 hmin[i]->SetMaximum(3.);
760 hmin[i]->Divide(hjets[i]); // normalize for jet spectrum
761
762
763 //Multiplicity MAX ------------------------------------------
764 hmultmax[i] = (TH1F*) ldata[i]->FindObject("hRegionMultMaxVsEt");
765 hmultmax[i]->GetXaxis()->SetTitle("Leading Jet P_{T} (GeV/c)");
766 hmultmax[i]->GetYaxis()->SetTitle("N_{ch}^{90,max}");
767 hmultmax[i]->SetMarkerStyle(20);
768 hmultmax[i]->SetMarkerColor(1+i);
769 hmultmax[i]->SetMaximum(10.);
770 hmultmax[i]->Divide(hjets[i]); // normalize for jet spectrum
771
772
773 //Multiplicity MIN ------------------------------------------
774 hmultmin[i] = (TH1F*) ldata[i]->FindObject("hRegionMultMinVsEt");
775 hmultmin[i]->GetXaxis()->SetTitle("Leading Jet P_{T} (GeV/c)");
776 hmultmin[i]->GetYaxis()->SetTitle("N_{ch}^{90,min}");
777 hmultmin[i]->SetMarkerStyle(20);
778 hmultmin[i]->SetMarkerColor(1+i);
779 hmultmin[i]->SetMaximum(3.);
780 hmultmin[i]->Divide(hjets[i]); // normalize for jet spectrum
781
782
783 // Phi-correlation with leading jet --------------------------
784 hetaphi[i] = (TH2F*) ldata[i]->FindObject("hdNdEtaPhiDist");
785 hetaphi[i]->GetXaxis()->SetTitle("#Delta #phi (w.r.t. leading jet)");
786 hetaphi[i]->SetMarkerStyle(20);
787
788
789 // pT distribution in full region vs. jet pT --------------------------
790 hptfull[i] = (TH2F*) ldata[i]->FindObject("hFullRegPartPtDistVsEt");
791 hptfull[i]->GetYaxis()->SetTitle("Leading Jet P_{T} (GeV/c)");
792 hptfull[i]->GetXaxis()->SetTitle("Track P_{T} (GeV/c)");
793
794
795 // pT distribution in transv region vs. jet pT --------------------------
796 hpttransv[i] = (TH2F*) ldata[i]->FindObject("hTransRegPartPtDistVsEt");
797 hpttransv[i]->GetYaxis()->SetTitle("Leading Jet P_{T} (GeV/c)");
798 hpttransv[i]->GetXaxis()->SetTitle("Track P_{T} (GeV/c)");
799
800
801 // Return Histos
802 arr->Add(hnjets[i]); //at 0 number of jets
803 arr->Add(hjets[i]); //at 1 leading jets
804 arr->Add(hmax[i]); //at 2 sum pT MAX
805 arr->Add(hmin[i]); //at 3 sumpT MIN
806 arr->Add(hmultmax[i]); //at 4 multiplicity MAX
807 arr->Add(hmultmin[i]); //at 5 multiplicity MIN
808 arr->Add(hetaphi[i]); //at 6 phi correlation
809 arr->Add(hptfull[i]); //at 7 pT distr in full region
810 arr->Add(hpttransv[i]); //at 8 pT distr in transv region
811
812 }
813 return arr;
814}
815
816//_______________________________________________________________________________________
817void AliHistogramsUE::SetStyle(){
818
819 // Set plotting style
820 gPad->SetFrameFillColor(0);
821 gPad->SetFillColor(0);
822 gPad->SetBorderSize(2);
823 gPad->SetGridy();
824 gPad->SetFrameBorderMode(0);
825 //gStyle->SetOptStat(0);
826 gStyle->SetOptTitle(0);
827
828}
829
830
831//____________________________________________________________________
832TList* AliHistogramsUE::GetListOfHistos(){
833
834 // Return list of relevant histograms
835 return fListOfHistos;
836
837}
838
839//____________________________________________________________________
840void AliHistogramsUE::PlotBranchesUE(TString file, TString branches, Double_t minJetProjection){
841
842 // Function to be called by external macro to plot analysis from different jet branches
843
844 //Count the branches
845 TObjArray *arr=branches.Tokenize(";");
846 TIter next(arr);
847 TObject *o=0;
848 Int_t br=0;
849 while ( (o=next()) ){
850 br++;
851 }
852
853
854 //Canvas
855 TObjArray *arrC=CreateCanvas(9);
856
857 //Histograms
858 TObjArray *arrH=GetHistosForPlotting(file.Data(),branches.Data());
859 TH1F *hnjets[20];
860 TH1F *hjets[20];
861 TH1F *hmax[20];
862 TH1F *hmin[20];
863 TH1F *hmultmax[20];
864 TH1F *hmultmin[20];
865 TH2F *hetaphi[20];
866 TH2F *hptfull[20];
867 TH2F *hpttransv[20];
868
869 for (Int_t i= 0; i<br; i++){
870 hnjets[i]=((TH1F*)arrH->At(9*i));
871 hjets[i]=((TH1F*)arrH->At(1+(9*i)));
872 hmax[i]=((TH1F*)arrH->At(2+(9*i)));
873 hmin[i]=((TH1F*)arrH->At(3+(9*i)));
874 hmultmax[i]=((TH1F*)arrH->At(4+(9*i)));
875 hmultmin[i]=((TH1F*)arrH->At(5+(9*i)));
876 hetaphi[i]=((TH2F*)arrH->At(6+(9*i)));
877 hptfull[i]=((TH2F*)arrH->At(7+(9*i)));
878 hpttransv[i]=((TH2F*)arrH->At(8+(9*i)));
879 }
880
881 //Define jet-pT range in projections
882 Int_t binstartprojection = hetaphi[0]->GetYaxis()->FindFixBin(minJetProjection); //be careful...
883 Int_t binstopprojection = hetaphi[0]->GetYaxis()->GetNbins();
884
885 Double_t normphi[20];
886
887 for (Int_t i= 0; i<br; i++){
888 normphi[i]=hjets[i]->Integral(binstartprojection,binstopprojection);
889 hnjets[i]->Scale(1./(hnjets[i]->GetBinWidth(1)*hnjets[i]->GetEntries()));
890 hjets[i]->Scale(1./(hjets[i]->GetBinWidth(1)*hnjets[i]->GetEntries()));
891 }
892
893 //LEGENDS
894 //Legend jets
895 TLegend *leg=new TLegend(0.5,0.6,0.89,0.89); // for jet spectrum
896 leg->SetFillColor(0);
897 leg->SetHeader("Jet Finders:");
898 //Legend density
899 TLegend *legd = new TLegend(0.1077586,0.6016949,0.4971264,0.8919492,NULL,"brNDC");
900 legd->SetFillColor(0);
901 legd->SetHeader("Jet Finders:");
902 //Legend pT distributions
903 TLegend *legp = new TLegend(0.1364943,0.1292373,0.5258621,0.4194915,NULL,"brNDC");
904 legp->SetFillColor(0);
905 legp->SetHeader("Jet Finders:");
906
907 arr=branches.Tokenize(";");
908 TIter next1(arr);
909 o=0;
910 Int_t brleg=0;
911 while ( (o=next1()) ){
912 leg->AddEntry(hjets[brleg],Form("UE%s",o->GetName()),"p");
913 legd->AddEntry(hjets[brleg],Form("UE%s",o->GetName()),"p");
914 legp->AddEntry(hjets[brleg],Form("UE%s",o->GetName()),"p");
915 brleg++;
916 }
917
918 //1) NUMBER OF CLUSTERS
919 TCanvas *c0=((TCanvas*)arrC->At(0));
920 c0->cd();
921 SetStyle();
922 gPad->SetLogy();
923 hnjets[0]->SetMaximum(1.4);
924 hnjets[0]->Draw("E1");
925 for (Int_t i= 1; i<br; i++){
926 hnjets[i]->Draw("same");
927 }
928 leg->Draw("same");
929
930 //2) LEADING CLUSTERS pT
931 TCanvas *c1=((TCanvas*)arrC->At(1));
932 c1->cd();
933 SetStyle();
934 gPad->SetLogy();
935 hjets[0]->Draw("E1");
936 hjets[0]->SetMaximum(1.4);
937 for (Int_t i= 1; i<br; i++){
938 hjets[i]->Draw("same");
939 }
940 leg->Draw("same");
941
942 //3) SUM-pT IN MAX REGION
943 TCanvas *c2=((TCanvas*)arrC->At(2));
944 c2->cd();
945 SetStyle();
946 hmax[0]->Draw("E1");
947 for (Int_t i= 1; i<br; i++){
948 hmax[i]->Draw("same");
949 }
950 leg->Draw("same");
951
952 //4) SUM-pT IN MIN REGION
953 TCanvas *c3=((TCanvas*)arrC->At(3));
954 c3->cd();
955 SetStyle();
956 hmin[0]->GetYaxis()->SetRangeUser(0.,2.5);
957 hmin[0]->Draw("E1");
958 for (Int_t i= 1; i<br; i++){
959 hmin[i]->Draw("same");
960 }
961 leg->Draw("same");
962
963 //5) MULTIPLICITY IN MAX REGION
964 TCanvas *c4=((TCanvas*)arrC->At(4));
965 c4->cd();
966 SetStyle();
967 hmultmax[0]->GetYaxis()->SetRangeUser(0.,5.);
968 hmultmax[0]->Draw("E1");
969 for (Int_t i= 1; i<br; i++){
970 hmultmax[i]->Draw("same");
971 }
972 leg->Draw("same");
973
974 //6) MULTIPLICITY IN MIN REGION
975 TCanvas *c5=((TCanvas*)arrC->At(5));
976 c5->cd();
977 SetStyle();
978 hmultmin[0]->GetYaxis()->SetRangeUser(0.,2.5);
979 hmultmin[0]->Draw("E1");
980 for (Int_t i= 1; i<br; i++){
981 hmultmin[i]->Draw("same");
982 }
983 leg->Draw("same");
984
985 //7) JET-TRACK CORRELATION
986 TCanvas *c6=((TCanvas*)arrC->At(6));
987 c6->cd();
988 SetStyle();
989 gPad->SetLogy();
990 TH1D *tmpetaphi[20];
991 for (Int_t i= 0; i<br; i++){
992 tmpetaphi[i]=new TH1D();
993 tmpetaphi[i]=hetaphi[i]->ProjectionX(Form("data%d",i),binstartprojection);
994 tmpetaphi[i]->GetYaxis()->SetTitle("1/n_{lj} dN/d#Delta #phi");
995 tmpetaphi[i]->Scale(1./(hetaphi[i]->GetBinWidth(1)*normphi[i]));
996 tmpetaphi[i]->SetMarkerColor(i+1);
997 tmpetaphi[i]->GetXaxis()->SetLimits(-3.*TMath::Pi()/2.,TMath::Pi()/2.);
998 tmpetaphi[i]->GetYaxis()->SetRangeUser(0.5,20.);
999 if (i==0) tmpetaphi[i]->Draw("E1");
1000 else tmpetaphi[i]->Draw("same");
1001 // evaluate mean multiplicity in transverse regions
1002 Double_t err1=0.;
1003 Double_t err2=0.;
1004 Double_t err3=0.;
1005 Double_t mean=(tmpetaphi[i]->IntegralAndError(1,5,err1)+tmpetaphi[i]->IntegralAndError(27,36,err2)+tmpetaphi[i]->IntegralAndError(58,62,err3))/(20.);
1006 err1=TMath::Sqrt(err1*err1+err2*err2+err3*err3)/20.;
1007 Printf("Branch: %d MeanTransvMult: %f err: %f",i,mean,err1);
1008 }
1009 legd->Draw("same");
1010
1011 //8) TRACK-pT DISTRIBUTION IN FULL REGION
1012 TCanvas *c7=((TCanvas*)arrC->At(7));
1013 c7->cd();
1014 SetStyle();
1015 gPad->SetLogy();
1016 gPad->SetLogx();
1017 gPad->SetTitle("Full region (projection X)");
1018 TH1D *tmpfull[20];
1019 for (Int_t i= 0; i<br; i++){
1020 tmpfull[i]=new TH1D();
1021 tmpfull[i]=hptfull[i]->ProjectionX(Form("full%d",i),binstartprojection);
1022 tmpfull[i]->GetYaxis()->SetTitle("1/n_{lj} dN/dp_{T}");
1023 tmpfull[i]->Scale(1./(tmpfull[i]->GetBinWidth(1)*normphi[i]));
1024 tmpfull[i]->SetMarkerStyle(20);
1025 tmpfull[i]->SetMarkerColor(i+1);
1026 tmpfull[i]->SetMaximum(100.);
1027 if (i==0)tmpfull[i]->Draw("E1");
1028 else tmpfull[i]->Draw("same");
1029 }
1030 legp->Draw("same");
1031
1032 //9) TRACK-pT DISTRIBUTION IN TRANSVERSE (MIN+MAX) REGION
1033 TCanvas *c8=((TCanvas*)arrC->At(8));
1034 c8->cd();
1035 SetStyle();
1036 gPad->SetLogy();
1037 gPad->SetLogx();
1038 gPad->SetTitle("Transverse regions (projection X)");
1039 TH1D *tmptransv[20];
1040 for (Int_t i= 0; i<br; i++){
1041 tmptransv[i]=new TH1D();
1042 tmptransv[i]=hpttransv[i]->ProjectionX(Form("transv%d",i),binstartprojection);
1043 tmptransv[i]->GetYaxis()->SetTitle("1/n_{lj} dN/dp_{T}");
1044 tmptransv[i]->Scale(1./(tmptransv[i]->GetBinWidth(1)*normphi[i]));
1045 tmptransv[i]->SetMarkerStyle(20);
1046 tmptransv[i]->SetMarkerColor(i+1);
1047 tmptransv[i]->SetMaximum(10.);
1048 if (i==0)tmptransv[i]->Draw("E1");
1049 else tmptransv[i]->Draw("same");
1050 }
1051 legp->Draw("same");
1052}