+
+ const Int_t nRefMultBins = 100;
+ const Double_t refMultUp = 100.;
+ const Double_t refMultDown = 0.;
+
+ const Int_t nJetPtBins = 20;
+ const Double_t jetPtUp = 100.;
+ const Double_t jetPtDown = 0.;
+
+ const Int_t nCentBins = 12;
+ const Double_t binsCentpp[nCentBins+1] = { 0, 0.01, 0.1, 1, 5, 10, 15, 20, 30, 40, 50, 70, 100};
+
+ fhJetPtRefMultEta5 = new TH2F("fhJetPtRefMultEta5",
+ "Correlation between jet energy and event multiplicity (|#eta| < 0.5);Ref. mult. (|#eta| < 0.5);#it{p}_{T, jet}^{ch} (GeV/#it{c})",
+ nRefMultBins, refMultDown, refMultUp, nJetPtBins, jetPtDown, jetPtUp);
+ fhJetPtRefMultEta8 = new TH2F("fhJetPtRefMultEta8",
+ "Correlation between jet energy and event multiplicity (|#eta| < 0.8);Ref. mult. (|#eta| < 0.8);#it{p}_{T, jet}^{ch} (GeV/#it{c})",
+ nRefMultBins, refMultDown, refMultUp, nJetPtBins, jetPtDown, jetPtUp);
+ fhJetPtMultPercent = new TH2F("fhJetPtMultPercent",
+ "Correlation between jet energy and event multiplicity percentile (V0M);Multiplicity Percentile (V0M);#it{p}_{T, jet}^{ch} (GeV/#it{c})",
+ nCentBins, binsCentpp, nJetPtBins, jetPtDown, jetPtUp);
+ fCommonHistList->Add(fhJetPtRefMultEta5);
+ fCommonHistList->Add(fhJetPtRefMultEta8);
+ fCommonHistList->Add(fhJetPtMultPercent);
+
+ if (fUseJetPIDtask) {
+ const Int_t nPtBins = 68;
+ Double_t binsPt[nPtBins+1] = {0. , 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45,
+ 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95,
+ 1.0, 1.1 , 1.2, 1.3 , 1.4, 1.5 , 1.6, 1.7 , 1.8, 1.9 ,
+ 2.0, 2.2 , 2.4, 2.6 , 2.8, 3.0 , 3.2, 3.4 , 3.6, 3.8 ,
+ 4.0, 4.5 , 5.0, 5.5 , 6.0, 6.5 , 7.0, 8.0 , 9.0, 10.0,
+ 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 18.0, 20.0, 22.0, 24.0,
+ 26.0, 28.0, 30.0, 32.0, 34.0, 36.0, 40.0, 45.0, 50.0 };
+
+ const Int_t DCAbins = 320;
+ const Double_t DCA_Z_max = 3.5;
+ const Double_t DCA_XY_max = 2.5;
+
+ fhDCA_XY = new TH2F("fhDCA_XY", "All rec. tracks;#it{p}_{T} (GeV/#it{c});DCA_{XY}", nPtBins, binsPt, DCAbins, -DCA_XY_max, DCA_XY_max);
+ fhDCA_Z = new TH2F("fhDCA_Z", "All rec. tracks;#it{p}_{T} (GeV/#it{c});DCA_{Z}", nPtBins, binsPt, DCAbins, -DCA_Z_max, DCA_Z_max);
+ fCommonHistList->Add(fhDCA_XY);
+ fCommonHistList->Add(fhDCA_Z);
+
+ for (Int_t i = 0; i < AliPID::kSPECIES; i++) {
+ fhDCA_XY_prim_MCID[i] = new TH2F(Form("fhDCA_XY_prim_MCID_%s", AliPID::ParticleShortName(i)),
+ Form("Rec. %s (prim.);#it{p}_{T} (GeV/#it{c});DCA_{XY}", AliPID::ParticleLatexName(i)),
+ nPtBins, binsPt, DCAbins, -DCA_XY_max, DCA_XY_max);
+ fhDCA_Z_prim_MCID[i] = new TH2F(Form("fhDCA_Z_prim_MCID_%s", AliPID::ParticleShortName(i)),
+ Form("Rec. %s (prim.);#it{p}_{T} (GeV/#it{c});DCA_{Z}", AliPID::ParticleLatexName(i)),
+ nPtBins, binsPt, DCAbins, -DCA_Z_max, DCA_Z_max);
+ fCommonHistList->Add(fhDCA_XY_prim_MCID[i]);
+ fCommonHistList->Add(fhDCA_Z_prim_MCID[i]);
+
+ fhDCA_XY_sec_MCID[i] = new TH2F(Form("fhDCA_XY_sec_MCID_%s", AliPID::ParticleShortName(i)),
+ Form("Rec. %s (sec.);#it{p}_{T} (GeV/#it{c});DCA_{XY}", AliPID::ParticleLatexName(i)),
+ nPtBins, binsPt, DCAbins, -DCA_XY_max, DCA_XY_max);
+ fhDCA_Z_sec_MCID[i] = new TH2F(Form("fhDCA_Z_sec_MCID_%s", AliPID::ParticleShortName(i)),
+ Form("Rec. %s (sec.);#it{p}_{T} (GeV/#it{c});DCA_{Z}", AliPID::ParticleLatexName(i)),
+ nPtBins, binsPt, DCAbins, -DCA_Z_max, DCA_Z_max);
+ fCommonHistList->Add(fhDCA_XY_sec_MCID[i]);
+ fCommonHistList->Add(fhDCA_Z_sec_MCID[i]);
+ }
+ }
+
+
+ // =========== Switch on Sumw2 for all histos ===========
+ for (Int_t i=0; i<fCommonHistList->GetEntries(); ++i){
+ TH1 *h1 = dynamic_cast<TH1*>(fCommonHistList->At(i));
+ if (h1) h1->Sumw2();
+ else{
+ THnSparse *hnSparse = dynamic_cast<THnSparse*>(fCommonHistList->At(i));
+ if(hnSparse) hnSparse->Sumw2();
+ }
+ }
+
+ TH1::AddDirectory(oldStatus);