Escape special characters for latex (Diego)
[u/mrichter/AliRoot.git] / PWGPP / Centrality / AliAnalysisTaskHIMultCorr.cxx
CommitLineData
44a84424 1//-*- Mode: C++ -*-
2
3#if __GNUS__ >= 3
4using namespace std;
5#endif
6
7#include "TFile.h"
8#include "TChain.h"
9#include "TTree.h"
10#include "TH1F.h"
11#include "TH2F.h"
12#include "TCanvas.h"
13#include "TStopwatch.h"
14#include "TMath.h"
15
16#include "AliAnalysisTask.h"
17#include "AliAnalysisManager.h"
18#include "AliTracker.h"
19#include "AliESDEvent.h"
20#include "AliESDInputHandler.h"
21#include "AliStack.h"
22#include "AliMCEvent.h"
23#include "AliMCEventHandler.h"
24#include "AliESDtrackCuts.h"
25#include "AliKineTrackCuts.h"
26#include "AliMCParticle.h"
27#include "AliESDVZERO.h"
28#include "AliAnalysisTaskHIMultCorr.h"
5da23c19 29#include "AliCentrality.h"
44a84424 30
31#include "AliESDVZERO.h"
32#include "AliTriggerAnalysis.h"
33#include "AliMultiplicity.h"
34
35#include "AliMultiplicityCorrelations.h"
36
37// Task for HI Multiplicity correlation checks
38// Author: Jochen Thaeder <jochen@thaeder.de>
39
40ClassImp(AliAnalysisTaskHIMultCorr)
41
42#define USE_STREAMER 0
43
44/*
45 * ---------------------------------------------------------------------------------
46 * Constructor / Destructor
47 * ---------------------------------------------------------------------------------
48 */
49
50//________________________________________________________________________
51AliAnalysisTaskHIMultCorr::AliAnalysisTaskHIMultCorr(const char *name) :
52 AliAnalysisTaskSE(name),
53 fpcstream(NULL),
54 fIsMC(kFALSE),
55 fHStat(NULL),
56 fOutList(NULL),
57 fESD(NULL), fESDTrackCuts(NULL), fESDTrackCuts2(NULL),
58 fUseCentralitySel(0),
59 fCentralityBin(-1),
60 fCentralitySPDBin(-1),
61 fCentralityVZEROBin(-1),
62 fCentralitySPD(-1.),
63 fCentralityVZERO(-1.),
64 fMaxVertexZ(30.),
65 fTriggerAnalysis(NULL),
66 fCorrObj(NULL),
67 fCorrObjCent0(NULL), fCorrObjCent1(NULL), fCorrObjCent2(NULL) {
68 // Constructor
69
70#if USE_STREAMER
71 fpcstream = new TTreeSRedirector("eventInfoCorr.root");
72#endif
73
74 // Output slot #1 writes into a TList container
75 DefineOutput(1, TList::Class());
76}
77
78//________________________________________________________________________
79AliAnalysisTaskHIMultCorr::~AliAnalysisTaskHIMultCorr() {
80 // Destructor
81
82 if (fTriggerAnalysis)
83 delete fTriggerAnalysis;
84 fTriggerAnalysis = NULL;
85
86#if USE_STREAMER
87 if (fpcstream)
88 delete fpcstream;
89 fpcstream = NULL;
90#endif
91
92}
93
94/*
95 * ---------------------------------------------------------------------------------
96 * Methods
97 * ---------------------------------------------------------------------------------
98 */
99
100//________________________________________________________________________
101void AliAnalysisTaskHIMultCorr::UserCreateOutputObjects() {
102 // Create histograms
103 // Called once
104
105 Bool_t oldStatus = TH1::AddDirectoryStatus();
106 TH1::AddDirectory(kFALSE);
107
108 fOutList = new TList;
109 fOutList->SetName(GetName()) ;
110 fOutList->SetOwner(kTRUE) ;
111
112 fTriggerAnalysis = new AliTriggerAnalysis;
113
114 // ------------------------------------------------------------------
115 // -- Correlations
116 // ------------------------------------------------------------------
117
118 fCorrObj = new AliMultiplicityCorrelations;
119 fCorrObj->Initialize("All");
120 fOutList->Add(fCorrObj->GetHistList());
121
122 fCorrObjCent0 = new AliMultiplicityCorrelations;
123 fCorrObjCent0->Initialize("Cent_0-5");
124 // fCorrObjCent0->SetCleanSample(1700.,2900.); // Cleaning of Sample
125 fOutList->Add(fCorrObjCent0->GetHistList());
126
127 fCorrObjCent1 = new AliMultiplicityCorrelations;
128 fCorrObjCent1->Initialize("Cent_80-90");
129 fOutList->Add(fCorrObjCent1->GetHistList());
130
131 fCorrObjCent2 = new AliMultiplicityCorrelations;
132 fCorrObjCent2->Initialize("Cent_70-80");
133 // fCorrObjCent2->SetCleanSample(0.,180.); // Cleaning of Sample
134 fOutList->Add(fCorrObjCent2->GetHistList());
135
136 if (fIsMC) {
0c81f7b3 137 fCorrObj->SetIsMC();
138 fCorrObjCent0->SetIsMC();
139 fCorrObjCent1->SetIsMC();
140 fCorrObjCent2->SetIsMC();
44a84424 141 }
142
143 if (fESDTrackCuts) {
0c81f7b3 144 fCorrObj->SetESDTrackCuts(fESDTrackCuts);
145 fCorrObjCent0->SetESDTrackCuts(fESDTrackCuts);
146 fCorrObjCent1->SetESDTrackCuts(fESDTrackCuts);
147 fCorrObjCent2->SetESDTrackCuts(fESDTrackCuts);
44a84424 148 }
149 else
150 AliError("No ESD trackCuts!!");
151
152 if (fESDTrackCuts2) {
0c81f7b3 153 fCorrObj->SetESDTrackCuts2(fESDTrackCuts2);
154 fCorrObjCent0->SetESDTrackCuts2(fESDTrackCuts2);
155 fCorrObjCent1->SetESDTrackCuts2(fESDTrackCuts2);
156 fCorrObjCent2->SetESDTrackCuts2(fESDTrackCuts2);
44a84424 157 }
158
159 // ------------------------------------------------------------------
160 // -- Centrality Distributions
161 // ------------------------------------------------------------------
162
163 fOutList->Add(new TH1F("fVZEROCentrality", "VZERO Centrality Distribution;Centrality;N_{Events}", 102,0.,101.));
164 fOutList->Add(new TH1F("fVZEROCentralityBin","VZERO Centrality Bin Distribution;Centrality;N_{Events}", 102,0.,101.));
165
166 fOutList->Add(new TH1F("fSPDCentrality", "SPD Centrality Distribution;Centrality;N_{Events}", 102,0.,101.));
167 fOutList->Add(new TH1F("fSPDCentralityBin","SPD Centrality Bin Distribution;Centrality;N_{Events}", 102,0.,101.));
168
169 fOutList->Add(new TH2F("fCorrCentrality",
170 "Centrality - SPD vs VZERO;Centrality_{SPD};Centrality_{VZERO}", 102,0.,101., 102,0.,101.));
171 fOutList->Add(new TH2F("fCorrCentralityBin",
172 "Centrality Bin - SPD vs VZERO;Centrality_{SPD};Centrality_{VZERO}", 102,0.,101., 102,0.,101.));
173
174 // ------------------------------------------------------------------
175 // -- Vertex : TPC - Tracks
176 // ------------------------------------------------------------------
177
178 fOutList->Add(new TH2F("fCorrVertexTracks",
179 "Tracks Vertex - Status vs Contributers;Status;N_{contributers}", 2,-1.,3., 500,0.,3000.));
180 fOutList->Add(new TH2F("fCrossCorrVertexTracks",
181 "Tracks Vertex - TPC Status vs Contributers;Status;N_{contributers}", 2,-1.,3., 500,0.,3000.));
182
183 fOutList->Add(new TH2F("fCorrVertexTPC",
184 "TPC Vertex - Status vs Contributers;Status;N_{contributers}", 2,-1.,3., 500,0.,3000.));
185 fOutList->Add(new TH2F("fCrossCorrVertexTPC",
186 "TPC Vertex - Tracks Status vs Contributers;Status;N_{contributers}", 2,-1.,3., 500,0.,3000.));
187
188 fOutList->Add(new TH2F("fCorrVertexNCont",
189 "N Contributers - Vertex Tracks vs Vertex TPC;N_{contributers}^{Tracks};N_{contributers}^{TPC}",
190 500,0.,3000., 500,0.,3000.));
191
192 fOutList->Add(new TH1F("fDeltaVx","Vertex Tracks - Vertex TPC #Deltax;#Delta X [cm]", 201,-10.,10.));
193 fOutList->Add(new TH1F("fDeltaVy","Vertex Tracks - Vertex TPC #Deltay;#Delta Y [cm]", 201,-10.,10.));
194 fOutList->Add(new TH1F("fDeltaVz","Vertex Tracks - Vertex TPC #Deltaz;#Delta z [cm]", 201,-10.,10.));
195
196 // ------------------------------------------------------------------
197 // -- Vertex Z : TPC - Tracks - SPD
198 // ------------------------------------------------------------------
199
200 fOutList->Add(new TH2F("fVzTPCvsSPD",
201 "Z - Vertex TPC vs Vertex SPD;Vz_{TPC} [cm];Vz_{SPD} [cm]", 401,-20.,20., 401,-20.,20.));
202 fOutList->Add(new TH2F("fVzTPCvsTracks",
203 "Z - Vertex TPC vs Vertex Tracks;Vz_{TPC} [cm];Vz_{Tracks} [cm]", 401,-20.,20., 401,-20.,20.));
204 fOutList->Add(new TH2F("fVzSPDvsTracks",
205 "Z - Vertex SPD vs Vertex Tracks;Vz_{SPD} [cm];Vz_{Tracks} [cm]", 401,-20.,20., 401,-20.,20.));
206
207 // ------------------------------------------------------------------
208 // -- Cuts
209 // ------------------------------------------------------------------
210
211 fOutList->Add(new TH1F("fHStat","Cut statistics", 20,0.,19));
212 fHStat = dynamic_cast<TH1F*>(fOutList->Last());
213
214 TH1::AddDirectory(oldStatus);
215}
216
217//________________________________________________________________________
218void AliAnalysisTaskHIMultCorr::UserExec(Option_t *) {
219 // Called for each event
220
221#if USE_STREAMER
222 Int_t i1 = -1;
223 Int_t i2 = -1;
224 Float_t f1 = -1.;
225 Float_t f2 = -1.;
226 Float_t f3 = -1.;
227#endif
228
229 // -- Setup Event
230 // ----------------
231 if ( !SetupEvent() ) {
232#if USE_STREAMER
233 (*fpcstream) <<"event"<<
234 "nTracksTPC=" << i1 <<
235 "nTracks=" << i2 <<
236 "v0A=" << f1 <<
237 "v0C=" << f2 <<
238 "v0Corr=" << f3 << "\n";
239#endif
240 return;
241 }
242
243 // -- Centrality dependence
244 // --------------------------
245 if (fCentralityBin < 0) {
246#if USE_STREAMER
247 (*fpcstream) <<"event"<<
248 "nTracksTPC=" << i1 <<
249 "nTracks=" << i2 <<
250 "v0A=" << f1 <<
251 "v0C=" << f2 <<
252 "v0Corr=" << f3 << "\n";
253#endif
254 return;
255 }
256
257 // -- Process Event
258 // ------------------
259 Int_t iResult = 0;
260 if (fCorrObj) iResult = fCorrObj->ProcessEvent(fESD);
261 if (iResult == -3) {
262#if USE_STREAMER
263 (*fpcstream) <<"event"<<
264 "nTracksTPC=" << i1 <<
265 "nTracks=" << i2 <<
266 "v0A=" << f1 <<
267 "v0C=" << f2 <<
268 "v0Corr=" << f3 << "\n";
269#endif
270 return;
271 }
272
273 // -- After "Jochen's cut"
274 fHStat->Fill(7);
275 fHStat->Fill(17);
276
277 if (fCentralityBin == 0 && fCorrObjCent0) fCorrObjCent0->ProcessEvent(fESD);
278 else if (fCentralityBin == 80 && fCorrObjCent1) fCorrObjCent1->ProcessEvent(fESD);
279 else if (fCentralityBin == 70 && fCorrObjCent2) fCorrObjCent2->ProcessEvent(fESD);
280
281 // -- Marian's Streamer
282 // ----------------------
283#if USE_STREAMER
284 if (fCorrObj) {
285 i1 = fCorrObj->GetNTracksTPC();
286 i2 = fCorrObj->GetNTracks();
287 f1 = fCorrObj->GetVZEROA();
288 f2 = fCorrObj->GetVZEROC();
289 f3 = fCorrObj->GetVZEROCorr();
290 }
291
292 (*fpcstream) <<"event"<<
293 "nTracksTPC=" << i1 <<
294 "nTracks=" << i2 <<
295 "v0A=" << f1 <<
296 "v0C=" << f2 <<
297 "v0Corr=" << f3 << "\n";
298#endif
299
300 // -- Post output data
301 // ---------------------
302 PostData(1,fOutList);
303
304 return;
305}
306
307/*
308 * ---------------------------------------------------------------------------------
309 * Setup Methods - private
310 * ---------------------------------------------------------------------------------
311 */
312
313//________________________________________________________________________
314Bool_t AliAnalysisTaskHIMultCorr::SetupEvent() {
315 // Setup Reading of event
316
317 Bool_t aCuts[] = {kFALSE, kFALSE, kFALSE, kFALSE, kFALSE, kFALSE, kFALSE};
318
319 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
320 // -- ESD Event
321 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
322
323 AliESDInputHandler *esdH = dynamic_cast<AliESDInputHandler*>
324 (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
325 if (!esdH) {
326 AliError("Could not get ESDInputHandler");
327 return kFALSE;
328 }
329
330 fESD = esdH->GetEvent();
331 if (!fESD) {
332 AliError("fESD not available");
333 return kFALSE;
334 }
335
336 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
337 // -- Check object existence
338 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
339
340 const AliESDVertex* vtxESD = fESD->GetPrimaryVertexTracks();
341 const AliESDVertex* vtxESDTPC = fESD->GetPrimaryVertexTPC();
342 const AliESDVertex* vtxESDSPD = fESD->GetPrimaryVertexSPD();
343 const AliMultiplicity* multESD = fESD->GetMultiplicity();
344
345 if ( !vtxESD ){
346 AliError("No Tracks Vertex");
347 return kFALSE;
348 }
349
350 if ( !vtxESDTPC ){
351 AliError("No TPC Vertex");
352 return kFALSE;
353 }
354
355 if ( !vtxESDSPD ){
356 AliError("No SPD Vertex");
357 return kFALSE;
358 }
359
360 if ( !multESD ){
361 AliError("No Multiplicity");
362 return kFALSE;
363 }
364
365 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
366 // -- After Physics Selection;
367 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
368
369 aCuts[0] = kFALSE;
370
371 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
372 // -- ZDC cut
373 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
374
375 if ( !fIsMC ) {
376 if (!fTriggerAnalysis->ZDCTimeTrigger(fESD))
377 aCuts[1] = kTRUE;
378 }
379
380 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
381 // -- Francesco Prino cut's
382 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
383
384 if (!vtxESDTPC->GetStatus())
385 aCuts[2] = kTRUE;
386
387 if( vtxESDTPC->GetNContributors() < (-10.+0.25*multESD->GetNumberOfITSClusters(0)) )
388 aCuts[3] = kTRUE;
389
390 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
391 // -- Vertex Cuts - Tracks
392 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
393
394 // if (!vtxESD->GetStatus())
395 // aCuts[4] = kTRUE;
396
397 if(vtxESDTPC->GetZv() > fMaxVertexZ || vtxESDTPC->GetZv() < (-1.*fMaxVertexZ))
398 aCuts[5] = kTRUE;
399
400 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
401 // -- CentralityBin
402 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
403
404 if (GetCentralityBin() < 0)
405 aCuts[6] = kTRUE;
406
407 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
408 // -- Marian's debug streamer
409 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
410#if USE_STREAMER
411 // TFile *inFile=(dynamic_cast<TChain*>(GetInputData(0)))->GetCurrentFile();
412 // TString fileName(inFile->GetName());
413 static Int_t runNo = fESD->GetRunNumber();
414 static UInt_t tt = fESD->GetTimeStamp();
415 static Float_t zvTPC = vtxESDTPC->GetZv();
416 static Float_t zvSPD = vtxESDSPD->GetZv();
417 static Float_t zvTra = vtxESD->GetZv();
418 static Int_t ncontTPCV = vtxESDTPC->GetNContributors();
419 static Int_t ncontSPDV = vtxESDSPD->GetNContributors();
420 static Int_t ncontTraV = vtxESD->GetNContributors();
421 static Bool_t cutPrino = aCuts[3];
422 static Bool_t cutZDCTi = aCuts[1];
423 static Int_t spd1 = multESD->GetNumberOfITSClusters(0);
424 static Int_t spd2 = multESD->GetNumberOfITSClusters(1);
425
426 (*fpcstream) << "event" <<
427 "run=" << runNo << //runNo
428 "time=" << tt <<
429 // "fname=" << fileName <<
430 // "eventNr=" << eventNumber<<
431 "zvTPC=" << zvTPC << //zvertex TPC
432 "zvSPD=" << zvSPD << //zvertex SPD
433 "zvTra=" << zvTra << //zvertex Tracks
434 "ncontTPCV"<< ncontTPCV << // N contributors to TPC vtx
435 "ncontSPDV"<< ncontSPDV << // N contributors to SPD vtx
436 "ncontTraV"<< ncontTraV << // N contributors to Tracks vtx
437 "cutPrino="<< cutPrino << //francecos cut
438 "cutZDCTi="<< cutZDCTi << //zdctiming cut
439 "spd1=" << spd1 << // N clus in SPD0
440 "spd2=" << spd2 << // N clus in SPD1
441 "centSPD=" << fCentralitySPD <<
442 "centSPDB="<< fCentralitySPDBin <<
443 "centV0=" << fCentralityVZERO <<
444 "centV0B=" << fCentralityVZEROBin;
445#endif
446 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
447 // -- Fill statistics / reject event
448 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
449
450 Bool_t isRejected = kFALSE;
451
452 for (Int_t idx = 0; idx < 7; ++idx) {
453 if (aCuts[idx])
454 isRejected = kTRUE;
455 else
456 fHStat->Fill(idx);
457 }
458
459 for (Int_t idx = 0; idx < 7; ++idx) {
460 if (aCuts[idx])
461 break;
462 fHStat->Fill(10+idx);
463 }
464
465 if (isRejected)
466 return kFALSE;
467
468 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
469 // -- Fill Centrality Histograms
470 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
471
472 (static_cast<TH1F*>(fOutList->FindObject("fVZEROCentrality")))->Fill(fCentralityVZERO);
473 (static_cast<TH1F*>(fOutList->FindObject("fVZEROCentralityBin")))->Fill(fCentralityVZEROBin);
474
475 (static_cast<TH1F*>(fOutList->FindObject("fSPDCentrality")))->Fill(fCentralitySPD);
476 (static_cast<TH1F*>(fOutList->FindObject("fSPDCentralityBin")))->Fill(fCentralitySPDBin);
477
478 (static_cast<TH2F*>(fOutList->FindObject("fCorrCentrality")))->Fill(fCentralitySPD,fCentralityVZERO);
479 (static_cast<TH2F*>(fOutList->FindObject("fCorrCentralityBin")))->Fill(fCentralitySPDBin,fCentralityVZEROBin);
480
481 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
482 // -- Fill Vertex Histograms
483 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
484
485 (static_cast<TH2F*>(fOutList->FindObject("fCorrVertexTracks")))->Fill(vtxESD->GetStatus(),vtxESD->GetNContributors());
486 (static_cast<TH2F*>(fOutList->FindObject("fCorrVertexTPC")))->Fill(vtxESDTPC->GetStatus(),vtxESDTPC->GetNContributors());
487 (static_cast<TH2F*>(fOutList->FindObject("fCrossCorrVertexTracks")))->Fill(vtxESDTPC->GetStatus(),vtxESD->GetNContributors());
488 (static_cast<TH2F*>(fOutList->FindObject("fCrossCorrVertexTPC")))->Fill(vtxESD->GetStatus(),vtxESDTPC->GetNContributors());
489 (static_cast<TH2F*>(fOutList->FindObject("fCorrVertexNCont")))->Fill(vtxESD->GetNContributors(),vtxESDTPC->GetNContributors());
490
491 (static_cast<TH1F*>(fOutList->FindObject("fDeltaVx")))->Fill(vtxESD->GetXv() - vtxESDTPC->GetXv());
492 (static_cast<TH1F*>(fOutList->FindObject("fDeltaVy")))->Fill(vtxESD->GetYv() - vtxESDTPC->GetYv());
493 (static_cast<TH1F*>(fOutList->FindObject("fDeltaVz")))->Fill(vtxESD->GetZv() - vtxESDTPC->GetZv());
494
495 (static_cast<TH2F*>(fOutList->FindObject("fVzTPCvsSPD")))->Fill(vtxESDTPC->GetZv(),vtxESDSPD->GetZv());
496 (static_cast<TH2F*>(fOutList->FindObject("fVzTPCvsTracks")))->Fill(vtxESDTPC->GetZv(),vtxESD->GetZv());
497 (static_cast<TH2F*>(fOutList->FindObject("fVzSPDvsTracks")))->Fill(vtxESDSPD->GetZv(),vtxESD->GetZv());
498
499 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
500
501 return kTRUE;
502}
503
504//________________________________________________________________________
505Int_t AliAnalysisTaskHIMultCorr::GetCentralityBin(){
506 // Get Centrality bin
507
508 fCentralityBin = -1;
509
510 if (fUseCentralitySel == 0)
511 return fCentralityBin;
512
5da23c19 513 AliCentrality *esdCentrality = fESD->GetCentrality();
44a84424 514 fCentralityVZERO = esdCentrality->GetCentralityPercentile("V0M");
515 fCentralitySPD = esdCentrality->GetCentralityPercentile("CL1");
516
517 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
518
519 fCentralityVZEROBin = esdCentrality->GetCentralityClass10("V0M")*10;
520 if (fCentralityVZEROBin == 0 && esdCentrality->GetCentralityClass5("V0M") == 1)
521 fCentralityVZEROBin = 5;
522
523 fCentralitySPDBin = esdCentrality->GetCentralityClass10("CL1")*10;
524 if (fCentralitySPDBin == 0 && esdCentrality->GetCentralityClass5("CL1") == 1)
525 fCentralitySPDBin = 5;
526
527 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
528
529 if ( fUseCentralitySel == 1 )
530 fCentralityBin = fCentralityVZEROBin;
531 else if ( fUseCentralitySel == 2 )
532 fCentralityBin = fCentralitySPDBin;
533
534 if ( fCentralityBin == 100 )
535 fCentralityBin = -1;
536
537 // -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --
538
539 return fCentralityBin;
540}
541
542//________________________________________________________________________
543void AliAnalysisTaskHIMultCorr::Terminate(Option_t *){
544 // Terminate
545#if USE_STREAMER
546 if (fpcstream) fpcstream->GetFile()->Write();
547#endif
548
549}