Skip to content

Commit 63ec9e1

Browse files
ddobrigkalibuild
andauthored
[Common] legacy centrality QA: add vertex Z cut (AliceO2Group#16317)
Co-authored-by: ALICE Builder <alibuild@users.noreply.github.com>
1 parent 592883d commit 63ec9e1

1 file changed

Lines changed: 31 additions & 6 deletions

File tree

Common/Tasks/centralityQa.cxx

Lines changed: 31 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,7 @@ using namespace o2::framework;
2828
struct CentralityQa {
2929
Configurable<int> nBins{"nBins", 1050, "number of bins"};
3030
Configurable<bool> INELgtZERO{"INELgtZERO", 1, "0 - no, 1 - yes"};
31+
Configurable<float> vertexZcut{"vertexZcut", 10, "vertex-Z position cut (absolute value in cm)"};
3132
OutputObj<TH1F> hCentRun2V0M{TH1F("hCentRun2V0M", "V0M", nBins, 0, 105.)};
3233
OutputObj<TH1F> hCentRun2V0A{TH1F("hCentRun2V0A", "V0A", nBins, 0, 105.)};
3334
OutputObj<TH1F> hCentRun2SPDTks{TH1F("hCentRun2SPDTks", "SPD Tracklets", nBins, 0, 105.)};
@@ -115,8 +116,12 @@ struct CentralityQa {
115116
if (INELgtZERO && col.multNTracksPVeta1() < 1) {
116117
return;
117118
}
118-
if (!col.sel8())
119+
if (!col.sel8()) {
119120
return;
121+
}
122+
if (std::abs(col.posZ()) > vertexZcut) {
123+
return;
124+
}
120125
LOGF(debug, "centFV0A=%.0f", col.centFV0A());
121126
hCentFV0A->Fill(col.centFV0A());
122127
hCentProfileFV0A->Fill(col.centFV0A(), col.multNTracksPVetaHalf());
@@ -128,8 +133,12 @@ struct CentralityQa {
128133
if (INELgtZERO && col.multNTracksPVeta1() < 1) {
129134
return;
130135
}
131-
if (!col.sel8())
136+
if (!col.sel8()) {
137+
return;
138+
}
139+
if (std::abs(col.posZ()) > vertexZcut) {
132140
return;
141+
}
133142
LOGF(debug, "centFT0M=%.0f", col.centFT0M());
134143
hCentFT0M->Fill(col.centFT0M());
135144
hCentProfileFT0M->Fill(col.centFT0M(), col.multNTracksPVetaHalf());
@@ -141,8 +150,12 @@ struct CentralityQa {
141150
if (INELgtZERO && col.multNTracksPVeta1() < 1) {
142151
return;
143152
}
144-
if (!col.sel8())
153+
if (!col.sel8()) {
154+
return;
155+
}
156+
if (std::abs(col.posZ()) > vertexZcut) {
145157
return;
158+
}
146159
hCentFT0A->Fill(col.centFT0A());
147160
hCentProfileFT0A->Fill(col.centFT0A(), col.multNTracksPVetaHalf());
148161
}
@@ -153,8 +166,12 @@ struct CentralityQa {
153166
if (INELgtZERO && col.multNTracksPVeta1() < 1) {
154167
return;
155168
}
156-
if (!col.sel8())
169+
if (!col.sel8()) {
157170
return;
171+
}
172+
if (std::abs(col.posZ()) > vertexZcut) {
173+
return;
174+
}
158175
hCentFT0C->Fill(col.centFT0C());
159176
hCentProfileFT0C->Fill(col.centFT0C(), col.multNTracksPVetaHalf());
160177
}
@@ -165,8 +182,12 @@ struct CentralityQa {
165182
if (INELgtZERO && col.multNTracksPVeta1() < 1) {
166183
return;
167184
}
168-
if (!col.sel8())
185+
if (!col.sel8()) {
186+
return;
187+
}
188+
if (std::abs(col.posZ()) > vertexZcut) {
169189
return;
190+
}
170191
hCentFDDM->Fill(col.centFDDM());
171192
hCentProfileFDDM->Fill(col.centFDDM(), col.multNTracksPVetaHalf());
172193
}
@@ -177,8 +198,12 @@ struct CentralityQa {
177198
if (INELgtZERO && col.multNTracksPVeta1() < 1) {
178199
return;
179200
}
180-
if (!col.sel8())
201+
if (!col.sel8()) {
202+
return;
203+
}
204+
if (std::abs(col.posZ()) > vertexZcut) {
181205
return;
206+
}
182207
hCentNTPV->Fill(col.centNTPV());
183208
hCentProfileNTPV->Fill(col.centNTPV(), col.multNTracksPVetaHalf());
184209
}

0 commit comments

Comments
 (0)