]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG4/JetTasks/AliUEHistograms.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Added fit macro from M. Putis
[u/mrichter/AliRoot.git] / PWG4 / JetTasks / AliUEHistograms.cxx
CommitLineData
a75aacd6 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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: AliUEHistograms.cxx 20164 2007-08-14 15:31:50Z morsch $ */
17
18//
19//
20// encapsulates several AliUEHist objects for a full UE analysis plus additional control histograms
21//
22//
23// Author: Jan Fiete Grosse-Oetringhaus, Sara Vallero
24
25#include "AliUEHistograms.h"
26
27#include "AliCFContainer.h"
28#include "AliVParticle.h"
29
30#include "TList.h"
31#include "TH2F.h"
32#include "TH1F.h"
33#include "TH3F.h"
34#include "TMath.h"
35
36ClassImp(AliUEHistograms)
37
bf58cbde 38const Int_t AliUEHistograms::fgkUEHists = 3;
39
e0331fd9 40AliUEHistograms::AliUEHistograms(const char* name, const char* histograms) :
41 TNamed(name, name),
a75aacd6 42 fNumberDensitypT(0),
43 fSumpT(0),
44 fNumberDensityPhi(0),
45 fCorrelationpT(0),
46 fCorrelationEta(0),
47 fCorrelationPhi(0),
48 fCorrelationR(0),
49 fCorrelationLeading2Phi(0),
50 fCorrelationMultiplicity(0),
51 fEventCount(0),
52 fEventCountDifferential(0),
bf58cbde 53 fVertexContributors(0),
54 fCentralityDistribution(0)
a75aacd6 55{
56 // Constructor
bf58cbde 57 //
58 // the string histograms defines which histograms are created:
59 // 1 = NumberDensitypT
60 // 2 = SumpT
61 // 3 = NumberDensityPhi
62 // 4 = NumberDensityPhiCentrality (other multiplicity for Pb)
63
64 TString histogramsStr(histograms);
65
66 if (histogramsStr.Contains("1"))
67 fNumberDensitypT = new AliUEHist("NumberDensitypT");
68 if (histogramsStr.Contains("2"))
69 fSumpT = new AliUEHist("SumpT");
a75aacd6 70
bf58cbde 71 if (histogramsStr.Contains("3"))
72 fNumberDensityPhi = new AliUEHist("NumberDensityPhi");
73 else if (histogramsStr.Contains("4"))
74 fNumberDensityPhi = new AliUEHist("NumberDensityPhiCentrality");
a75aacd6 75
76 // do not add this hists to the directory
77 Bool_t oldStatus = TH1::AddDirectoryStatus();
78 TH1::AddDirectory(kFALSE);
79
80 fCorrelationpT = new TH2F("fCorrelationpT", ";p_{T,lead} (MC);p_{T,lead} (RECO)", 200, 0, 50, 200, 0, 50);
81 fCorrelationEta = new TH2F("fCorrelationEta", ";#eta_{T,lead} (MC);#eta_{T,lead} (RECO)", 200, -1, 1, 200, -1, 1);
82 fCorrelationPhi = new TH2F("fCorrelationPhi", ";#phi_{T,lead} (MC);#phi_{T,lead} (RECO)", 200, 0, TMath::TwoPi(), 200, 0, TMath::TwoPi());
83 fCorrelationR = new TH2F("fCorrelationR", ";R;p_{T,lead} (MC)", 200, 0, 2, 200, 0, 50);
84 fCorrelationLeading2Phi = new TH2F("fCorrelationLeading2Phi", ";#Delta #phi;p_{T,lead} (MC)", 200, -TMath::Pi(), TMath::Pi(), 200, 0, 50);
85 fCorrelationMultiplicity = new TH2F("fCorrelationMultiplicity", ";MC tracks;Reco tracks", 100, -0.5, 99.5, 100, -0.5, 99.5);
86
144bd037 87 fEventCount = new TH2F("fEventCount", ";step;event type;count", AliUEHist::fgkCFSteps+2, -2.5, -0.5 + AliUEHist::fgkCFSteps, 3, -0.5, 2.5);
a75aacd6 88 fEventCount->GetYaxis()->SetBinLabel(1, "ND");
89 fEventCount->GetYaxis()->SetBinLabel(2, "SD");
90 fEventCount->GetYaxis()->SetBinLabel(3, "DD");
91
92 fEventCountDifferential = new TH3F("fEventCountDifferential", ";p_{T,lead};step;event type", 100, 0, 50, AliUEHist::fgkCFSteps, -0.5, -0.5 + AliUEHist::fgkCFSteps, 3, -0.5, 2.5);
93 fEventCountDifferential->GetZaxis()->SetBinLabel(1, "ND");
94 fEventCountDifferential->GetZaxis()->SetBinLabel(2, "SD");
95 fEventCountDifferential->GetZaxis()->SetBinLabel(3, "DD");
96
97 fVertexContributors = new TH1F("fVertexContributors", ";contributors;count", 100, -0.5, 99.5);
98
bf58cbde 99 fCentralityDistribution = new TH1F("fCentralityDistribution", ";;count", fNumberDensityPhi->GetEventHist()->GetNBins(1), fNumberDensityPhi->GetEventHist()->GetAxis(1, 0)->GetBinLowEdge(1), fNumberDensityPhi->GetEventHist()->GetAxis(1, 0)->GetBinUpEdge(fNumberDensityPhi->GetEventHist()->GetNBins(1)));
100
a75aacd6 101 TH1::AddDirectory(oldStatus);
102}
103
d1c75d06 104//_____________________________________________________________________________
105AliUEHistograms::AliUEHistograms(const AliUEHistograms &c) :
e0331fd9 106 TNamed(fName, fTitle),
d1c75d06 107 fNumberDensitypT(0),
108 fSumpT(0),
109 fNumberDensityPhi(0),
110 fCorrelationpT(0),
111 fCorrelationEta(0),
112 fCorrelationPhi(0),
113 fCorrelationR(0),
114 fCorrelationLeading2Phi(0),
115 fCorrelationMultiplicity(0),
116 fEventCount(0),
117 fEventCountDifferential(0),
bf58cbde 118 fVertexContributors(0),
119 fCentralityDistribution(0)
d1c75d06 120{
121 //
122 // AliUEHistograms copy constructor
123 //
124
125 ((AliUEHistograms &) c).Copy(*this);
126}
127
a75aacd6 128//____________________________________________________________________
129AliUEHistograms::~AliUEHistograms()
130{
131 // Destructor
132
133 if (fNumberDensitypT)
134 {
135 delete fNumberDensitypT;
136 fNumberDensitypT = 0;
137 }
138
139 if (fSumpT)
140 {
141 delete fSumpT;
142 fSumpT = 0;
143 }
144
145 if (fNumberDensityPhi)
146 {
147 delete fNumberDensityPhi;
148 fNumberDensityPhi = 0;
149 }
150
151 if (fCorrelationpT)
152 {
153 delete fCorrelationpT;
154 fCorrelationpT = 0;
155 }
156
157 if (fCorrelationEta)
158 {
159 delete fCorrelationEta;
160 fCorrelationEta = 0;
161 }
162
163 if (fCorrelationPhi)
164 {
165 delete fCorrelationPhi;
166 fCorrelationPhi = 0;
167 }
168
169 if (fCorrelationR)
170 {
171 delete fCorrelationR;
172 fCorrelationR = 0;
173 }
174
175 if (fCorrelationLeading2Phi)
176 {
177 delete fCorrelationLeading2Phi;
178 fCorrelationLeading2Phi = 0;
179 }
180
181 if (fCorrelationMultiplicity)
182 {
183 delete fCorrelationMultiplicity;
184 fCorrelationMultiplicity = 0;
185 }
186
187 if (fEventCount)
188 {
189 delete fEventCount;
190 fEventCount = 0;
191 }
192
193 if (fEventCountDifferential)
194 {
195 delete fEventCountDifferential;
196 fEventCountDifferential = 0;
197 }
198
199 if (fVertexContributors)
200 {
201 delete fVertexContributors;
202 fVertexContributors = 0;
203 }
bf58cbde 204
205 if (fCentralityDistribution)
206 {
207 delete fCentralityDistribution;
208 fCentralityDistribution = 0;
209 }
a75aacd6 210}
211
ada1a03f 212AliUEHist* AliUEHistograms::GetUEHist(Int_t id)
213{
214 // returns AliUEHist object, useful for loops
215
216 switch (id)
217 {
218 case 0: return fNumberDensitypT; break;
219 case 1: return fSumpT; break;
220 case 2: return fNumberDensityPhi; break;
221 }
222
223 return 0;
224}
225
a75aacd6 226//____________________________________________________________________
227Int_t AliUEHistograms::CountParticles(TList* list, Float_t ptMin)
228{
229 // counts the number of particles in the list with a pT above ptMin
230 // TODO eta cut needed here?
231
232 Int_t count = 0;
233 for (Int_t j=0; j<list->GetEntries(); j++)
234 if (((AliVParticle*) list->At(j))->Pt() > ptMin)
235 count++;
236
237 return count;
238}
239
240//____________________________________________________________________
241void AliUEHistograms::Fill(Int_t eventType, AliUEHist::CFStep step, AliVParticle* leading, TList* toward, TList* away, TList* min, TList* max)
242{
243 // fills the UE event histograms
244 //
245 // this function needs the leading (track or jet or ...) and four lists of AliVParticles which contain the particles/tracks to be filled in the four regions
246
247 // if leading is not set, just fill event statistics
248 if (leading)
249 {
250 Int_t multiplicity = 0;
251
252 // TODO configurable?
253 Float_t ptMin = 0.15;
254 if (leading->Pt() > ptMin)
255 multiplicity++;
256
257 multiplicity += CountParticles(toward, ptMin);
258 multiplicity += CountParticles(away, ptMin);
259 multiplicity += CountParticles(min, ptMin);
260 multiplicity += CountParticles(max, ptMin);
261
262 FillRegion(AliUEHist::kToward, step, leading, toward, multiplicity);
263 FillRegion(AliUEHist::kAway, step, leading, away, multiplicity);
264 FillRegion(AliUEHist::kMin, step, leading, min, multiplicity);
265 FillRegion(AliUEHist::kMax, step, leading, max, multiplicity);
b1831bcb 266
a75aacd6 267 Double_t vars[2];
268 vars[0] = leading->Pt();
269 vars[1] = multiplicity;
bf58cbde 270 for (Int_t i=0; i<fgkUEHists; i++)
271 if (GetUEHist(i))
272 GetUEHist(i)->GetEventHist()->Fill(vars, step);
a75aacd6 273
274 fEventCountDifferential->Fill(leading->Pt(), step, eventType);
275 }
276
277 FillEvent(eventType, step);
278}
279
280//____________________________________________________________________
281void AliUEHistograms::FillRegion(AliUEHist::Region region, AliUEHist::CFStep step, AliVParticle* leading, TList* list, Int_t multiplicity)
282{
283 // loops over AliVParticles in list and fills the given region at the given step
284 //
285 // See also Fill(...)
286
287 for (Int_t i=0; i<list->GetEntries(); i++)
288 {
289 AliVParticle* particle = (AliVParticle*) list->At(i);
290
291 Double_t vars[5];
292 vars[0] = particle->Eta();
293 vars[1] = particle->Pt();
294 vars[2] = leading->Pt();
295 vars[3] = multiplicity;
296 vars[4] = leading->Phi() - particle->Phi();
2ac8dc5c 297 if (vars[4] > 1.5 * TMath::Pi())
298 vars[4] -= TMath::TwoPi();
299 if (vars[4] < -0.5 * TMath::Pi())
300 vars[4] += TMath::TwoPi();
d4abbc90 301
bf58cbde 302 if (fNumberDensitypT)
303 fNumberDensitypT->GetTrackHist(region)->Fill(vars, step);
304
305 if (fSumpT)
306 fSumpT->GetTrackHist(region)->Fill(vars, step, particle->Pt());
a75aacd6 307
308 // fill all in toward region (is anyway as function of delta phi!)
bf58cbde 309 if (fNumberDensityPhi)
310 fNumberDensityPhi->GetTrackHist(AliUEHist::kToward)->Fill(vars, step);
a75aacd6 311 }
312}
313
314//____________________________________________________________________
315void AliUEHistograms::Fill(AliVParticle* leadingMC, AliVParticle* leadingReco)
316{
317 // fills the correlation histograms
318
319 if (leadingMC)
320 {
321 fCorrelationpT->Fill(leadingMC->Pt(), leadingReco->Pt());
322 if (leadingMC->Pt() > 0.5)
323 {
324 fCorrelationEta->Fill(leadingMC->Eta(), leadingReco->Eta());
325 fCorrelationPhi->Fill(leadingMC->Phi(), leadingReco->Phi());
326 }
327
328 Float_t phiDiff = leadingMC->Phi() - leadingReco->Phi();
329 if (phiDiff > TMath::Pi())
330 phiDiff -= TMath::TwoPi();
331 if (phiDiff < -TMath::Pi())
332 phiDiff += TMath::TwoPi();
333
334 Float_t etaDiff = leadingMC->Eta() - leadingReco->Eta();
335
336 fCorrelationR->Fill(TMath::Sqrt(phiDiff * phiDiff + etaDiff * etaDiff), leadingMC->Pt());
337 fCorrelationLeading2Phi->Fill(phiDiff, leadingMC->Pt());
338 }
339 else
340 {
341 fCorrelationpT->Fill(1.0, leadingReco->Pt());
342 if (leadingReco->Pt() > 0.5)
343 {
344 fCorrelationEta->Fill(0.0, leadingReco->Eta());
345 fCorrelationPhi->Fill(0.0, leadingReco->Phi());
346 }
347 }
348}
bf58cbde 349
350//____________________________________________________________________
e0331fd9 351void AliUEHistograms::FillCorrelations(Int_t eventType, Int_t centrality, AliUEHist::CFStep step, TSeqCollection* particles, TSeqCollection* mixed)
bf58cbde 352{
353 // fills the fNumberDensityPhi histogram
354 //
355 // this function need a list of AliVParticles which contain the particles/tracks to be filled
e0331fd9 356 //
357 // if mixed is non-0, mixed events are filled, the trigger particle is from particles, the associated from mixed
bf58cbde 358
359 // if particles is not set, just fill event statistics
360 if (particles)
361 {
362 for (Int_t i=0; i<particles->GetEntries(); i++)
363 {
364 AliVParticle* triggerParticle = (AliVParticle*) particles->At(i);
e0331fd9 365 Int_t jMax = particles->GetEntries();
366 if (mixed)
367 jMax = mixed->GetEntries();
368
369 for (Int_t j=0; j<jMax; j++)
bf58cbde 370 {
e0331fd9 371 if (!mixed && i == j)
bf58cbde 372 continue;
373
e0331fd9 374 AliVParticle* particle = 0;
375 if (!mixed)
376 particle = (AliVParticle*) particles->At(j);
377 else
378 particle = (AliVParticle*) mixed->At(j);
bf58cbde 379
380 Double_t vars[5];
e0331fd9 381 vars[0] = triggerParticle->Eta() - particle->Eta();
bf58cbde 382 vars[1] = particle->Pt();
383 vars[2] = triggerParticle->Pt();
384 vars[3] = centrality;
385 vars[4] = triggerParticle->Phi() - particle->Phi();
386 if (vars[4] > 1.5 * TMath::Pi())
387 vars[4] -= TMath::TwoPi();
388 if (vars[4] < -0.5 * TMath::Pi())
389 vars[4] += TMath::TwoPi();
390
391 // fill all in toward region and no not use the other regions
392 fNumberDensityPhi->GetTrackHist(AliUEHist::kToward)->Fill(vars, step);
393 }
394
395 // once per trigger particle
396 Double_t vars[2];
397 vars[0] = triggerParticle->Pt();
398 vars[1] = centrality;
399 fNumberDensityPhi->GetEventHist()->Fill(vars, step);
400
401 fEventCountDifferential->Fill(triggerParticle->Pt(), step, eventType);
402 }
403 }
404
405 fCentralityDistribution->Fill(centrality);
406 FillEvent(eventType, step);
407}
a75aacd6 408
b1831bcb 409//____________________________________________________________________
410void AliUEHistograms::FillTrackingEfficiency(TObjArray* mc, TObjArray* recoPrim, TObjArray* recoAll, Int_t particleType)
411{
412 // fills the tracking efficiency objects
413 //
414 // mc: all primary MC particles
415 // recoPrim: reconstructed primaries (again MC particles)
416 // recoAll: reconstructed (again MC particles)
417 // particleType is: 0 for pion, 1 for kaon, 2 for proton, 3 for others
418
419 for (Int_t step=0; step<3; step++)
420 {
421 TObjArray* list = mc;
422 if (step == 1)
423 list = recoPrim;
424 else if (step == 2)
425 list = recoAll;
426
427 for (Int_t i=0; i<list->GetEntries(); i++)
428 {
429 AliVParticle* particle = (AliVParticle*) list->At(i);
a700ef3a 430 Double_t vars[3];
b1831bcb 431 vars[0] = particle->Eta();
432 vars[1] = particle->Pt();
433 vars[2] = particleType;
434
bf58cbde 435 for (Int_t j=0; j<fgkUEHists; j++)
436 if (GetUEHist(j))
437 GetUEHist(j)->GetTrackHistEfficiency()->Fill(vars, step);
b1831bcb 438 }
439 }
440}
441
a75aacd6 442//____________________________________________________________________
443void AliUEHistograms::FillEvent(Int_t eventType, Int_t step)
444{
445 // fills the number of events at the given step and the given enty type
446 //
447 // WARNING: This function is called from Fill, so only call it for steps where Fill is not called
448
449 fEventCount->Fill(step, eventType);
450}
451
452//____________________________________________________________________
453void AliUEHistograms::SetEtaRange(Float_t etaMin, Float_t etaMax)
454{
455 // sets eta min and max for all contained AliUEHist classes
456
bf58cbde 457 for (Int_t i=0; i<fgkUEHists; i++)
458 if (GetUEHist(i))
459 GetUEHist(i)->SetEtaRange(etaMin, etaMax);
a75aacd6 460}
461
462//____________________________________________________________________
463void AliUEHistograms::SetPtRange(Float_t ptMin, Float_t ptMax)
464{
465 // sets pT min and max for all contained AliUEHist classes
466
bf58cbde 467 for (Int_t i=0; i<fgkUEHists; i++)
468 if (GetUEHist(i))
469 GetUEHist(i)->SetPtRange(ptMin, ptMax);
a75aacd6 470}
471
144bd037 472//____________________________________________________________________
473void AliUEHistograms::SetContaminationEnhancement(TH1F* hist)
474{
475 // sets the contamination enhancement histogram in all contained AliUEHist classes
476
bf58cbde 477 for (Int_t i=0; i<fgkUEHists; i++)
478 if (GetUEHist(i))
479 GetUEHist(i)->SetContaminationEnhancement(hist);
144bd037 480}
481
a75aacd6 482//____________________________________________________________________
483void AliUEHistograms::SetCombineMinMax(Bool_t flag)
484{
485 // sets pT min and max for all contained AliUEHist classes
486
bf58cbde 487 for (Int_t i=0; i<fgkUEHists; i++)
488 if (GetUEHist(i))
489 GetUEHist(i)->SetCombineMinMax(flag);
a75aacd6 490}
491
492//____________________________________________________________________
493void AliUEHistograms::Correct(AliUEHistograms* corrections)
494{
495 // corrects the contained histograms by calling AliUEHist::Correct
496
bf58cbde 497 for (Int_t i=0; i<fgkUEHists; i++)
498 if (GetUEHist(i))
499 GetUEHist(i)->Correct(corrections->GetUEHist(i));
a75aacd6 500}
501
502//____________________________________________________________________
503AliUEHistograms &AliUEHistograms::operator=(const AliUEHistograms &c)
504{
505 // assigment operator
506
507 if (this != &c)
508 ((AliUEHistograms &) c).Copy(*this);
509
510 return *this;
511}
512
513//____________________________________________________________________
514void AliUEHistograms::Copy(TObject& c) const
515{
516 // copy function
517
518 AliUEHistograms& target = (AliUEHistograms &) c;
519
520 if (fNumberDensitypT)
521 target.fNumberDensitypT = dynamic_cast<AliUEHist*> (fNumberDensitypT->Clone());
522
523 if (fSumpT)
524 target.fSumpT = dynamic_cast<AliUEHist*> (fSumpT->Clone());
525
526 if (fNumberDensityPhi)
527 target.fNumberDensityPhi = dynamic_cast<AliUEHist*> (fNumberDensityPhi->Clone());
528
529 if (fCorrelationpT)
530 target.fCorrelationpT = dynamic_cast<TH2F*> (fCorrelationpT->Clone());
531
532 if (fCorrelationEta)
533 target.fCorrelationEta = dynamic_cast<TH2F*> (fCorrelationEta->Clone());
534
535 if (fCorrelationPhi)
536 target.fCorrelationPhi = dynamic_cast<TH2F*> (fCorrelationPhi->Clone());
537
538 if (fCorrelationR)
539 target.fCorrelationR = dynamic_cast<TH2F*> (fCorrelationR->Clone());
540
541 if (fCorrelationLeading2Phi)
542 target.fCorrelationLeading2Phi = dynamic_cast<TH2F*> (fCorrelationLeading2Phi->Clone());
543
544 if (fCorrelationMultiplicity)
545 target.fCorrelationMultiplicity = dynamic_cast<TH2F*> (fCorrelationMultiplicity->Clone());
546
547 if (fEventCount)
548 target.fEventCount = dynamic_cast<TH2F*> (fEventCount->Clone());
549
550 if (fEventCountDifferential)
551 target.fEventCountDifferential = dynamic_cast<TH3F*> (fEventCountDifferential->Clone());
552
553 if (fVertexContributors)
554 target.fVertexContributors = dynamic_cast<TH1F*> (fVertexContributors->Clone());
bf58cbde 555
556 if (fCentralityDistribution)
557 target.fCentralityDistribution = dynamic_cast<TH1F*> (fCentralityDistribution->Clone());
a75aacd6 558}
559
560//____________________________________________________________________
561Long64_t AliUEHistograms::Merge(TCollection* list)
562{
563 // Merge a list of AliUEHistograms objects with this (needed for
564 // PROOF).
565 // Returns the number of merged objects (including this).
566
567 if (!list)
568 return 0;
569
570 if (list->IsEmpty())
571 return 1;
572
573 TIterator* iter = list->MakeIterator();
574 TObject* obj;
575
576 // collections of objects
bf58cbde 577 const Int_t kMaxLists = 13;
a75aacd6 578 TList* lists[kMaxLists];
579
580 for (Int_t i=0; i<kMaxLists; i++)
581 lists[i] = new TList;
582
583 Int_t count = 0;
584 while ((obj = iter->Next())) {
585
586 AliUEHistograms* entry = dynamic_cast<AliUEHistograms*> (obj);
587 if (entry == 0)
588 continue;
589
bf58cbde 590 if (entry->fNumberDensitypT)
591 lists[0]->Add(entry->fNumberDensitypT);
592 if (entry->fSumpT)
593 lists[1]->Add(entry->fSumpT);
594 if (entry->fNumberDensityPhi)
595 lists[2]->Add(entry->fNumberDensityPhi);
a75aacd6 596 lists[3]->Add(entry->fCorrelationpT);
597 lists[4]->Add(entry->fCorrelationEta);
598 lists[5]->Add(entry->fCorrelationPhi);
599 lists[6]->Add(entry->fCorrelationR);
600 lists[7]->Add(entry->fCorrelationLeading2Phi);
601 lists[8]->Add(entry->fCorrelationMultiplicity);
602 lists[9]->Add(entry->fEventCount);
603 lists[10]->Add(entry->fEventCountDifferential);
604 lists[11]->Add(entry->fVertexContributors);
bf58cbde 605 lists[12]->Add(entry->fCentralityDistribution);
a75aacd6 606
607 count++;
608 }
609
bf58cbde 610 if (fNumberDensitypT)
611 fNumberDensitypT->Merge(lists[0]);
612 if (fSumpT)
613 fSumpT->Merge(lists[1]);
614 if (fNumberDensityPhi)
615 fNumberDensityPhi->Merge(lists[2]);
a75aacd6 616 fCorrelationpT->Merge(lists[3]);
617 fCorrelationEta->Merge(lists[4]);
618 fCorrelationPhi->Merge(lists[5]);
619 fCorrelationR->Merge(lists[6]);
620 fCorrelationLeading2Phi->Merge(lists[7]);
621 fCorrelationMultiplicity->Merge(lists[8]);
622 fEventCount->Merge(lists[9]);
623 fEventCountDifferential->Merge(lists[10]);
624 fVertexContributors->Merge(lists[11]);
bf58cbde 625 fCentralityDistribution->Merge(lists[12]);
a75aacd6 626
627 for (Int_t i=0; i<kMaxLists; i++)
628 delete lists[i];
629
630 return count+1;
631}
b1831bcb 632
633void AliUEHistograms::CopyReconstructedData(AliUEHistograms* from)
634{
635 // copies those histograms extracted from ESD to this object
636
bf58cbde 637 for (Int_t i=0; i<fgkUEHists; i++)
638 if (GetUEHist(i))
639 GetUEHist(i)->CopyReconstructedData(from->GetUEHist(i));
b1831bcb 640}
6f803f6c 641
642void AliUEHistograms::ExtendTrackingEfficiency()
643{
644 // delegates to AliUEHists
645
bf58cbde 646 for (Int_t i=0; i<fgkUEHists; i++)
647 if (GetUEHist(i))
648 GetUEHist(i)->ExtendTrackingEfficiency();
6f803f6c 649}
650