ANALYSIS/AliCentralitySelectionTask.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-2008, 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 //*****************************************************
  17 //   Class AliCentralitySelectionTask
  18 //   Class to analyze determine centrality
  19 //   author: Alberica Toia
  20 //*****************************************************
  21
  22 #include "AliCentralitySelectionTask.h"
  23
  24 #include <TTree.h>
  25 #include <TList.h>
  26 #include <TH1F.h>
  27 #include <TH2F.h>
  28 #include <TF1.h>
  29 #include <TProfile.h>
  30 #include <TFile.h>
  31 #include <TObjString.h>
  32 #include <TString.h>
  33 #include <TCanvas.h>
  34 #include <TROOT.h>
  35 #include <TDirectory.h>
  36 #include <TSystem.h>
  37 #include <iostream>
  38
  39 #include "AliAnalysisManager.h"
  40 #include "AliHeader.h"
  41 #include "AliVEvent.h"
  42 #include "AliESD.h"
  43 #include "AliESDEvent.h"
  44 #include "AliESDHeader.h"
  45 #include "AliESDInputHandler.h"
  46 #include "AliESDZDC.h"
  47 #include "AliESDFMD.h"
  48 #include "AliESDTZERO.h"
  49 #include "AliESDVZERO.h"
  50 #include "AliESDtrackCuts.h"
  51 #include "AliESDVertex.h"
  52 #include "AliCentrality.h"
  53 #include "AliOADBCentrality.h"
  54 #include "AliOADBContainer.h"
  55 #include "AliMultiplicity.h"
  56 #include "AliAODHandler.h"
  57 #include "AliAODHeader.h"
  58 #include "AliAODEvent.h"
  59 #include "AliAODVertex.h"
  60 #include "AliAODVZERO.h"
  61 #include "AliAODTracklets.h"
  62 #include "AliAODMCHeader.h"
  63 #include "AliMCEventHandler.h"
  64 #include "AliMCEvent.h"
  65 #include "AliAODMCParticle.h"
  66 #include "AliMCParticle.h"
  67 #include "AliStack.h"
  68 #include "AliAnalysisTaskSE.h"
  69 #include "AliGenEventHeader.h"
  70 #include "AliGenHijingEventHeader.h"
  71 #include "AliGenDPMjetEventHeader.h"
  72 #include "AliGenCocktailEventHeader.h"
  73 #include "AliPhysicsSelectionTask.h"
  74 #include "AliPhysicsSelection.h"
  75 #include "AliBackgroundSelection.h"
  76 #include "AliESDUtils.h"
  77
  78 ClassImp(AliCentralitySelectionTask)
  79
  80
  81 //________________________________________________________________________
  82 AliCentralitySelectionTask::AliCentralitySelectionTask():
  83 AliAnalysisTaskSE(),
  84   fAnalysisInput("ESD"),
  85   fIsMCInput(kFALSE),
  86   fCurrentRun(-1),
  87   fUseScaling(0),
  88   fUseCleaning(0),
  89   fFillHistos(0),
  90   fV0MScaleFactor(0),
  91   fSPDScaleFactor(0),
  92   fTPCScaleFactor(0),
  93   fV0MScaleFactorMC(0),
  94   fV0MSPDOutlierPar0(0),
  95   fV0MSPDOutlierPar1(0),
  96   fV0MTPCOutlierPar0(0),
  97   fV0MTPCOutlierPar1(0),
  98   fV0MSPDSigmaOutlierPar0(0),
  99   fV0MSPDSigmaOutlierPar1(0),
 100   fV0MSPDSigmaOutlierPar2(0),
 101   fV0MTPCSigmaOutlierPar0(0),
 102   fV0MTPCSigmaOutlierPar1(0),
 103   fV0MTPCSigmaOutlierPar2(0),
 104   fV0MZDCOutlierPar0(0),
 105   fV0MZDCOutlierPar1(0),
 106   fV0MZDCEcalOutlierPar0(0),
 107   fV0MZDCEcalOutlierPar1(0),
 108   fTrackCuts(0),
 109   fEsdTrackCuts(0),
 110   fEsdTrackCutsExtra1(0),
 111   fEsdTrackCutsExtra2(0),
 112   fZVCut(10),
 113   fOutliersCut(5),
 114   fQuality(999),
 115   fIsSelected(0),
 116   fMSL(0),
 117   fMSH(0),
 118   fMUL(0),
 119   fMLL(0),
 120   fEJE(0),
 121   fEGA(0),
 122   fPHS(0),
 123   fMB(0),
 124   fCVHN(0),
 125   fCVLN(0),
 126   fCVHNbit(0),
 127   fCVLNbit(0),
 128   fCCENT(0),
 129   fCSEMI(0),
 130   fCCENTbit(0),
 131   fCSEMIbit(0),
 132   fCentV0M(0),
 133   fCentV0A(0),
 134   fCentV0A0(0),
 135   fCentV0A123(0),
 136   fCentV0C(0),
 137   fCentV0MEq(0),
 138   fCentV0AEq(0),
 139   fCentV0CEq(0),
 140   fCentFMD(0),
 141   fCentTRK(0),
 142   fCentTKL(0),
 143   fCentCL0(0),
 144   fCentCL1(0),
 145   fCentCND(0),
 146   fCentZNA(0),
 147   fCentZNC(0),
 148   fCentZPA(0),
 149   fCentZPC(0),
 150   fCentNPA(0),
 151   fCentV0MvsFMD(0),
 152   fCentTKLvsV0M(0),
 153   fCentZEMvsZDC(0),
 154   fCentV0Mtrue(0),
 155   fCentV0Atrue(0),
 156   fCentV0Ctrue(0),
 157   fCentV0MEqtrue(0),
 158   fCentV0AEqtrue(0),
 159   fCentV0CEqtrue(0),
 160   fCentFMDtrue(0),
 161   fCentTRKtrue(0),
 162   fCentTKLtrue(0),
 163   fCentCL0true(0),
 164   fCentCL1true(0),
 165   fCentCNDtrue(0),
 166   fCentZNAtrue(0),
 167   fCentZNCtrue(0),
 168   fCentZPAtrue(0),
 169   fCentZPCtrue(0),
 170   fHtempV0M(0),
 171   fHtempV0A(0),
 172   fHtempV0A0(0),
 173   fHtempV0A123(0),
 174   fHtempV0C(0),
 175   fHtempV0MEq(0),
 176   fHtempV0AEq(0),
 177   fHtempV0CEq(0),
 178   fHtempFMD(0),
 179   fHtempTRK(0),
 180   fHtempTKL(0),
 181   fHtempCL0(0),
 182   fHtempCL1(0),
 183   fHtempCND(0),
 184   fHtempZNA(0),
 185   fHtempZNC(0),
 186   fHtempZPA(0),
 187   fHtempZPC(0),
 188   fHtempV0MvsFMD(0),
 189   fHtempTKLvsV0M(0),
 190   fHtempZEMvsZDC(0),
 191   fHtempNPA(0),
 192   fHtempV0Mtrue(0),
 193   fHtempV0Atrue(0),
 194   fHtempV0Ctrue(0),
 195   fHtempV0MEqtrue(0),
 196   fHtempV0AEqtrue(0),
 197   fHtempV0CEqtrue(0),
 198   fHtempFMDtrue(0),
 199   fHtempTRKtrue(0),
 200   fHtempTKLtrue(0),
 201   fHtempCL0true(0),
 202   fHtempCL1true(0),
 203   fHtempCNDtrue(0),
 204   fHtempZNAtrue(0),
 205   fHtempZNCtrue(0),
 206   fHtempZPAtrue(0),
 207   fHtempZPCtrue(0),
 208   fOutputList(0),
 209   fHOutCentV0M(0),
 210   fHOutCentV0A(0),
 211   fHOutCentV0A0(0),
 212   fHOutCentV0A123(0),
 213   fHOutCentV0C(0),
 214   fHOutCentV0MEq(0),
 215   fHOutCentV0AEq(0),
 216   fHOutCentV0CEq(0),
 217   fHOutCentV0MCVHN(0),
 218   fHOutCentV0MCVLN(0),
 219   fHOutCentV0MCVHNinMB(0),
 220   fHOutCentV0MCVLNinMB(0),
 221   fHOutCentV0MCCENT(0),
 222   fHOutCentV0MCSEMI(0),
 223   fHOutCentV0MCCENTinMB(0),
 224   fHOutCentV0MCSEMIinMB(0),
 225   fHOutCentV0MMSL(0),
 226   fHOutCentV0MMSH(0),
 227   fHOutCentV0MMUL(0),
 228   fHOutCentV0MMLL(0),
 229   fHOutCentV0MEJE(0),
 230   fHOutCentV0MEGA(0),
 231   fHOutCentV0MPHS(0),
 232   fHOutCentV0MMSLinMB(0),
 233   fHOutCentV0MMSHinMB(0),
 234   fHOutCentV0MMULinMB(0),
 235   fHOutCentV0MMLLinMB(0),
 236   fHOutCentV0MEJEinMB(0),
 237   fHOutCentV0MEGAinMB(0),
 238   fHOutCentV0MPHSinMB(0),
 239   fHOutCentFMD(0),
 240   fHOutCentTRK(0),
 241   fHOutCentTKL(0),
 242   fHOutCentCL0(0),
 243   fHOutCentCL1(0),
 244   fHOutCentCND(0),
 245   fHOutCentNPA(0),
 246   fHOutCentZNA(0),
 247   fHOutCentZNC(0),
 248   fHOutCentZPA(0),
 249   fHOutCentZPC(0),
 250   fHOutCentV0MvsFMD(0),
 251   fHOutCentTKLvsV0M(0),
 252   fHOutCentZEMvsZDC(0),
 253   fHOutCentV0MvsCentCL1(0),
 254   fHOutCentV0MvsCentTRK(0),
 255   fHOutCentTRKvsCentCL1(0),
 256   fHOutCentV0MvsCentZDC(0),
 257   fHOutCentV0AvsCentV0C(0),
 258   fHOutCentV0AvsCentTRK(0),
 259   fHOutCentV0AvsCentCND(0),
 260   fHOutCentV0AvsCentCL1(0),
 261   fHOutCentV0CvsCentTRK(0),
 262   fHOutCentV0CvsCentCND(0),
 263   fHOutCentV0CvsCentCL1(0),
 264   fHOutCentNPAvsCentV0A(0),
 265   fHOutCentNPAvsCentV0C(0),
 266   fHOutCentNPAvsCentTRK(0),
 267   fHOutCentNPAvsCentCND(0),
 268   fHOutCentNPAvsCentCL1(0),
 269   fHOutCentZNAvsCentV0A(0),
 270   fHOutCentZNAvsCentV0C(0),
 271   fHOutCentZNAvsCentTRK(0),
 272   fHOutCentZNAvsCentCND(0),
 273   fHOutCentZNAvsCentCL1(0),
 274   fHOutCentZNAvsCentZPA(0),
 275   fHOutMultV0AC(0),
 276   fHOutMultV0M(0),
 277   fHOutMultV0A(0),
 278   fHOutMultV0A0(0),
 279   fHOutMultV0A123(0),
 280   fHOutMultV0C(0),
 281   fHOutMultV0MEq(0),
 282   fHOutMultV0AEq(0),
 283   fHOutMultV0CEq(0),
 284   fHOutMultV0Mnc(0),
 285   fHOutMultV0Anc(0),
 286   fHOutMultV0Cnc(0),
 287   fHOutMultV0O(0),
 288   fHOutMultV0Cells(0),
 289   fHOutMultFMD(0),
 290   fHOutMultTRK(0),
 291   fHOutMultTKL(0),
 292   fHOutMultCL0(0),
 293   fHOutMultCL1(0),
 294   fHOutMultCND(0),
 295   fHOutMultNPA(0),
 296   fHOutMultZNA(0),
 297   fHOutMultZNC(0),
 298   fHOutMultZPA(0),
 299   fHOutMultZPC(0),
 300   fHOutMultV0MvsZDN(0),
 301   fHOutMultZEMvsZDN(0),
 302   fHOutMultV0MvsZDC(0),
 303   fHOutMultZEMvsZDC(0),
 304   fHOutMultZEMvsZDCw(0),
 305   fHOutMultV0MvsCL1(0),
 306   fHOutMultV0MvsTRK(0),
 307   fHOutMultTRKvsCL1(0),
 308   fHOutMultV0MvsV0O(0),
 309   fHOutMultV0OvsCL1(0),
 310   fHOutMultV0OvsTRK(0),
 311   fHOutMultCL1vsTKL(0),
 312   fHOutMultZNAvsZPA(0),
 313   fHOutCentV0Mqual1(0),
 314   fHOutCentTRKqual1(0),
 315   fHOutCentCL1qual1(0),
 316   fHOutMultV0MvsCL1qual1(0),
 317   fHOutMultV0MvsTRKqual1(0),
 318   fHOutMultTRKvsCL1qual1(0),
 319   fHOutCentV0Mqual2(0),
 320   fHOutCentTRKqual2(0),
 321   fHOutCentCL1qual2(0),
 322   fHOutMultV0MvsCL1qual2(0),
 323   fHOutMultV0MvsTRKqual2(0),
 324   fHOutMultTRKvsCL1qual2(0),
 325   fHOutQuality(0),
 326   fHOutVertex(0),
 327   fHOutVertexT0(0)
 328 {
 329   // Default constructor
 330   AliInfo("Centrality Selection enabled.");
 331
 332   fUseScaling=kTRUE;
 333   fUseCleaning=kTRUE;
 334   fFillHistos=kFALSE;
 335   fBranchNames="ESD:AliESDRun.,AliESDHeader.,AliESDZDC.,AliESDFMD.,AliESDVZERO.,AliESDTZERO."
 336     ",SPDVertex.,TPCVertex.,PrimaryVertex.,AliMultiplicity.,Tracks ";
 337 }
 338
 339 //________________________________________________________________________
 340 AliCentralitySelectionTask::AliCentralitySelectionTask(const char *name):
 341   AliAnalysisTaskSE(name),
 342   fAnalysisInput("ESD"),
 343   fIsMCInput(kFALSE),
 344   fCurrentRun(-1),
 345   fUseScaling(0),
 346   fUseCleaning(0),
 347   fFillHistos(0),
 348   fV0MScaleFactor(0),
 349   fSPDScaleFactor(0),
 350   fTPCScaleFactor(0),
 351   fV0MScaleFactorMC(0),
 352   fV0MSPDOutlierPar0(0),
 353   fV0MSPDOutlierPar1(0),
 354   fV0MTPCOutlierPar0(0),
 355   fV0MTPCOutlierPar1(0),
 356   fV0MSPDSigmaOutlierPar0(0),
 357   fV0MSPDSigmaOutlierPar1(0),
 358   fV0MSPDSigmaOutlierPar2(0),
 359   fV0MTPCSigmaOutlierPar0(0),
 360   fV0MTPCSigmaOutlierPar1(0),
 361   fV0MTPCSigmaOutlierPar2(0),
 362   fV0MZDCOutlierPar0(0),
 363   fV0MZDCOutlierPar1(0),
 364   fV0MZDCEcalOutlierPar0(0),
 365   fV0MZDCEcalOutlierPar1(0),
 366   fTrackCuts(0),
 367   fEsdTrackCuts(0),
 368   fEsdTrackCutsExtra1(0),
 369   fEsdTrackCutsExtra2(0),
 370   fZVCut(10),
 371   fOutliersCut(5),
 372   fQuality(999),
 373   fIsSelected(0),
 374   fMSL(0),
 375   fMSH(0),
 376   fMUL(0),
 377   fMLL(0),
 378   fEJE(0),
 379   fEGA(0),
 380   fPHS(0),
 381   fMB(0),
 382   fCVHN(0),
 383   fCVLN(0),
 384   fCVHNbit(0),
 385   fCVLNbit(0),
 386   fCCENT(0),
 387   fCSEMI(0),
 388   fCCENTbit(0),
 389   fCSEMIbit(0),
 390   fCentV0M(0),
 391   fCentV0A(0),
 392   fCentV0A0(0),
 393   fCentV0A123(0),
 394   fCentV0C(0),
 395   fCentV0MEq(0),
 396   fCentV0AEq(0),
 397   fCentV0CEq(0),
 398   fCentFMD(0),
 399   fCentTRK(0),
 400   fCentTKL(0),
 401   fCentCL0(0),
 402   fCentCL1(0),
 403   fCentCND(0),
 404   fCentZNA(0),
 405   fCentZNC(0),
 406   fCentZPA(0),
 407   fCentZPC(0),
 408   fCentNPA(0),
 409   fCentV0MvsFMD(0),
 410   fCentTKLvsV0M(0),
 411   fCentZEMvsZDC(0),
 412   fCentV0Mtrue(0),
 413   fCentV0Atrue(0),
 414   fCentV0Ctrue(0),
 415   fCentV0MEqtrue(0),
 416   fCentV0AEqtrue(0),
 417   fCentV0CEqtrue(0),
 418   fCentFMDtrue(0),
 419   fCentTRKtrue(0),
 420   fCentTKLtrue(0),
 421   fCentCL0true(0),
 422   fCentCL1true(0),
 423   fCentCNDtrue(0),
 424   fCentZNAtrue(0),
 425   fCentZNCtrue(0),
 426   fCentZPAtrue(0),
 427   fCentZPCtrue(0),
 428   fHtempV0M(0),
 429   fHtempV0A(0),
 430   fHtempV0A0(0),
 431   fHtempV0A123(0),
 432   fHtempV0C(0),
 433   fHtempV0MEq(0),
 434   fHtempV0AEq(0),
 435   fHtempV0CEq(0),
 436   fHtempFMD(0),
 437   fHtempTRK(0),
 438   fHtempTKL(0),
 439   fHtempCL0(0),
 440   fHtempCL1(0),
 441   fHtempCND(0),
 442   fHtempZNA(0),
 443   fHtempZNC(0),
 444   fHtempZPA(0),
 445   fHtempZPC(0),
 446   fHtempV0MvsFMD(0),
 447   fHtempTKLvsV0M(0),
 448   fHtempZEMvsZDC(0),
 449   fHtempNPA(0),
 450   fHtempV0Mtrue(0),
 451   fHtempV0Atrue(0),
 452   fHtempV0Ctrue(0),
 453   fHtempV0MEqtrue(0),
 454   fHtempV0AEqtrue(0),
 455   fHtempV0CEqtrue(0),
 456   fHtempFMDtrue(0),
 457   fHtempTRKtrue(0),
 458   fHtempTKLtrue(0),
 459   fHtempCL0true(0),
 460   fHtempCL1true(0),
 461   fHtempCNDtrue(0),
 462   fHtempZNAtrue(0),
 463   fHtempZNCtrue(0),
 464   fHtempZPAtrue(0),
 465   fHtempZPCtrue(0),
 466   fOutputList(0),
 467   fHOutCentV0M(0),
 468   fHOutCentV0A(0),
 469   fHOutCentV0A0(0),
 470   fHOutCentV0A123(0),
 471   fHOutCentV0C(0),
 472   fHOutCentV0MEq(0),
 473   fHOutCentV0AEq(0),
 474   fHOutCentV0CEq(0),
 475   fHOutCentV0MCVHN(0),
 476   fHOutCentV0MCVLN(0),
 477   fHOutCentV0MCVHNinMB(0),
 478   fHOutCentV0MCVLNinMB(0),
 479   fHOutCentV0MCCENT(0),
 480   fHOutCentV0MCSEMI(0),
 481   fHOutCentV0MCCENTinMB(0),
 482   fHOutCentV0MCSEMIinMB(0),
 483   fHOutCentV0MMSL(0),
 484   fHOutCentV0MMSH(0),
 485   fHOutCentV0MMUL(0),
 486   fHOutCentV0MMLL(0),
 487   fHOutCentV0MEJE(0),
 488   fHOutCentV0MEGA(0),
 489   fHOutCentV0MPHS(0),
 490   fHOutCentV0MMSLinMB(0),
 491   fHOutCentV0MMSHinMB(0),
 492   fHOutCentV0MMULinMB(0),
 493   fHOutCentV0MMLLinMB(0),
 494   fHOutCentV0MEJEinMB(0),
 495   fHOutCentV0MEGAinMB(0),
 496   fHOutCentV0MPHSinMB(0),
 497   fHOutCentFMD(0),
 498   fHOutCentTRK(0),
 499   fHOutCentTKL(0),
 500   fHOutCentCL0(0),
 501   fHOutCentCL1(0),
 502   fHOutCentCND(0),
 503   fHOutCentNPA(0),
 504   fHOutCentZNA(0),
 505   fHOutCentZNC(0),
 506   fHOutCentZPA(0),
 507   fHOutCentZPC(0),
 508   fHOutCentV0MvsFMD(0),
 509   fHOutCentTKLvsV0M(0),
 510   fHOutCentZEMvsZDC(0),
 511   fHOutCentV0MvsCentCL1(0),
 512   fHOutCentV0MvsCentTRK(0),
 513   fHOutCentTRKvsCentCL1(0),
 514   fHOutCentV0MvsCentZDC(0),
 515   fHOutCentV0AvsCentV0C(0),
 516   fHOutCentV0AvsCentTRK(0),
 517   fHOutCentV0AvsCentCND(0),
 518   fHOutCentV0AvsCentCL1(0),
 519   fHOutCentV0CvsCentTRK(0),
 520   fHOutCentV0CvsCentCND(0),
 521   fHOutCentV0CvsCentCL1(0),
 522   fHOutCentNPAvsCentV0A(0),
 523   fHOutCentNPAvsCentV0C(0),
 524   fHOutCentNPAvsCentTRK(0),
 525   fHOutCentNPAvsCentCND(0),
 526   fHOutCentNPAvsCentCL1(0),
 527   fHOutCentZNAvsCentV0A(0),
 528   fHOutCentZNAvsCentV0C(0),
 529   fHOutCentZNAvsCentTRK(0),
 530   fHOutCentZNAvsCentCND(0),
 531   fHOutCentZNAvsCentCL1(0),
 532   fHOutCentZNAvsCentZPA(0),
 533   fHOutMultV0AC(0),
 534   fHOutMultV0M(0),
 535   fHOutMultV0A(0),
 536   fHOutMultV0A0(0),
 537   fHOutMultV0A123(0),
 538   fHOutMultV0C(0),
 539   fHOutMultV0MEq(0),
 540   fHOutMultV0AEq(0),
 541   fHOutMultV0CEq(0),
 542   fHOutMultV0Mnc(0),
 543   fHOutMultV0Anc(0),
 544   fHOutMultV0Cnc(0),
 545   fHOutMultV0O(0),
 546   fHOutMultV0Cells(0),
 547   fHOutMultFMD(0),
 548   fHOutMultTRK(0),
 549   fHOutMultTKL(0),
 550   fHOutMultCL0(0),
 551   fHOutMultCL1(0),
 552   fHOutMultCND(0),
 553   fHOutMultNPA(0),
 554   fHOutMultZNA(0),
 555   fHOutMultZNC(0),
 556   fHOutMultZPA(0),
 557   fHOutMultZPC(0),
 558   fHOutMultV0MvsZDN(0),
 559   fHOutMultZEMvsZDN(0),
 560   fHOutMultV0MvsZDC(0),
 561   fHOutMultZEMvsZDC(0),
 562   fHOutMultZEMvsZDCw(0),
 563   fHOutMultV0MvsCL1(0),
 564   fHOutMultV0MvsTRK(0),
 565   fHOutMultTRKvsCL1(0),
 566   fHOutMultV0MvsV0O(0),
 567   fHOutMultV0OvsCL1(0),
 568   fHOutMultV0OvsTRK(0),
 569   fHOutMultCL1vsTKL(0),
 570   fHOutMultZNAvsZPA(0),
 571   fHOutCentV0Mqual1(0),
 572   fHOutCentTRKqual1(0),
 573   fHOutCentCL1qual1(0),
 574   fHOutMultV0MvsCL1qual1(0),
 575   fHOutMultV0MvsTRKqual1(0),
 576   fHOutMultTRKvsCL1qual1(0),
 577   fHOutCentV0Mqual2(0),
 578   fHOutCentTRKqual2(0),
 579   fHOutCentCL1qual2(0),
 580   fHOutMultV0MvsCL1qual2(0),
 581   fHOutMultV0MvsTRKqual2(0),
 582   fHOutMultTRKvsCL1qual2(0),
 583   fHOutQuality(0),
 584   fHOutVertex(0),
 585   fHOutVertexT0(0)
 586 {
 587   // Default constructor
 588   AliInfo("Centrality Selection enabled.");
 589   //DefineOutput(1, TList::Class());
 590   fUseScaling=kTRUE;
 591   fUseCleaning=kTRUE;
 592   fFillHistos=kFALSE;
 593   fBranchNames="ESD:AliESDRun.,AliESDHeader.,AliESDZDC.,AliESDFMD.,AliESDVZERO.,AliESDTZERO."
 594     ",SPDVertex.,TPCVertex.,PrimaryVertex.,AliMultiplicity.,Tracks ";
 595 }
 596
 597 //________________________________________________________________________
 598 AliCentralitySelectionTask& AliCentralitySelectionTask::operator=(const AliCentralitySelectionTask& c)
 599 {
 600   // Assignment operator
 601   if (this!=&c) {
 602     AliAnalysisTaskSE::operator=(c);
 603   }
 604   return *this;
 605 }
 606
 607 //________________________________________________________________________
 608 AliCentralitySelectionTask::AliCentralitySelectionTask(const AliCentralitySelectionTask& ana):
 609   AliAnalysisTaskSE(ana),
 610   fAnalysisInput(ana.fAnalysisInput),
 611   fIsMCInput(ana.fIsMCInput),
 612   fCurrentRun(ana.fCurrentRun),
 613   fUseScaling(ana.fUseScaling),
 614   fUseCleaning(ana.fUseCleaning),
 615   fFillHistos(ana.fFillHistos),
 616   fV0MScaleFactor(ana.fV0MScaleFactor),
 617   fSPDScaleFactor(ana.fSPDScaleFactor),
 618   fTPCScaleFactor(ana.fTPCScaleFactor),
 619   fV0MScaleFactorMC(ana.fV0MScaleFactorMC),
 620   fV0MSPDOutlierPar0(ana.fV0MSPDOutlierPar0),
 621   fV0MSPDOutlierPar1(ana.fV0MSPDOutlierPar1),
 622   fV0MTPCOutlierPar0(ana.fV0MTPCOutlierPar0),
 623   fV0MTPCOutlierPar1(ana.fV0MTPCOutlierPar1),
 624   fV0MSPDSigmaOutlierPar0(ana.fV0MSPDSigmaOutlierPar0),
 625   fV0MSPDSigmaOutlierPar1(ana.fV0MSPDSigmaOutlierPar1),
 626   fV0MSPDSigmaOutlierPar2(ana.fV0MSPDSigmaOutlierPar2),
 627   fV0MTPCSigmaOutlierPar0(ana.fV0MTPCSigmaOutlierPar0),
 628   fV0MTPCSigmaOutlierPar1(ana.fV0MTPCSigmaOutlierPar1),
 629   fV0MTPCSigmaOutlierPar2(ana.fV0MTPCSigmaOutlierPar2),
 630   fV0MZDCOutlierPar0(ana.fV0MZDCOutlierPar0),
 631   fV0MZDCOutlierPar1(ana.fV0MZDCOutlierPar1),
 632   fV0MZDCEcalOutlierPar0(ana.fV0MZDCEcalOutlierPar0),
 633   fV0MZDCEcalOutlierPar1(ana.fV0MZDCEcalOutlierPar1),
 634   fTrackCuts(ana.fTrackCuts),
 635   fEsdTrackCuts(ana.fEsdTrackCuts),
 636   fEsdTrackCutsExtra1(ana.fEsdTrackCutsExtra1),
 637   fEsdTrackCutsExtra2(ana.fEsdTrackCutsExtra2),
 638   fZVCut(ana.fZVCut),
 639   fOutliersCut(ana.fOutliersCut),
 640   fQuality(ana.fQuality),
 641   fIsSelected(ana.fIsSelected),
 642   fMSL(ana.fMSL),
 643   fMSH(ana.fMSH),
 644   fMUL(ana.fMUL),
 645   fMLL(ana.fMLL),
 646   fEJE(ana.fEJE),
 647   fEGA(ana.fEGA),
 648   fPHS(ana.fPHS),
 649   fMB(ana.fMB),
 650   fCVHN(ana.fCVHN),
 651   fCVLN(ana.fCVLN),
 652   fCVHNbit(ana.fCVHNbit),
 653   fCVLNbit(ana.fCVLNbit),
 654   fCCENT(ana.fCCENT),
 655   fCSEMI(ana.fCSEMI),
 656   fCCENTbit(ana.fCCENTbit),
 657   fCSEMIbit(ana.fCSEMIbit),
 658   fCentV0M(ana.fCentV0M),
 659   fCentV0A(ana.fCentV0A),
 660   fCentV0A0(ana.fCentV0A0),
 661   fCentV0A123(ana.fCentV0A123),
 662   fCentV0C(ana.fCentV0C),
 663   fCentV0MEq(ana.fCentV0MEq),
 664   fCentV0AEq(ana.fCentV0AEq),
 665   fCentV0CEq(ana.fCentV0CEq),
 666   fCentFMD(ana.fCentFMD),
 667   fCentTRK(ana.fCentTRK),
 668   fCentTKL(ana.fCentTKL),
 669   fCentCL0(ana.fCentCL0),
 670   fCentCL1(ana.fCentCL1),
 671   fCentCND(ana.fCentCND),
 672   fCentZNA(ana.fCentZNA),
 673   fCentZNC(ana.fCentZNC),
 674   fCentZPA(ana.fCentZPA),
 675   fCentZPC(ana.fCentZPC),
 676   fCentNPA(ana.fCentNPA),
 677   fCentV0MvsFMD(ana.fCentV0MvsFMD),
 678   fCentTKLvsV0M(ana.fCentTKLvsV0M),
 679   fCentZEMvsZDC(ana.fCentZEMvsZDC),
 680   fCentV0Mtrue(ana.fCentV0Mtrue),
 681   fCentV0Atrue(ana.fCentV0Atrue),
 682   fCentV0Ctrue(ana.fCentV0Ctrue),
 683   fCentV0MEqtrue(ana.fCentV0MEqtrue),
 684   fCentV0AEqtrue(ana.fCentV0AEqtrue),
 685   fCentV0CEqtrue(ana.fCentV0CEqtrue),
 686   fCentFMDtrue(ana.fCentFMDtrue),
 687   fCentTRKtrue(ana.fCentTRKtrue),
 688   fCentTKLtrue(ana.fCentTKLtrue),
 689   fCentCL0true(ana.fCentCL0true),
 690   fCentCL1true(ana.fCentCL1true),
 691   fCentCNDtrue(ana.fCentCNDtrue),
 692   fCentZNAtrue(ana.fCentZNAtrue),
 693   fCentZNCtrue(ana.fCentZNCtrue),
 694   fCentZPAtrue(ana.fCentZPAtrue),
 695   fCentZPCtrue(ana.fCentZPCtrue),
 696   fHtempV0M(ana.fHtempV0M),
 697   fHtempV0A(ana.fHtempV0A),
 698   fHtempV0A0(ana.fHtempV0A0),
 699   fHtempV0A123(ana.fHtempV0A123),
 700   fHtempV0C(ana.fHtempV0C),
 701   fHtempV0MEq(ana.fHtempV0MEq),
 702   fHtempV0AEq(ana.fHtempV0AEq),
 703   fHtempV0CEq(ana.fHtempV0CEq),
 704   fHtempFMD(ana.fHtempFMD),
 705   fHtempTRK(ana.fHtempTRK),
 706   fHtempTKL(ana.fHtempTKL),
 707   fHtempCL0(ana.fHtempCL0),
 708   fHtempCL1(ana.fHtempCL1),
 709   fHtempCND(ana.fHtempCND),
 710   fHtempZNA(ana.fHtempZNA),
 711   fHtempZNC(ana.fHtempZNC),
 712   fHtempZPA(ana.fHtempZPA),
 713   fHtempZPC(ana.fHtempZPC),
 714   fHtempV0MvsFMD(ana.fHtempV0MvsFMD),
 715   fHtempTKLvsV0M(ana.fHtempTKLvsV0M),
 716   fHtempZEMvsZDC(ana.fHtempZEMvsZDC),
 717   fHtempNPA(ana.fHtempNPA),
 718   fHtempV0Mtrue(ana.fHtempV0Mtrue),
 719   fHtempV0Atrue(ana.fHtempV0Atrue),
 720   fHtempV0Ctrue(ana.fHtempV0Ctrue),
 721   fHtempV0MEqtrue(ana.fHtempV0MEqtrue),
 722   fHtempV0AEqtrue(ana.fHtempV0AEqtrue),
 723   fHtempV0CEqtrue(ana.fHtempV0CEqtrue),
 724   fHtempFMDtrue(ana.fHtempFMDtrue),
 725   fHtempTRKtrue(ana.fHtempTRKtrue),
 726   fHtempTKLtrue(ana.fHtempTKLtrue),
 727   fHtempCL0true(ana.fHtempCL0true),
 728   fHtempCL1true(ana.fHtempCL1true),
 729   fHtempCNDtrue(ana.fHtempCNDtrue),
 730   fHtempZNAtrue(ana.fHtempZNAtrue),
 731   fHtempZNCtrue(ana.fHtempZNCtrue),
 732   fHtempZPAtrue(ana.fHtempZPAtrue),
 733   fHtempZPCtrue(ana.fHtempZPCtrue),
 734   fOutputList(ana.fOutputList),
 735   fHOutCentV0M(ana.fHOutCentV0M),
 736   fHOutCentV0A(ana.fHOutCentV0A),
 737   fHOutCentV0A0(ana.fHOutCentV0A0),
 738   fHOutCentV0A123(ana.fHOutCentV0A123),
 739   fHOutCentV0C(ana.fHOutCentV0C),
 740   fHOutCentV0MEq(ana.fHOutCentV0MEq),
 741   fHOutCentV0AEq(ana.fHOutCentV0AEq),
 742   fHOutCentV0CEq(ana.fHOutCentV0CEq),
 743   fHOutCentV0MCVHN(ana.fHOutCentV0MCVHN),
 744   fHOutCentV0MCVLN(ana.fHOutCentV0MCVLN),
 745   fHOutCentV0MCVHNinMB(ana.fHOutCentV0MCVHNinMB),
 746   fHOutCentV0MCVLNinMB(ana.fHOutCentV0MCVLNinMB),
 747   fHOutCentV0MCCENT(ana.fHOutCentV0MCCENT),
 748   fHOutCentV0MCSEMI(ana.fHOutCentV0MCSEMI),
 749   fHOutCentV0MCCENTinMB(ana.fHOutCentV0MCCENTinMB),
 750   fHOutCentV0MCSEMIinMB(ana.fHOutCentV0MCSEMIinMB),
 751   fHOutCentV0MMSL(ana.fHOutCentV0MMSL),
 752   fHOutCentV0MMSH(ana.fHOutCentV0MMSH),
 753   fHOutCentV0MMUL(ana.fHOutCentV0MMUL),
 754   fHOutCentV0MMLL(ana.fHOutCentV0MMLL),
 755   fHOutCentV0MEJE(ana.fHOutCentV0MEJE),
 756   fHOutCentV0MEGA(ana.fHOutCentV0MEGA),
 757   fHOutCentV0MPHS(ana.fHOutCentV0MPHS),
 758   fHOutCentV0MMSLinMB(ana.fHOutCentV0MMSLinMB),
 759   fHOutCentV0MMSHinMB(ana.fHOutCentV0MMSHinMB),
 760   fHOutCentV0MMULinMB(ana.fHOutCentV0MMULinMB),
 761   fHOutCentV0MMLLinMB(ana.fHOutCentV0MMLLinMB),
 762   fHOutCentV0MEJEinMB(ana.fHOutCentV0MEJEinMB),
 763   fHOutCentV0MEGAinMB(ana.fHOutCentV0MEGAinMB),
 764   fHOutCentV0MPHSinMB(ana.fHOutCentV0MPHSinMB),
 765   fHOutCentFMD(ana.fHOutCentFMD),
 766   fHOutCentTRK(ana.fHOutCentTRK),
 767   fHOutCentTKL(ana.fHOutCentTKL),
 768   fHOutCentCL0(ana.fHOutCentCL0),
 769   fHOutCentCL1(ana.fHOutCentCL1),
 770   fHOutCentCND(ana.fHOutCentCND),
 771   fHOutCentNPA(ana.fHOutCentNPA),
 772   fHOutCentZNA(ana.fHOutCentZNA),
 773   fHOutCentZNC(ana.fHOutCentZNC),
 774   fHOutCentZPA(ana.fHOutCentZPA),
 775   fHOutCentZPC(ana.fHOutCentZPC),
 776   fHOutCentV0MvsFMD(ana.fHOutCentV0MvsFMD),
 777   fHOutCentTKLvsV0M(ana.fHOutCentTKLvsV0M),
 778   fHOutCentZEMvsZDC(ana.fHOutCentZEMvsZDC),
 779   fHOutCentV0MvsCentCL1(ana.fHOutCentV0MvsCentCL1),
 780   fHOutCentV0MvsCentTRK(ana.fHOutCentV0MvsCentTRK),
 781   fHOutCentTRKvsCentCL1(ana.fHOutCentTRKvsCentCL1),
 782   fHOutCentV0MvsCentZDC(ana.fHOutCentV0MvsCentZDC),
 783   fHOutCentV0AvsCentV0C(ana.fHOutCentV0AvsCentV0C),
 784   fHOutCentV0AvsCentTRK(ana.fHOutCentV0AvsCentTRK),
 785   fHOutCentV0AvsCentCND(ana.fHOutCentV0AvsCentCND),
 786   fHOutCentV0AvsCentCL1(ana.fHOutCentV0AvsCentCL1),
 787   fHOutCentV0CvsCentTRK(ana.fHOutCentV0CvsCentTRK),
 788   fHOutCentV0CvsCentCND(ana.fHOutCentV0CvsCentCND),
 789   fHOutCentV0CvsCentCL1(ana.fHOutCentV0CvsCentCL1),
 790   fHOutCentNPAvsCentV0A(ana.fHOutCentNPAvsCentV0A),
 791   fHOutCentNPAvsCentV0C(ana.fHOutCentNPAvsCentV0C),
 792   fHOutCentNPAvsCentTRK(ana.fHOutCentNPAvsCentTRK),
 793   fHOutCentNPAvsCentCND(ana.fHOutCentNPAvsCentCND),
 794   fHOutCentNPAvsCentCL1(ana.fHOutCentNPAvsCentCL1),
 795   fHOutCentZNAvsCentV0A(ana.fHOutCentZNAvsCentV0A),
 796   fHOutCentZNAvsCentV0C(ana.fHOutCentZNAvsCentV0C),
 797   fHOutCentZNAvsCentTRK(ana.fHOutCentZNAvsCentTRK),
 798   fHOutCentZNAvsCentCND(ana.fHOutCentZNAvsCentCND),
 799   fHOutCentZNAvsCentCL1(ana.fHOutCentZNAvsCentCL1),
 800   fHOutCentZNAvsCentZPA(ana.fHOutCentZNAvsCentZPA),
 801   fHOutMultV0AC(ana.fHOutMultV0AC),
 802   fHOutMultV0M(ana.fHOutMultV0M),
 803   fHOutMultV0A(ana.fHOutMultV0A),
 804   fHOutMultV0A0(ana.fHOutMultV0A0),
 805   fHOutMultV0A123(ana.fHOutMultV0A123),
 806   fHOutMultV0C(ana.fHOutMultV0C),
 807   fHOutMultV0MEq(ana.fHOutMultV0MEq),
 808   fHOutMultV0AEq(ana.fHOutMultV0AEq),
 809   fHOutMultV0CEq(ana.fHOutMultV0CEq),
 810   fHOutMultV0Mnc(ana.fHOutMultV0Mnc),
 811   fHOutMultV0Anc(ana.fHOutMultV0Anc),
 812   fHOutMultV0Cnc(ana.fHOutMultV0Cnc),
 813   fHOutMultV0O(ana.fHOutMultV0O),
 814   fHOutMultV0Cells(ana.fHOutMultV0Cells),
 815   fHOutMultFMD(ana.fHOutMultFMD),
 816   fHOutMultTRK(ana.fHOutMultTRK),
 817   fHOutMultTKL(ana.fHOutMultTKL),
 818   fHOutMultCL0(ana.fHOutMultCL0),
 819   fHOutMultCL1(ana.fHOutMultCL1),
 820   fHOutMultCND(ana.fHOutMultCND),
 821   fHOutMultNPA(ana.fHOutMultNPA),
 822   fHOutMultZNA(ana.fHOutMultZNA),
 823   fHOutMultZNC(ana.fHOutMultZNC),
 824   fHOutMultZPA(ana.fHOutMultZPA),
 825   fHOutMultZPC(ana.fHOutMultZPC),
 826   fHOutMultV0MvsZDN(ana.fHOutMultV0MvsZDN),
 827   fHOutMultZEMvsZDN(ana.fHOutMultZEMvsZDN),
 828   fHOutMultV0MvsZDC(ana.fHOutMultV0MvsZDC),
 829   fHOutMultZEMvsZDC(ana.fHOutMultZEMvsZDC),
 830   fHOutMultZEMvsZDCw(ana.fHOutMultZEMvsZDCw),
 831   fHOutMultV0MvsCL1(ana.fHOutMultV0MvsCL1),
 832   fHOutMultV0MvsTRK(ana.fHOutMultV0MvsTRK),
 833   fHOutMultTRKvsCL1(ana.fHOutMultTRKvsCL1),
 834   fHOutMultV0MvsV0O(ana.fHOutMultV0MvsV0O),
 835   fHOutMultV0OvsCL1(ana.fHOutMultV0OvsCL1),
 836   fHOutMultV0OvsTRK(ana.fHOutMultV0OvsTRK),
 837   fHOutMultCL1vsTKL(ana.fHOutMultCL1vsTKL),
 838   fHOutMultZNAvsZPA(ana.fHOutMultZNAvsZPA),
 839   fHOutCentV0Mqual1(ana.fHOutCentV0Mqual1),
 840   fHOutCentTRKqual1(ana.fHOutCentTRKqual1),
 841   fHOutCentCL1qual1(ana.fHOutCentCL1qual1),
 842   fHOutMultV0MvsCL1qual1(ana.fHOutMultV0MvsCL1qual1),
 843   fHOutMultV0MvsTRKqual1(ana.fHOutMultV0MvsTRKqual1),
 844   fHOutMultTRKvsCL1qual1(ana.fHOutMultTRKvsCL1qual1),
 845   fHOutCentV0Mqual2(ana.fHOutCentV0Mqual2),
 846   fHOutCentTRKqual2(ana.fHOutCentTRKqual2),
 847   fHOutCentCL1qual2(ana.fHOutCentCL1qual2),
 848   fHOutMultV0MvsCL1qual2(ana.fHOutMultV0MvsCL1qual2),
 849   fHOutMultV0MvsTRKqual2(ana.fHOutMultV0MvsTRKqual2),
 850   fHOutMultTRKvsCL1qual2(ana.fHOutMultTRKvsCL1qual2),
 851   fHOutQuality(ana.fHOutQuality),
 852   fHOutVertex(ana.fHOutVertex),
 853   fHOutVertexT0(ana.fHOutVertexT0)
 854 {
 855   // Copy Constructor
 856
 857 }
 858
 859 //________________________________________________________________________
 860 AliCentralitySelectionTask::~AliCentralitySelectionTask()
 861 {
 862   // Destructor
 863   if (fOutputList && !AliAnalysisManager::GetAnalysisManager()->IsProofMode()) delete fOutputList;
 864   if (fTrackCuts) delete fTrackCuts;
 865   if (fEsdTrackCuts) delete fEsdTrackCuts;
 866   if (fEsdTrackCutsExtra1) delete fEsdTrackCutsExtra1;
 867   if (fEsdTrackCutsExtra2) delete fEsdTrackCutsExtra2;
 868 }
 869
 870 //________________________________________________________________________
 871 void AliCentralitySelectionTask::UserCreateOutputObjects()
 872 {
 873   // Create the output containers
 874   if(fDebug>1) printf("AnalysisCentralitySelectionTask::UserCreateOutputObjects() \n");
 875   AliLog::SetClassDebugLevel("AliCentralitySelectionTask", AliLog::kInfo);
 876
 877   if (fFillHistos) {
 878     fOutputList = new TList();
 879     fOutputList->SetOwner();
 880     fHOutCentV0M     = new TH1F("fHOutCentV0M","fHOutCentV0M; Centrality V0",505,0,101);
 881     fHOutCentV0A    = new TH1F("fHOutCentV0A","fHOutCentV0A; Centrality V0A",505,0,101);
 882     fHOutCentV0A0 = new TH1F("fHOutCentV0A0","fHOutCentV0A0; Centrality V0A-0",505,0,101);
 883     fHOutCentV0A123 = new TH1F("fHOutCentV0A123","fHOutCentV0A123; Centrality V0A-123",505,0,101);
 884     fHOutCentV0C    = new TH1F("fHOutCentV0C","fHOutCentV0C; Centrality V0C",505,0,101);
 885     fHOutCentV0MEq     = new TH1F("fHOutCentV0MEq","fHOutCentV0MEq; Centrality V0 equalized",505,0,101);
 886     fHOutCentV0AEq    = new TH1F("fHOutCentV0AEq","fHOutCentV0AEq; Centrality V0A equalized",505,0,101);
 887     fHOutCentV0CEq    = new TH1F("fHOutCentV0CEq","fHOutCentV0CEq; Centrality V0C equalized",505,0,101);
 888     fHOutCentV0MCVHN= new TH1F("fHOutCentV0M_CVHN","fHOutCentV0M_CVHN; Centrality V0",505,0,101);
 889     fHOutCentV0MCVLN= new TH1F("fHOutCentV0M_CVLN","fHOutCentV0M_CVLN; Centrality V0",505,0,101);
 890     fHOutCentV0MCVHNinMB= new TH1F("fHOutCentV0M_CVHNinMB","fHOutCentV0M_CVHN; Centrality V0",505,0,101);
 891     fHOutCentV0MCVLNinMB= new TH1F("fHOutCentV0M_CVLNinMB","fHOutCentV0M_CVLN; Centrality V0",505,0,101);
 892     fHOutCentV0MCCENT= new TH1F("fHOutCentV0M_CCENT","fHOutCentV0M_CCENT; Centrality V0",505,0,101);
 893     fHOutCentV0MCSEMI= new TH1F("fHOutCentV0M_CSEMI","fHOutCentV0M_CSEMI; Centrality V0",505,0,101);
 894     fHOutCentV0MCCENTinMB= new TH1F("fHOutCentV0M_CCENTinMB","fHOutCentV0M_CCENT; Centrality V0",505,0,101);
 895     fHOutCentV0MCSEMIinMB= new TH1F("fHOutCentV0M_CSEMIinMB","fHOutCentV0M_CSEMI; Centrality V0",505,0,101);
 896     fHOutCentV0MMSL= new TH1F("fHOutCentV0M_MSL","fHOutCentV0M_MSL; Centrality V0",505,0,101);
 897     fHOutCentV0MMSH= new TH1F("fHOutCentV0M_MSH","fHOutCentV0M_MSH; Centrality V0",505,0,101);
 898     fHOutCentV0MMUL= new TH1F("fHOutCentV0M_MUL","fHOutCentV0M_MUL; Centrality V0",505,0,101);
 899     fHOutCentV0MMLL= new TH1F("fHOutCentV0M_MLL","fHOutCentV0M_MLL; Centrality V0",505,0,101);
 900     fHOutCentV0MEJE= new TH1F("fHOutCentV0M_EJE","fHOutCentV0M_EJE; Centrality V0",505,0,101);
 901     fHOutCentV0MEGA= new TH1F("fHOutCentV0M_EGA","fHOutCentV0M_EGA; Centrality V0",505,0,101);
 902     fHOutCentV0MPHS= new TH1F("fHOutCentV0M_PHS","fHOutCentV0M_PHS; Centrality V0",505,0,101);
 903     fHOutCentV0MMSLinMB= new TH1F("fHOutCentV0M_MSLinMB","fHOutCentV0M_MSLinMB; Centrality V0",505,0,101);
 904     fHOutCentV0MMSHinMB= new TH1F("fHOutCentV0M_MSHinMB","fHOutCentV0M_MSHinMB; Centrality V0",505,0,101);
 905     fHOutCentV0MMULinMB= new TH1F("fHOutCentV0M_MULinMB","fHOutCentV0M_MULinMB; Centrality V0",505,0,101);
 906     fHOutCentV0MMLLinMB= new TH1F("fHOutCentV0M_MLLinMB","fHOutCentV0M_MLLinMB; Centrality V0",505,0,101);
 907     fHOutCentV0MEJEinMB= new TH1F("fHOutCentV0M_EJEinMB","fHOutCentV0M_EJEinMB; Centrality V0",505,0,101);
 908     fHOutCentV0MEGAinMB= new TH1F("fHOutCentV0M_EGAinMB","fHOutCentV0M_EGAinMB; Centrality V0",505,0,101);
 909     fHOutCentV0MPHSinMB= new TH1F("fHOutCentV0M_PHSinMB","fHOutCentV0M_PHSinMB; Centrality V0",505,0,101);
 910     fHOutCentFMD     = new TH1F("fHOutCentFMD","fHOutCentFMD; Centrality FMD",505,0,101);
 911     fHOutCentTRK     = new TH1F("fHOutCentTRK","fHOutCentTRK; Centrality TPC",505,0,101);
 912     fHOutCentTKL     = new TH1F("fHOutCentTKL","fHOutCentTKL; Centrality tracklets",505,0,101);
 913     fHOutCentCL0     = new TH1F("fHOutCentCL0","fHOutCentCL0; Centrality SPD inner",505,0,101);
 914     fHOutCentCL1     = new TH1F("fHOutCentCL1","fHOutCentCL1; Centrality SPD outer",505,0,101);
 915     fHOutCentCND     = new TH1F("fHOutCentCND","fHOutCentCND; Centrality candle",505,0,101);
 916     fHOutCentNPA     = new TH1F("fHOutCentNPA","fHOutCentNPA; Centrality Npart",505,0,101);
 917     fHOutCentZNA     = new TH1F("fHOutCentZNA","fHOutCentZNA; Centrality ZNA",505,0,101);
 918     fHOutCentZNC     = new TH1F("fHOutCentZNC","fHOutCentZNC; Centrality ZNC",505,0,101);
 919     fHOutCentZPA     = new TH1F("fHOutCentZPA","fHOutCentZPA; Centrality ZPA",505,0,101);
 920     fHOutCentZPC     = new TH1F("fHOutCentZPC","fHOutCentZPC; Centrality ZPC",505,0,101);
 921     fHOutCentV0MvsFMD= new TH1F("fHOutCentV0MvsFMD","fHOutCentV0MvsFMD; Centrality V0 vs FMD",505,0,101);
 922     fHOutCentTKLvsV0M= new TH1F("fHOutCentTKLvsV0M","fHOutCentTKLvsV0M; Centrality tracklets vs V0",505,0,101);
 923     fHOutCentZEMvsZDC= new TH1F("fHOutCentZEMvsZDC","fHOutCentZEMvsZDC; Centrality ZEM vs ZDC",505,0,101);
 924     fHOutCentV0MvsCentCL1= new TH2F("fHOutCentV0MvsCentCL1","fHOutCentV0MvsCentCL1; Cent V0; Cent SPD",505,0,101,505,0,101);
 925     fHOutCentV0MvsCentTRK= new TH2F("fHOutCentV0MvsCentTRK","fHOutCentV0MvsCentTRK; Cent V0; Cent TPC",505,0,101,505,0,101);
 926     fHOutCentTRKvsCentCL1= new TH2F("fHOutCentTRKvsCentCL1","fHOutCentTRKvsCentCL1; Cent TPC; Cent SPD",505,0,101,505,0,101);
 927     fHOutCentV0MvsCentZDC= new TH2F("fHOutCentV0MvsCentZDC","fHOutCentV0MvsCentZDC; Cent V0; Cent ZDC",505,0,101,505,0,101);
 928     fHOutCentV0AvsCentV0C= new TH2F("fHOutCentV0AvsCentV0C","fHOutCentV0AvsCentV0C; Cent V0A; Cent V0C;", 505,0,101,505,0,101);
 929     fHOutCentV0AvsCentTRK= new TH2F("fHOutCentV0AvsCentTRK","fHOutCentV0AvsCentTRK; Cent V0A; Cent TRK;", 505,0,101,505,0,101);
 930     fHOutCentV0AvsCentCND= new TH2F("fHOutCentV0AvsCentCND","fHOutCentV0AvsCentCND; Cent V0A; Cent CND;", 505,0,101,505,0,101);
 931     fHOutCentV0AvsCentCL1= new TH2F("fHOutCentV0AvsCentCL1","fHOutCentV0AvsCentCL1; Cent V0A; Cent CL1;", 505,0,101,505,0,101);
 932     fHOutCentV0CvsCentTRK= new TH2F("fHOutCentV0CvsCentTRK","fHOutCentV0CvsCentTRK; Cent V0C; Cent TRK;", 505,0,101,505,0,101);
 933     fHOutCentV0CvsCentCND= new TH2F("fHOutCentV0CvsCentCND","fHOutCentV0CvsCentCND; Cent V0C; Cent CND;", 505,0,101,505,0,101);
 934     fHOutCentV0CvsCentCL1= new TH2F("fHOutCentV0CvsCentCL1","fHOutCentV0CvsCentCL1; Cent V0C; Cent CL1;", 505,0,101,505,0,101);
 935     fHOutCentNPAvsCentV0A= new TH2F("fHOutCentNPAvsCentV0A","fHOutCentNPAvsCentV0A; Cent NPA; Cent V0A;", 505,0,101,505,0,101);
 936     fHOutCentNPAvsCentV0C= new TH2F("fHOutCentNPAvsCentV0C","fHOutCentNPAvsCentV0C; Cent NPA; Cent V0C;", 505,0,101,505,0,101);
 937     fHOutCentNPAvsCentTRK= new TH2F("fHOutCentNPAvsCentTRK","fHOutCentNPAvsCentTRK; Cent NPA; Cent TRK;", 505,0,101,505,0,101);
 938     fHOutCentNPAvsCentCND= new TH2F("fHOutCentNPAvsCentCND","fHOutCentNPAvsCentCND; Cent NPA; Cent CND;", 505,0,101,505,0,101);
 939     fHOutCentNPAvsCentCL1= new TH2F("fHOutCentNPAvsCentCL1","fHOutCentNPAvsCentCL1; Cent NPA; Cent CL1;", 505,0,101,505,0,101);
 940     fHOutCentZNAvsCentV0A= new TH2F("fHOutCentZNAvsCentV0A","fHOutCentZNAvsCentV0A; Cent ZNA; Cent V0A;", 505,0,101,505,0,101);
 941     fHOutCentZNAvsCentV0C= new TH2F("fHOutCentZNAvsCentV0C","fHOutCentZNAvsCentV0C; Cent ZNA; Cent V0C;", 505,0,101,505,0,101);
 942     fHOutCentZNAvsCentTRK= new TH2F("fHOutCentZNAvsCentTRK","fHOutCentZNAvsCentTRK; Cent ZNA; Cent TRK;", 505,0,101,505,0,101);
 943     fHOutCentZNAvsCentCND= new TH2F("fHOutCentZNAvsCentCND","fHOutCentZNAvsCentCND; Cent ZNA; Cent CND;", 505,0,101,505,0,101);
 944     fHOutCentZNAvsCentCL1= new TH2F("fHOutCentZNAvsCentCL1","fHOutCentZNAvsCentCL1; Cent ZNA; Cent CL1;", 505,0,101,505,0,101);
 945     fHOutCentZNAvsCentZPA= new TH2F("fHOutCentZNAvsCentZPA","fHOutCentZNAvsCentZPA; Cent ZNA; Cent ZPA;", 505,0,101,505,0,101);
 946
 947     fHOutMultV0AC = new TH2F("fHOutMultV0AC","fHOutMultV0AC; Multiplicity V0A; Multiplicity V0C",1000,0,1000,1000,0,1000);
 948     fHOutMultV0M  = new TH1F("fHOutMultV0M","fHOutMultV0M; Multiplicity V0",25000,0,25000);
 949     fHOutMultV0A  = new TH1F("fHOutMultV0A","fHOutMultV0A; Multiplicity V0",25000,0,25000);
 950     fHOutMultV0A0  = new TH1F("fHOutMultV0A0","fHOutMultV0A0; Multiplicity V0A-0",25000,0,25000);
 951     fHOutMultV0A123  = new TH1F("fHOutMultV0A123","fHOutMultV0A123; Multiplicity V0A-123",25000,0,25000);
 952     fHOutMultV0C  = new TH1F("fHOutMultV0C","fHOutMultV0C; Multiplicity V0",25000,0,25000);
 953     fHOutMultV0MEq  = new TH1F("fHOutMultV0MEq","fHOutMultV0MEq; Multiplicity V0",25000,0,25000);
 954     fHOutMultV0AEq  = new TH1F("fHOutMultV0AEq","fHOutMultV0AEq; Multiplicity V0",25000,0,25000);
 955     fHOutMultV0CEq  = new TH1F("fHOutMultV0CEq","fHOutMultV0CEq; Multiplicity V0",25000,0,25000);
 956     fHOutMultV0Mnc= new TH1F("fHOutMultV0Mnc","fHOutMultV0Mnc; Multiplicity V0",25000,0,25000);
 957     fHOutMultV0Anc= new TH1F("fHOutMultV0Anc","fHOutMultV0Anc; Multiplicity V0",25000,0,25000);
 958     fHOutMultV0Cnc= new TH1F("fHOutMultV0Cnc","fHOutMultV0Cnc; Multiplicity V0",25000,0,25000);
 959     fHOutMultV0O  = new TH1F("fHOutMultV0O","fHOutMultV0O; Multiplicity V0",40000,0,40000);
 960     fHOutMultV0Cells = new TH2F("fHOutMultV0Cells","fHOutMultV0Cells",33,-0.5,32.5,33,-0.5,32.5);
 961     fHOutMultFMD = new TH1F("fHOutMultFMD","fHOutMultFMD; Multiplicity FMD",24000,0,24000);
 962     fHOutMultTRK = new TH1F("fHOutMultTRK","fHOutMultTRK; Multiplicity TPC",4000,0,4000);
 963     fHOutMultTKL = new TH1F("fHOutMultTKL","fHOutMultTKL; Multiplicity tracklets",5000,0,5000);
 964     fHOutMultCL0 = new TH1F("fHOutMultCL0","fHOutMultCL0; Multiplicity SPD inner",7000,0,7000);
 965     fHOutMultCL1 = new TH1F("fHOutMultCL1","fHOutMultCL1; Multiplicity SPD outer",7000,0,7000);
 966     fHOutMultCND = new TH1F("fHOutMultCND","fHOutMultCND; Multiplicity candle",4000,0,4000);
 967     fHOutMultNPA = new TH1F("fHOutMultNPA","fHOutMultNPA; Nparticipants",450,0,450);
 968     fHOutMultZNA = new TH1F("fHOutMultZNA","fHOutMultZNA; ZNA Energy",2000,0,2000);
 969     fHOutMultZNC = new TH1F("fHOutMultZNC","fHOutMultZNC; ZNC Energy",2000,0,2000);
 970     fHOutMultZPA = new TH1F("fHOutMultZPA","fHOutMultZPA; ZPA Energy",2000,0,2000);
 971     fHOutMultZPC = new TH1F("fHOutMultZPC","fHOutMultZPC; ZPC Energy",2000,0,2000);
 972
 973     fHOutMultV0MvsZDN = new TH2F("fHOutMultV0MvsZDN","fHOutMultV0MvsZDN; Multiplicity V0; Energy ZDC-N",500,0,1200,500,0,2000);
 974     fHOutMultZEMvsZDN = new TH2F("fHOutMultZEMvsZDN","fHOutMultZEMvsZDN; Energy ZEM; Energy ZDC-N",500,0,2500,500,0,180000);
 975     fHOutMultV0MvsZDC = new TH2F("fHOutMultV0MvsZDC","fHOutMultV0MvsZDC; Multiplicity V0; Energy ZDC",500,0,30000,500,0,200000);
 976     fHOutMultZEMvsZDC = new TH2F("fHOutMultZEMvsZDC","fHOutMultZEMvsZDC; Energy ZEM; Energy ZDC",500,0,2500,500,0,200000);
 977     fHOutMultZEMvsZDCw = new TH2F("fHOutMultZEMvsZDCw","fHOutMultZEMvsZDCw; Energy ZEM; Energy ZDC (weigthed with V0 percentile)",500,0,2500,500,0,200000);
 978     fHOutMultV0MvsCL1 = new TH2F("fHOutMultV0MvsCL1","fHOutMultV0MvsCL1; Multiplicity V0; Multiplicity SPD outer",2500,0,30000,700,0,7000);
 979     fHOutMultV0MvsTRK = new TH2F("fHOutMultV0MvsTRK","fHOutMultV0MvsTRK; Multiplicity V0; Multiplicity TPC",2500,0,30000,400,0,4000);
 980     fHOutMultTRKvsCL1 = new TH2F("fHOutMultTRKvsCL1","fHOutMultTRKvsCL1; Multiplicity TPC; Multiplicity SPD outer",400,0,4000,700,0,7000);
 981     fHOutMultV0MvsV0O = new TH2F("fHOutMultV0MvsV0O","fHOutMultV0MvsV0O; Multiplicity V0; Multiplicity V0 Online",500,0,30000,500,0,40000);
 982     fHOutMultV0OvsCL1 = new TH2F("fHOutMultV0OvsCL1","fHOutMultV0OvsCL1; Multiplicity V0; Multiplicity SPD outer",500,0,40000,700,0,7000);
 983     fHOutMultV0OvsTRK = new TH2F("fHOutMultV0OvsTRK","fHOutMultV0OvsTRK; Multiplicity V0; Multiplicity TPC",500,0,40000,400,0,4000);
 984     fHOutMultV0MvsV0O = new TH2F("fHOutMultV0MvsV0O","fHOutMultV0MvsV0O; Multiplicity V0; Multiplicity V0 Online",500,0,30000,500,0,30000);
 985     fHOutMultV0OvsCL1 = new TH2F("fHOutMultV0OvsCL1","fHOutMultV0OvsCL1; Multiplicity V0; Multiplicity SPD outer",2500,0,30000,700,0,7000);
 986     fHOutMultV0OvsTRK = new TH2F("fHOutMultV0OvsTRK","fHOutMultV0OvsTRK; Multiplicity V0; Multiplicity TPC",2500,0,30000,400,0,4000);
 987     fHOutMultCL1vsTKL = new TH2F ("fHOutMultCL1vsTKL","fHOutMultCL1vsTKL; Multiplicity SPD outer; Multiplicity tracklets",700,0,7000,700,0,7000);
 988     fHOutMultZNAvsZPA = new TH2F ("fHOutMultZNAvsZPA","fHOutMultZNAvsZPA; Energy ZNA; Energy ZPA",500,0,2000,500,0,2000);
 989
 990     fHOutCentV0Mqual1 = new TH1F("fHOutCentV0M_qual1","fHOutCentV0M_qual1; Centrality V0",505,0,101);
 991     fHOutCentTRKqual1 = new TH1F("fHOutCentTRK_qual1","fHOutCentTRK_qual1; Centrality TPC",505,0,101);
 992     fHOutCentCL1qual1 = new TH1F("fHOutCentCL1_qual1","fHOutCentCL1_qual1; Centrality SPD outer",505,0,101);
 993     fHOutMultV0MvsCL1qual1 = new TH2F("fHOutMultV0MvsCL1_qual1","fHOutMultV0MvsCL1_qual1; Multiplicity V0; Multiplicity SPD outer",2500,0,25000,700,0,7000);
 994     fHOutMultV0MvsTRKqual1 = new TH2F("fHOutMultV0MvsTRK_qual1","fHOutMultV0MvsTRK_qual1; Multiplicity V0; Multiplicity TPC",2500,0,25000,400,0,4000);
 995     fHOutMultTRKvsCL1qual1 = new TH2F("fHOutMultTRKvsCL1_qual1","fHOutMultTRKvsCL1_qual1; Multiplicity TPC; Multiplicity SPD outer",400,0,4000,700,0,7000);
 996
 997     fHOutCentV0Mqual2 = new TH1F("fHOutCentV0M_qual2","fHOutCentV0M_qual2; Centrality V0",505,0,101);
 998     fHOutCentTRKqual2 = new TH1F("fHOutCentTRK_qual2","fHOutCentTRK_qual2; Centrality TPC",505,0,101);
 999     fHOutCentCL1qual2 = new TH1F("fHOutCentCL1_qual2","fHOutCentCL1_qual2; Centrality SPD outer",505,0,101);
1000     fHOutMultV0MvsCL1qual2 = new TH2F("fHOutMultV0MvsCL1_qual2","fHOutMultV0MvsCL1_qual2; Multiplicity V0; Multiplicity SPD outer",2500,0,25000,700,0,7000);
1001     fHOutMultV0MvsTRKqual2 = new TH2F("fHOutMultV0MvsTRK_qual2","fHOutMultV0MvsTRK_qual2; Multiplicity V0; Multiplicity TPC",2500,0,25000,400,0,4000);
1002     fHOutMultTRKvsCL1qual2 = new TH2F("fHOutMultTRKvsCL1_qual2","fHOutMultTRKvsCL1_qual2; Multiplicity TPC; Multiplicity SPD outer",400,0,4000,700,0,7000);
1003
1004     fHOutQuality = new TH1F("fHOutQuality", "fHOutQuality", 100,-0.5,99.5);
1005     fHOutVertex  = new TH1F("fHOutVertex", "fHOutVertex", 100,-20,20);
1006     fHOutVertexT0  = new TH1F("fHOutVertexT0", "fHOutVertexT0", 100,-20,20);
1007
1008     fOutputList->Add(fHOutCentV0M);
1009     fOutputList->Add(fHOutCentV0A);
1010     fOutputList->Add(fHOutCentV0A0);
1011     fOutputList->Add(fHOutCentV0A123);
1012     fOutputList->Add(fHOutCentV0C);
1013     fOutputList->Add(fHOutCentV0MEq);
1014     fOutputList->Add(fHOutCentV0AEq);
1015     fOutputList->Add(fHOutCentV0CEq);
1016     fOutputList->Add(fHOutCentV0MCVHN);
1017     fOutputList->Add(fHOutCentV0MCVLN);
1018     fOutputList->Add(fHOutCentV0MCVHNinMB);
1019     fOutputList->Add(fHOutCentV0MCVLNinMB);
1020     fOutputList->Add(fHOutCentV0MCCENT);
1021     fOutputList->Add(fHOutCentV0MCSEMI);
1022     fOutputList->Add(fHOutCentV0MCCENTinMB);
1023     fOutputList->Add(fHOutCentV0MCSEMIinMB);
1024     fOutputList->Add(fHOutCentV0MMSL);
1025     fOutputList->Add(fHOutCentV0MMSH);
1026     fOutputList->Add(fHOutCentV0MMUL);
1027     fOutputList->Add(fHOutCentV0MMLL);
1028     fOutputList->Add(fHOutCentV0MEJE);
1029     fOutputList->Add(fHOutCentV0MEGA);
1030     fOutputList->Add(fHOutCentV0MPHS);
1031     fOutputList->Add(fHOutCentV0MMSLinMB);
1032     fOutputList->Add(fHOutCentV0MMSHinMB);
1033     fOutputList->Add(fHOutCentV0MMULinMB);
1034     fOutputList->Add(fHOutCentV0MMLLinMB);
1035     fOutputList->Add(fHOutCentV0MEJEinMB);
1036     fOutputList->Add(fHOutCentV0MEGAinMB);
1037     fOutputList->Add(fHOutCentV0MPHSinMB);
1038     fOutputList->Add(fHOutCentFMD);
1039     fOutputList->Add(fHOutCentTRK);
1040     fOutputList->Add(fHOutCentTKL);
1041     fOutputList->Add(fHOutCentCL0);
1042     fOutputList->Add(fHOutCentCL1);
1043     fOutputList->Add(fHOutCentCND);
1044     fOutputList->Add(fHOutCentNPA);
1045     fOutputList->Add(fHOutCentZNA);
1046     fOutputList->Add(fHOutCentZNC);
1047     fOutputList->Add(fHOutCentZPA);
1048     fOutputList->Add(fHOutCentZPC);
1049     fOutputList->Add(fHOutCentV0MvsFMD);
1050     fOutputList->Add(fHOutCentTKLvsV0M);
1051     fOutputList->Add(fHOutCentZEMvsZDC);
1052     fOutputList->Add(fHOutCentV0MvsCentCL1);
1053     fOutputList->Add(fHOutCentV0MvsCentTRK);
1054     fOutputList->Add(fHOutCentTRKvsCentCL1);
1055     fOutputList->Add(fHOutCentV0MvsCentZDC);
1056     fOutputList->Add(fHOutCentV0AvsCentV0C);
1057     fOutputList->Add(fHOutCentV0AvsCentTRK);
1058     fOutputList->Add(fHOutCentV0AvsCentCND);
1059     fOutputList->Add(fHOutCentV0AvsCentCL1);
1060     fOutputList->Add(fHOutCentV0CvsCentTRK);
1061     fOutputList->Add(fHOutCentV0CvsCentCND);
1062     fOutputList->Add(fHOutCentV0CvsCentCL1);
1063     fOutputList->Add(fHOutCentNPAvsCentV0A);
1064     fOutputList->Add(fHOutCentNPAvsCentV0C);
1065     fOutputList->Add(fHOutCentNPAvsCentTRK);
1066     fOutputList->Add(fHOutCentNPAvsCentCND);
1067     fOutputList->Add(fHOutCentNPAvsCentCL1);
1068     fOutputList->Add(fHOutCentZNAvsCentV0A);
1069     fOutputList->Add(fHOutCentZNAvsCentV0C);
1070     fOutputList->Add(fHOutCentZNAvsCentTRK);
1071     fOutputList->Add(fHOutCentZNAvsCentCND);
1072     fOutputList->Add(fHOutCentZNAvsCentCL1);
1073     fOutputList->Add(fHOutCentZNAvsCentZPA);
1074
1075     fOutputList->Add(fHOutMultV0AC);
1076     fOutputList->Add(fHOutMultV0M);
1077     fOutputList->Add(fHOutMultV0A);
1078     fOutputList->Add(fHOutMultV0A0);
1079     fOutputList->Add(fHOutMultV0A123);
1080     fOutputList->Add(fHOutMultV0C);
1081     fOutputList->Add(fHOutMultV0MEq);
1082     fOutputList->Add(fHOutMultV0AEq);
1083     fOutputList->Add(fHOutMultV0CEq);
1084     fOutputList->Add(fHOutMultV0Mnc);
1085     fOutputList->Add(fHOutMultV0Anc);
1086     fOutputList->Add(fHOutMultV0Cnc);
1087     fOutputList->Add(fHOutMultV0O);
1088     fOutputList->Add(fHOutMultV0Cells) ;
1089     fOutputList->Add(fHOutMultFMD);
1090     fOutputList->Add(fHOutMultTRK);
1091     fOutputList->Add(fHOutMultTKL);
1092     fOutputList->Add(fHOutMultCL0);
1093     fOutputList->Add(fHOutMultCL1);
1094     fOutputList->Add(fHOutMultCND);
1095     fOutputList->Add(fHOutMultNPA);
1096     fOutputList->Add(fHOutMultZNA);
1097     fOutputList->Add(fHOutMultZNC);
1098     fOutputList->Add(fHOutMultZPA);
1099     fOutputList->Add(fHOutMultZPC);
1100     fOutputList->Add(fHOutMultV0MvsZDN);
1101     fOutputList->Add(fHOutMultZEMvsZDN);
1102     fOutputList->Add(fHOutMultV0MvsZDC);
1103     fOutputList->Add(fHOutMultZEMvsZDC);
1104     fOutputList->Add(fHOutMultZEMvsZDCw);
1105     fOutputList->Add(fHOutMultV0MvsCL1);
1106     fOutputList->Add(fHOutMultV0MvsTRK);
1107     fOutputList->Add(fHOutMultTRKvsCL1);
1108     fOutputList->Add(fHOutMultV0MvsV0O);
1109     fOutputList->Add(fHOutMultV0OvsCL1);
1110     fOutputList->Add(fHOutMultV0OvsTRK);
1111     fOutputList->Add(fHOutMultCL1vsTKL);
1112     fOutputList->Add(fHOutMultZNAvsZPA);
1113     fOutputList->Add(fHOutCentV0Mqual1);
1114     fOutputList->Add(fHOutCentTRKqual1);
1115     fOutputList->Add(fHOutCentCL1qual1);
1116     fOutputList->Add(fHOutMultV0MvsCL1qual1);
1117     fOutputList->Add(fHOutMultV0MvsTRKqual1);
1118     fOutputList->Add(fHOutMultTRKvsCL1qual1);
1119     fOutputList->Add(fHOutCentV0Mqual2);
1120     fOutputList->Add(fHOutCentTRKqual2);
1121     fOutputList->Add(fHOutCentCL1qual2);
1122     fOutputList->Add(fHOutMultV0MvsCL1qual2);
1123     fOutputList->Add(fHOutMultV0MvsTRKqual2);
1124     fOutputList->Add(fHOutMultTRKvsCL1qual2);
1125     fOutputList->Add(fHOutQuality);
1126     fOutputList->Add(fHOutVertex);
1127     fOutputList->Add(fHOutVertexT0);
1128
1129     PostData(1, fOutputList);
1130   }
1131
1132   fTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
1133   fEsdTrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011();
1134   fEsdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kOff);
1135   // Add SPD requirement
1136   fEsdTrackCutsExtra1 = new AliESDtrackCuts("SPD", "Require 1 cluster in SPD");
1137   fEsdTrackCutsExtra1->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1138   // Add SDD requirement
1139   fEsdTrackCutsExtra2 = new AliESDtrackCuts("SDD", "Require 1 cluster in first layer SDD");
1140   fEsdTrackCutsExtra2->SetClusterRequirementITS(AliESDtrackCuts::kSDD,AliESDtrackCuts::kFirst);
1141 }
1142
1143 //________________________________________________________________________
1144 void AliCentralitySelectionTask::UserExec(Option_t */*option*/)
1145 {
1146   // Execute analysis for current event:
1147   if(fDebug>1) printf(" **** AliCentralitySelectionTask::UserExec() \n");
1148
1149   Int_t    runType = 0;             // 0:PbPb, 1:pPb or Pbp
1150
1151   Float_t  zncEnergy = 0.;          //  ZNC Energy
1152   Float_t  zpcEnergy = 0.;          //  ZPC Energy
1153   Float_t  znaEnergy = 0.;          //  ZNA Energy
1154   Float_t  zpaEnergy = 0.;          //  ZPA Energy
1155   Float_t  zem1Energy = 0.;         //  ZEM1 Energy
1156   Float_t  zem2Energy = 0.;         //  ZEM2 Energy
1157   Bool_t   zdcEnergyCal = kFALSE;   // if zdc is calibrated (in pass2)
1158   Double_t znaTower = 0.;           // common PMT of ZNA
1159   Double_t zncTower = 0.;           // common PMT of ZNC
1160   Double_t zpaTower = 0.;           // common PMT of ZPA
1161   Double_t zpcTower = 0.;           // common PMT of ZPC
1162   Bool_t   znaFired = kFALSE;
1163   Bool_t   zncFired = kFALSE;
1164   Bool_t   zpaFired = kFALSE;
1165   Bool_t   zpcFired = kFALSE;
1166
1167   Int_t    nTracks = 0;             //  no. tracks
1168   Int_t    nTracklets = 0;          //  no. tracklets
1169   Int_t    nClusters[6] = {0};      //  no. clusters on 6 ITS layers
1170   Int_t    nChips[2] = {0,0};       //  no. chips on 2 SPD layers
1171   Float_t  spdCorr =0;              //  corrected spd2 multiplicity
1172   Int_t    multCND = 0;             //  no. tracks (candle condition)
1173
1174   Float_t  multV0A0 = 0;            //  multiplicity from V0 reco side A (ring 0)
1175   Float_t  multV0A123 = 0;          //  multiplicity from V0 reco side A (ring 1-2-3)
1176   Float_t  multV0A  = 0;            //  multiplicity from V0 reco side A
1177   Float_t  multV0C  = 0;            //  multiplicity from V0 reco side C
1178   Float_t  multV0AEq  = 0;          //  multiplicity from V0 reco side A
1179   Float_t  multV0CEq  = 0;          //  multiplicity from V0 reco side C
1180   Float_t  multV0A0Corr  = 0;       //  multiplicity from V0 reco side A (ring 0)
1181   Float_t  multV0A123Corr  = 0;     //  multiplicity from V0 reco side A (ring 1-2-3)
1182   Float_t  multV0ACorr  = 0;        //  multiplicity from V0 reco side A
1183   Float_t  multV0CCorr  = 0;        //  multiplicity from V0 reco side C
1184   Short_t  multV0AOnline  = 0;      //  multiplicity from V0 reco side A
1185   Short_t  multV0COnline  = 0;      //  multiplicity from V0 reco side C
1186   Float_t  v0Corr = 0;              // corrected V0 multiplicity (used for MC)
1187   Int_t nV0A = 0;
1188   Int_t nV0C = 0;
1189
1190   Float_t  multFMDA = 0;            //  multiplicity from FMD on detector A
1191   Float_t  multFMDC = 0;            //  multiplicity from FMD on detector C
1192
1193   Float_t zvtx =0;                  // z-vertex SPD
1194   Int_t zvtxNcont =0;               // contributors to z-vertex SPD
1195
1196   Float_t zvtxT0 =0;                // z-vertex T0
1197
1198   Int_t Npart =0;                   // N. of participants (true MC)
1199
1200   AliCentrality *esdCent = 0;
1201
1202   AliVEvent *event = InputEvent();
1203   AliESDEvent *esd = dynamic_cast<AliESDEvent*>(event);
1204   AliAODEvent *aod =  dynamic_cast<AliAODEvent*>(event);
1205   if(fAnalysisInput.CompareTo("ESD")==0){
1206     if (!esd) {
1207       AliError("No ESD Event");
1208       return;
1209     }
1210   } else if(fAnalysisInput.CompareTo("AOD")==0){
1211     if (!aod) {
1212       AliError("No AOD Event");
1213       return;
1214     }
1215   }
1216   LoadBranches();
1217
1218   if (SetupRun(event)<0) {
1219     AliError("Centrality File not available for this run");
1220     return;
1221   }
1222
1223   if (esd) {
1224     if (strcmp(esd->GetESDRun()->GetBeamType(), "A-A") == 0) runType=0;
1225     else runType=1;
1226     if (strcmp(esd->GetESDRun()->GetBeamType(), "p-p") == 0) runType=2;
1227   } else {
1228     Int_t runNumber = event->GetRunNumber();
1229     if ((runNumber >= 136851 && runNumber <= 139517) ||  // LHC10h
1230         (runNumber >= 166529 && runNumber <= 170593))    // LHC11h
1231       runType=0;
1232   }
1233
1234   esdCent = event->GetCentrality();
1235
1236   // ***** Vertex Info
1237   if (esd) {
1238     const AliESDVertex* vtxESD = esd->GetPrimaryVertexSPD();
1239     zvtx        = vtxESD->GetZ();
1240     zvtxNcont   = vtxESD->GetNContributors();
1241   } else {
1242     const AliAODVertex* spdVtx = aod->GetPrimaryVertexSPD();
1243     zvtx        = spdVtx->GetZ();
1244     zvtxNcont   = spdVtx->GetNContributors();
1245   }
1246
1247     // ***** V0 info
1248   AliVVZERO* esdV0 = event->GetVZEROData();
1249   if (!esdV0) {
1250     AliError("AliVVZERO not available");
1251     return;
1252   }
1253
1254   multV0A=esdV0->GetMTotV0A();
1255   multV0C=esdV0->GetMTotV0C();
1256
1257   // excluding innermost V0A ring
1258   for (int i=32; i<40; i++)
1259     multV0A0 += esdV0->GetMultiplicity(i);
1260   for (int i=40; i<64; i++)
1261     multV0A123 += esdV0->GetMultiplicity(i);
1262
1263   multV0A0Corr   = AliESDUtils::GetCorrV0A0(multV0A0,zvtx);
1264   multV0A123Corr = AliESDUtils::GetCorrV0A(multV0A123,zvtx);
1265   multV0ACorr    = AliESDUtils::GetCorrV0A(multV0A,zvtx);
1266   multV0CCorr    = AliESDUtils::GetCorrV0C(multV0C,zvtx);
1267
1268   v0Corr = multV0A+multV0C; // Todo: C.L. not clear why here we do not use the sum of the corrected values?
1269
1270   multV0AOnline=esdV0->GetTriggerChargeA();
1271   multV0COnline=esdV0->GetTriggerChargeC();
1272
1273   // Count V0 flags
1274   for(Int_t i = 0; i < 32; ++i) {
1275     if (esdV0->GetBBFlag(i)) nV0C++;
1276     if (esdV0->GetBBFlag(i+32)) nV0A++;
1277   }
1278
1279   // Equalized signals
1280   multV0AEq=0.;
1281   multV0CEq=0.;
1282   if (esd) { // only for ESD?
1283     for(Int_t iCh = 32; iCh < 64; ++iCh) {
1284       Double_t mult = esd->GetVZEROEqMultiplicity(iCh);
1285       multV0AEq += mult;
1286     }
1287     for(Int_t iCh = 0; iCh < 32; ++iCh) {
1288       Double_t mult = esd->GetVZEROEqMultiplicity(iCh);
1289       multV0CEq += mult;
1290     }
1291   }
1292   else {
1293     for(Int_t iCh = 32; iCh < 64; ++iCh) {
1294       Double_t mult = aod->GetVZEROEqMultiplicity(iCh);
1295       multV0AEq += mult;
1296     }
1297     for(Int_t iCh = 0; iCh < 32; ++iCh) {
1298       Double_t mult = aod->GetVZEROEqMultiplicity(iCh);
1299       multV0CEq += mult;
1300     }
1301   }
1302   Bool_t kT0BB = kFALSE;
1303   if (esd) {
1304     // ***** T0 info
1305     const AliESDTZERO* esdT0 = esd->GetESDTZERO();
1306     if (!esdT0)
1307       {
1308         AliError("AliESDTZERO not available");
1309         return;
1310       }
1311     Int_t trig=esdT0->GetT0Trig();
1312     if(trig&1) kT0BB=kTRUE;
1313     zvtxT0=esdT0->GetT0zVertex();
1314   } else {
1315     const AliAODTZERO* esdT0 = aod->GetTZEROData();
1316     if (!esdT0)
1317       {
1318         AliError("AliAODTZERO not available");
1319         return;
1320       }
1321     Int_t trig=1;//esdT0->GetT0Trig(); //* Todo: C.L. This info is not in AOD? */
1322     if(trig&1) kT0BB=kTRUE;
1323     zvtxT0=esdT0->GetT0zVertex();
1324   }
1325
1326   // ***** Trigger info
1327   fIsSelected = ((esdV0->GetV0ADecision()==1) && (esdV0->GetV0CDecision()==1));
1328   TString trigStr;
1329   if (esd)
1330     trigStr = esd->GetFiredTriggerClasses();
1331   else
1332     trigStr = aod->GetFiredTriggerClasses();
1333
1334   fMB=kFALSE;
1335   fCVHN=kFALSE; fCVLN=kFALSE; fCCENT=kFALSE; fCSEMI=kFALSE;
1336   fMSL=kFALSE;  fMSH=kFALSE;  fMUL=kFALSE;   fMLL=kFALSE;
1337   fEJE=kFALSE;  fEGA=kFALSE;  fPHS=kFALSE;
1338
1339   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CPBI")) && (fIsSelected))
1340     fMB=kTRUE;
1341   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CVHN")) && (fIsSelected))
1342     fCVHN=kTRUE;
1343   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CVLN")) && (fIsSelected))
1344     fCVLN=kTRUE;
1345   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CCENT")) && (fIsSelected))
1346     fCCENT=kTRUE;
1347   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CSEMI")) && (fIsSelected))
1348     fCSEMI=kTRUE;
1349
1350   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CPBI1MSL")) && (fIsSelected))
1351     fMSL=kTRUE;
1352   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CPBI1MSH")) && (fIsSelected))
1353     fMSH=kTRUE;
1354   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CPBI1MUL")) && (fIsSelected))
1355     fMUL=kTRUE;
1356   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CPBI1MLL")) && (fIsSelected))
1357     fMLL=kTRUE;
1358   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CPBI1EJE")) && (fIsSelected))
1359     fEJE=kTRUE;
1360   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CPBI1EGA")) && (fIsSelected))
1361     fEGA=kTRUE;
1362   if ( (trigStr.Contains("-B-")) &&  (trigStr.Contains("CPBI1PHS")) && (fIsSelected))
1363     fPHS=kTRUE;
1364
1365   fCVHNbit=kFALSE;    fCVLNbit=kFALSE;       fCCENTbit=kFALSE;    fCSEMIbit=kFALSE;
1366   if (esdV0->GetTriggerBits() & (1<<8))
1367     fCVHNbit=kTRUE;
1368   if (esdV0->GetTriggerBits() & (1<<6))
1369     fCVLNbit=kTRUE;
1370
1371   if (kT0BB && fCVHNbit)
1372     fCCENTbit=kTRUE;
1373   if (kT0BB && fCVLNbit)
1374     fCSEMIbit=kTRUE;
1375
1376   if (esd) {
1377     // ***** CB info (tracklets, clusters, chips)
1378     //nTracks    = event->GetNumberOfTracks();
1379     nTracks    = fTrackCuts ? (Short_t)fTrackCuts->GetReferenceMultiplicity(esd,kTRUE):-1;
1380   } else {
1381     AliAODHeader *h = aod->GetHeader();
1382     nTracks    = h!=0 ? (Short_t)h->GetTPConlyRefMultiplicity():-1;
1383   }
1384
1385   if (esd) {
1386     Short_t nTrTPCcandle = 0;
1387     for (Int_t iTracks = 0; iTracks < esd->GetNumberOfTracks(); iTracks++) {
1388
1389       AliESDtrack* track = esd->GetTrack(iTracks);
1390       if (!track) continue;
1391
1392       if (! fEsdTrackCuts->IsSelected(track) )continue;
1393
1394       if (fEsdTrackCutsExtra1 && fEsdTrackCutsExtra2 &&
1395           !fEsdTrackCutsExtra1->IsSelected(track) &&
1396           !fEsdTrackCutsExtra2->IsSelected(track)) continue;
1397
1398       if (track->Pt() > 0.4 && TMath::Abs(track->Eta()) < 0.9)  nTrTPCcandle++;
1399     }
1400     multCND = nTrTPCcandle;
1401   } else {
1402     Short_t nTrTPCcandle = 0;
1403     for (Int_t iTracks = 0; iTracks < aod->GetNumberOfTracks(); iTracks++) {
1404
1405       AliAODTrack* track = aod->GetTrack(iTracks);
1406
1407       if (!track) continue;
1408       if (!track->TestFilterBit(1<<5) &&
1409           !track->TestFilterBit(1<<6)) continue;
1410
1411       if (track->Pt() > 0.4 && TMath::Abs(track->Eta()) < 0.9)  nTrTPCcandle++;
1412     }
1413     multCND = nTrTPCcandle;
1414   }
1415
1416   if (esd) {
1417     const AliMultiplicity *mult = esd->GetMultiplicity();
1418     nTracklets = mult->GetNumberOfTracklets();
1419
1420     for(Int_t ilay=0; ilay<6; ilay++){
1421       nClusters[ilay] = mult->GetNumberOfITSClusters(ilay);
1422     }
1423
1424     for(Int_t ilay=0; ilay<2; ilay++){
1425       nChips[ilay] = mult->GetNumberOfFiredChips(ilay);
1426     }
1427   } else {
1428     AliAODTracklets *mult = aod->GetTracklets();
1429     nTracklets = mult->GetNumberOfTracklets();
1430     AliAODHeader *h = aod->GetHeader();
1431     for(Int_t ilay=0; ilay<6; ilay++){
1432       nClusters[ilay] = h->GetNumberOfITSClusters(ilay);
1433     }
1434   }
1435   spdCorr = AliESDUtils::GetCorrSPD2(nClusters[1],zvtx);
1436
1437   if (esd) {
1438     // ***** FMD info
1439     AliESDFMD *fmd = esd->GetFMDData();
1440     Float_t totalMultA = 0;
1441     Float_t totalMultC = 0;
1442     const Float_t fFMDLowCut = 0.4;
1443
1444     for(UShort_t det=1;det<=3;det++) {
1445       Int_t nRings = (det==1 ? 1 : 2);
1446       for (UShort_t ir = 0; ir < nRings; ir++) {
1447         Char_t   ring = (ir == 0 ? 'I' : 'O');
1448         UShort_t nsec = (ir == 0 ? 20  : 40);
1449         UShort_t nstr = (ir == 0 ? 512 : 256);
1450         for(UShort_t sec =0; sec < nsec;  sec++)  {
1451           for(UShort_t strip = 0; strip < nstr; strip++) {
1452
1453             Float_t fmdMult = fmd->Multiplicity(det,ring,sec,strip);
1454             if(fmdMult == 0 || fmdMult == AliESDFMD::kInvalidMult) continue;
1455
1456             Float_t nParticles=0;
1457
1458             if(fmdMult > fFMDLowCut) {
1459               nParticles = 1.;
1460             }
1461
1462             if (det<3) totalMultA = totalMultA + nParticles;
1463             else totalMultC = totalMultC + nParticles;
1464
1465           }
1466         }
1467       }
1468     }
1469     multFMDA = totalMultA;
1470     multFMDC = totalMultC;
1471   }
1472
1473   if (esd) {
1474     // ***** ZDC info
1475     AliESDZDC *esdZDC = esd->GetESDZDC();
1476     zdcEnergyCal = esdZDC->AliESDZDC::TestBit(AliESDZDC::kEnergyCalibratedSignal);
1477     if (zdcEnergyCal) {
1478       zncEnergy = (Float_t) (esdZDC->GetZDCN1Energy());
1479       zpcEnergy = (Float_t) (esdZDC->GetZDCP1Energy());
1480       znaEnergy = (Float_t) (esdZDC->GetZDCN2Energy());
1481       zpaEnergy = (Float_t) (esdZDC->GetZDCP2Energy());
1482     } else {
1483       zncEnergy = (Float_t) (esdZDC->GetZDCN1Energy())/8.;
1484       zpcEnergy = (Float_t) (esdZDC->GetZDCP1Energy())/8.;
1485       znaEnergy = (Float_t) (esdZDC->GetZDCN2Energy())/8.;
1486       zpaEnergy = (Float_t) (esdZDC->GetZDCP2Energy())/8.;
1487     }
1488     zem1Energy = (Float_t) (esdZDC->GetZDCEMEnergy(0))/8.;
1489     zem2Energy = (Float_t) (esdZDC->GetZDCEMEnergy(1))/8.;
1490
1491     for (Int_t j = 0; j < 4; ++j)
1492       if (esdZDC->GetZDCTDCData(12,j) != 0)
1493         znaFired = kTRUE;
1494
1495     for (Int_t j = 0; j < 4; ++j)
1496       if (esdZDC->GetZDCTDCData(10,j) != 0)
1497         zncFired = kTRUE;
1498
1499     for (Int_t j = 0; j < 4; ++j)
1500       if (esdZDC->GetZDCTDCData(13,j) != 0)
1501         zpaFired = kTRUE;
1502
1503     for (Int_t j = 0; j < 4; ++j)
1504       if (esdZDC->GetZDCTDCData(11,j) != 0)
1505         zpcFired = kTRUE;
1506
1507     const Double_t *ZNAtower = esdZDC->GetZN2TowerEnergy();
1508     const Double_t *ZNCtower = esdZDC->GetZN1TowerEnergy();
1509     const Double_t *ZPAtower = esdZDC->GetZP2TowerEnergy();
1510     const Double_t *ZPCtower = esdZDC->GetZP1TowerEnergy();
1511     if (znaFired) znaTower = ZNAtower[0];
1512     if (zncFired) zncTower = ZNCtower[0];
1513     if (zpaFired) zpaTower = ZPAtower[0];
1514     if (zpcFired) zpcTower = ZPCtower[0];
1515
1516   } else {
1517     AliAODHeader *h = aod->GetHeader();
1518     zncEnergy  = (Float_t) (h->GetZDCN1Energy());
1519     zpcEnergy  = (Float_t) (h->GetZDCP1Energy());
1520     znaEnergy  = (Float_t) (h->GetZDCN2Energy());
1521     zpaEnergy  = (Float_t) (h->GetZDCP2Energy());
1522     zem1Energy = (Float_t) (h->GetZDCEMEnergy(0))/8.; //Todo: C.L. Should we devide here by 8? It is done in the ESD case!
1523     zem2Energy = (Float_t) (h->GetZDCEMEnergy(1))/8.;
1524
1525     AliAODZDC *aodZDC = aod->GetZDCData();
1526     const Double_t *ZNAtower = aodZDC->GetZNATowerEnergy();
1527     const Double_t *ZNCtower = aodZDC->GetZNCTowerEnergy();
1528     const Double_t *ZPAtower = aodZDC->GetZPATowerEnergy();
1529     const Double_t *ZPCtower = aodZDC->GetZPCTowerEnergy();
1530     znaTower = ZNAtower[0];
1531     zncTower = ZNCtower[0];
1532     zpaTower = ZPAtower[0];
1533     zpcTower = ZPCtower[0];
1534
1535     znaFired = kFALSE; // trick because info is not stored in AOD
1536     zncFired = kFALSE; // trick because info is not stored in AOD
1537     zpaFired = kFALSE; // trick because info is not stored in AOD
1538     zpcFired = kFALSE; // trick because info is not stored in AOD
1539     if (esdCent->GetCentralityPercentile("ZNA") != 101) znaFired = kTRUE;
1540     if (esdCent->GetCentralityPercentile("ZNC") != 101) zncFired = kTRUE;
1541     if (esdCent->GetCentralityPercentile("ZPA") != 101) zpaFired = kTRUE;
1542     if (esdCent->GetCentralityPercentile("ZPC") != 101) zpcFired = kTRUE;
1543   }
1544
1545   if (esd) {
1546     // ***** MC info
1547     AliAnalysisManager* anMan = AliAnalysisManager::GetAnalysisManager();
1548     AliMCEventHandler* eventHandler = (AliMCEventHandler*)anMan->GetMCtruthEventHandler();
1549     AliStack*    stack=0;
1550     AliMCEvent*  mcEvent=0;
1551     if (fIsMCInput && eventHandler && (mcEvent=eventHandler->MCEvent()) && (stack=mcEvent->Stack())) {
1552       AliGenHijingEventHeader* hHijing=0;
1553       AliGenDPMjetEventHeader* dpmHeader=0;
1554
1555       AliGenEventHeader* mcGenH = mcEvent->GenEventHeader();
1556       if (mcGenH->InheritsFrom(AliGenHijingEventHeader::Class()))
1557         hHijing = (AliGenHijingEventHeader*)mcGenH;
1558       else if (mcGenH->InheritsFrom(AliGenCocktailEventHeader::Class())) {
1559         TList* headers = ((AliGenCocktailEventHeader*)mcGenH)->GetHeaders();
1560         //headers->Print();
1561         hHijing = dynamic_cast<AliGenHijingEventHeader*>(headers->FindObject("Hijing"));
1562         if (!hHijing) hHijing = dynamic_cast<AliGenHijingEventHeader*>(headers->FindObject("Hijing pPb_0"));
1563         if (!hHijing) hHijing = dynamic_cast<AliGenHijingEventHeader*>(headers->FindObject("Hijing_0"));
1564       }
1565       else if (mcGenH->InheritsFrom(AliGenDPMjetEventHeader::Class())) {
1566         dpmHeader = (AliGenDPMjetEventHeader*)mcGenH;
1567       }
1568       if(hHijing)   Npart = hHijing->ProjectileParticipants()+hHijing->TargetParticipants();
1569       if(dpmHeader) Npart = dpmHeader->ProjectileParticipants()+ dpmHeader->TargetParticipants();
1570     }
1571   }
1572
1573   // ***** Scaling for MC
1574   if (fIsMCInput) {
1575     fUseScaling=kFALSE;
1576     v0Corr   = Short_t((multV0A+multV0C)  * fV0MScaleFactorMC);
1577     multV0A0 = multV0A0 * fV0MScaleFactorMC;
1578     multV0A  = multV0A * fV0MScaleFactorMC;
1579     multV0C  = multV0C * fV0MScaleFactorMC;
1580   }
1581   // ***** Scaling for Data
1582   if (fUseScaling) {
1583     v0Corr  = Short_t(v0Corr /   fV0MScaleFactor);
1584     spdCorr = spdCorr / fSPDScaleFactor;
1585     nTracks = Int_t(nTracks /   fTPCScaleFactor);
1586   }
1587
1588   // ***** Centrality Selection
1589   if(fHtempV0M) fCentV0M = fHtempV0M->GetBinContent(fHtempV0M->FindBin((v0Corr)));
1590   if(fHtempV0A) fCentV0A = fHtempV0A->GetBinContent(fHtempV0A->FindBin((multV0ACorr)));
1591   if(fHtempV0A0) fCentV0A0 = fHtempV0A0->GetBinContent(fHtempV0A0->FindBin((multV0A0Corr)));
1592   if(fHtempV0A123) fCentV0A123 = fHtempV0A123->GetBinContent(fHtempV0A123->FindBin((multV0A123Corr)));
1593   if(fHtempV0C) fCentV0C = fHtempV0C->GetBinContent(fHtempV0C->FindBin((multV0CCorr)));
1594   if(fHtempV0MEq) fCentV0MEq = fHtempV0MEq->GetBinContent(fHtempV0MEq->FindBin((multV0AEq+multV0CEq)));
1595   if(fHtempV0AEq) fCentV0AEq = fHtempV0AEq->GetBinContent(fHtempV0AEq->FindBin((multV0AEq)));
1596   if(fHtempV0CEq) fCentV0CEq = fHtempV0CEq->GetBinContent(fHtempV0CEq->FindBin((multV0CEq)));
1597   if(fHtempFMD) fCentFMD = fHtempFMD->GetBinContent(fHtempFMD->FindBin((multFMDA+multFMDC)));
1598   if(fHtempTRK) fCentTRK = fHtempTRK->GetBinContent(fHtempTRK->FindBin(nTracks));
1599   if(fHtempTKL) fCentTKL = fHtempTKL->GetBinContent(fHtempTKL->FindBin(nTracklets));
1600   if(fHtempCL0) fCentCL0 = fHtempCL0->GetBinContent(fHtempCL0->FindBin(nClusters[0]));
1601   if(fHtempCL1) fCentCL1 = fHtempCL1->GetBinContent(fHtempCL1->FindBin(spdCorr));
1602   if(fHtempCND) fCentCND = fHtempCND->GetBinContent(fHtempCND->FindBin(multCND));
1603   if(fHtempZNA) {
1604     if(znaFired) fCentZNA = fHtempZNA->GetBinContent(fHtempZNA->FindBin(znaTower));
1605     else fCentZNA = 101;
1606   }
1607   if(fHtempZNC) {
1608     if(zncFired) fCentZNC = fHtempZNC->GetBinContent(fHtempZNC->FindBin(zncTower));
1609     else fCentZNC = 101;
1610   }
1611   if(fHtempZPA) {
1612     if(znaFired) fCentZPA = fHtempZPA->GetBinContent(fHtempZPA->FindBin(zpaTower));
1613     else fCentZPA = 101;
1614   }
1615   if(fHtempZPC) {
1616     if(zpcFired) fCentZPC = fHtempZPC->GetBinContent(fHtempZPC->FindBin(zpcTower));
1617     else fCentZPC = 101;
1618   }
1619
1620
1621   if(fHtempV0MvsFMD) fCentV0MvsFMD = fHtempV0MvsFMD->GetBinContent(fHtempV0MvsFMD->FindBin((multV0A+multV0C)));
1622   if(fHtempTKLvsV0M) fCentTKLvsV0M = fHtempTKLvsV0M->GetBinContent(fHtempTKLvsV0M->FindBin(nTracklets));
1623   if(fHtempZEMvsZDC) fCentZEMvsZDC = fHtempZEMvsZDC->GetBinContent(fHtempZEMvsZDC->FindBin(zem1Energy+zem2Energy,zncEnergy+znaEnergy+zpcEnergy+zpaEnergy));
1624
1625   if(fHtempNPA) fCentNPA = fHtempNPA->GetBinContent(fHtempNPA->FindBin(Npart));
1626   if(fHtempV0Mtrue) fCentV0Mtrue = fHtempV0Mtrue->GetBinContent(fHtempV0Mtrue->FindBin((multV0ACorr+multV0CCorr)));
1627   if(fHtempV0Atrue) fCentV0Atrue = fHtempV0Atrue->GetBinContent(fHtempV0Atrue->FindBin((multV0ACorr)));
1628   if(fHtempV0Ctrue) fCentV0Ctrue = fHtempV0Ctrue->GetBinContent(fHtempV0Ctrue->FindBin((multV0CCorr)));
1629   if(fHtempV0MEqtrue) fCentV0MEqtrue = fHtempV0MEqtrue->GetBinContent(fHtempV0MEqtrue->FindBin((multV0AEq+multV0CEq)));
1630   if(fHtempV0AEqtrue) fCentV0AEqtrue = fHtempV0AEqtrue->GetBinContent(fHtempV0AEqtrue->FindBin((multV0AEq)));
1631   if(fHtempV0CEqtrue) fCentV0CEqtrue = fHtempV0CEqtrue->GetBinContent(fHtempV0CEqtrue->FindBin((multV0CEq)));
1632   if(fHtempFMDtrue) fCentFMDtrue = fHtempFMDtrue->GetBinContent(fHtempFMDtrue->FindBin((multFMDA+multFMDC)));
1633   if(fHtempTRKtrue) fCentTRKtrue = fHtempTRKtrue->GetBinContent(fHtempTRKtrue->FindBin(nTracks));
1634   if(fHtempTKLtrue) fCentTKLtrue = fHtempTKLtrue->GetBinContent(fHtempTKLtrue->FindBin(nTracklets));
1635   if(fHtempCL0true) fCentCL0true = fHtempCL0true->GetBinContent(fHtempCL0true->FindBin(nClusters[0]));
1636   if(fHtempCL1true) fCentCL1true = fHtempCL1true->GetBinContent(fHtempCL1true->FindBin(spdCorr));
1637   if(fHtempCNDtrue) fCentCNDtrue = fHtempCNDtrue->GetBinContent(fHtempCNDtrue->FindBin(multCND));
1638   if(fHtempZNAtrue) fCentZNAtrue = fHtempZNAtrue->GetBinContent(fHtempZNAtrue->FindBin(znaTower));
1639   if(fHtempZNCtrue) fCentZNCtrue = fHtempZNCtrue->GetBinContent(fHtempZNCtrue->FindBin(zncTower));
1640
1641
1642   // ***** Cleaning
1643   if (fUseCleaning) {
1644     fQuality=0;
1645
1646     // ***** vertex
1647     if (TMath::Abs(zvtx)>fZVCut || zvtxNcont<1) fQuality += 1;
1648
1649     // ***** outliers, skip in case of MC input
1650     if (!fIsMCInput) {
1651       // **** V0 vs SPD
1652       if (IsOutlierV0MSPD(spdCorr, v0Corr, int(fCentV0M))) fQuality  += 2;
1653       // ***** V0 vs TPC
1654       if (IsOutlierV0MTPC(nTracks, v0Corr, int(fCentV0M))) fQuality  += 4;
1655       // ***** V0 vs ZDC
1656       if (IsOutlierV0MZDC((zncEnergy+znaEnergy+zpcEnergy+zpaEnergy), v0Corr) &&
1657           (zdcEnergyCal==kFALSE) ) fQuality  += 8;
1658       if (IsOutlierV0MZDCECal((zncEnergy+znaEnergy+zpcEnergy+zpaEnergy), v0Corr) &&
1659           (zdcEnergyCal==kTRUE) ) fQuality  += 8;
1660     }
1661   } else {
1662     fQuality = 0;
1663   }
1664
1665
1666   if (esdCent) {
1667     if (aod&&(fDebug>1)) {
1668       Double_t v0m = esdCent->GetCentralityPercentile("V0M");
1669       Double_t cl1 = esdCent->GetCentralityPercentile("CL1");
1670       Double_t trk = esdCent->GetCentralityPercentile("TRK");
1671       Double_t cnd = esdCent->GetCentralityPercentile("CND");
1672       Double_t zna = esdCent->GetCentralityPercentile("ZNA");
1673       printf("AOD: v0m %.2f %.2f (%.2f) cl1 %.2f %.2f (%.2f) trk %.2f %.2f (%.2f) cnd %.2f %.2f (%.2f) zna %.2f %.2f (%.2f)\n",
1674              v0m, fCentV0M, fCentV0M!=0?v0m/fCentV0M:1, cl1, fCentCL1, fCentCL1!=0?cl1/fCentCL1:1, trk, fCentTRK,
1675              fCentTRK!=0?trk/fCentTRK:1, cnd, fCentCND, fCentCND!=0?cnd/fCentCND:1, zna, fCentZNA, fCentZNA!=0?zna/fCentZNA:1);
1676     }
1677     esdCent->SetQuality(fQuality);
1678     esdCent->SetCentralityV0M(fCentV0M);
1679     esdCent->SetCentralityV0A(fCentV0A);
1680     esdCent->SetCentralityV0A0(fCentV0A0);
1681     esdCent->SetCentralityV0A123(fCentV0A123);
1682     esdCent->SetCentralityV0C(fCentV0C);
1683     esdCent->SetCentralityV0MEq(fCentV0MEq);
1684     esdCent->SetCentralityV0AEq(fCentV0AEq);
1685     esdCent->SetCentralityV0CEq(fCentV0CEq);
1686     esdCent->SetCentralityFMD(fCentFMD);
1687     esdCent->SetCentralityTRK(fCentTRK);
1688     esdCent->SetCentralityTKL(fCentTKL);
1689     esdCent->SetCentralityCL0(fCentCL0);
1690     esdCent->SetCentralityCL1(fCentCL1);
1691     esdCent->SetCentralityCND(fCentCND);
1692     esdCent->SetCentralityNPA(fCentNPA);
1693     esdCent->SetCentralityZNA(fCentZNA);
1694     esdCent->SetCentralityZNC(fCentZNC);
1695     esdCent->SetCentralityZPA(fCentZPA);
1696     esdCent->SetCentralityZPC(fCentZPC);
1697     esdCent->SetCentralityV0MvsFMD(fCentV0MvsFMD);
1698     esdCent->SetCentralityTKLvsV0M(fCentTKLvsV0M);
1699     esdCent->SetCentralityZEMvsZDC(fCentZEMvsZDC);
1700   }
1701
1702   // filling QA histograms
1703   if (fFillHistos) {
1704
1705     if (fIsMCInput) { // fill histo with true centrality for simulations
1706       fCentV0M = fCentV0Mtrue;
1707       fCentV0A = fCentV0Atrue;
1708       fCentV0C = fCentV0Ctrue;
1709       fCentV0MEq = fCentV0MEqtrue;
1710       fCentV0AEq = fCentV0AEqtrue;
1711       fCentV0CEq = fCentV0CEqtrue;
1712       fCentFMD = fCentFMDtrue;
1713       fCentTRK = fCentTRKtrue;
1714       fCentTKL = fCentTKLtrue;
1715       fCentCL0 = fCentCL0true;
1716       fCentCL1 = fCentCL1true;
1717       fCentCND = fCentCNDtrue;
1718       fCentZNA = fCentZNAtrue;
1719       fCentZNC = fCentZNCtrue;
1720       fCentZPA = fCentZPAtrue;
1721       fCentZPC = fCentZPCtrue;
1722     }
1723
1724
1725     if ((fMB) && (abs(zvtx)<10))        fHOutMultCL1vsTKL->Fill(spdCorr,nTracklets);
1726
1727     if (fCVHN)   fHOutCentV0MCVHN->Fill(fCentV0M);
1728     if (fCVLN)   fHOutCentV0MCVLN->Fill(fCentV0M);
1729     if (fCCENT)  fHOutCentV0MCCENT->Fill(fCentV0M);
1730     if (fCSEMI)  fHOutCentV0MCSEMI->Fill(fCentV0M);
1731     if (fMSL) fHOutCentV0MMSL->Fill(fCentV0M);
1732     if (fMSH) fHOutCentV0MMSH->Fill(fCentV0M);
1733     if (fMUL) fHOutCentV0MMUL->Fill(fCentV0M);
1734     if (fMLL) fHOutCentV0MMLL->Fill(fCentV0M);
1735     if (fEJE) fHOutCentV0MEJE->Fill(fCentV0M);
1736     if (fEGA) fHOutCentV0MEGA->Fill(fCentV0M);
1737     if (fPHS) fHOutCentV0MPHS->Fill(fCentV0M);
1738
1739     if (((((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & AliVEvent::kMB) && (runType==0)) ||
1740         ((((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & AliVEvent::kINT7) && (runType==1)) ||
1741         ((((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & AliVEvent::kMB) && (runType==2)) ) { // fill the QA histograms only for MB events!
1742
1743       fHOutQuality->Fill(fQuality);
1744       fHOutVertex->Fill(zvtx);
1745       fHOutVertexT0->Fill(zvtxT0);
1746
1747       if (fQuality==0) {
1748         fHOutCentV0M->Fill(fCentV0M);
1749         fHOutCentV0A->Fill(fCentV0A);
1750         fHOutCentV0A0->Fill(fCentV0A0);
1751         fHOutCentV0A123->Fill(fCentV0A123);
1752         fHOutCentV0C->Fill(fCentV0C);
1753         fHOutCentV0MEq->Fill(fCentV0MEq);
1754         fHOutCentV0AEq->Fill(fCentV0AEq);
1755         fHOutCentV0CEq->Fill(fCentV0CEq);
1756
1757         if (fCVHNbit)  fHOutCentV0MCVHNinMB->Fill(fCentV0M);
1758         if (fCVLNbit)  fHOutCentV0MCVLNinMB->Fill(fCentV0M);
1759         if (fCCENTbit) fHOutCentV0MCCENTinMB->Fill(fCentV0M);
1760         if (fCSEMIbit) fHOutCentV0MCSEMIinMB->Fill(fCentV0M);
1761         if (fMSL) fHOutCentV0MMSLinMB->Fill(fCentV0M);
1762         if (fMSH) fHOutCentV0MMSHinMB->Fill(fCentV0M);
1763         if (fMUL) fHOutCentV0MMULinMB->Fill(fCentV0M);
1764         if (fMLL) fHOutCentV0MMLLinMB->Fill(fCentV0M);
1765         if (fEJE) fHOutCentV0MEJEinMB->Fill(fCentV0M);
1766         if (fEGA) fHOutCentV0MEGAinMB->Fill(fCentV0M);
1767         if (fPHS) fHOutCentV0MPHSinMB->Fill(fCentV0M);
1768
1769         fHOutCentFMD->Fill(fCentFMD);
1770         fHOutCentTRK->Fill(fCentTRK);
1771         fHOutCentTKL->Fill(fCentTKL);
1772         fHOutCentCL0->Fill(fCentCL0);
1773         fHOutCentCL1->Fill(fCentCL1);
1774         fHOutCentCND->Fill(fCentCND);
1775         fHOutCentNPA->Fill(fCentNPA);
1776         fHOutCentZNA->Fill(fCentZNA);
1777         fHOutCentZNC->Fill(fCentZNC);
1778         fHOutCentZPA->Fill(fCentZPA);
1779         fHOutCentZPC->Fill(fCentZPC);
1780         fHOutCentV0MvsFMD->Fill(fCentV0MvsFMD);
1781         fHOutCentTKLvsV0M->Fill(fCentTKLvsV0M);
1782         fHOutCentZEMvsZDC->Fill(fCentZEMvsZDC);
1783         fHOutCentV0MvsCentCL1->Fill(fCentV0M,fCentCL1);
1784         fHOutCentV0MvsCentTRK->Fill(fCentV0M,fCentTRK);
1785         fHOutCentTRKvsCentCL1->Fill(fCentTRK,fCentCL1);
1786         fHOutCentV0MvsCentZDC->Fill(fCentV0M,fCentZEMvsZDC);
1787         fHOutCentV0AvsCentV0C->Fill(fCentV0A,fCentV0C);
1788         fHOutCentV0AvsCentTRK->Fill(fCentV0A,fCentTRK);
1789         fHOutCentV0AvsCentCND->Fill(fCentV0A,fCentCND);
1790         fHOutCentV0AvsCentCL1->Fill(fCentV0A,fCentCL1);
1791         fHOutCentV0CvsCentTRK->Fill(fCentV0C,fCentTRK);
1792         fHOutCentV0CvsCentCND->Fill(fCentV0C,fCentCND);
1793         fHOutCentV0CvsCentCL1->Fill(fCentV0C,fCentCL1);
1794         fHOutCentNPAvsCentV0A->Fill(fCentNPA,fCentV0A);
1795         fHOutCentNPAvsCentV0C->Fill(fCentNPA,fCentV0C);
1796         fHOutCentNPAvsCentTRK->Fill(fCentNPA,fCentTRK);
1797         fHOutCentNPAvsCentCND->Fill(fCentNPA,fCentCND);
1798         fHOutCentNPAvsCentCL1->Fill(fCentNPA,fCentCL1);
1799         fHOutCentZNAvsCentV0A->Fill(fCentZNA,fCentV0A);
1800         fHOutCentZNAvsCentV0C->Fill(fCentZNA,fCentV0C);
1801         fHOutCentZNAvsCentTRK->Fill(fCentZNA,fCentTRK);
1802         fHOutCentZNAvsCentCND->Fill(fCentZNA,fCentCND);
1803         fHOutCentZNAvsCentCL1->Fill(fCentZNA,fCentCL1);
1804         fHOutCentZNAvsCentZPA->Fill(fCentZNA,fCentZPA);
1805
1806         fHOutMultV0AC->Fill(multV0A,multV0C);
1807         fHOutMultV0M->Fill(multV0ACorr+multV0CCorr);
1808         fHOutMultV0A->Fill(multV0ACorr);
1809         fHOutMultV0A0->Fill(multV0A0Corr);
1810         fHOutMultV0A123->Fill(multV0A123Corr);
1811         fHOutMultV0C->Fill(multV0CCorr);
1812         fHOutMultV0MEq->Fill(multV0AEq+multV0CEq);
1813         fHOutMultV0AEq->Fill(multV0AEq);
1814         fHOutMultV0CEq->Fill(multV0CEq);
1815         fHOutMultV0Mnc->Fill(multV0A+multV0C);
1816         fHOutMultV0Anc->Fill(multV0A);
1817         fHOutMultV0Cnc->Fill(multV0C);
1818         fHOutMultV0O->Fill(multV0AOnline+multV0COnline);
1819         fHOutMultV0Cells->Fill(nV0A,nV0C);
1820         fHOutMultFMD->Fill(multFMDA+multFMDC);
1821         fHOutMultTRK->Fill(nTracks);
1822         fHOutMultTKL->Fill(nTracklets);
1823         fHOutMultCL0->Fill(nClusters[0]);
1824         fHOutMultCL1->Fill(spdCorr);
1825         fHOutMultCND->Fill(multCND);
1826         fHOutMultNPA->Fill(Npart);
1827         if(znaFired)fHOutMultZNA->Fill(znaTower);
1828         if(zncFired)fHOutMultZNC->Fill(zncTower);
1829         if(zpaFired)fHOutMultZPA->Fill(zpaTower);
1830         if(zpcFired)fHOutMultZPC->Fill(zpcTower);
1831         fHOutMultZNAvsZPA->Fill(znaTower,zpaTower);
1832
1833         //fHOutMultV0MvsZDN->Fill(v0Corr,(zncEnergy+znaEnergy));
1834         fHOutMultV0MvsZDN->Fill(v0Corr,znaTower);
1835         fHOutMultZEMvsZDN->Fill((zem1Energy+zem2Energy),(zncEnergy+znaEnergy));
1836         fHOutMultV0MvsZDC->Fill(v0Corr,(zncEnergy+znaEnergy+zpcEnergy+zpaEnergy));
1837         fHOutMultZEMvsZDC->Fill((zem1Energy+zem2Energy),(zncEnergy+znaEnergy+zpcEnergy+zpaEnergy));
1838         fHOutMultZEMvsZDCw->Fill((zem1Energy+zem2Energy),(zncEnergy+znaEnergy+zpcEnergy+zpaEnergy),fCentV0M);
1839         fHOutMultV0MvsCL1->Fill(v0Corr,spdCorr);
1840         fHOutMultV0MvsTRK->Fill(v0Corr,nTracks);
1841         fHOutMultTRKvsCL1->Fill(nTracks,spdCorr);
1842         fHOutMultV0MvsV0O->Fill(v0Corr,(multV0AOnline+multV0COnline));
1843         fHOutMultV0OvsCL1->Fill((multV0AOnline+multV0COnline),spdCorr);
1844         fHOutMultV0OvsTRK->Fill((multV0AOnline+multV0COnline),nTracks);
1845       } else if (fQuality%2 == 0) {
1846         fHOutCentV0Mqual1->Fill(fCentV0M);
1847         fHOutCentTRKqual1->Fill(fCentTRK);
1848         fHOutCentCL1qual1->Fill(fCentCL1);
1849         fHOutMultV0MvsCL1qual1->Fill(v0Corr,spdCorr);
1850         fHOutMultV0MvsTRKqual1->Fill(v0Corr,nTracks);
1851         fHOutMultTRKvsCL1qual1->Fill(nTracks,spdCorr);
1852       } else {
1853         fHOutCentV0Mqual2->Fill(fCentV0M);
1854         fHOutCentTRKqual2->Fill(fCentTRK);
1855         fHOutCentCL1qual2->Fill(fCentCL1);
1856         fHOutMultV0MvsCL1qual2->Fill(v0Corr,spdCorr);
1857         fHOutMultV0MvsTRKqual2->Fill(v0Corr,nTracks);
1858         fHOutMultTRKvsCL1qual2->Fill(nTracks,spdCorr);
1859       }
1860     }
1861     PostData(1, fOutputList);
1862   }
1863 }
1864 //________________________________________________________________________
1865 void AliCentralitySelectionTask::Terminate(Option_t */*option*/)
1866 {
1867   // Terminate analysis
1868 }
1869 //________________________________________________________________________
1870 Int_t AliCentralitySelectionTask::SetupRun(const AliVEvent* const esd)
1871 {
1872   // Setup files for run
1873
1874   if (!esd)
1875     return -1;
1876
1877   // check if something to be done
1878   if (fCurrentRun == esd->GetRunNumber())
1879     return 0;
1880   else
1881     fCurrentRun = esd->GetRunNumber();
1882
1883   TString fileName =(Form("%s/COMMON/CENTRALITY/data/centrality.root", AliAnalysisManager::GetOADBPath()));
1884   AliInfo(Form("Setup Centrality Selection for run %d with file %s\n",fCurrentRun,fileName.Data()));
1885
1886   AliOADBContainer *con = new AliOADBContainer("OADB");
1887   con->InitFromFile(fileName,"Centrality");
1888
1889   AliOADBCentrality*  centOADB = 0;
1890   centOADB = (AliOADBCentrality*)(con->GetObject(fCurrentRun));
1891   if (!centOADB) {
1892     AliWarning(Form("Centrality OADB does not exist for run %d, using Default \n",fCurrentRun ));
1893     centOADB  = (AliOADBCentrality*)(con->GetDefaultObject("oadbDefault"));
1894   }
1895
1896   Bool_t isHijing=kFALSE;
1897   Bool_t isDpmjet=kFALSE;
1898   AliAnalysisManager* anMan = AliAnalysisManager::GetAnalysisManager();
1899   AliMCEventHandler* eventHandler = (AliMCEventHandler*)anMan->GetMCtruthEventHandler();
1900   AliMCEvent*  mcEvent=0;
1901   if (fIsMCInput && eventHandler && (mcEvent=eventHandler->MCEvent()) ) {
1902     AliGenEventHeader* mcGenH = mcEvent->GenEventHeader();
1903     if (mcGenH->InheritsFrom(AliGenHijingEventHeader::Class())) isHijing=kTRUE;
1904     else if (mcGenH->InheritsFrom(AliGenCocktailEventHeader::Class())) isHijing=kTRUE;
1905     else if (mcGenH->InheritsFrom(AliGenDPMjetEventHeader::Class())) isDpmjet=kTRUE;
1906   }
1907
1908
1909   // modes
1910   fUseScaling   = centOADB->UseScaling();
1911   fUseCleaning  = centOADB->UseCleaning();
1912
1913   // cuts
1914   fZVCut        = centOADB->ZVCut();
1915   fOutliersCut  = centOADB->OutliersCut();
1916
1917   // centrality histos
1918   fHtempV0M       = centOADB->V0hist();
1919   fHtempV0A       = centOADB->V0Ahist();
1920   fHtempV0A0      = centOADB->V0A0hist();
1921   fHtempV0A123    = centOADB->V0A123hist();
1922   fHtempV0C       = centOADB->V0Chist();
1923   fHtempV0MEq     = centOADB->V0Eqhist();
1924   fHtempV0AEq     = centOADB->V0AEqhist();
1925   fHtempV0CEq     = centOADB->V0CEqhist();
1926   fHtempTRK       = centOADB->TPChist();
1927   fHtempCL1       = centOADB->SPDhist();
1928   fHtempCND       = centOADB->CNDhist();
1929   fHtempFMD       = centOADB->FMDhist();
1930   fHtempZNA       = centOADB->ZNAhist();
1931   fHtempZNC       = centOADB->ZNChist();
1932   fHtempZPA       = centOADB->ZPAhist();
1933   fHtempZPC       = centOADB->ZPChist();
1934   fHtempZEMvsZDC  = centOADB->ZEMvsZDChist();
1935
1936    if (isHijing) {
1937      fHtempNPA           = centOADB->NPAhist();
1938      fHtempV0Mtrue       = centOADB->V0histtrue();
1939      fHtempV0Atrue       = centOADB->V0Ahisttrue();
1940      fHtempV0Ctrue       = centOADB->V0Chisttrue();
1941      fHtempV0MEqtrue     = centOADB->V0Eqhisttrue();
1942      fHtempV0AEqtrue     = centOADB->V0AEqhisttrue();
1943      fHtempV0CEqtrue     = centOADB->V0CEqhisttrue();
1944      fHtempTRKtrue       = centOADB->TPChisttrue();
1945      fHtempCL1true       = centOADB->SPDhisttrue();
1946      fHtempCNDtrue       = centOADB->CNDhisttrue();
1947      fHtempFMDtrue       = centOADB->FMDhisttrue();
1948      fHtempZNAtrue       = centOADB->ZNAhisttrue();
1949      fHtempZNCtrue       = centOADB->ZNChisttrue();
1950      fHtempZPAtrue       = centOADB->ZPAhisttrue();
1951      fHtempZPCtrue       = centOADB->ZPChisttrue();
1952    }   else if (isDpmjet)   {
1953      fHtempNPA           = centOADB->NPAhistDPM();
1954      fHtempV0Mtrue       = centOADB->V0histtrueDPM();
1955      fHtempV0Atrue       = centOADB->V0AhisttrueDPM();
1956      fHtempV0Ctrue       = centOADB->V0ChisttrueDPM();
1957      fHtempV0MEqtrue     = centOADB->V0EqhisttrueDPM();
1958      fHtempV0AEqtrue     = centOADB->V0AEqhisttrueDPM();
1959      fHtempV0CEqtrue     = centOADB->V0CEqhisttrueDPM();
1960      fHtempTRKtrue       = centOADB->TPChisttrueDPM();
1961      fHtempCL1true       = centOADB->SPDhisttrueDPM();
1962      fHtempCNDtrue       = centOADB->CNDhisttrueDPM();
1963      fHtempFMDtrue       = centOADB->FMDhisttrueDPM();
1964      fHtempZNAtrue       = centOADB->ZNAhisttrueDPM();
1965      fHtempZNCtrue       = centOADB->ZNChisttrueDPM();
1966      fHtempZPAtrue       = centOADB->ZPAhisttrueDPM();
1967      fHtempZPCtrue       = centOADB->ZPChisttrueDPM();
1968    }
1969
1970
1971     TString path = gSystem->ExpandPathName(fileName.Data());
1972   if (!fHtempV0M) AliWarning(Form("Calibration for V0M does not exist in %s", path.Data()));
1973   if (!fHtempV0A) AliWarning(Form("Calibration for V0A does not exist in %s", path.Data()));
1974   if (!fHtempV0C) AliWarning(Form("Calibration for V0C does not exist in %s", path.Data()));
1975   if (!fHtempV0MEq) AliWarning(Form("Calibration for V0MEq does not exist in %s", path.Data()));
1976   if (!fHtempV0AEq) AliWarning(Form("Calibration for V0AEq does not exist in %s", path.Data()));
1977   if (!fHtempV0CEq) AliWarning(Form("Calibration for V0CEq does not exist in %s", path.Data()));
1978   if (!fHtempTRK) AliWarning(Form("Calibration for TRK does not exist in %s", path.Data()));
1979   if (!fHtempCL1) AliWarning(Form("Calibration for CL1 does not exist in %s", path.Data()));
1980   if (!fHtempCND) AliWarning(Form("Calibration for CND does not exist in %s", path.Data()));
1981   if (!fHtempZNA) AliWarning(Form("Calibration for ZNA does not exist in %s", path.Data()));
1982   if (!fHtempZNC) AliWarning(Form("Calibration for ZNC does not exist in %s", path.Data()));
1983   if (!fHtempZPA) AliWarning(Form("Calibration for ZPA does not exist in %s", path.Data()));
1984   if (!fHtempZPC) AliWarning(Form("Calibration for ZPC does not exist in %s", path.Data()));
1985   if (!fHtempFMD) AliWarning(Form("Calibration for FMD does not exist in %s", path.Data()));
1986   if (!fHtempZEMvsZDC) AliWarning(Form("Calibration for ZEMvsZDC does not exist in %s", path.Data()));
1987   if (!fHtempNPA) AliWarning(Form("Calibration for NPA does not exist in %s", path.Data()));
1988
1989   if (!fHtempV0Mtrue) AliWarning(Form("Calibration for V0Mtrue does not exist in %s", path.Data()));
1990   if (!fHtempV0Atrue) AliWarning(Form("Calibration for V0Atrue does not exist in %s", path.Data()));
1991   if (!fHtempV0Ctrue) AliWarning(Form("Calibration for V0Ctrue does not exist in %s", path.Data()));
1992   if (!fHtempV0MEqtrue) AliWarning(Form("Calibration for V0MEqtrue does not exist in %s", path.Data()));
1993   if (!fHtempV0AEqtrue) AliWarning(Form("Calibration for V0AEqtrue does not exist in %s", path.Data()));
1994   if (!fHtempV0CEqtrue) AliWarning(Form("Calibration for V0CEqtrue does not exist in %s", path.Data()));
1995   if (!fHtempTRKtrue) AliWarning(Form("Calibration for TRKtrue does not exist in %s", path.Data()));
1996   if (!fHtempCL1true) AliWarning(Form("Calibration for CL1true does not exist in %s", path.Data()));
1997   if (!fHtempCNDtrue) AliWarning(Form("Calibration for CNDtrue does not exist in %s", path.Data()));
1998   if (!fHtempZNAtrue) AliWarning(Form("Calibration for ZNAtrue does not exist in %s", path.Data()));
1999   if (!fHtempZNCtrue) AliWarning(Form("Calibration for ZNCtrue does not exist in %s", path.Data()));
2000   if (!fHtempZPAtrue) AliWarning(Form("Calibration for ZPAtrue does not exist in %s", path.Data()));
2001   if (!fHtempZPCtrue) AliWarning(Form("Calibration for ZPCtrue does not exist in %s", path.Data()));
2002   if (!fHtempFMDtrue) AliWarning(Form("Calibration for FMDtrue does not exist in %s", path.Data()));
2003
2004
2005   // scale factors
2006   fV0MScaleFactor    = centOADB->V0MScaleFactor();
2007   fSPDScaleFactor    = centOADB->SPDScaleFactor();
2008   fTPCScaleFactor    = centOADB->TPCScaleFactor();
2009   fV0MScaleFactorMC  = centOADB->V0MScaleFactorMC();
2010
2011   // outliers parameters
2012   fV0MSPDOutlierPar0 = centOADB->V0MSPDOutlierPar0();
2013   fV0MSPDOutlierPar1 = centOADB->V0MSPDOutlierPar1();
2014   fV0MTPCOutlierPar0 = centOADB->V0MTPCOutlierPar0();
2015   fV0MTPCOutlierPar1 = centOADB->V0MTPCOutlierPar1();
2016
2017   fV0MSPDSigmaOutlierPar0 = centOADB->V0MSPDSigmaOutlierPar0();
2018   fV0MSPDSigmaOutlierPar1 = centOADB->V0MSPDSigmaOutlierPar1();
2019   fV0MSPDSigmaOutlierPar2 = centOADB->V0MSPDSigmaOutlierPar2();
2020   fV0MTPCSigmaOutlierPar0 = centOADB->V0MTPCSigmaOutlierPar0();
2021   fV0MTPCSigmaOutlierPar1 = centOADB->V0MTPCSigmaOutlierPar1();
2022   fV0MTPCSigmaOutlierPar2 = centOADB->V0MTPCSigmaOutlierPar2();
2023
2024   fV0MZDCOutlierPar0 =      centOADB->V0MZDCOutlierPar0();
2025   fV0MZDCOutlierPar1 =      centOADB->V0MZDCOutlierPar1();
2026   fV0MZDCEcalOutlierPar0 =  centOADB->V0MZDCEcalOutlierPar0();
2027   fV0MZDCEcalOutlierPar1 =  centOADB->V0MZDCEcalOutlierPar1();
2028
2029   return 0;
2030 }
2031
2032
2033
2034 //________________________________________________________________________
2035 Bool_t AliCentralitySelectionTask::IsOutlierV0MSPD(Float_t spd, Float_t v0, Int_t cent) const
2036 {
2037   // Clean outliers
2038   Float_t val = fV0MSPDOutlierPar0 +  fV0MSPDOutlierPar1 * v0;
2039   Float_t spdSigma = fV0MSPDSigmaOutlierPar0 + fV0MSPDSigmaOutlierPar1*cent + fV0MSPDSigmaOutlierPar2*cent*cent;
2040   if ( TMath::Abs(spd-val) > fOutliersCut*spdSigma )
2041     return kTRUE;
2042   else
2043     return kFALSE;
2044 }
2045
2046 //________________________________________________________________________
2047 Bool_t AliCentralitySelectionTask::IsOutlierV0MTPC(Int_t tracks, Float_t v0, Int_t cent) const
2048 {
2049   // Clean outliers
2050   Float_t val = fV0MTPCOutlierPar0 +  fV0MTPCOutlierPar1 * v0;
2051   Float_t tpcSigma = fV0MTPCSigmaOutlierPar0 + fV0MTPCSigmaOutlierPar1*cent + fV0MTPCSigmaOutlierPar2*cent*cent;
2052   if ( TMath::Abs(tracks-val) > fOutliersCut*tpcSigma )
2053     return kTRUE;
2054   else
2055     return kFALSE;
2056 }
2057
2058 //________________________________________________________________________
2059 Bool_t AliCentralitySelectionTask::IsOutlierV0MZDC(Float_t zdc, Float_t v0) const
2060 {
2061   // Clean outliers
2062   Float_t val = fV0MZDCOutlierPar0 + fV0MZDCOutlierPar1 * v0;
2063   if (zdc >  val)
2064     return kTRUE;
2065   else
2066   return kFALSE;
2067 }
2068
2069 //________________________________________________________________________
2070 Bool_t AliCentralitySelectionTask::IsOutlierV0MZDCECal(Float_t zdc, Float_t v0) const
2071 {
2072   // Clean outliers
2073   Float_t val = fV0MZDCEcalOutlierPar0 + fV0MZDCEcalOutlierPar1 * v0;
2074   if (zdc >  val)
2075     return kTRUE;
2076   else
2077     return kFALSE;
2078 }
2079