Skip to content

Commit 8ab692a

Browse files
Preet-BhanjanPreet Patialibuild
authored
[PWGCF] Addition of Nch analysis (#13047)
Co-authored-by: Preet Pati <preet@preet-6.local> Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent bdfab25 commit 8ab692a

File tree

1 file changed

+25
-34
lines changed

1 file changed

+25
-34
lines changed

‎PWGCF/Flow/Tasks/flowPbpbPikp.cxx‎

Lines changed: 25 additions & 34 deletions
Original file line numberDiff line numberDiff line change
@@ -107,6 +107,7 @@ struct FlowPbpbPikp {
107107
O2_DEFINE_CONFIGURABLE(cfgV0AT0Acut, int, 5, "V0AT0A cut")
108108
O2_DEFINE_CONFIGURABLE(cfgUseAsymmetricPID, bool, false, "Use asymmetric PID cuts")
109109
O2_DEFINE_CONFIGURABLE(cfgUseItsPID, bool, true, "Use ITS PID for particle identification")
110+
O2_DEFINE_CONFIGURABLE(cfgUseNch, bool, true, "Use Nch for multiplicity selection instead of centrality")
110111

111112
Configurable<std::vector<double>> cfgTrackDensityP0{"cfgTrackDensityP0", std::vector<double>{0.7217476707, 0.7384792571, 0.7542625668, 0.7640680200, 0.7701951667, 0.7755299053, 0.7805901710, 0.7849446786, 0.7957356586, 0.8113039262, 0.8211968966, 0.8280558878, 0.8329342135}, "parameter 0 for track density efficiency correction"};
112113
Configurable<std::vector<double>> cfgTrackDensityP1{"cfgTrackDensityP1", std::vector<double>{-2.169488e-05, -2.191913e-05, -2.295484e-05, -2.556538e-05, -2.754463e-05, -2.816832e-05, -2.846502e-05, -2.843857e-05, -2.705974e-05, -2.477018e-05, -2.321730e-05, -2.203315e-05, -2.109474e-05}, "parameter 1 for track density efficiency correction"};
@@ -129,6 +130,7 @@ struct FlowPbpbPikp {
129130
ConfigurableAxis axisParticles{"axisParticles", {3, 0, 3}, "axis for different hadrons"};
130131
ConfigurableAxis axisTPCsignal{"axisTPCsignal", {10000, 0, 1000}, "axis for TPC signal"};
131132
ConfigurableAxis axisTOFbeta{"axisTOFbeta", {200, 0, 2}, "axis for TOF beta"};
133+
ConfigurableAxis axisNch{"axisNch", {200, 2000, 4000}, "N_{ch}"};
132134

133135
std::vector<double> tofNsigmaCut;
134136
std::vector<double> itsNsigmaCut;
@@ -232,6 +234,8 @@ struct FlowPbpbPikp {
232234
regions.Print();
233235
configs.Print();
234236

237+
const AxisSpec axisCentForQA{100, 0, 100, "centrality (%)"};
238+
235239
histos.add("hVtxZ", "", {HistType::kTH1D, {axisVertex}});
236240
histos.add("hMult", "", {HistType::kTH1D, {{3000, 0.5, 3000.5}}});
237241
histos.add("hCent", "", {HistType::kTH1D, {{90, 0, 90}}});
@@ -240,25 +244,13 @@ struct FlowPbpbPikp {
240244
histos.add("hEta", "", {HistType::kTH1D, {axisEta}});
241245
histos.add("hPt", "", {HistType::kTH1D, {axisPt}});
242246
histos.add("c22_full_ch", "", {HistType::kTProfile, {axisMultiplicity}});
243-
histos.add("c22_full_pi", "", {HistType::kTProfile, {axisMultiplicity}});
244-
histos.add("c22_full_ka", "", {HistType::kTProfile, {axisMultiplicity}});
245-
histos.add("c22_full_pr", "", {HistType::kTProfile, {axisMultiplicity}});
246-
histos.add("c22_gap08F_ch", "", {HistType::kTProfile, {axisMultiplicity}});
247-
histos.add("c22_gap08F_pi", "", {HistType::kTProfile, {axisMultiplicity}});
248-
histos.add("c22_gap08F_ka", "", {HistType::kTProfile, {axisMultiplicity}});
249-
histos.add("c22_gap08F_pr", "", {HistType::kTProfile, {axisMultiplicity}});
250-
histos.add("c22_gap08B_ch", "", {HistType::kTProfile, {axisMultiplicity}});
251-
histos.add("c22_gap08B_pi", "", {HistType::kTProfile, {axisMultiplicity}});
252-
histos.add("c22_gap08B_ka", "", {HistType::kTProfile, {axisMultiplicity}});
253-
histos.add("c22_gap08B_pr", "", {HistType::kTProfile, {axisMultiplicity}});
254-
histos.add("c24_full_ch", "", {HistType::kTProfile, {axisMultiplicity}});
255-
histos.add("c24_full_pi", "", {HistType::kTProfile, {axisMultiplicity}});
256-
histos.add("c24_full_ka", "", {HistType::kTProfile, {axisMultiplicity}});
257-
histos.add("c24_full_pr", "", {HistType::kTProfile, {axisMultiplicity}});
247+
histos.add("c22_full_ch_Nch", "", {HistType::kTProfile, {axisNch}});
258248

259249
histos.add("TpcdEdx", "", {HistType::kTH2D, {axisPt, axisTPCsignal}});
260250
histos.add("TofBeta", "", {HistType::kTH2D, {axisPt, axisTOFbeta}});
261251

252+
histos.add("globalTracks_centT0C", "after cut;Centrality T0C;mulplicity global tracks", {HistType::kTH2D, {axisCentForQA, axisNch}});
253+
262254
histos.add("TofTpcNsigma_before", "", {HistType::kTHnSparseD, {{axisParticles, axisNsigmaTPC, axisNsigmaTOF, axisPt}}});
263255
if (!cfgUseItsPID)
264256
histos.add("TofTpcNsigma_after", "", {HistType::kTHnSparseD, {{axisParticles, axisNsigmaTPC, axisNsigmaTOF, axisPt}}});
@@ -360,7 +352,13 @@ struct FlowPbpbPikp {
360352

361353
fFC->SetName("FlowContainer");
362354
fFC->SetXAxis(fPtAxis);
363-
fFC->Initialize(oba, axisMultiplicity, cfgNbootstrap);
355+
356+
if (!cfgUseNch) {
357+
fFC->Initialize(oba, axisMultiplicity, cfgNbootstrap);
358+
} else {
359+
fFC->Initialize(oba, axisNch, cfgNbootstrap);
360+
}
361+
364362
delete oba;
365363

366364
if (eventCuts[kUseMultCorrCut]) {
@@ -898,6 +896,7 @@ struct FlowPbpbPikp {
898896
histos.fill(HIST("hVtxZ"), vtxz);
899897
histos.fill(HIST("hMult"), nTot);
900898
histos.fill(HIST("hCent"), cent);
899+
histos.fill(HIST("globalTracks_centT0C"), cent, nTot);
901900
fGFW->Clear();
902901

903902
float weff = 1;
@@ -1008,24 +1007,16 @@ struct FlowPbpbPikp {
10081007

10091008
// Filling cumulants with ROOT TProfile
10101009
fillProfile(corrconfigs.at(0), HIST("c22_full_ch"), cent);
1011-
fillProfile(corrconfigs.at(1), HIST("c22_full_pi"), cent);
1012-
fillProfile(corrconfigs.at(2), HIST("c22_full_ka"), cent);
1013-
fillProfile(corrconfigs.at(3), HIST("c22_full_pr"), cent);
1014-
fillProfile(corrconfigs.at(4), HIST("c22_gap08F_ch"), cent);
1015-
fillProfile(corrconfigs.at(5), HIST("c22_gap08F_pi"), cent);
1016-
fillProfile(corrconfigs.at(6), HIST("c22_gap08F_ka"), cent);
1017-
fillProfile(corrconfigs.at(7), HIST("c22_gap08F_pr"), cent);
1018-
fillProfile(corrconfigs.at(8), HIST("c22_gap08B_ch"), cent);
1019-
fillProfile(corrconfigs.at(9), HIST("c22_gap08B_pi"), cent);
1020-
fillProfile(corrconfigs.at(10), HIST("c22_gap08B_ka"), cent);
1021-
fillProfile(corrconfigs.at(11), HIST("c22_gap08B_pr"), cent);
1022-
fillProfile(corrconfigs.at(12), HIST("c24_full_ch"), cent);
1023-
fillProfile(corrconfigs.at(13), HIST("c24_full_pi"), cent);
1024-
fillProfile(corrconfigs.at(14), HIST("c24_full_ka"), cent);
1025-
fillProfile(corrconfigs.at(15), HIST("c24_full_pr"), cent);
1026-
1027-
for (uint l_ind = 0; l_ind < corrconfigs.size(); l_ind++) {
1028-
fillFC(corrconfigs.at(l_ind), cent, lRandom);
1010+
fillProfile(corrconfigs.at(0), HIST("c22_full_ch_Nch"), nTot);
1011+
1012+
if (!cfgUseNch) {
1013+
for (uint l_ind = 0; l_ind < corrconfigs.size(); l_ind++) {
1014+
fillFC(corrconfigs.at(l_ind), cent, lRandom);
1015+
}
1016+
} else {
1017+
for (uint l_ind = 0; l_ind < corrconfigs.size(); l_ind++) {
1018+
fillFC(corrconfigs.at(l_ind), nTot, lRandom);
1019+
}
10291020
}
10301021

10311022
} // end of process

0 commit comments

Comments
 (0)