STEER/STEERBase/AliCentrality.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 AliCentrality
  18 //   author: Alberica Toia
  19 //*****************************************************
  20 /// A container for the centrality stored in AOD in ESD
  21
  22 #include "AliCentrality.h"
  23
  24 ClassImp(AliCentrality)
  25
  26 AliCentrality::AliCentrality() : TNamed("Centrality", "Centrality"),
  27   fQuality(999),
  28   fCentralityV0M(0),
  29   fCentralityV0A(0),
  30   fCentralityV0A0(0),
  31   fCentralityV0A123(0),
  32   fCentralityV0C(0),
  33   fCentralityV0A23(0),
  34   fCentralityV0C01(0),
  35   fCentralityV0S(0),
  36   fCentralityV0MEq(0),
  37   fCentralityV0AEq(0),
  38   fCentralityV0CEq(0),
  39   fCentralityFMD(0),
  40   fCentralityTRK(0),
  41   fCentralityTKL(0),
  42   fCentralityCL0(0),
  43   fCentralityCL1(0),
  44   fCentralityCND(0),
  45   fCentralityZNA(0),
  46   fCentralityZNC(0),
  47   fCentralityZPA(0),
  48   fCentralityZPC(0),
  49   fCentralityNPA(0),
  50   fCentralityV0MvsFMD(0),
  51   fCentralityTKLvsV0M(0),
  52   fCentralityZEMvsZDC(0),
  53   fCentralityV0Mtrue(0),
  54   fCentralityV0Atrue(0),
  55   fCentralityV0Ctrue(0),
  56   fCentralityV0MEqtrue(0),
  57   fCentralityV0AEqtrue(0),
  58   fCentralityV0CEqtrue(0),
  59   fCentralityFMDtrue(0),
  60   fCentralityTRKtrue(0),
  61   fCentralityTKLtrue(0),
  62   fCentralityCL0true(0),
  63   fCentralityCL1true(0),
  64   fCentralityCNDtrue(0),
  65   fCentralityZNAtrue(0),
  66   fCentralityZNCtrue(0),
  67   fCentralityZPAtrue(0),
  68   fCentralityZPCtrue(0)
  69 {
  70   /// constructor
  71 }
  72
  73 AliCentrality::AliCentrality(const AliCentrality& cnt) :
  74   TNamed(cnt),
  75   fQuality(cnt.fQuality),
  76   fCentralityV0M(cnt.fCentralityV0M),
  77   fCentralityV0A(cnt.fCentralityV0A),
  78   fCentralityV0A0(cnt.fCentralityV0A0),
  79   fCentralityV0A123(cnt.fCentralityV0A123),
  80   fCentralityV0C(cnt.fCentralityV0C),
  81   fCentralityV0A23(cnt.fCentralityV0A23),
  82   fCentralityV0C01(cnt.fCentralityV0C01),
  83   fCentralityV0S(cnt.fCentralityV0S),
  84   fCentralityV0MEq(cnt.fCentralityV0MEq),
  85   fCentralityV0AEq(cnt.fCentralityV0AEq),
  86   fCentralityV0CEq(cnt.fCentralityV0CEq),
  87   fCentralityFMD(cnt.fCentralityFMD),
  88   fCentralityTRK(cnt.fCentralityTRK),
  89   fCentralityTKL(cnt.fCentralityTKL),
  90   fCentralityCL0(cnt.fCentralityCL0),
  91   fCentralityCL1(cnt.fCentralityCL1),
  92   fCentralityCND(cnt.fCentralityCND),
  93   fCentralityZNA(cnt.fCentralityZNA),
  94   fCentralityZNC(cnt.fCentralityZNC),
  95   fCentralityZPA(cnt.fCentralityZPA),
  96   fCentralityZPC(cnt.fCentralityZPC),
  97   fCentralityNPA(cnt.fCentralityNPA),
  98   fCentralityV0MvsFMD(cnt.fCentralityV0MvsFMD),
  99   fCentralityTKLvsV0M(cnt.fCentralityTKLvsV0M),
 100   fCentralityZEMvsZDC(cnt.fCentralityZEMvsZDC),
 101   fCentralityV0Mtrue(cnt.fCentralityV0Mtrue),
 102   fCentralityV0Atrue(cnt.fCentralityV0Atrue),
 103   fCentralityV0Ctrue(cnt.fCentralityV0Ctrue),
 104   fCentralityV0MEqtrue(cnt.fCentralityV0MEqtrue),
 105   fCentralityV0AEqtrue(cnt.fCentralityV0AEqtrue),
 106   fCentralityV0CEqtrue(cnt.fCentralityV0CEqtrue),
 107   fCentralityFMDtrue(cnt.fCentralityFMDtrue),
 108   fCentralityTRKtrue(cnt.fCentralityTRKtrue),
 109   fCentralityTKLtrue(cnt.fCentralityTKLtrue),
 110   fCentralityCL0true(cnt.fCentralityCL0true),
 111   fCentralityCL1true(cnt.fCentralityCL1true),
 112   fCentralityCNDtrue(cnt.fCentralityCNDtrue),
 113   fCentralityZNAtrue(cnt.fCentralityZNAtrue),
 114   fCentralityZNCtrue(cnt.fCentralityZNCtrue),
 115   fCentralityZPAtrue(cnt.fCentralityZPAtrue),
 116   fCentralityZPCtrue(cnt.fCentralityZPCtrue)
 117 {
 118   /// Copy constructor
 119 }
 120
 121 AliCentrality& AliCentrality::operator=(const AliCentrality& c)
 122 {
 123   /// Assignment operator
 124   if (this!=&c) {
 125     TNamed::operator=(c);
 126     fQuality = c.fQuality;
 127     fCentralityV0M = c.fCentralityV0M;
 128     fCentralityV0A = c.fCentralityV0A;
 129     fCentralityV0A0 = c.fCentralityV0A0;
 130     fCentralityV0A123 = c.fCentralityV0A123;
 131     fCentralityV0C = c.fCentralityV0C;
 132     fCentralityV0A23 = c.fCentralityV0A23;
 133     fCentralityV0C01 = c.fCentralityV0C01;
 134     fCentralityV0S = c.fCentralityV0S;
 135     fCentralityV0MEq = c.fCentralityV0MEq;
 136     fCentralityV0AEq = c.fCentralityV0AEq;
 137     fCentralityV0CEq = c.fCentralityV0CEq;
 138     fCentralityFMD = c.fCentralityFMD;
 139     fCentralityTRK = c.fCentralityTRK;
 140     fCentralityTKL = c.fCentralityTKL;
 141     fCentralityCL0 = c.fCentralityCL0;
 142     fCentralityCL1 = c.fCentralityCL1;
 143     fCentralityCND = c.fCentralityCND;
 144     fCentralityZNA = c.fCentralityZNA;
 145     fCentralityZNC = c.fCentralityZNC;
 146     fCentralityZPA = c.fCentralityZPA;
 147     fCentralityZPC = c.fCentralityZPC;
 148     fCentralityNPA = c.fCentralityNPA;
 149     fCentralityV0MvsFMD = c.fCentralityV0MvsFMD;
 150     fCentralityTKLvsV0M = c.fCentralityTKLvsV0M;
 151     fCentralityZEMvsZDC = c.fCentralityZEMvsZDC;
 152     fCentralityV0Mtrue = c.fCentralityV0Mtrue;
 153     fCentralityV0Atrue = c.fCentralityV0Atrue;
 154     fCentralityV0Ctrue = c.fCentralityV0Ctrue;
 155     fCentralityV0MEqtrue = c.fCentralityV0MEqtrue;
 156     fCentralityV0AEqtrue = c.fCentralityV0AEqtrue;
 157     fCentralityV0CEqtrue = c.fCentralityV0CEqtrue;
 158     fCentralityFMDtrue = c.fCentralityFMDtrue;
 159     fCentralityTRKtrue = c.fCentralityTRKtrue;
 160     fCentralityTKLtrue = c.fCentralityTKLtrue;
 161     fCentralityCL0true = c.fCentralityCL0true;
 162     fCentralityCL1true = c.fCentralityCL1true;
 163     fCentralityCNDtrue = c.fCentralityCNDtrue;
 164     fCentralityZNAtrue = c.fCentralityZNAtrue;
 165     fCentralityZNCtrue = c.fCentralityZNCtrue;
 166     fCentralityZPAtrue = c.fCentralityZPAtrue;
 167     fCentralityZPCtrue = c.fCentralityZPCtrue;
 168   }
 169
 170   return *this;
 171 }
 172
 173 AliCentrality::~AliCentrality()
 174 {
 175   /// destructor
 176 }
 177
 178 Int_t AliCentrality::GetQuality() const
 179 {
 180   return fQuality;
 181 }
 182
 183 Float_t AliCentrality::GetCentralityPercentile(const char *x) const
 184 {
 185 // Return the centrality percentile
 186   if (fQuality == 0) {
 187     TString method = x;
 188     if(method.CompareTo("V0M")==0)      return fCentralityV0M;
 189     if(method.CompareTo("V0A")==0)      return fCentralityV0A;
 190     if(method.CompareTo("V0A0")==0)   return fCentralityV0A0;
 191     if(method.CompareTo("V0A123")==0)   return fCentralityV0A123;
 192     if(method.CompareTo("V0C")==0)      return fCentralityV0C;
 193     if(method.CompareTo("V0A23")==0)    return fCentralityV0A23;
 194     if(method.CompareTo("V0C01")==0)    return fCentralityV0C01;
 195     if(method.CompareTo("V0S")==0)      return fCentralityV0S;
 196     if(method.CompareTo("V0MEq")==0)    return fCentralityV0MEq;
 197     if(method.CompareTo("V0AEq")==0)    return fCentralityV0AEq;
 198     if(method.CompareTo("V0CEq")==0)    return fCentralityV0CEq;
 199     if(method.CompareTo("FMD")==0)      return fCentralityFMD;
 200     if(method.CompareTo("TRK")==0)      return fCentralityTRK;
 201     if(method.CompareTo("TKL")==0)      return fCentralityTKL;
 202     if(method.CompareTo("CL0")==0)      return fCentralityCL0;
 203     if(method.CompareTo("CL1")==0)      return fCentralityCL1;
 204     if(method.CompareTo("CND")==0)      return fCentralityCND;
 205     if(method.CompareTo("ZNA")==0)      return fCentralityZNA;
 206     if(method.CompareTo("ZNC")==0)      return fCentralityZNC;
 207     if(method.CompareTo("ZPA")==0)      return fCentralityZPA;
 208     if(method.CompareTo("ZPC")==0)      return fCentralityZPC;
 209     if(method.CompareTo("NPA")==0)      return fCentralityNPA;
 210     if(method.CompareTo("V0MvsFMD")==0) return fCentralityV0MvsFMD;
 211     if(method.CompareTo("TKLvsV0M")==0) return fCentralityTKLvsV0M;
 212     if(method.CompareTo("ZEMvsZDC")==0) return fCentralityZEMvsZDC;
 213     if(method.CompareTo("V0Mtrue")==0)      return fCentralityV0Mtrue;
 214     if(method.CompareTo("V0Atrue")==0)      return fCentralityV0Atrue;
 215     if(method.CompareTo("V0Ctrue")==0)      return fCentralityV0Ctrue;
 216     if(method.CompareTo("V0MEqtrue")==0)    return fCentralityV0MEqtrue;
 217     if(method.CompareTo("V0AEqtrue")==0)    return fCentralityV0AEqtrue;
 218     if(method.CompareTo("V0CEqtrue")==0)    return fCentralityV0CEqtrue;
 219     if(method.CompareTo("FMDtrue")==0)      return fCentralityFMDtrue;
 220     if(method.CompareTo("TRKtrue")==0)      return fCentralityTRKtrue;
 221     if(method.CompareTo("TKLtrue")==0)      return fCentralityTKLtrue;
 222     if(method.CompareTo("CL0true")==0)      return fCentralityCL0true;
 223     if(method.CompareTo("CL1true")==0)      return fCentralityCL1true;
 224     if(method.CompareTo("CNDtrue")==0)      return fCentralityCNDtrue;
 225     if(method.CompareTo("ZNAtrue")==0)      return fCentralityZNAtrue;
 226     if(method.CompareTo("ZNCtrue")==0)      return fCentralityZNCtrue;
 227     if(method.CompareTo("ZPAtrue")==0)      return fCentralityZPAtrue;
 228     if(method.CompareTo("ZPCtrue")==0)      return fCentralityZPCtrue;
 229     return -1;
 230   } else {
 231     return -1;
 232   }
 233 }
 234
 235 Int_t AliCentrality::GetCentralityClass10(const char *x) const
 236 {
 237 // Return the centrality class
 238   if (fQuality == 0) {
 239     TString method = x;
 240     if(method.CompareTo("V0M")==0)      return (Int_t) (fCentralityV0M / 10.0);
 241     if(method.CompareTo("V0A")==0)      return (Int_t) (fCentralityV0A / 10.0);
 242     if(method.CompareTo("V0A0")==0)     return (Int_t) (fCentralityV0A0 / 10.0);
 243     if(method.CompareTo("V0A123")==0)   return (Int_t) (fCentralityV0A123 / 10.0);
 244     if(method.CompareTo("V0C")==0)      return (Int_t) (fCentralityV0C / 10.0);
 245     if(method.CompareTo("V0A23")==0)    return (Int_t) (fCentralityV0A23 / 10.0);
 246     if(method.CompareTo("V0C01")==0)    return (Int_t) (fCentralityV0C01 / 10.0);
 247     if(method.CompareTo("V0S")==0)      return (Int_t) (fCentralityV0S / 10.0);
 248     if(method.CompareTo("V0MEq")==0)    return (Int_t) (fCentralityV0MEq / 10.0);
 249     if(method.CompareTo("V0AEq")==0)    return (Int_t) (fCentralityV0AEq / 10.0);
 250     if(method.CompareTo("V0CEq")==0)    return (Int_t) (fCentralityV0CEq / 10.0);
 251     if(method.CompareTo("FMD")==0)      return (Int_t) (fCentralityFMD / 10.0);
 252     if(method.CompareTo("TRK")==0)      return (Int_t) (fCentralityTRK / 10.0);
 253     if(method.CompareTo("TKL")==0)      return (Int_t) (fCentralityTKL / 10.0);
 254     if(method.CompareTo("CL0")==0)      return (Int_t) (fCentralityCL0 / 10.0);
 255     if(method.CompareTo("CL1")==0)      return (Int_t) (fCentralityCL1 / 10.0);
 256     if(method.CompareTo("CND")==0)      return (Int_t) (fCentralityCND / 10.0);
 257     if(method.CompareTo("ZNA")==0)      return (Int_t) (fCentralityZNA / 10.0);
 258     if(method.CompareTo("ZNC")==0)      return (Int_t) (fCentralityZNC / 10.0);
 259     if(method.CompareTo("ZPA")==0)      return (Int_t) (fCentralityZPA / 10.0);
 260     if(method.CompareTo("ZPC")==0)      return (Int_t) (fCentralityZPC / 10.0);
 261     if(method.CompareTo("NPA")==0)      return (Int_t) (fCentralityNPA / 10.0);
 262     if(method.CompareTo("V0MvsFMD")==0) return (Int_t) (fCentralityV0MvsFMD / 10.0);
 263     if(method.CompareTo("TKLvsV0M")==0) return (Int_t) (fCentralityTKLvsV0M / 10.0);
 264     if(method.CompareTo("ZEMvsZDC")==0) return (Int_t) (fCentralityZEMvsZDC / 10.0);
 265     if(method.CompareTo("V0Mtrue")==0)  return (Int_t) (fCentralityV0Mtrue / 10.0);
 266     if(method.CompareTo("V0Atrue")==0)  return (Int_t) (fCentralityV0Atrue / 10.0);
 267     if(method.CompareTo("V0Ctrue")==0)  return (Int_t) (fCentralityV0Ctrue / 10.0);
 268     if(method.CompareTo("V0MEqtrue")==0)return (Int_t) (fCentralityV0MEqtrue / 10.0);
 269     if(method.CompareTo("V0AEqtrue")==0)return (Int_t) (fCentralityV0AEqtrue / 10.0);
 270     if(method.CompareTo("V0CEqtrue")==0)return (Int_t) (fCentralityV0CEqtrue / 10.0);
 271     if(method.CompareTo("FMDtrue")==0)  return (Int_t) (fCentralityFMDtrue / 10.0);
 272     if(method.CompareTo("TRKtrue")==0)  return (Int_t) (fCentralityTRKtrue / 10.0);
 273     if(method.CompareTo("TKLtrue")==0)  return (Int_t) (fCentralityTKLtrue / 10.0);
 274     if(method.CompareTo("CL0true")==0)  return (Int_t) (fCentralityCL0true / 10.0);
 275     if(method.CompareTo("CL1true")==0)  return (Int_t) (fCentralityCL1true / 10.0);
 276     if(method.CompareTo("CNDtrue")==0)  return (Int_t) (fCentralityCNDtrue / 10.0);
 277     if(method.CompareTo("ZNAtrue")==0)  return (Int_t) (fCentralityZNAtrue / 10.0);
 278     if(method.CompareTo("ZNCtrue")==0)  return (Int_t) (fCentralityZNCtrue / 10.0);
 279     if(method.CompareTo("ZPAtrue")==0)  return (Int_t) (fCentralityZPAtrue / 10.0);
 280     if(method.CompareTo("ZPCtrue")==0)  return (Int_t) (fCentralityZPCtrue / 10.0);
 281     return -1;
 282   } else {
 283     return -1;
 284   }
 285 }
 286
 287 Int_t AliCentrality::GetCentralityClass5(const char *x) const
 288 {
 289 // Return the centrality class
 290   if (fQuality == 0) {
 291     TString method = x;
 292     if(method.CompareTo("V0M")==0)      return (Int_t) (fCentralityV0M / 5.0);
 293     if(method.CompareTo("V0A")==0)      return (Int_t) (fCentralityV0A / 5.0);
 294     if(method.CompareTo("V0A0")==0)     return (Int_t) (fCentralityV0A0 / 5.0);
 295     if(method.CompareTo("V0A123")==0)   return (Int_t) (fCentralityV0A123 / 5.0);
 296     if(method.CompareTo("V0C")==0)      return (Int_t) (fCentralityV0C / 5.0);
 297     if(method.CompareTo("V0A23")==0)    return (Int_t) (fCentralityV0A23 / 5.0);
 298     if(method.CompareTo("V0C01")==0)    return (Int_t) (fCentralityV0C01 / 5.0);
 299     if(method.CompareTo("V0S")==0)      return (Int_t) (fCentralityV0S / 5.0);
 300     if(method.CompareTo("V0MEq")==0)    return (Int_t) (fCentralityV0MEq / 5.0);
 301     if(method.CompareTo("V0AEq")==0)    return (Int_t) (fCentralityV0AEq / 5.0);
 302     if(method.CompareTo("V0CEq")==0)    return (Int_t) (fCentralityV0CEq / 5.0);
 303     if(method.CompareTo("FMD")==0)      return (Int_t) (fCentralityFMD / 5.0);
 304     if(method.CompareTo("TRK")==0)      return (Int_t) (fCentralityTRK / 5.0);
 305     if(method.CompareTo("TKL")==0)      return (Int_t) (fCentralityTKL / 5.0);
 306     if(method.CompareTo("CL0")==0)      return (Int_t) (fCentralityCL0 / 5.0);
 307     if(method.CompareTo("CL1")==0)      return (Int_t) (fCentralityCL1 / 5.0);
 308     if(method.CompareTo("CND")==0)      return (Int_t) (fCentralityCND / 5.0);
 309     if(method.CompareTo("ZNA")==0)      return (Int_t) (fCentralityZNA / 5.0);
 310     if(method.CompareTo("ZNC")==0)      return (Int_t) (fCentralityZNC / 5.0);
 311     if(method.CompareTo("ZPA")==0)      return (Int_t) (fCentralityZPA / 5.0);
 312     if(method.CompareTo("ZPC")==0)      return (Int_t) (fCentralityZPC / 5.0);
 313     if(method.CompareTo("NPA")==0)      return (Int_t) (fCentralityNPA / 5.0);
 314     if(method.CompareTo("V0MvsFMD")==0) return (Int_t) (fCentralityV0MvsFMD / 5.0);
 315     if(method.CompareTo("TKLvsV0M")==0) return (Int_t) (fCentralityTKLvsV0M / 5.0);
 316     if(method.CompareTo("ZEMvsZDC")==0) return (Int_t) (fCentralityZEMvsZDC / 5.0);
 317     if(method.CompareTo("V0Mtrue")==0)  return (Int_t) (fCentralityV0Mtrue / 5.0);
 318     if(method.CompareTo("V0Atrue")==0)  return (Int_t) (fCentralityV0Atrue / 5.0);
 319     if(method.CompareTo("V0Ctrue")==0)  return (Int_t) (fCentralityV0Ctrue / 5.0);
 320     if(method.CompareTo("V0MEqtrue")==0)return (Int_t) (fCentralityV0MEqtrue / 5.0);
 321     if(method.CompareTo("V0AEqtrue")==0)return (Int_t) (fCentralityV0AEqtrue / 5.0);
 322     if(method.CompareTo("V0CEqtrue")==0)return (Int_t) (fCentralityV0CEqtrue / 5.0);
 323     if(method.CompareTo("FMDtrue")==0)  return (Int_t) (fCentralityFMDtrue / 5.0);
 324     if(method.CompareTo("TRKtrue")==0)  return (Int_t) (fCentralityTRKtrue / 5.0);
 325     if(method.CompareTo("TKLtrue")==0)  return (Int_t) (fCentralityTKLtrue / 5.0);
 326     if(method.CompareTo("CL0true")==0)  return (Int_t) (fCentralityCL0true / 5.0);
 327     if(method.CompareTo("CL1true")==0)  return (Int_t) (fCentralityCL1true / 5.0);
 328     if(method.CompareTo("CNDtrue")==0)  return (Int_t) (fCentralityCNDtrue / 5.0);
 329     if(method.CompareTo("ZNAtrue")==0)  return (Int_t) (fCentralityZNAtrue / 5.0);
 330     if(method.CompareTo("ZNCtrue")==0)  return (Int_t) (fCentralityZNCtrue / 5.0);
 331     if(method.CompareTo("ZPAtrue")==0)  return (Int_t) (fCentralityZPAtrue / 5.0);
 332     if(method.CompareTo("ZPCtrue")==0)  return (Int_t) (fCentralityZPCtrue / 5.0);
 333     return -1;
 334   } else {
 335     return -1;
 336   }
 337 }
 338
 339
 340 Bool_t AliCentrality::IsEventInCentralityClass(Float_t a, Float_t b, const char *x) const
 341 {
 342 // True if event is inside a given class
 343   if (fQuality == 0) {
 344     TString method = x;
 345     if ((method.CompareTo("V0M")==0) && (fCentralityV0M >=a && fCentralityV0M < b)) return kTRUE;
 346     if ((method.CompareTo("V0A")==0) && (fCentralityV0A >=a && fCentralityV0A < b)) return kTRUE;
 347     if ((method.CompareTo("V0A0")==0) && (fCentralityV0A0 >=a && fCentralityV0A0 < b)) return kTRUE;
 348     if ((method.CompareTo("V0A123")==0) && (fCentralityV0A123 >=a && fCentralityV0A123 < b)) return kTRUE;
 349     if ((method.CompareTo("V0C")==0) && (fCentralityV0C >=a && fCentralityV0C < b)) return kTRUE;
 350     if ((method.CompareTo("V0A23")==0) && (fCentralityV0A23 >=a && fCentralityV0A23 < b)) return kTRUE;
 351     if ((method.CompareTo("V0C01")==0) && (fCentralityV0C01 >=a && fCentralityV0C01 < b)) return kTRUE;
 352     if ((method.CompareTo("V0S")==0) && (fCentralityV0S >=a && fCentralityV0S < b)) return kTRUE;
 353     if ((method.CompareTo("V0MEq")==0) && (fCentralityV0MEq >=a && fCentralityV0MEq < b)) return kTRUE;
 354     if ((method.CompareTo("V0AEq")==0) && (fCentralityV0AEq >=a && fCentralityV0AEq < b)) return kTRUE;
 355     if ((method.CompareTo("V0CEq")==0) && (fCentralityV0CEq >=a && fCentralityV0CEq < b)) return kTRUE;
 356     if ((method.CompareTo("FMD")==0) && (fCentralityFMD >=a && fCentralityFMD < b)) return kTRUE;
 357     if ((method.CompareTo("TRK")==0) && (fCentralityTRK >=a && fCentralityTRK < b)) return kTRUE;
 358     if ((method.CompareTo("TKL")==0) && (fCentralityTKL >=a && fCentralityTKL < b)) return kTRUE;
 359     if ((method.CompareTo("CL0")==0) && (fCentralityCL0 >=a && fCentralityCL0 < b)) return kTRUE;
 360     if ((method.CompareTo("CL1")==0) && (fCentralityCL1 >=a && fCentralityCL1 < b)) return kTRUE;
 361     if ((method.CompareTo("CND")==0) && (fCentralityCND >=a && fCentralityCND < b)) return kTRUE;
 362     if ((method.CompareTo("ZNA")==0) && (fCentralityZNA >=a && fCentralityZNA < b)) return kTRUE;
 363     if ((method.CompareTo("ZNC")==0) && (fCentralityZNC >=a && fCentralityZNC < b)) return kTRUE;
 364     if ((method.CompareTo("ZPA")==0) && (fCentralityZPA >=a && fCentralityZPA < b)) return kTRUE;
 365     if ((method.CompareTo("ZPC")==0) && (fCentralityZPC >=a && fCentralityZPC < b)) return kTRUE;
 366     if ((method.CompareTo("NPA")==0) && (fCentralityNPA >=a && fCentralityNPA < b)) return kTRUE;
 367     if ((method.CompareTo("V0MvsFMD")==0) && (fCentralityV0MvsFMD >=a && fCentralityV0MvsFMD < b)) return kTRUE;
 368     if ((method.CompareTo("TKLvsV0M")==0) && (fCentralityTKLvsV0M >=a && fCentralityTKLvsV0M < b)) return kTRUE;
 369     if ((method.CompareTo("ZEMvsZDC")==0) && (fCentralityZEMvsZDC >=a && fCentralityZEMvsZDC < b)) return kTRUE;
 370     if ((method.CompareTo("V0Mtrue")==0) && (fCentralityV0Mtrue >=a && fCentralityV0Mtrue < b)) return kTRUE;
 371     if ((method.CompareTo("V0Atrue")==0) && (fCentralityV0Atrue >=a && fCentralityV0Atrue < b)) return kTRUE;
 372     if ((method.CompareTo("V0Ctrue")==0) && (fCentralityV0Ctrue >=a && fCentralityV0Ctrue < b)) return kTRUE;
 373     if ((method.CompareTo("V0MEqtrue")==0) && (fCentralityV0MEqtrue >=a && fCentralityV0MEqtrue < b)) return kTRUE;
 374     if ((method.CompareTo("V0AEqtrue")==0) && (fCentralityV0AEqtrue >=a && fCentralityV0AEqtrue < b)) return kTRUE;
 375     if ((method.CompareTo("V0CEqtrue")==0) && (fCentralityV0CEqtrue >=a && fCentralityV0CEqtrue < b)) return kTRUE;
 376     if ((method.CompareTo("FMDtrue")==0) && (fCentralityFMDtrue >=a && fCentralityFMDtrue < b)) return kTRUE;
 377     if ((method.CompareTo("TRKtrue")==0) && (fCentralityTRKtrue >=a && fCentralityTRKtrue < b)) return kTRUE;
 378     if ((method.CompareTo("TKLtrue")==0) && (fCentralityTKLtrue >=a && fCentralityTKLtrue < b)) return kTRUE;
 379     if ((method.CompareTo("CL0true")==0) && (fCentralityCL0true >=a && fCentralityCL0true < b)) return kTRUE;
 380     if ((method.CompareTo("CL1true")==0) && (fCentralityCL1true >=a && fCentralityCL1true < b)) return kTRUE;
 381     if ((method.CompareTo("CNDtrue")==0) && (fCentralityCNDtrue >=a && fCentralityCNDtrue < b)) return kTRUE;
 382     if ((method.CompareTo("ZNAtrue")==0) && (fCentralityZNAtrue >=a && fCentralityZNAtrue < b)) return kTRUE;
 383     if ((method.CompareTo("ZNCtrue")==0) && (fCentralityZNCtrue >=a && fCentralityZNCtrue < b)) return kTRUE;
 384     if ((method.CompareTo("ZPAtrue")==0) && (fCentralityZPAtrue >=a && fCentralityZPAtrue < b)) return kTRUE;
 385     if ((method.CompareTo("ZPCtrue")==0) && (fCentralityZPCtrue >=a && fCentralityZPCtrue < b)) return kTRUE;
 386     else return kFALSE;
 387   } else {
 388     return kFALSE;
 389   }
 390 }
 391
 392 Float_t AliCentrality::GetCentralityPercentileUnchecked(const char *x) const
 393 {
 394 // Return the centrality percentile
 395   TString method = x;
 396   if(method.CompareTo("V0M")==0)      return fCentralityV0M;
 397   if(method.CompareTo("V0A")==0)      return fCentralityV0A;
 398   if(method.CompareTo("V0A0")==0)     return fCentralityV0A0;
 399   if(method.CompareTo("V0A123")==0)   return fCentralityV0A123;
 400   if(method.CompareTo("V0C")==0)      return fCentralityV0C;
 401   if(method.CompareTo("V0A23")==0)    return fCentralityV0A23;
 402   if(method.CompareTo("V0C01")==0)    return fCentralityV0C01;
 403   if(method.CompareTo("V0S")==0)      return fCentralityV0S;
 404   if(method.CompareTo("V0MEq")==0)    return fCentralityV0MEq;
 405   if(method.CompareTo("V0AEq")==0)    return fCentralityV0AEq;
 406   if(method.CompareTo("V0CEq")==0)    return fCentralityV0CEq;
 407   if(method.CompareTo("FMD")==0)      return fCentralityFMD;
 408   if(method.CompareTo("TRK")==0)      return fCentralityTRK;
 409   if(method.CompareTo("TKL")==0)      return fCentralityTKL;
 410   if(method.CompareTo("CL0")==0)      return fCentralityCL0;
 411   if(method.CompareTo("CL1")==0)      return fCentralityCL1;
 412   if(method.CompareTo("CND")==0)      return fCentralityCND;
 413   if(method.CompareTo("ZNA")==0)      return fCentralityZNA;
 414   if(method.CompareTo("ZNC")==0)      return fCentralityZNC;
 415   if(method.CompareTo("ZPA")==0)      return fCentralityZPA;
 416   if(method.CompareTo("ZPC")==0)      return fCentralityZPC;
 417   if(method.CompareTo("NPA")==0)      return fCentralityNPA;
 418   if(method.CompareTo("V0MvsFMD")==0) return fCentralityV0MvsFMD;
 419   if(method.CompareTo("TKLvsV0M")==0) return fCentralityTKLvsV0M;
 420   if(method.CompareTo("ZEMvsZDC")==0) return fCentralityZEMvsZDC;
 421   if(method.CompareTo("V0Mtrue")==0)  return fCentralityV0Mtrue;
 422   if(method.CompareTo("V0Atrue")==0)  return fCentralityV0Atrue;
 423   if(method.CompareTo("V0Ctrue")==0)  return fCentralityV0Ctrue;
 424   if(method.CompareTo("V0MEqtrue")==0)    return fCentralityV0MEqtrue;
 425   if(method.CompareTo("V0AEqtrue")==0)    return fCentralityV0AEqtrue;
 426   if(method.CompareTo("V0CEqtrue")==0)    return fCentralityV0CEqtrue;
 427   if(method.CompareTo("FMDtrue")==0)  return fCentralityFMDtrue;
 428   if(method.CompareTo("TRKtrue")==0)  return fCentralityTRKtrue;
 429   if(method.CompareTo("TKLtrue")==0)  return fCentralityTKLtrue;
 430   if(method.CompareTo("CL0true")==0)  return fCentralityCL0true;
 431   if(method.CompareTo("CL1true")==0)  return fCentralityCL1true;
 432   if(method.CompareTo("CNDtrue")==0)  return fCentralityCNDtrue;
 433   if(method.CompareTo("ZNAtrue")==0)  return fCentralityZNAtrue;
 434   if(method.CompareTo("ZNCtrue")==0)  return fCentralityZNCtrue;
 435   if(method.CompareTo("ZPAtrue")==0)  return fCentralityZPAtrue;
 436   if(method.CompareTo("ZPCtrue")==0)  return fCentralityZPCtrue;
 437   return -1;
 438 }
 439
 440 Int_t AliCentrality::GetCentralityClass10Unchecked(const char *x) const
 441 {
 442 // Return the centrality class
 443   TString method = x;
 444     if(method.CompareTo("V0M")==0)      return (Int_t) (fCentralityV0M / 10.0);
 445     if(method.CompareTo("V0A")==0)      return (Int_t) (fCentralityV0A / 10.0);
 446     if(method.CompareTo("V0A0")==0)     return (Int_t) (fCentralityV0A0 / 10.0);
 447     if(method.CompareTo("V0C")==0)      return (Int_t) (fCentralityV0C / 10.0);
 448     if(method.CompareTo("V0A23")==0)    return (Int_t) (fCentralityV0A23 / 10.0);
 449     if(method.CompareTo("V0C01")==0)    return (Int_t) (fCentralityV0C01 / 10.0);
 450     if(method.CompareTo("V0S")==0)      return (Int_t) (fCentralityV0S / 10.0);
 451     if(method.CompareTo("V0MEq")==0)    return (Int_t) (fCentralityV0MEq / 10.0);
 452     if(method.CompareTo("V0AEq")==0)    return (Int_t) (fCentralityV0AEq / 10.0);
 453     if(method.CompareTo("V0CEq")==0)    return (Int_t) (fCentralityV0CEq / 10.0);
 454     if(method.CompareTo("FMD")==0)      return (Int_t) (fCentralityFMD / 10.0);
 455     if(method.CompareTo("TRK")==0)      return (Int_t) (fCentralityTRK / 10.0);
 456     if(method.CompareTo("TKL")==0)      return (Int_t) (fCentralityTKL / 10.0);
 457     if(method.CompareTo("CL0")==0)      return (Int_t) (fCentralityCL0 / 10.0);
 458     if(method.CompareTo("CL1")==0)      return (Int_t) (fCentralityCL1 / 10.0);
 459     if(method.CompareTo("CND")==0)      return (Int_t) (fCentralityCND / 10.0);
 460     if(method.CompareTo("ZNA")==0)      return (Int_t) (fCentralityZNA / 10.0);
 461     if(method.CompareTo("ZNC")==0)      return (Int_t) (fCentralityZNC / 10.0);
 462     if(method.CompareTo("ZPA")==0)      return (Int_t) (fCentralityZPA / 10.0);
 463     if(method.CompareTo("ZPC")==0)      return (Int_t) (fCentralityZPC / 10.0);
 464     if(method.CompareTo("NPA")==0)      return (Int_t) (fCentralityNPA / 10.0);
 465     if(method.CompareTo("V0MvsFMD")==0) return (Int_t) (fCentralityV0MvsFMD / 10.0);
 466     if(method.CompareTo("TKLvsV0M")==0) return (Int_t) (fCentralityTKLvsV0M / 10.0);
 467     if(method.CompareTo("ZEMvsZDC")==0) return (Int_t) (fCentralityZEMvsZDC / 10.0);
 468     if(method.CompareTo("V0Mtrue")==0)  return (Int_t) (fCentralityV0Mtrue / 10.0);
 469     if(method.CompareTo("V0Atrue")==0)  return (Int_t) (fCentralityV0Atrue / 10.0);
 470     if(method.CompareTo("V0Ctrue")==0)  return (Int_t) (fCentralityV0Ctrue / 10.0);
 471     if(method.CompareTo("V0MEqtrue")==0)return (Int_t) (fCentralityV0MEqtrue / 10.0);
 472     if(method.CompareTo("V0AEqtrue")==0)return (Int_t) (fCentralityV0AEqtrue / 10.0);
 473     if(method.CompareTo("V0CEqtrue")==0)return (Int_t) (fCentralityV0CEqtrue / 10.0);
 474     if(method.CompareTo("FMDtrue")==0)  return (Int_t) (fCentralityFMDtrue / 10.0);
 475     if(method.CompareTo("TRKtrue")==0)  return (Int_t) (fCentralityTRKtrue / 10.0);
 476     if(method.CompareTo("TKLtrue")==0)  return (Int_t) (fCentralityTKLtrue / 10.0);
 477     if(method.CompareTo("CL0true")==0)  return (Int_t) (fCentralityCL0true / 10.0);
 478     if(method.CompareTo("CL1true")==0)  return (Int_t) (fCentralityCL1true / 10.0);
 479     if(method.CompareTo("CNDtrue")==0)  return (Int_t) (fCentralityCNDtrue / 10.0);
 480     if(method.CompareTo("ZNAtrue")==0)  return (Int_t) (fCentralityZNAtrue / 10.0);
 481     if(method.CompareTo("ZNCtrue")==0)  return (Int_t) (fCentralityZNCtrue / 10.0);
 482     if(method.CompareTo("ZPAtrue")==0)  return (Int_t) (fCentralityZPAtrue / 10.0);
 483     if(method.CompareTo("ZPCtrue")==0)  return (Int_t) (fCentralityZPCtrue / 10.0);
 484   return -1;
 485 }
 486
 487 Int_t AliCentrality::GetCentralityClass5Unchecked(const char *x) const
 488 {
 489 // Return the centrality class
 490   TString method = x;
 491     if(method.CompareTo("V0M")==0)      return (Int_t) (fCentralityV0M / 5.0);
 492     if(method.CompareTo("V0A")==0)      return (Int_t) (fCentralityV0A / 5.0);
 493     if(method.CompareTo("V0A0")==0)     return (Int_t) (fCentralityV0A0 / 5.0);
 494     if(method.CompareTo("V0A123")==0)   return (Int_t) (fCentralityV0A123 / 5.0);
 495     if(method.CompareTo("V0C")==0)      return (Int_t) (fCentralityV0C / 5.0);
 496     if(method.CompareTo("V0A23")==0)    return (Int_t) (fCentralityV0A23 / 5.0);
 497     if(method.CompareTo("V0C01")==0)    return (Int_t) (fCentralityV0C01 / 5.0);
 498     if(method.CompareTo("V0S")==0)      return (Int_t) (fCentralityV0S / 5.0);
 499     if(method.CompareTo("V0MEq")==0)    return (Int_t) (fCentralityV0MEq / 5.0);
 500     if(method.CompareTo("V0AEq")==0)    return (Int_t) (fCentralityV0AEq / 5.0);
 501     if(method.CompareTo("V0CEq")==0)    return (Int_t) (fCentralityV0CEq / 5.0);
 502     if(method.CompareTo("FMD")==0)      return (Int_t) (fCentralityFMD / 5.0);
 503     if(method.CompareTo("TRK")==0)      return (Int_t) (fCentralityTRK / 5.0);
 504     if(method.CompareTo("TKL")==0)      return (Int_t) (fCentralityTKL / 5.0);
 505     if(method.CompareTo("CL0")==0)      return (Int_t) (fCentralityCL0 / 5.0);
 506     if(method.CompareTo("CL1")==0)      return (Int_t) (fCentralityCL1 / 5.0);
 507     if(method.CompareTo("CND")==0)      return (Int_t) (fCentralityCND / 5.0);
 508     if(method.CompareTo("ZNA")==0)      return (Int_t) (fCentralityZNA / 5.0);
 509     if(method.CompareTo("ZNC")==0)      return (Int_t) (fCentralityZNC / 5.0);
 510     if(method.CompareTo("ZPA")==0)      return (Int_t) (fCentralityZPA / 5.0);
 511     if(method.CompareTo("ZPC")==0)      return (Int_t) (fCentralityZPC / 5.0);
 512     if(method.CompareTo("NPA")==0)      return (Int_t) (fCentralityNPA / 5.0);
 513     if(method.CompareTo("V0MvsFMD")==0) return (Int_t) (fCentralityV0MvsFMD / 5.0);
 514     if(method.CompareTo("TKLvsV0M")==0) return (Int_t) (fCentralityTKLvsV0M / 5.0);
 515     if(method.CompareTo("ZEMvsZDC")==0) return (Int_t) (fCentralityZEMvsZDC / 5.0);
 516     if(method.CompareTo("V0Mtrue")==0)  return (Int_t) (fCentralityV0Mtrue / 5.0);
 517     if(method.CompareTo("V0Atrue")==0)  return (Int_t) (fCentralityV0Atrue / 5.0);
 518     if(method.CompareTo("V0Ctrue")==0)  return (Int_t) (fCentralityV0Ctrue / 5.0);
 519     if(method.CompareTo("V0MEqtrue")==0)return (Int_t) (fCentralityV0MEqtrue / 5.0);
 520     if(method.CompareTo("V0AEqtrue")==0)return (Int_t) (fCentralityV0AEqtrue / 5.0);
 521     if(method.CompareTo("V0CEqtrue")==0)return (Int_t) (fCentralityV0CEqtrue / 5.0);
 522     if(method.CompareTo("FMDtrue")==0)  return (Int_t) (fCentralityFMDtrue / 5.0);
 523     if(method.CompareTo("TRKtrue")==0)  return (Int_t) (fCentralityTRKtrue / 5.0);
 524     if(method.CompareTo("TKLtrue")==0)  return (Int_t) (fCentralityTKLtrue / 5.0);
 525     if(method.CompareTo("CL0true")==0)  return (Int_t) (fCentralityCL0true / 5.0);
 526     if(method.CompareTo("CL1true")==0)  return (Int_t) (fCentralityCL1true / 5.0);
 527     if(method.CompareTo("CNDtrue")==0)  return (Int_t) (fCentralityCNDtrue / 5.0);
 528     if(method.CompareTo("ZNAtrue")==0)  return (Int_t) (fCentralityZNAtrue / 5.0);
 529     if(method.CompareTo("ZNCtrue")==0)  return (Int_t) (fCentralityZNCtrue / 5.0);
 530     if(method.CompareTo("ZPAtrue")==0)  return (Int_t) (fCentralityZPAtrue / 5.0);
 531     if(method.CompareTo("ZPCtrue")==0)  return (Int_t) (fCentralityZPCtrue / 5.0);
 532   return -1;
 533 }
 534
 535 Bool_t AliCentrality::IsEventInCentralityClassUnchecked(Float_t a, Float_t b, const char *x) const
 536 {
 537 // True if event inside given centrality class
 538   TString method = x;
 539     if ((method.CompareTo("V0M")==0) && (fCentralityV0M >=a && fCentralityV0M < b)) return kTRUE;
 540     if ((method.CompareTo("V0A")==0) && (fCentralityV0A >=a && fCentralityV0A < b)) return kTRUE;
 541     if ((method.CompareTo("V0A0")==0) && (fCentralityV0A0 >=a && fCentralityV0A0 < b)) return kTRUE;
 542     if ((method.CompareTo("V0A123")==0) && (fCentralityV0A123 >=a && fCentralityV0A123 < b)) return kTRUE;
 543     if ((method.CompareTo("V0C")==0) && (fCentralityV0C >=a && fCentralityV0C < b)) return kTRUE;
 544     if ((method.CompareTo("V0A23")==0) && (fCentralityV0A23 >=a && fCentralityV0A23 < b)) return kTRUE;
 545     if ((method.CompareTo("V0C01")==0) && (fCentralityV0C01 >=a && fCentralityV0C01 < b)) return kTRUE;
 546     if ((method.CompareTo("V0S")==0) && (fCentralityV0S >=a && fCentralityV0S < b)) return kTRUE;
 547     if ((method.CompareTo("V0MEq")==0) && (fCentralityV0MEq >=a && fCentralityV0MEq < b)) return kTRUE;
 548     if ((method.CompareTo("V0AEq")==0) && (fCentralityV0AEq >=a && fCentralityV0AEq < b)) return kTRUE;
 549     if ((method.CompareTo("V0CEq")==0) && (fCentralityV0CEq >=a && fCentralityV0CEq < b)) return kTRUE;
 550     if ((method.CompareTo("FMD")==0) && (fCentralityFMD >=a && fCentralityFMD < b)) return kTRUE;
 551     if ((method.CompareTo("TRK")==0) && (fCentralityTRK >=a && fCentralityTRK < b)) return kTRUE;
 552     if ((method.CompareTo("TKL")==0) && (fCentralityTKL >=a && fCentralityTKL < b)) return kTRUE;
 553     if ((method.CompareTo("CL0")==0) && (fCentralityCL0 >=a && fCentralityCL0 < b)) return kTRUE;
 554     if ((method.CompareTo("CL1")==0) && (fCentralityCL1 >=a && fCentralityCL1 < b)) return kTRUE;
 555     if ((method.CompareTo("CND")==0) && (fCentralityCND >=a && fCentralityCND < b)) return kTRUE;
 556     if ((method.CompareTo("ZNA")==0) && (fCentralityZNA >=a && fCentralityZNA < b)) return kTRUE;
 557     if ((method.CompareTo("ZNC")==0) && (fCentralityZNC >=a && fCentralityZNC < b)) return kTRUE;
 558     if ((method.CompareTo("ZPA")==0) && (fCentralityZPA >=a && fCentralityZPA < b)) return kTRUE;
 559     if ((method.CompareTo("ZPC")==0) && (fCentralityZPC >=a && fCentralityZPC < b)) return kTRUE;
 560     if ((method.CompareTo("NPA")==0) && (fCentralityNPA >=a && fCentralityNPA < b)) return kTRUE;
 561     if ((method.CompareTo("V0MvsFMD")==0) && (fCentralityV0MvsFMD >=a && fCentralityV0MvsFMD < b)) return kTRUE;
 562     if ((method.CompareTo("TKLvsV0M")==0) && (fCentralityTKLvsV0M >=a && fCentralityTKLvsV0M < b)) return kTRUE;
 563     if ((method.CompareTo("ZEMvsZDC")==0) && (fCentralityZEMvsZDC >=a && fCentralityZEMvsZDC < b)) return kTRUE;
 564     if ((method.CompareTo("V0Mtrue")==0) && (fCentralityV0Mtrue >=a && fCentralityV0Mtrue < b)) return kTRUE;
 565     if ((method.CompareTo("V0Atrue")==0) && (fCentralityV0Atrue >=a && fCentralityV0Atrue < b)) return kTRUE;
 566     if ((method.CompareTo("V0Ctrue")==0) && (fCentralityV0Ctrue >=a && fCentralityV0Ctrue < b)) return kTRUE;
 567     if ((method.CompareTo("V0MEqtrue")==0) && (fCentralityV0MEqtrue >=a && fCentralityV0MEqtrue < b)) return kTRUE;
 568     if ((method.CompareTo("V0AEqtrue")==0) && (fCentralityV0AEqtrue >=a && fCentralityV0AEqtrue < b)) return kTRUE;
 569     if ((method.CompareTo("V0CEqtrue")==0) && (fCentralityV0CEqtrue >=a && fCentralityV0CEqtrue < b)) return kTRUE;
 570     if ((method.CompareTo("FMDtrue")==0) && (fCentralityFMDtrue >=a && fCentralityFMDtrue < b)) return kTRUE;
 571     if ((method.CompareTo("TRKtrue")==0) && (fCentralityTRKtrue >=a && fCentralityTRKtrue < b)) return kTRUE;
 572     if ((method.CompareTo("TKLtrue")==0) && (fCentralityTKLtrue >=a && fCentralityTKLtrue < b)) return kTRUE;
 573     if ((method.CompareTo("CL0true")==0) && (fCentralityCL0true >=a && fCentralityCL0true < b)) return kTRUE;
 574     if ((method.CompareTo("CL1true")==0) && (fCentralityCL1true >=a && fCentralityCL1true < b)) return kTRUE;
 575     if ((method.CompareTo("CNDtrue")==0) && (fCentralityCNDtrue >=a && fCentralityCNDtrue < b)) return kTRUE;
 576     if ((method.CompareTo("ZNAtrue")==0) && (fCentralityZNAtrue >=a && fCentralityZNAtrue < b)) return kTRUE;
 577     if ((method.CompareTo("ZNCtrue")==0) && (fCentralityZNCtrue >=a && fCentralityZNCtrue < b)) return kTRUE;
 578     if ((method.CompareTo("ZPAtrue")==0) && (fCentralityZPAtrue >=a && fCentralityZPAtrue < b)) return kTRUE;
 579     if ((method.CompareTo("ZPCtrue")==0) && (fCentralityZPCtrue >=a && fCentralityZPCtrue < b)) return kTRUE;
 580   else return kFALSE;
 581 }
 582
 583 void AliCentrality::Reset()
 584 {
 585 // Reset.
 586
 587   fQuality            =  999;
 588   fCentralityV0M      =  0;
 589   fCentralityV0A      =  0;
 590   fCentralityV0A0     =  0;
 591   fCentralityV0A123   =  0;
 592   fCentralityV0C      =  0;
 593   fCentralityV0A23    =  0;
 594   fCentralityV0C01    =  0;
 595   fCentralityV0S      =  0;
 596   fCentralityV0MEq    =  0;
 597   fCentralityV0AEq    =  0;
 598   fCentralityV0CEq    =  0;
 599   fCentralityFMD      =  0;
 600   fCentralityTRK      =  0;
 601   fCentralityTKL      =  0;
 602   fCentralityCL0      =  0;
 603   fCentralityCL1      =  0;
 604   fCentralityCND      =  0;
 605   fCentralityZNA      =  0;
 606   fCentralityZNC      =  0;
 607   fCentralityZPA      =  0;
 608   fCentralityZPC      =  0;
 609   fCentralityNPA      =  0;
 610   fCentralityV0MvsFMD =  0;
 611   fCentralityTKLvsV0M =  0;
 612   fCentralityZEMvsZDC =  0;
 613   fCentralityV0Mtrue  =  0;
 614   fCentralityV0Atrue  =  0;
 615   fCentralityV0Ctrue  =  0;
 616   fCentralityV0MEqtrue  =  0;
 617   fCentralityV0AEqtrue  =  0;
 618   fCentralityV0CEqtrue  =  0;
 619   fCentralityFMDtrue  =  0;
 620   fCentralityTRKtrue  =  0;
 621   fCentralityTKLtrue  =  0;
 622   fCentralityCL0true  =  0;
 623   fCentralityCL1true  =  0;
 624   fCentralityCNDtrue  =  0;
 625   fCentralityZNAtrue  =  0;
 626   fCentralityZNCtrue  =  0;
 627   fCentralityZPAtrue  =  0;
 628   fCentralityZPCtrue  =  0;
 629 }