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