Skip to content

Commit f7ef450

Browse files
author
Preet Pati
committed
Addition of Nch analysis
1 parent 1016b75 commit f7ef450

File tree

1 file changed

+26
-35
lines changed

1 file changed

+26
-35
lines changed

‎PWGCF/Flow/Tasks/flowPbpbPikp.cxx‎

Lines changed: 26 additions & 35 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,14 +130,15 @@ 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;
135137
std::vector<double> tpcNsigmaCut;
136138
std::vector<int> eventCuts;
137139

138140
Filter collisionFilter = nabs(aod::collision::posZ) < cfgCutVertex;
139-
Filter trackFilter = (nabs(aod::track::dcaXY) < cfgCutDCAxy) && (nabs(aod::track::dcaZ) < cfgCutDCAz) && (nabs(aod::track::eta) < cfgCutEta) && (aod::track::pt > cfgCutPtPOIMin) && (aod::track::pt < cfgCutPtPOIMax) && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t) true)) && (aod::track::tpcChi2NCl < cfgCutChi2prTPCcls);
141+
Filter trackFilter = (nabs(aod::track::dcaXY) < cfgCutDCAxy) && (nabs(aod::track::dcaZ) < cfgCutDCAz) && (nabs(aod::track::eta) < cfgCutEta) && (aod::track::pt > cfgCutPtPOIMin) && (aod::track::pt < cfgCutPtPOIMax) && ((requireGlobalTrackInFilter()) || (aod::track::isGlobalTrackSDD == (uint8_t)true)) && (aod::track::tpcChi2NCl < cfgCutChi2prTPCcls);
140142

141143
using AodCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::FV0Mults, aod::TPCMults, aod::CentFV0As, aod::CentFT0Ms, aod::CentFT0Cs, aod::CentFT0As, aod::Mults>>;
142144
using AodTracksWithoutBayes = soa::Filtered<soa::Join<aod::Tracks, aod::TrackSelection, aod::TracksExtra, aod::TracksDCA, aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr, aod::pidTOFbeta, aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr>>;
@@ -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)