fCentTRK(0),
fCentTKL(0),
fCentCL0(0),
+ fCentCL1(0),
fCentV0MvsFMD(0),
fCentTKLvsV0M(0),
fCentZEMvsZDC(0),
fHtempTRK(0),
fHtempTKL(0),
fHtempCL0(0),
+ fHtempCL1(0),
fHtempV0MvsFMD(0),
fHtempTKLvsV0M(0),
fHtempZEMvsZDC(0)
fCentTRK(0),
fCentTKL(0),
fCentCL0(0),
+ fCentCL1(0),
fCentV0MvsFMD(0),
fCentTKLvsV0M(0),
fCentZEMvsZDC(0),
fHtempTRK(0),
fHtempTKL(0),
fHtempCL0(0),
+ fHtempCL1(0),
fHtempV0MvsFMD(0),
fHtempTKLvsV0M(0),
fHtempZEMvsZDC(0)
fCentTRK(ana.fCentTRK),
fCentTKL(ana.fCentTKL),
fCentCL0(ana.fCentCL0),
+ fCentCL1(ana.fCentCL1),
fCentV0MvsFMD(ana.fCentV0MvsFMD),
fCentTKLvsV0M(ana.fCentTKLvsV0M),
fCentZEMvsZDC(ana.fCentZEMvsZDC),
fHtempTRK(ana.fHtempTRK),
fHtempTKL(ana.fHtempTKL),
fHtempCL0(ana.fHtempCL0),
+ fHtempCL1(ana.fHtempCL1),
fHtempV0MvsFMD(ana.fHtempV0MvsFMD),
fHtempTKLvsV0M(ana.fHtempTKLvsV0M),
fHtempZEMvsZDC(ana.fHtempZEMvsZDC)
// Execute analysis for current event:
if(fDebug>1) printf(" **** AliCentralitySelectionTask::UserExec() \n");
- Float_t zncEnergy; // ZNC Energy
- Float_t zpcEnergy; // ZPC Energy
- Float_t znaEnergy; // ZNA Energy
- Float_t zpaEnergy; // ZPA Energy
+ Float_t zncEnergy; // ZNC Energy
+ Float_t zpcEnergy; // ZPC Energy
+ Float_t znaEnergy; // ZNA Energy
+ Float_t zpaEnergy; // ZPA Energy
Float_t zem1Energy = 0.; // ZEM1 Energy
Float_t zem2Energy = 0.; // ZEM2 Energy
- Int_t nTracks = 0; // no. tracks
+ Int_t nTracks = 0; // no. tracks
Int_t nTracklets = 0; // no. tracklets
Int_t nClusters[6]; // no. clusters on 6 ITS layers
Int_t nChips[2]; // no. chips on 2 SPD layers
fCentTRK = fHtempTRK->GetBinContent(fHtempTRK->FindBin(nTracks));
fCentTKL = fHtempTKL->GetBinContent(fHtempTKL->FindBin(nTracklets));
fCentCL0 = fHtempCL0->GetBinContent(fHtempCL0->FindBin(nClusters[0]));
+ fCentCL1 = fHtempCL1->GetBinContent(fHtempCL1->FindBin(nClusters[1]));
fCentV0MvsFMD = fHtempV0MvsFMD->GetBinContent(fHtempV0MvsFMD->FindBin((multV0A+multV0C)));
fCentTKLvsV0M = fHtempTKLvsV0M->GetBinContent(fHtempTKLvsV0M->FindBin(nTracklets));
esdCent->SetCentralityTRK(fCentTRK);
esdCent->SetCentralityTKL(fCentTKL);
esdCent->SetCentralityCL0(fCentCL0);
+ esdCent->SetCentralityCL1(fCentCL1);
esdCent->SetCentralityV0MvsFMD(fCentV0MvsFMD);
esdCent->SetCentralityTKLvsV0M(fCentTKLvsV0M);
esdCent->SetCentralityZEMvsZDC(fCentZEMvsZDC);
fHtempTRK = (TH1D*) (fFile->Get("hNtracks_percentile"));
fHtempTKL = (TH1D*) (fFile->Get("hNtracklets_percentile"));
fHtempCL0 = (TH1D*) (fFile->Get("hNclusters0_percentile"));
+ fHtempCL1 = (TH1D*) (fFile->Get("hNclusters1_percentile"));
}
void AliCentralitySelectionTask::ReadCentralityHistos2()
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 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
TH1D *fHtempTRK; // histogram with centrality vs multiplicity using tracks
TH1D *fHtempTKL; // histogram with centrality vs multiplicity using tracklets
TH1D *fHtempCL0; // histogram with centrality vs multiplicity using clusters in layer 0
+ TH1D *fHtempCL1; // histogram with centrality vs multiplicity using clusters in layer 0
TH1D *fHtempV0MvsFMD; // histogram with centrality vs multiplicity using V0 vs FMD
TH1D *fHtempTKLvsV0M; // histogram with centrality vs multiplicity using tracklets vs V0
TH1D *fHtempZEMvsZDC; // histogram with centrality vs multiplicity using ZEM vs ZDC
fCentralityTRK(0),
fCentralityTKL(0),
fCentralityCL0(0),
+ fCentralityCL1(0),
fCentralityV0MvsFMD(0),
fCentralityTKLvsV0M(0),
fCentralityZEMvsZDC(0)
fCentralityTRK(cnt.fCentralityTRK),
fCentralityTKL(cnt.fCentralityTKL),
fCentralityCL0(cnt.fCentralityCL0),
+ fCentralityCL1(cnt.fCentralityCL1),
fCentralityV0MvsFMD(cnt.fCentralityV0MvsFMD),
fCentralityTKLvsV0M(cnt.fCentralityTKLvsV0M),
fCentralityZEMvsZDC(cnt.fCentralityZEMvsZDC)
fCentralityTRK = c.fCentralityTRK;
fCentralityTKL = c.fCentralityTKL;
fCentralityCL0 = c.fCentralityCL0;
+ fCentralityCL1 = c.fCentralityCL1;
fCentralityV0MvsFMD = c.fCentralityV0MvsFMD;
fCentralityTKLvsV0M = c.fCentralityTKLvsV0M;
fCentralityZEMvsZDC = c.fCentralityZEMvsZDC;
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("V0MvsFMD")==0) return fCentralityV0MvsFMD;
if(method.CompareTo("TKLvsV0M")==0) return fCentralityTKLvsV0M;
if(method.CompareTo("ZENvsZDC")==0) return fCentralityZEMvsZDC;
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("V0MvsFMD")==0) return (Int_t) (fCentralityV0MvsFMD / 10.0);
if(method.CompareTo("TKLvsV0M")==0) return (Int_t) (fCentralityTKLvsV0M / 10.0);
if(method.CompareTo("ZENvsZDC")==0) return (Int_t) (fCentralityZEMvsZDC / 10.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("V0MvsFMD")==0) return (Int_t) (fCentralityV0MvsFMD / 5.0);
if(method.CompareTo("TKLvsV0M")==0) return (Int_t) (fCentralityTKLvsV0M / 5.0);
if(method.CompareTo("ZENvsZDC")==0) return (Int_t) (fCentralityZEMvsZDC / 5.0);
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("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;
void SetCentralityTRK(Float_t cent) {fCentralityTRK = cent;}
void SetCentralityTKL(Float_t cent) {fCentralityTKL = cent;}
void SetCentralityCL0(Float_t cent) {fCentralityCL0 = cent;}
+ void SetCentralityCL1(Float_t cent) {fCentralityCL1 = cent;}
void SetCentralityV0MvsFMD(Float_t cent) {fCentralityV0MvsFMD = cent;}
void SetCentralityTKLvsV0M(Float_t cent) {fCentralityTKLvsV0M = cent;}
void SetCentralityZEMvsZDC(Float_t cent) {fCentralityZEMvsZDC = cent;}
Float_t fCentralityTRK; // Centrality from tracks
Float_t fCentralityTKL; // Centrality from tracklets
Float_t fCentralityCL0; // Centrality from Clusters in layer 0
+ Float_t fCentralityCL1; // Centrality from Clusters in layer 0
Float_t fCentralityV0MvsFMD; // Centrality from V0 vs FMD
Float_t fCentralityTKLvsV0M; // Centrality from tracklets vs V0
Float_t fCentralityZEMvsZDC; // Centrality from ZEM vs ZDC