fV0MSPDSigmaOutlierPar2(0),
fV0MTPCSigmaOutlierPar0(0),
fV0MTPCSigmaOutlierPar1(0),
- fV0MTPCSigmaOutlierPar2(0),
+ fV0MTPCSigmaOutlierPar2(0),
fV0MZDCOutlierPar0(0),
fV0MZDCOutlierPar1(0),
fV0MZDCEcalOutlierPar0(0),
fEJE(0),
fEGA(0),
fPHS(0),
+ fMB(0),
fCVHN(0),
fCVLN(0),
fCVHNbit(0),
fCCENTbit(0),
fCSEMIbit(0),
fCentV0M(0),
+ fCentV0A(0),
+ fCentV0C(0),
fCentFMD(0),
fCentTRK(0),
fCentTKL(0),
fCentCL0(0),
fCentCL1(0),
+ fCentCND(0),
fCentV0MvsFMD(0),
fCentTKLvsV0M(0),
fCentZEMvsZDC(0),
fHtempV0M(0),
+ fHtempV0A(0),
+ fHtempV0C(0),
fHtempFMD(0),
fHtempTRK(0),
fHtempTKL(0),
fHtempCL0(0),
fHtempCL1(0),
+ fHtempCND(0),
fHtempV0MvsFMD(0),
fHtempTKLvsV0M(0),
fHtempZEMvsZDC(0),
fOutputList(0),
- fHOutCentV0M (0),
+ fHOutCentV0M (0),
+ fHOutCentV0A (0),
+ fHOutCentV0C (0),
fHOutCentV0MCVHN(0),
fHOutCentV0MCVLN(0),
fHOutCentV0MCVHNinMB(0),
fHOutCentTKL (0),
fHOutCentCL0 (0),
fHOutCentCL1 (0),
+ fHOutCentCND (0),
fHOutCentV0MvsFMD(0),
fHOutCentTKLvsV0M(0),
fHOutCentZEMvsZDC(0),
fHOutCentV0MvsCentTRK(0),
fHOutCentTRKvsCentCL1(0),
fHOutCentV0MvsCentZDC(0),
+ fHOutCentV0AvsCentV0C(0),
+ fHOutCentV0AvsCentTRK(0),
+ fHOutCentV0AvsCentCND(0),
+ fHOutCentV0AvsCentCL1(0),
+ fHOutCentV0CvsCentTRK(0),
+ fHOutCentV0CvsCentCND(0),
+ fHOutCentV0CvsCentCL1(0),
+ fHOutMultV0AC(0),
fHOutMultV0M(0),
+ fHOutMultV0A(0),
+ fHOutMultV0C(0),
fHOutMultV0O(0),
fHOutMultFMD(0),
fHOutMultTRK(0),
fHOutMultTKL(0),
fHOutMultCL0(0),
fHOutMultCL1(0),
+ fHOutMultCND(0),
fHOutMultV0MvsZDN(0),
fHOutMultZEMvsZDN(0),
fHOutMultV0MvsZDC(0),
fHOutMultV0MvsV0O(0),
fHOutMultV0OvsCL1(0),
fHOutMultV0OvsTRK(0),
+ fHOutMultCL1vsTKL(0),
fHOutCentV0Mqual1(0),
fHOutCentTRKqual1(0),
fHOutCentCL1qual1(0),
fV0MSPDSigmaOutlierPar2(0),
fV0MTPCSigmaOutlierPar0(0),
fV0MTPCSigmaOutlierPar1(0),
- fV0MTPCSigmaOutlierPar2(0),
+ fV0MTPCSigmaOutlierPar2(0),
fV0MZDCOutlierPar0(0),
fV0MZDCOutlierPar1(0),
fV0MZDCEcalOutlierPar0(0),
fEJE(0),
fEGA(0),
fPHS(0),
+ fMB(0),
fCVHN(0),
fCVLN(0),
fCVHNbit(0),
fCCENTbit(0),
fCSEMIbit(0),
fCentV0M(0),
+ fCentV0A(0),
+ fCentV0C(0),
fCentFMD(0),
fCentTRK(0),
fCentTKL(0),
fCentCL0(0),
fCentCL1(0),
+ fCentCND(0),
fCentV0MvsFMD(0),
fCentTKLvsV0M(0),
fCentZEMvsZDC(0),
fHtempV0M(0),
+ fHtempV0A(0),
+ fHtempV0C(0),
fHtempFMD(0),
fHtempTRK(0),
fHtempTKL(0),
fHtempCL0(0),
fHtempCL1(0),
+ fHtempCND(0),
fHtempV0MvsFMD(0),
fHtempTKLvsV0M(0),
fHtempZEMvsZDC(0),
fOutputList(0),
- fHOutCentV0M (0),
+ fHOutCentV0M (0),
+ fHOutCentV0A (0),
+ fHOutCentV0C (0),
fHOutCentV0MCVHN(0),
fHOutCentV0MCVLN(0),
fHOutCentV0MCVHNinMB(0),
fHOutCentTKL (0),
fHOutCentCL0 (0),
fHOutCentCL1 (0),
+ fHOutCentCND (0),
fHOutCentV0MvsFMD(0),
fHOutCentTKLvsV0M(0),
fHOutCentZEMvsZDC(0),
fHOutCentV0MvsCentTRK(0),
fHOutCentTRKvsCentCL1(0),
fHOutCentV0MvsCentZDC(0),
+ fHOutCentV0AvsCentV0C(0),
+ fHOutCentV0AvsCentTRK(0),
+ fHOutCentV0AvsCentCND(0),
+ fHOutCentV0AvsCentCL1(0),
+ fHOutCentV0CvsCentTRK(0),
+ fHOutCentV0CvsCentCND(0),
+ fHOutCentV0CvsCentCL1(0),
+ fHOutMultV0AC(0),
fHOutMultV0M(0),
+ fHOutMultV0A(0),
+ fHOutMultV0C(0),
fHOutMultV0O(0),
fHOutMultFMD(0),
fHOutMultTRK(0),
fHOutMultTKL(0),
fHOutMultCL0(0),
fHOutMultCL1(0),
+ fHOutMultCND(0),
fHOutMultV0MvsZDN(0),
fHOutMultZEMvsZDN(0),
fHOutMultV0MvsZDC(0),
fHOutMultV0MvsV0O(0),
fHOutMultV0OvsCL1(0),
fHOutMultV0OvsTRK(0),
+ fHOutMultCL1vsTKL(0),
fHOutCentV0Mqual1(0),
fHOutCentTRKqual1(0),
fHOutCentCL1qual1(0),
fEJE(ana.fEJE),
fEGA(ana.fEGA),
fPHS(ana.fPHS),
+ fMB(ana.fMB),
fCVHN(ana.fCVHN),
fCVLN(ana.fCVLN),
fCVHNbit(ana.fCVHNbit),
fCCENTbit(ana.fCCENTbit),
fCSEMIbit(ana.fCSEMIbit),
fCentV0M(ana.fCentV0M),
+ fCentV0A(ana.fCentV0A),
+ fCentV0C(ana.fCentV0C),
fCentFMD(ana.fCentFMD),
fCentTRK(ana.fCentTRK),
fCentTKL(ana.fCentTKL),
fCentCL0(ana.fCentCL0),
fCentCL1(ana.fCentCL1),
+ fCentCND(ana.fCentCND),
fCentV0MvsFMD(ana.fCentV0MvsFMD),
fCentTKLvsV0M(ana.fCentTKLvsV0M),
fCentZEMvsZDC(ana.fCentZEMvsZDC),
fHtempV0M(ana.fHtempV0M),
+ fHtempV0A(ana.fHtempV0A),
+ fHtempV0C(ana.fHtempV0C),
fHtempFMD(ana.fHtempFMD),
fHtempTRK(ana.fHtempTRK),
fHtempTKL(ana.fHtempTKL),
fHtempCL0(ana.fHtempCL0),
fHtempCL1(ana.fHtempCL1),
+ fHtempCND(ana.fHtempCND),
fHtempV0MvsFMD(ana.fHtempV0MvsFMD),
fHtempTKLvsV0M(ana.fHtempTKLvsV0M),
fHtempZEMvsZDC(ana.fHtempZEMvsZDC),
fOutputList(ana.fOutputList),
- fHOutCentV0M (ana.fHOutCentV0M ),
+ fHOutCentV0M (ana.fHOutCentV0M ),
+ fHOutCentV0A (ana.fHOutCentV0A ),
+ fHOutCentV0C (ana.fHOutCentV0C ),
fHOutCentV0MCVHN(ana.fHOutCentV0MCVHN),
fHOutCentV0MCVLN(ana.fHOutCentV0MCVLN),
fHOutCentV0MCVHNinMB(ana.fHOutCentV0MCVHNinMB),
fHOutCentTKL (ana.fHOutCentTKL ),
fHOutCentCL0 (ana.fHOutCentCL0 ),
fHOutCentCL1 (ana.fHOutCentCL1 ),
+ fHOutCentCND (ana.fHOutCentCND ),
fHOutCentV0MvsFMD(ana.fHOutCentV0MvsFMD),
fHOutCentTKLvsV0M(ana.fHOutCentTKLvsV0M),
fHOutCentZEMvsZDC(ana.fHOutCentZEMvsZDC),
fHOutCentV0MvsCentTRK(ana.fHOutCentV0MvsCentTRK),
fHOutCentTRKvsCentCL1(ana.fHOutCentTRKvsCentCL1),
fHOutCentV0MvsCentZDC(ana.fHOutCentV0MvsCentZDC),
+ fHOutCentV0AvsCentV0C(ana.fHOutCentV0AvsCentV0C),
+ fHOutCentV0AvsCentTRK(ana.fHOutCentV0AvsCentTRK),
+ fHOutCentV0AvsCentCND(ana.fHOutCentV0AvsCentCND),
+ fHOutCentV0AvsCentCL1(ana.fHOutCentV0AvsCentCL1),
+ fHOutCentV0CvsCentTRK(ana.fHOutCentV0CvsCentTRK),
+ fHOutCentV0CvsCentCND(ana.fHOutCentV0CvsCentCND),
+ fHOutCentV0CvsCentCL1(ana.fHOutCentV0CvsCentCL1),
+ fHOutMultV0AC(ana.fHOutMultV0AC),
fHOutMultV0M(ana.fHOutMultV0M),
+ fHOutMultV0A(ana.fHOutMultV0A),
+ fHOutMultV0C(ana.fHOutMultV0C),
fHOutMultV0O(ana.fHOutMultV0O),
fHOutMultFMD(ana.fHOutMultFMD),
fHOutMultTRK(ana.fHOutMultTRK),
fHOutMultTKL(ana.fHOutMultTKL),
fHOutMultCL0(ana.fHOutMultCL0),
fHOutMultCL1(ana.fHOutMultCL1),
+ fHOutMultCND(ana.fHOutMultCND),
fHOutMultV0MvsZDN(ana.fHOutMultV0MvsZDN),
fHOutMultZEMvsZDN(ana.fHOutMultZEMvsZDN),
fHOutMultV0MvsZDC(ana.fHOutMultV0MvsZDC),
fHOutMultV0MvsV0O(ana.fHOutMultV0MvsV0O),
fHOutMultV0OvsCL1(ana.fHOutMultV0OvsCL1),
fHOutMultV0OvsTRK(ana.fHOutMultV0OvsTRK),
+ fHOutMultCL1vsTKL(ana.fHOutMultCL1vsTKL),
fHOutCentV0Mqual1(ana.fHOutCentV0Mqual1),
fHOutCentTRKqual1(ana.fHOutCentTRKqual1),
fHOutCentCL1qual1(ana.fHOutCentCL1qual1),
fOutputList = new TList();
fOutputList->SetOwner();
fHOutCentV0M = new TH1F("fHOutCentV0M","fHOutCentV0M; Centrality V0",505,0,101);
+ fHOutCentV0A = new TH1F("fHOutCentV0A","fHOutCentV0A; Centrality V0A",505,0,101);
+ fHOutCentV0C = new TH1F("fHOutCentV0C","fHOutCentV0C; Centrality V0C",505,0,101);
fHOutCentV0MCVHN= new TH1F("fHOutCentV0M_CVHN","fHOutCentV0M_CVHN; Centrality V0",505,0,101);
fHOutCentV0MCVLN= new TH1F("fHOutCentV0M_CVLN","fHOutCentV0M_CVLN; Centrality V0",505,0,101);
fHOutCentV0MCVHNinMB= new TH1F("fHOutCentV0M_CVHNinMB","fHOutCentV0M_CVHN; Centrality V0",505,0,101);
fHOutCentTKL = new TH1F("fHOutCentTKL","fHOutCentTKL; Centrality tracklets",505,0,101);
fHOutCentCL0 = new TH1F("fHOutCentCL0","fHOutCentCL0; Centrality SPD inner",505,0,101);
fHOutCentCL1 = new TH1F("fHOutCentCL1","fHOutCentCL1; Centrality SPD outer",505,0,101);
+ fHOutCentCND = new TH1F("fHOutCentCND","fHOutCentCND; Centrality candle",505,0,101);
fHOutCentV0MvsFMD= new TH1F("fHOutCentV0MvsFMD","fHOutCentV0MvsFMD; Centrality V0 vs FMD",505,0,101);
fHOutCentTKLvsV0M= new TH1F("fHOutCentTKLvsV0M","fHOutCentTKLvsV0M; Centrality tracklets vs V0",505,0,101);
fHOutCentZEMvsZDC= new TH1F("fHOutCentZEMvsZDC","fHOutCentZEMvsZDC; Centrality ZEM vs ZDC",505,0,101);
fHOutCentV0MvsCentTRK= new TH2F("fHOutCentV0MvsCentTRK","fHOutCentV0MvsCentTRK; Cent V0; Cent TPC",505,0,101,505,0,101);
fHOutCentTRKvsCentCL1= new TH2F("fHOutCentTRKvsCentCL1","fHOutCentTRKvsCentCL1; Cent TPC; Cent SPD",505,0,101,505,0,101);
fHOutCentV0MvsCentZDC= new TH2F("fHOutCentV0MvsCentZDC","fHOutCentV0MvsCentZDC; Cent V0; Cent ZDC",505,0,101,505,0,101);
- fHOutMultV0M = new TH1F("fHOutMultV0M","fHOutMultV0M; Multiplicity V0",25000,0,30000);
- fHOutMultV0O = new TH1F("fHOutMultV0O","fHOutMultV0O; Multiplicity V0",40000,0,40000);
+ fHOutCentV0AvsCentV0C= new TH2F("fHOutCentV0AvsCentV0C","fHOutCentV0AvsCentV0C; Cent V0A; Cent V0C;", 505,0,101,505,0,101);
+ fHOutCentV0AvsCentTRK= new TH2F("fHOutCentV0AvsCentTRK","fHOutCentV0AvsCentTRK; Cent V0A; Cent TRK;", 505,0,101,505,0,101);
+ fHOutCentV0AvsCentCND= new TH2F("fHOutCentV0AvsCentCND","fHOutCentV0AvsCentCND; Cent V0A; Cent CND;", 505,0,101,505,0,101);
+ fHOutCentV0AvsCentCL1= new TH2F("fHOutCentV0AvsCentCL1","fHOutCentV0AvsCentCL1; Cent V0A; Cent CL1;", 505,0,101,505,0,101);
+ fHOutCentV0CvsCentTRK= new TH2F("fHOutCentV0CvsCentTRK","fHOutCentV0CvsCentTRK; Cent V0C; Cent TRK;", 505,0,101,505,0,101);
+ fHOutCentV0CvsCentCND= new TH2F("fHOutCentV0CvsCentCND","fHOutCentV0CvsCentCND; Cent V0C; Cent CND;", 505,0,101,505,0,101);
+ fHOutCentV0CvsCentCL1= new TH2F("fHOutCentV0CvsCentCL1","fHOutCentV0CvsCentCL1; Cent V0C; Cent CL1;", 505,0,101,505,0,101);
+
+ fHOutMultV0AC = new TH2F("fHOutMultV0AC","fHOutMultV0AC; Multiplicity V0A; Multiplicity V0C",1000,0,1000,1000,0,1000);
+ fHOutMultV0M = new TH1F("fHOutMultV0M","fHOutMultV0M; Multiplicity V0",25000,0,25000);
+ fHOutMultV0A = new TH1F("fHOutMultV0A","fHOutMultV0A; Multiplicity V0",25000,0,25000);
+ fHOutMultV0C = new TH1F("fHOutMultV0C","fHOutMultV0C; Multiplicity V0",25000,0,25000);
+ fHOutMultV0O = new TH1F("fHOutMultV0O","fHOutMultV0O; Multiplicity V0",40000,0,40000);
fHOutMultFMD = new TH1F("fHOutMultFMD","fHOutMultFMD; Multiplicity FMD",24000,0,24000);
fHOutMultTRK = new TH1F("fHOutMultTRK","fHOutMultTRK; Multiplicity TPC",4000,0,4000);
fHOutMultTKL = new TH1F("fHOutMultTKL","fHOutMultTKL; Multiplicity tracklets",5000,0,5000);
fHOutMultCL0 = new TH1F("fHOutMultCL0","fHOutMultCL0; Multiplicity SPD inner",7000,0,7000);
fHOutMultCL1 = new TH1F("fHOutMultCL1","fHOutMultCL1; Multiplicity SPD outer",7000,0,7000);
+ fHOutMultCND = new TH1F("fHOutMultCND","fHOutMultCND; Multiplicity candle",4000,0,4000);
+
fHOutMultV0MvsZDN = new TH2F("fHOutMultV0MvsZDN","fHOutMultV0MvsZDN; Multiplicity V0; Energy ZDC-N",500,0,30000,500,0,180000);
fHOutMultZEMvsZDN = new TH2F("fHOutMultZEMvsZDN","fHOutMultZEMvsZDN; Energy ZEM; Energy ZDC-N",500,0,2500,500,0,180000);
fHOutMultV0MvsZDC = new TH2F("fHOutMultV0MvsZDC","fHOutMultV0MvsZDC; Multiplicity V0; Energy ZDC",500,0,30000,500,0,200000);
fHOutMultV0MvsV0O = new TH2F("fHOutMultV0MvsV0O","fHOutMultV0MvsV0O; Multiplicity V0; Multiplicity V0 Online",500,0,30000,500,0,30000);
fHOutMultV0OvsCL1 = new TH2F("fHOutMultV0OvsCL1","fHOutMultV0OvsCL1; Multiplicity V0; Multiplicity SPD outer",2500,0,30000,700,0,7000);
fHOutMultV0OvsTRK = new TH2F("fHOutMultV0OvsTRK","fHOutMultV0OvsTRK; Multiplicity V0; Multiplicity TPC",2500,0,30000,400,0,4000);
+ fHOutMultCL1vsTKL = new TH2F ("fHOutMultCL1vsTKL","fHOutMultCL1vsTKL; Multiplicity SPD outer; Multiplicity tracklets",700,0,7000,700,0,7000);
fHOutCentV0Mqual1 = new TH1F("fHOutCentV0M_qual1","fHOutCentV0M_qual1; Centrality V0",505,0,101);
fHOutCentTRKqual1 = new TH1F("fHOutCentTRK_qual1","fHOutCentTRK_qual1; Centrality TPC",505,0,101);
fHOutVertex = new TH1F("fHOutVertex", "fHOutVertex", 100,-20,20);
fHOutVertexT0 = new TH1F("fHOutVertexT0", "fHOutVertexT0", 100,-20,20);
- fOutputList->Add( fHOutCentV0M );
+ fOutputList->Add( fHOutCentV0M );
+ fOutputList->Add( fHOutCentV0A );
+ fOutputList->Add( fHOutCentV0C );
fOutputList->Add( fHOutCentV0MCVHN);
fOutputList->Add( fHOutCentV0MCVLN);
fOutputList->Add( fHOutCentV0MCVHNinMB);
fOutputList->Add( fHOutCentTKL );
fOutputList->Add( fHOutCentCL0 );
fOutputList->Add( fHOutCentCL1 );
+ fOutputList->Add( fHOutCentCND );
fOutputList->Add( fHOutCentV0MvsFMD);
fOutputList->Add( fHOutCentTKLvsV0M);
fOutputList->Add( fHOutCentZEMvsZDC);
fOutputList->Add( fHOutCentV0MvsCentTRK);
fOutputList->Add( fHOutCentTRKvsCentCL1);
fOutputList->Add( fHOutCentV0MvsCentZDC);
+ fOutputList->Add( fHOutCentV0AvsCentV0C);
+ fOutputList->Add( fHOutCentV0AvsCentTRK);
+ fOutputList->Add( fHOutCentV0AvsCentCND);
+ fOutputList->Add( fHOutCentV0AvsCentCL1);
+ fOutputList->Add( fHOutCentV0CvsCentTRK);
+ fOutputList->Add( fHOutCentV0CvsCentCND);
+ fOutputList->Add( fHOutCentV0CvsCentCL1);
+ fOutputList->Add( fHOutMultV0AC);
fOutputList->Add( fHOutMultV0M);
+ fOutputList->Add( fHOutMultV0A);
+ fOutputList->Add( fHOutMultV0C);
fOutputList->Add( fHOutMultV0O);
fOutputList->Add( fHOutMultFMD);
fOutputList->Add( fHOutMultTRK);
fOutputList->Add( fHOutMultTKL);
fOutputList->Add( fHOutMultCL0);
fOutputList->Add( fHOutMultCL1);
+ fOutputList->Add( fHOutMultCND);
fOutputList->Add( fHOutMultV0MvsZDN);
fOutputList->Add( fHOutMultZEMvsZDN);
fOutputList->Add( fHOutMultV0MvsZDC);
fOutputList->Add( fHOutMultV0MvsV0O);
fOutputList->Add( fHOutMultV0OvsCL1);
fOutputList->Add( fHOutMultV0OvsTRK);
+ fOutputList->Add( fHOutMultCL1vsTKL);
fOutputList->Add( fHOutCentV0Mqual1 );
fOutputList->Add( fHOutCentTRKqual1 );
fOutputList->Add( fHOutCentCL1qual1 );
Int_t nClusters[6] = {0}; // no. clusters on 6 ITS layers
Int_t nChips[2]; // no. chips on 2 SPD layers
Float_t spdCorr =0; // corrected spd2 multiplicity
+ Int_t multCND = 0; // no. tracks (candle condition)
Float_t multV0A = 0; // multiplicity from V0 reco side A
Float_t multV0C = 0; // multiplicity from V0 reco side C
fIsSelected = ((esdV0->GetV0ADecision()==1) && (esdV0->GetV0CDecision()==1));
TString trigStr(esd->GetFiredTriggerClasses());
+ fMB=kFALSE;
fCVHN=kFALSE; fCVLN=kFALSE; fCCENT=kFALSE; fCSEMI=kFALSE;
fMSL=kFALSE; fMSH=kFALSE; fMUL=kFALSE; fMLL=kFALSE;
fEJE=kFALSE; fEGA=kFALSE; fPHS=kFALSE;
+ if ( (trigStr.Contains("-B-")) && (trigStr.Contains("CPBI")) && (fIsSelected))
+ fMB=kTRUE;
if ( (trigStr.Contains("-B-")) && (trigStr.Contains("CVHN")) && (fIsSelected))
fCVHN=kTRUE;
if ( (trigStr.Contains("-B-")) && (trigStr.Contains("CVLN")) && (fIsSelected))
// ***** CB info (tracklets, clusters, chips)
//nTracks = event->GetNumberOfTracks();
nTracks = fTrackCuts ? (Short_t)fTrackCuts->GetReferenceMultiplicity(esd,kTRUE):-1;
+ Short_t nTrTPCcandle = 0;
+ for (Int_t iTracks = 0; iTracks < esd->GetNumberOfTracks(); iTracks++) {
+ AliESDtrack* track = esd->GetTrack(iTracks);
+ if (!track) continue;
+ if (track->Pt() > 0.5 && TMath::Abs(track->Eta()) < 0.8) nTrTPCcandle++;
+ }
+ multCND = nTrTPCcandle;
const AliMultiplicity *mult = esd->GetMultiplicity();
-
nTracklets = mult->GetNumberOfTracklets();
for(Int_t ilay=0; ilay<6; ilay++){
// ***** Centrality Selection
if(fHtempV0M) fCentV0M = fHtempV0M->GetBinContent(fHtempV0M->FindBin((v0Corr)));
+ if(fHtempV0A) fCentV0A = fHtempV0A->GetBinContent(fHtempV0A->FindBin((multV0A)));
+ if(fHtempV0C) fCentV0C = fHtempV0C->GetBinContent(fHtempV0C->FindBin((multV0C)));
if(fHtempFMD) fCentFMD = fHtempFMD->GetBinContent(fHtempFMD->FindBin((multFMDA+multFMDC)));
if(fHtempTRK) fCentTRK = fHtempTRK->GetBinContent(fHtempTRK->FindBin(nTracks));
if(fHtempTKL) fCentTKL = fHtempTKL->GetBinContent(fHtempTKL->FindBin(nTracklets));
if(fHtempCL0) fCentCL0 = fHtempCL0->GetBinContent(fHtempCL0->FindBin(nClusters[0]));
if(fHtempCL1) fCentCL1 = fHtempCL1->GetBinContent(fHtempCL1->FindBin(spdCorr));
+ if(fHtempCND) fCentCND = fHtempCND->GetBinContent(fHtempCND->FindBin(multCND));
if(fHtempV0MvsFMD) fCentV0MvsFMD = fHtempV0MvsFMD->GetBinContent(fHtempV0MvsFMD->FindBin((multV0A+multV0C)));
if(fHtempTKLvsV0M) fCentTKLvsV0M = fHtempTKLvsV0M->GetBinContent(fHtempTKLvsV0M->FindBin(nTracklets));
if (esdCent) {
esdCent->SetQuality(fQuality);
esdCent->SetCentralityV0M(fCentV0M);
+ esdCent->SetCentralityV0A(fCentV0A);
+ esdCent->SetCentralityV0C(fCentV0C);
esdCent->SetCentralityFMD(fCentFMD);
esdCent->SetCentralityTRK(fCentTRK);
esdCent->SetCentralityTKL(fCentTKL);
esdCent->SetCentralityCL0(fCentCL0);
esdCent->SetCentralityCL1(fCentCL1);
+ esdCent->SetCentralityCND(fCentCND);
esdCent->SetCentralityV0MvsFMD(fCentV0MvsFMD);
esdCent->SetCentralityTKLvsV0M(fCentTKLvsV0M);
esdCent->SetCentralityZEMvsZDC(fCentZEMvsZDC);
// filling QA histograms
if (fFillHistos) {
+ if ((fMB) && (abs(zvtx)<10)) fHOutMultCL1vsTKL->Fill(spdCorr,nTracklets);
+
if (fCVHN) fHOutCentV0MCVHN->Fill(fCentV0M);
if (fCVLN) fHOutCentV0MCVLN->Fill(fCentV0M);
if (fCCENT) fHOutCentV0MCCENT->Fill(fCentV0M);
if (fQuality==0) {
fHOutCentV0M->Fill(fCentV0M);
+ fHOutCentV0A->Fill(fCentV0A);
+ fHOutCentV0C->Fill(fCentV0C);
if (fCVHNbit) fHOutCentV0MCVHNinMB->Fill(fCentV0M);
if (fCVLNbit) fHOutCentV0MCVLNinMB->Fill(fCentV0M);
if (fEJE) fHOutCentV0MEJEinMB->Fill(fCentV0M);
if (fEGA) fHOutCentV0MEGAinMB->Fill(fCentV0M);
if (fPHS) fHOutCentV0MPHSinMB->Fill(fCentV0M);
-
-
+
fHOutCentFMD->Fill(fCentFMD);
fHOutCentTRK->Fill(fCentTRK);
fHOutCentTKL->Fill(fCentTKL);
fHOutCentCL0->Fill(fCentCL0);
fHOutCentCL1->Fill(fCentCL1);
+ fHOutCentCND->Fill(fCentCND);
fHOutCentV0MvsFMD->Fill(fCentV0MvsFMD);
fHOutCentTKLvsV0M->Fill(fCentTKLvsV0M);
fHOutCentZEMvsZDC->Fill(fCentZEMvsZDC);
fHOutCentV0MvsCentTRK->Fill(fCentV0M,fCentTRK);
fHOutCentTRKvsCentCL1->Fill(fCentTRK,fCentCL1);
fHOutCentV0MvsCentZDC->Fill(fCentV0M,fCentZEMvsZDC);
+ fHOutCentV0AvsCentV0C->Fill(fCentV0A,fCentV0C);
+ fHOutCentV0AvsCentTRK->Fill(fCentV0A,fCentTRK);
+ fHOutCentV0AvsCentCND->Fill(fCentV0A,fCentCND);
+ fHOutCentV0AvsCentCL1->Fill(fCentV0A,fCentCL1);
+ fHOutCentV0CvsCentTRK->Fill(fCentV0C,fCentTRK);
+ fHOutCentV0CvsCentCND->Fill(fCentV0C,fCentCND);
+ fHOutCentV0CvsCentCL1->Fill(fCentV0C,fCentCL1);
+
+ fHOutMultV0AC->Fill(multV0A,multV0C);
fHOutMultV0M->Fill(multV0A+multV0C);
+ fHOutMultV0A->Fill(multV0A);
+ fHOutMultV0C->Fill(multV0C);
fHOutMultV0O->Fill(multV0AOnline+multV0COnline);
fHOutMultFMD->Fill(multFMDA+multFMDC);
fHOutMultTRK->Fill(nTracks);
fHOutMultTKL->Fill(nTracklets);
fHOutMultCL0->Fill(nClusters[0]);
fHOutMultCL1->Fill(spdCorr);
+ fHOutMultCND->Fill(multCND);
+
fHOutMultV0MvsZDN->Fill(v0Corr,(zncEnergy+znaEnergy));
fHOutMultZEMvsZDN->Fill((zem1Energy+zem2Energy),(zncEnergy+znaEnergy));
fHOutMultV0MvsZDC->Fill(v0Corr,(zncEnergy+znaEnergy+zpcEnergy+zpaEnergy));
// centrality histos
fHtempV0M = centOADB->V0hist();
+ fHtempV0A = centOADB->V0Ahist();
+ fHtempV0C = centOADB->V0Chist();
fHtempTRK = centOADB->TPChist();
fHtempCL1 = centOADB->SPDhist();
+ fHtempCND = centOADB->CNDhist();
fHtempZEMvsZDC = centOADB->ZEMvsZDChist();
TString path = gSystem->ExpandPathName(fileName.Data());
if (!fHtempV0M) AliWarning(Form("Calibration for V0M does not exist in %s", path.Data()));
+ if (!fHtempV0A) AliWarning(Form("Calibration for V0A does not exist in %s", path.Data()));
+ if (!fHtempV0C) AliWarning(Form("Calibration for V0C does not exist in %s", path.Data()));
if (!fHtempTRK) AliWarning(Form("Calibration for TRK does not exist in %s", path.Data()));
if (!fHtempCL1) AliWarning(Form("Calibration for CL1 does not exist in %s", path.Data()));
+ if (!fHtempCND) AliWarning(Form("Calibration for CND does not exist in %s", path.Data()));
if (!fHtempZEMvsZDC) AliWarning(Form("Calibration for ZEMvsZDC does not exist in %s", path.Data()));
// scale factors
Bool_t fEGA; //!
Bool_t fPHS; //!
+ Bool_t fMB; //! if the event is MB
Bool_t fCVHN; //! if the event is central trigger
Bool_t fCVLN; //! if the event is semicentral trigger
- Bool_t fCVHNbit; //! if the event is central trigger
- Bool_t fCVLNbit; //! if the event is semicentral trigger
+ Bool_t fCVHNbit; //! if the event is central trigger
+ Bool_t fCVLNbit; //! if the event is semicentral trigger
Bool_t fCCENT; //! if the event is central trigger
Bool_t fCSEMI; //! if the event is semicentral trigger
Bool_t fCCENTbit; //! if the event is central trigger
Bool_t fCSEMIbit; //! if the event is semicentral trigger
Float_t fCentV0M; // percentile centrality from V0
+ Float_t fCentV0A; // percentile centrality from V0A
+ Float_t fCentV0C; // percentile centrality from V0C
Float_t fCentFMD; // percentile centrality from FMD
Float_t fCentTRK; // percentile centrality from tracks
Float_t fCentTKL; // percentile centrality from tracklets
Float_t fCentCL0; // percentile centrality from clusters in layer 0
- Float_t fCentCL1; // percentile centrality from clusters in layer 0
+ Float_t fCentCL1; // percentile centrality from clusters in layer 1
+ Float_t fCentCND; // percentile centrality from candle
Float_t fCentV0MvsFMD; // percentile centrality from V0 vs FMD
Float_t fCentTKLvsV0M; // percentile centrality from tracklets vs V0
Float_t fCentZEMvsZDC; // percentile centrality from ZEM vs ZDC
TH1F *fHtempV0M; // histogram with centrality vs multiplicity using V0
+ TH1F *fHtempV0A; // histogram with centrality vs multiplicity using V0A
+ TH1F *fHtempV0C; // histogram with centrality vs multiplicity using V0C
TH1F *fHtempFMD; // histogram with centrality vs multiplicity using FMD
TH1F *fHtempTRK; // histogram with centrality vs multiplicity using tracks
TH1F *fHtempTKL; // histogram with centrality vs multiplicity using tracklets
TH1F *fHtempCL0; // histogram with centrality vs multiplicity using clusters in layer 0
- TH1F *fHtempCL1; // histogram with centrality vs multiplicity using clusters in layer 0
- TH1F *fHtempV0MvsFMD; // histogram with centrality vs multiplicity using V0 vs FMD
- TH1F *fHtempTKLvsV0M; // histogram with centrality vs multiplicity using tracklets vs V0
- TH2F *fHtempZEMvsZDC; // histogram with centrality vs multiplicity using ZEM vs ZDC
+ TH1F *fHtempCL1; // histogram with centrality vs multiplicity using clusters in layer 1
+ TH1F *fHtempCND; // histogram with centrality vs multiplicity using candle
+ TH1F *fHtempV0MvsFMD; // histogram with centrality vs multiplicity using V0 vs FMD
+ TH1F *fHtempTKLvsV0M; // histogram with centrality vs multiplicity using tracklets vs V0
+ TH2F *fHtempZEMvsZDC; // histogram with centrality vs multiplicity using ZEM vs ZDC
TList *fOutputList; // output list
- TH1F *fHOutCentV0M ; //control histogram for centrality
- TH1F *fHOutCentV0MCVHN; //control histogram for centrality
- TH1F *fHOutCentV0MCVLN; //control histogram for centrality
- TH1F *fHOutCentV0MCVHNinMB; //control histogram for centrality
- TH1F *fHOutCentV0MCVLNinMB; //control histogram for centrality
- TH1F *fHOutCentV0MCCENT; //control histogram for centrality
- TH1F *fHOutCentV0MCSEMI; //control histogram for centrality
- TH1F *fHOutCentV0MCCENTinMB; //control histogram for centrality
- TH1F *fHOutCentV0MCSEMIinMB; //control histogram for centrality
- TH1F *fHOutCentV0MMSL; //control histogram for centrality
- TH1F *fHOutCentV0MMSH; //control histogram for centrality
- TH1F *fHOutCentV0MMUL; //control histogram for centrality
- TH1F *fHOutCentV0MMLL; //control histogram for centrality
- TH1F *fHOutCentV0MEJE; //control histogram for centrality
- TH1F *fHOutCentV0MEGA; //control histogram for centrality
- TH1F *fHOutCentV0MPHS; //control histogram for centrality
+ TH1F *fHOutCentV0M ; //control histogram for centrality
+ TH1F *fHOutCentV0A ; //control histogram for centrality
+ TH1F *fHOutCentV0C ; //control histogram for centrality
+ TH1F *fHOutCentV0MCVHN; //control histogram for centrality
+ TH1F *fHOutCentV0MCVLN; //control histogram for centrality
+ TH1F *fHOutCentV0MCVHNinMB; //control histogram for centrality
+ TH1F *fHOutCentV0MCVLNinMB; //control histogram for centrality
+ TH1F *fHOutCentV0MCCENT; //control histogram for centrality
+ TH1F *fHOutCentV0MCSEMI; //control histogram for centrality
+ TH1F *fHOutCentV0MCCENTinMB; //control histogram for centrality
+ TH1F *fHOutCentV0MCSEMIinMB; //control histogram for centrality
+ TH1F *fHOutCentV0MMSL; //control histogram for centrality
+ TH1F *fHOutCentV0MMSH; //control histogram for centrality
+ TH1F *fHOutCentV0MMUL; //control histogram for centrality
+ TH1F *fHOutCentV0MMLL; //control histogram for centrality
+ TH1F *fHOutCentV0MEJE; //control histogram for centrality
+ TH1F *fHOutCentV0MEGA; //control histogram for centrality
+ TH1F *fHOutCentV0MPHS; //control histogram for centrality
TH1F *fHOutCentV0MMSLinMB; //control histogram for centrality
TH1F *fHOutCentV0MMSHinMB; //control histogram for centrality
TH1F *fHOutCentV0MMULinMB; //control histogram for centrality
TH1F *fHOutCentV0MEJEinMB; //control histogram for centrality
TH1F *fHOutCentV0MEGAinMB; //control histogram for centrality
TH1F *fHOutCentV0MPHSinMB; //control histogram for centrality
- TH1F *fHOutCentFMD ; //control histogram for centrality
- TH1F *fHOutCentTRK ; //control histogram for centrality
- TH1F *fHOutCentTKL ; //control histogram for centrality
- TH1F *fHOutCentCL0 ; //control histogram for centrality
- TH1F *fHOutCentCL1 ; //control histogram for centrality
- TH1F *fHOutCentV0MvsFMD; //control histogram for centrality
- TH1F *fHOutCentTKLvsV0M; //control histogram for centrality
- TH1F *fHOutCentZEMvsZDC; //control histogram for centrality
- TH2F *fHOutCentV0MvsCentCL1; //control histogram for centrality
- TH2F *fHOutCentV0MvsCentTRK; //control histogram for centrality
- TH2F *fHOutCentTRKvsCentCL1; //control histogram for centrality
- TH2F *fHOutCentV0MvsCentZDC; //control histogram for centrality
-
- TH1F *fHOutMultV0M ; //control histogram for multiplicity
- TH1F *fHOutMultV0O ; //control histogram for multiplicity
- TH1F *fHOutMultFMD ; //control histogram for multiplicity
- TH1F *fHOutMultTRK ; //control histogram for multiplicity
- TH1F *fHOutMultTKL ; //control histogram for multiplicity
- TH1F *fHOutMultCL0 ; //control histogram for multiplicity
- TH1F *fHOutMultCL1 ; //control histogram for multiplicity
-
- TH2F *fHOutMultV0MvsZDN; //control histogram for multiplicity
- TH2F *fHOutMultZEMvsZDN; //control histogram for multiplicity
- TH2F *fHOutMultV0MvsZDC; //control histogram for multiplicity
- TH2F *fHOutMultZEMvsZDC; //control histogram for multiplicity
- TH2F *fHOutMultZEMvsZDCw; //control histogram for multiplicity
- TH2F *fHOutMultV0MvsCL1; //control histogram for multiplicity
- TH2F *fHOutMultV0MvsTRK; //control histogram for multiplicity
- TH2F *fHOutMultTRKvsCL1; //control histogram for multiplicity
- TH2F *fHOutMultV0MvsV0O; //control histogram for multiplicity
- TH2F *fHOutMultV0OvsCL1; //control histogram for multiplicity
- TH2F *fHOutMultV0OvsTRK; //control histogram for multiplicity
-
- TH1F *fHOutCentV0Mqual1 ; //control histogram for centrality quality 1
- TH1F *fHOutCentTRKqual1 ; //control histogram for centrality quality 1
- TH1F *fHOutCentCL1qual1 ; //control histogram for centrality quality 1
- TH2F *fHOutMultV0MvsCL1qual1; //control histogram for multiplicity quality 1
- TH2F *fHOutMultV0MvsTRKqual1; //control histogram for multiplicity quality 1
- TH2F *fHOutMultTRKvsCL1qual1; //control histogram for multiplicity quality 1
-
- TH1F *fHOutCentV0Mqual2 ; //control histogram for centrality quality 2
- TH1F *fHOutCentTRKqual2 ; //control histogram for centrality quality 2
- TH1F *fHOutCentCL1qual2 ; //control histogram for centrality quality 2
- TH2F *fHOutMultV0MvsCL1qual2; //control histogram for multiplicity quality 2
- TH2F *fHOutMultV0MvsTRKqual2; //control histogram for multiplicity quality 2
- TH2F *fHOutMultTRKvsCL1qual2; //control histogram for multiplicity quality 2
-
- TH1F *fHOutQuality ; //control histogram for quality
- TH1F *fHOutVertex ; //control histogram for vertex
- TH1F *fHOutVertexT0 ; //control histogram for vertex
-
- ClassDef(AliCentralitySelectionTask, 18);
+ TH1F *fHOutCentFMD ; //control histogram for centrality
+ TH1F *fHOutCentTRK ; //control histogram for centrality
+ TH1F *fHOutCentTKL ; //control histogram for centrality
+ TH1F *fHOutCentCL0 ; //control histogram for centrality
+ TH1F *fHOutCentCL1 ; //control histogram for centrality
+ TH1F *fHOutCentCND ; //control histogram for centrality
+ TH1F *fHOutCentV0MvsFMD; //control histogram for centrality
+ TH1F *fHOutCentTKLvsV0M; //control histogram for centrality
+ TH1F *fHOutCentZEMvsZDC; //control histogram for centrality
+ TH2F *fHOutCentV0MvsCentCL1; //control histogram for centrality
+ TH2F *fHOutCentV0MvsCentTRK; //control histogram for centrality
+ TH2F *fHOutCentTRKvsCentCL1; //control histogram for centrality
+ TH2F *fHOutCentV0MvsCentZDC; //control histogram for centrality
+ TH2F *fHOutCentV0AvsCentV0C; //control histogram for centrality
+ TH2F *fHOutCentV0AvsCentTRK; //control histogram for centrality
+ TH2F *fHOutCentV0AvsCentCND; //control histogram for centrality
+ TH2F *fHOutCentV0AvsCentCL1; //control histogram for centrality
+ TH2F *fHOutCentV0CvsCentTRK; //control histogram for centrality
+ TH2F *fHOutCentV0CvsCentCND; //control histogram for centrality
+ TH2F *fHOutCentV0CvsCentCL1; //control histogram for centrality
+
+ TH2F *fHOutMultV0AC; //control histogram for multiplicity
+ TH1F *fHOutMultV0M ; //control histogram for multiplicity
+ TH1F *fHOutMultV0A ; //control histogram for multiplicity
+ TH1F *fHOutMultV0C ; //control histogram for multiplicity
+ TH1F *fHOutMultV0O ; //control histogram for multiplicity
+ TH1F *fHOutMultFMD ; //control histogram for multiplicity
+ TH1F *fHOutMultTRK ; //control histogram for multiplicity
+ TH1F *fHOutMultTKL ; //control histogram for multiplicity
+ TH1F *fHOutMultCL0 ; //control histogram for multiplicity
+ TH1F *fHOutMultCL1 ; //control histogram for multiplicity
+ TH1F *fHOutMultCND ; //control histogram for multiplicity
+
+ TH2F *fHOutMultV0MvsZDN; //control histogram for multiplicity
+ TH2F *fHOutMultZEMvsZDN; //control histogram for multiplicity
+ TH2F *fHOutMultV0MvsZDC; //control histogram for multiplicity
+ TH2F *fHOutMultZEMvsZDC; //control histogram for multiplicity
+ TH2F *fHOutMultZEMvsZDCw; //control histogram for multiplicity
+ TH2F *fHOutMultV0MvsCL1; //control histogram for multiplicity
+ TH2F *fHOutMultV0MvsTRK; //control histogram for multiplicity
+ TH2F *fHOutMultTRKvsCL1; //control histogram for multiplicity
+ TH2F *fHOutMultV0MvsV0O; //control histogram for multiplicity
+ TH2F *fHOutMultV0OvsCL1; //control histogram for multiplicity
+ TH2F *fHOutMultV0OvsTRK; //control histogram for multiplicity
+ TH2F *fHOutMultCL1vsTKL; //control histogram for multiplicity
+
+ TH1F *fHOutCentV0Mqual1 ; //control histogram for centrality quality 1
+ TH1F *fHOutCentTRKqual1 ; //control histogram for centrality quality 1
+ TH1F *fHOutCentCL1qual1 ; //control histogram for centrality quality 1
+ TH2F *fHOutMultV0MvsCL1qual1; //control histogram for multiplicity quality 1
+ TH2F *fHOutMultV0MvsTRKqual1; //control histogram for multiplicity quality 1
+ TH2F *fHOutMultTRKvsCL1qual1; //control histogram for multiplicity quality 1
+
+ TH1F *fHOutCentV0Mqual2 ; //control histogram for centrality quality 2
+ TH1F *fHOutCentTRKqual2 ; //control histogram for centrality quality 2
+ TH1F *fHOutCentCL1qual2 ; //control histogram for centrality quality 2
+ TH2F *fHOutMultV0MvsCL1qual2; //control histogram for multiplicity quality 2
+ TH2F *fHOutMultV0MvsTRKqual2; //control histogram for multiplicity quality 2
+ TH2F *fHOutMultTRKvsCL1qual2; //control histogram for multiplicity quality 2
+
+ TH1F *fHOutQuality ; //control histogram for quality
+ TH1F *fHOutVertex ; //control histogram for vertex SPD
+ TH1F *fHOutVertexT0 ; //control histogram for vertex T0
+
+ ClassDef(AliCentralitySelectionTask, 19);
};
#endif
AliCentrality::AliCentrality() : TNamed("Centrality", "Centrality"),
fQuality(999),
fCentralityV0M(0),
+ fCentralityV0A(0),
+ fCentralityV0C(0),
fCentralityFMD(0),
fCentralityTRK(0),
fCentralityTKL(0),
fCentralityCL0(0),
fCentralityCL1(0),
+ fCentralityCND(0),
fCentralityV0MvsFMD(0),
fCentralityTKLvsV0M(0),
fCentralityZEMvsZDC(0)
TNamed(cnt),
fQuality(cnt.fQuality),
fCentralityV0M(cnt.fCentralityV0M),
+ fCentralityV0A(cnt.fCentralityV0A),
+ fCentralityV0C(cnt.fCentralityV0C),
fCentralityFMD(cnt.fCentralityFMD),
fCentralityTRK(cnt.fCentralityTRK),
fCentralityTKL(cnt.fCentralityTKL),
fCentralityCL0(cnt.fCentralityCL0),
fCentralityCL1(cnt.fCentralityCL1),
+ fCentralityCND(cnt.fCentralityCND),
fCentralityV0MvsFMD(cnt.fCentralityV0MvsFMD),
fCentralityTKLvsV0M(cnt.fCentralityTKLvsV0M),
fCentralityZEMvsZDC(cnt.fCentralityZEMvsZDC)
TNamed::operator=(c);
fQuality = c.fQuality;
fCentralityV0M = c.fCentralityV0M;
+ fCentralityV0A = c.fCentralityV0A;
+ fCentralityV0C = c.fCentralityV0C;
fCentralityFMD = c.fCentralityFMD;
fCentralityTRK = c.fCentralityTRK;
fCentralityTKL = c.fCentralityTKL;
fCentralityCL0 = c.fCentralityCL0;
fCentralityCL1 = c.fCentralityCL1;
+ fCentralityCND = c.fCentralityCND;
fCentralityV0MvsFMD = c.fCentralityV0MvsFMD;
fCentralityTKLvsV0M = c.fCentralityTKLvsV0M;
fCentralityZEMvsZDC = c.fCentralityZEMvsZDC;
if (fQuality == 0) {
TString method = x;
if(method.CompareTo("V0M")==0) return fCentralityV0M;
+ if(method.CompareTo("V0A")==0) return fCentralityV0A;
+ if(method.CompareTo("V0C")==0) return fCentralityV0C;
if(method.CompareTo("FMD")==0) return fCentralityFMD;
if(method.CompareTo("TRK")==0) return fCentralityTRK;
if(method.CompareTo("TKL")==0) return fCentralityTKL;
if(method.CompareTo("CL0")==0) return fCentralityCL0;
if(method.CompareTo("CL1")==0) return fCentralityCL1;
+ if(method.CompareTo("CND")==0) return fCentralityCND;
if(method.CompareTo("V0MvsFMD")==0) return fCentralityV0MvsFMD;
if(method.CompareTo("TKLvsV0M")==0) return fCentralityTKLvsV0M;
if(method.CompareTo("ZEMvsZDC")==0) return fCentralityZEMvsZDC;
if (fQuality == 0) {
TString method = x;
if(method.CompareTo("V0M")==0) return (Int_t) (fCentralityV0M / 10.0);
+ if(method.CompareTo("V0A")==0) return (Int_t) (fCentralityV0A / 10.0);
+ if(method.CompareTo("V0C")==0) return (Int_t) (fCentralityV0C / 10.0);
if(method.CompareTo("FMD")==0) return (Int_t) (fCentralityFMD / 10.0);
if(method.CompareTo("TRK")==0) return (Int_t) (fCentralityTRK / 10.0);
if(method.CompareTo("TKL")==0) return (Int_t) (fCentralityTKL / 10.0);
if(method.CompareTo("CL0")==0) return (Int_t) (fCentralityCL0 / 10.0);
if(method.CompareTo("CL1")==0) return (Int_t) (fCentralityCL1 / 10.0);
+ if(method.CompareTo("CND")==0) return (Int_t) (fCentralityCND / 10.0);
if(method.CompareTo("V0MvsFMD")==0) return (Int_t) (fCentralityV0MvsFMD / 10.0);
if(method.CompareTo("TKLvsV0M")==0) return (Int_t) (fCentralityTKLvsV0M / 10.0);
if(method.CompareTo("ZEMvsZDC")==0) return (Int_t) (fCentralityZEMvsZDC / 10.0);
if (fQuality == 0) {
TString method = x;
if(method.CompareTo("V0M")==0) return (Int_t) (fCentralityV0M / 5.0);
+ if(method.CompareTo("V0A")==0) return (Int_t) (fCentralityV0A / 5.0);
+ if(method.CompareTo("V0C")==0) return (Int_t) (fCentralityV0C / 5.0);
if(method.CompareTo("FMD")==0) return (Int_t) (fCentralityFMD / 5.0);
if(method.CompareTo("TRK")==0) return (Int_t) (fCentralityTRK / 5.0);
if(method.CompareTo("TKL")==0) return (Int_t) (fCentralityTKL / 5.0);
if(method.CompareTo("CL0")==0) return (Int_t) (fCentralityCL0 / 5.0);
if(method.CompareTo("CL1")==0) return (Int_t) (fCentralityCL1 / 5.0);
+ if(method.CompareTo("CND")==0) return (Int_t) (fCentralityCND / 5.0);
if(method.CompareTo("V0MvsFMD")==0) return (Int_t) (fCentralityV0MvsFMD / 5.0);
if(method.CompareTo("TKLvsV0M")==0) return (Int_t) (fCentralityTKLvsV0M / 5.0);
if(method.CompareTo("ZEMvsZDC")==0) return (Int_t) (fCentralityZEMvsZDC / 5.0);
if (fQuality == 0) {
TString method = x;
if ((method.CompareTo("V0M")==0) && (fCentralityV0M >=a && fCentralityV0M < b)) return kTRUE;
+ if ((method.CompareTo("V0A")==0) && (fCentralityV0A >=a && fCentralityV0A < b)) return kTRUE;
+ if ((method.CompareTo("V0C")==0) && (fCentralityV0C >=a && fCentralityV0C < b)) return kTRUE;
if ((method.CompareTo("FMD")==0) && (fCentralityFMD >=a && fCentralityFMD < b)) return kTRUE;
if ((method.CompareTo("TRK")==0) && (fCentralityTRK >=a && fCentralityTRK < b)) return kTRUE;
if ((method.CompareTo("TKL")==0) && (fCentralityTKL >=a && fCentralityTKL < b)) return kTRUE;
if ((method.CompareTo("CL0")==0) && (fCentralityCL0 >=a && fCentralityCL0 < b)) return kTRUE;
if ((method.CompareTo("CL1")==0) && (fCentralityCL1 >=a && fCentralityCL1 < b)) return kTRUE;
+ if ((method.CompareTo("CND")==0) && (fCentralityCND >=a && fCentralityCND < b)) return kTRUE;
if ((method.CompareTo("V0MvsFMD")==0) && (fCentralityV0MvsFMD >=a && fCentralityV0MvsFMD < b)) return kTRUE;
if ((method.CompareTo("TKLvsV0M")==0) && (fCentralityTKLvsV0M >=a && fCentralityTKLvsV0M < b)) return kTRUE;
if ((method.CompareTo("ZEMvsZDC")==0) && (fCentralityZEMvsZDC >=a && fCentralityZEMvsZDC < b)) return kTRUE;
// Return the centrality percentile
TString method = x;
if(method.CompareTo("V0M")==0) return fCentralityV0M;
+ if(method.CompareTo("V0A")==0) return fCentralityV0A;
+ if(method.CompareTo("V0C")==0) return fCentralityV0C;
if(method.CompareTo("FMD")==0) return fCentralityFMD;
if(method.CompareTo("TRK")==0) return fCentralityTRK;
if(method.CompareTo("TKL")==0) return fCentralityTKL;
if(method.CompareTo("CL0")==0) return fCentralityCL0;
if(method.CompareTo("CL1")==0) return fCentralityCL1;
+ if(method.CompareTo("CND")==0) return fCentralityCND;
if(method.CompareTo("V0MvsFMD")==0) return fCentralityV0MvsFMD;
if(method.CompareTo("TKLvsV0M")==0) return fCentralityTKLvsV0M;
if(method.CompareTo("ZEMvsZDC")==0) return fCentralityZEMvsZDC;
// Return the centrality class
TString method = x;
if(method.CompareTo("V0M")==0) return (Int_t) (fCentralityV0M / 10.0);
+ if(method.CompareTo("V0A")==0) return (Int_t) (fCentralityV0A / 10.0);
+ if(method.CompareTo("V0C")==0) return (Int_t) (fCentralityV0C / 10.0);
if(method.CompareTo("FMD")==0) return (Int_t) (fCentralityFMD / 10.0);
if(method.CompareTo("TRK")==0) return (Int_t) (fCentralityTRK / 10.0);
if(method.CompareTo("TKL")==0) return (Int_t) (fCentralityTKL / 10.0);
if(method.CompareTo("CL0")==0) return (Int_t) (fCentralityCL0 / 10.0);
if(method.CompareTo("CL1")==0) return (Int_t) (fCentralityCL1 / 10.0);
+ if(method.CompareTo("CND")==0) return (Int_t) (fCentralityCND / 10.0);
if(method.CompareTo("V0MvsFMD")==0) return (Int_t) (fCentralityV0MvsFMD / 10.0);
if(method.CompareTo("TKLvsV0M")==0) return (Int_t) (fCentralityTKLvsV0M / 10.0);
if(method.CompareTo("ZEMvsZDC")==0) return (Int_t) (fCentralityZEMvsZDC / 10.0);
// Return the centrality class
TString method = x;
if(method.CompareTo("V0M")==0) return (Int_t) (fCentralityV0M / 5.0);
+ if(method.CompareTo("V0A")==0) return (Int_t) (fCentralityV0A / 5.0);
+ if(method.CompareTo("V0C")==0) return (Int_t) (fCentralityV0C / 5.0);
if(method.CompareTo("FMD")==0) return (Int_t) (fCentralityFMD / 5.0);
if(method.CompareTo("TRK")==0) return (Int_t) (fCentralityTRK / 5.0);
if(method.CompareTo("TKL")==0) return (Int_t) (fCentralityTKL / 5.0);
if(method.CompareTo("CL0")==0) return (Int_t) (fCentralityCL0 / 5.0);
if(method.CompareTo("CL1")==0) return (Int_t) (fCentralityCL1 / 5.0);
+ if(method.CompareTo("CND")==0) return (Int_t) (fCentralityCND / 5.0);
if(method.CompareTo("V0MvsFMD")==0) return (Int_t) (fCentralityV0MvsFMD / 5.0);
if(method.CompareTo("TKLvsV0M")==0) return (Int_t) (fCentralityTKLvsV0M / 5.0);
if(method.CompareTo("ZEMvsZDC")==0) return (Int_t) (fCentralityZEMvsZDC / 5.0);
// True if event inside given centrality class
TString method = x;
if ((method.CompareTo("V0M")==0) && (fCentralityV0M >=a && fCentralityV0M < b)) return kTRUE;
+ if ((method.CompareTo("V0A")==0) && (fCentralityV0A >=a && fCentralityV0A < b)) return kTRUE;
+ if ((method.CompareTo("V0C")==0) && (fCentralityV0C >=a && fCentralityV0C < b)) return kTRUE;
if ((method.CompareTo("FMD")==0) && (fCentralityFMD >=a && fCentralityFMD < b)) return kTRUE;
if ((method.CompareTo("TRK")==0) && (fCentralityTRK >=a && fCentralityTRK < b)) return kTRUE;
if ((method.CompareTo("TKL")==0) && (fCentralityTKL >=a && fCentralityTKL < b)) return kTRUE;
if ((method.CompareTo("CL0")==0) && (fCentralityCL0 >=a && fCentralityCL0 < b)) return kTRUE;
if ((method.CompareTo("CL1")==0) && (fCentralityCL1 >=a && fCentralityCL1 < b)) return kTRUE;
+ if ((method.CompareTo("CND")==0) && (fCentralityCND >=a && fCentralityCND < b)) return kTRUE;
if ((method.CompareTo("V0MvsFMD")==0) && (fCentralityV0MvsFMD >=a && fCentralityV0MvsFMD < b)) return kTRUE;
if ((method.CompareTo("TKLvsV0M")==0) && (fCentralityTKLvsV0M >=a && fCentralityTKLvsV0M < b)) return kTRUE;
if ((method.CompareTo("ZEMvsZDC")==0) && (fCentralityZEMvsZDC >=a && fCentralityZEMvsZDC < b)) return kTRUE;
fQuality = 999;
fCentralityV0M = 0;
+ fCentralityV0A = 0;
+ fCentralityV0C = 0;
fCentralityFMD = 0;
fCentralityTRK = 0;
fCentralityTKL = 0;
fCentralityCL0 = 0;
fCentralityCL1 = 0;
+ fCentralityCND = 0;
fCentralityV0MvsFMD = 0;
fCentralityTKLvsV0M = 0;
fCentralityZEMvsZDC = 0;