Skip to content

Commit 0dc9d91

Browse files
allesgraualibuild
allesgrau
and
alibuild
authored
[PWGCF] Add TOF nsigma PID selections for V0 and cascade daughters and bachelors (#17302)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent bef8c65 commit 0dc9d91

2 files changed

Lines changed: 63 additions & 3 deletions

File tree

‎PWGCF/Femto/Core/cascadeBuilder.h‎

Lines changed: 33 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -76,18 +76,22 @@ struct ConfCascadeFilters : o2::framework::ConfigurableGroup {
7676
o2::framework::Configurable<std::vector<float>> dauAbsDcaxyMin{"dauAbsDcaxyMin", {0.05f}, "Minimum |DCAxy| of the daughters and bachelor from primary vertex (cm)"}; \
7777
o2::framework::Configurable<std::vector<float>> dauTpcClustersMin{"dauTpcClustersMin", {80.f}, "Minimum number of TPC clusters for daughter tracks"}; \
7878
o2::framework::Configurable<std::vector<float>> posDauTpc{"posDauTpc", {5.f}, "Maximum |nsimga_Pion/Proton| TPC for positive daughter tracks"}; \
79-
o2::framework::Configurable<std::vector<float>> negDauTpc{"negDauTpc", {5.f}, "Maximum |nsimga_Pion/Proton| TPC for negative daughter tracks"};
79+
o2::framework::Configurable<std::vector<float>> negDauTpc{"negDauTpc", {5.f}, "Maximum |nsimga_Pion/Proton| TPC for negative daughter tracks"}; \
80+
o2::framework::Configurable<std::vector<float>> posDauTof{"posDauTof", {}, "Maximum |nsimga_Pion/Proton| TOF for positive daughter tracks"}; \
81+
o2::framework::Configurable<std::vector<float>> negDauTof{"negDauTof", {}, "Maximum |nsigma_Pion/Proton| TOF for negative daughter tracks"};
8082

8183
struct ConfXiBits : o2::framework::ConfigurableGroup {
8284
std::string prefix = std::string("XiBits");
8385
CASCADE_DEFAULT_BITS
8486
o2::framework::Configurable<std::vector<float>> bachelorTpcPion{"bachelorTpcPion", {5.f}, "Maximum |nsimga_Pion| TPC for bachelor tracks"};
87+
o2::framework::Configurable<std::vector<float>> bachelorTofPion{"bachelorTofPion", {}, "Maximum |nsimga_Pion| TOF for bachelor tracks"};
8588
};
8689

8790
struct ConfOmegaBits : o2::framework::ConfigurableGroup {
8891
std::string prefix = std::string("OmegaBits");
8992
CASCADE_DEFAULT_BITS
9093
o2::framework::Configurable<std::vector<float>> bachelorTpcKaon{"bachelorTpcKaon", {5.f}, "Maximum |nsimga_Kaon| TPC for bachelor tracks"};
94+
o2::framework::Configurable<std::vector<float>> bachelorTofKaon{"bachelorTofKaon", {}, "Maximum |nsimga_Kaon| TOF for bachelor tracks"};
9195
};
9296

9397
#undef CASCADE_DEFAULT_BITS
@@ -139,10 +143,14 @@ enum CascadeSels {
139143
// PID selection for cascade bachelor
140144
kBachelorTpcPion, ///< TPC Pion PID for bachelor
141145
kBachelorTpcKaon, ///< TPC Kaon PID for bachelor
146+
kBachelorTofPion, ///< TOF Pion PID for bachelor
147+
kBachelorTofKaon, ///< TOF Kaon PID for bachelor
142148
///
143149
// PID selection for lambda daughers
144150
kPosDauTpc, ///< TPC PID for positive daughter
145151
kNegDauTpc, ///< TPC PID for negative daughter
152+
kPosDauTof, ///< TOF PID for positive daughter
153+
kNegDauTof, ///< TOF PID for negative daughter
146154

147155
kCascadeSelsMax
148156
};
@@ -166,9 +174,12 @@ const std::unordered_map<CascadeSels, std::string> cascadeSelectionNames = {
166174

167175
{kBachelorTpcPion, "Bachelor TPC Pion PID"},
168176
{kBachelorTpcKaon, "Bachelor TPC Kaon PID"},
169-
177+
{kBachelorTofPion, "Bachelor TOF Pion PID"},
178+
{kBachelorTofKaon, "Bachelor TOF Kaon PID"},
170179
{kPosDauTpc, "Positive Daughter TPC PID"},
171180
{kNegDauTpc, "Negative Daughter TPC PID"},
181+
{kPosDauTof, "Positive Daughter TOF PID"},
182+
{kNegDauTof, "Negative Daughter TOF PID"},
172183

173184
{kCascadeSelsMax, "Cascade Selections Max"}};
174185

@@ -231,6 +242,7 @@ class CascadeSelection : public baseselection::BaseSelection<float, o2::analysis
231242
mOmegaMassLowerLimit = filter.rejectMassOmegaMin.value;
232243
mOmegaMassUpperLimit = filter.rejectMassOmegaMax.value;
233244
this->addSelection(kBachelorTpcPion, cascadeSelectionNames.at(kBachelorTpcPion), config.bachelorTpcPion.value, limits::kAbsUpperLimit, true, true, false);
245+
this->addSelection(kBachelorTofPion, cascadeSelectionNames.at(kBachelorTofPion), config.bachelorTofPion.value, limits::kAbsUpperLimit, true, true, false);
234246
}
235247
if constexpr (modes::isEqual(cascadeType, modes::Cascade::kOmega)) {
236248
mOmegaMassLowerLimit = filter.massOmegaMin.value;
@@ -239,6 +251,7 @@ class CascadeSelection : public baseselection::BaseSelection<float, o2::analysis
239251
mXiMassLowerLimit = filter.rejectMassXiMin.value;
240252
mXiMassUpperLimit = filter.rejectMassXiMax.value;
241253
this->addSelection(kBachelorTpcKaon, cascadeSelectionNames.at(kBachelorTpcKaon), config.bachelorTpcKaon.value, limits::kAbsUpperLimit, true, true, false);
254+
this->addSelection(kBachelorTofKaon, cascadeSelectionNames.at(kBachelorTofKaon), config.bachelorTofKaon.value, limits::kAbsUpperLimit, true, true, false);
242255
}
243256

244257
mPtMin = filter.ptMin.value;
@@ -252,6 +265,8 @@ class CascadeSelection : public baseselection::BaseSelection<float, o2::analysis
252265

253266
this->addSelection(kPosDauTpc, cascadeSelectionNames.at(kPosDauTpc), config.posDauTpc.value, limits::kAbsUpperLimit, true, true, false);
254267
this->addSelection(kNegDauTpc, cascadeSelectionNames.at(kNegDauTpc), config.negDauTpc.value, limits::kAbsUpperLimit, true, true, false);
268+
this->addSelection(kPosDauTof, cascadeSelectionNames.at(kPosDauTof), config.posDauTof.value, limits::kAbsUpperLimit, true, true, false);
269+
this->addSelection(kNegDauTof, cascadeSelectionNames.at(kNegDauTof), config.negDauTof.value, limits::kAbsUpperLimit, true, true, false);
255270

256271
this->addSelection(kCascadeCpaMin, cascadeSelectionNames.at(kCascadeCpaMin), config.cascadeCpaMin.value, limits::kLowerLimit, true, true, false);
257272
this->addSelection(kCascadeTransRadMin, cascadeSelectionNames.at(kCascadeTransRadMin), config.cascadeTransRadMin.value, limits::kLowerLimit, true, true, false);
@@ -318,14 +333,30 @@ class CascadeSelection : public baseselection::BaseSelection<float, o2::analysis
318333
// check both pion and kaon PID for xi and omega
319334
this->evaluateObservable(kBachelorTpcPion, bachelor.tpcNSigmaPi());
320335
this->evaluateObservable(kBachelorTpcKaon, bachelor.tpcNSigmaKa());
336+
if (bachelor.hasTOF()) {
337+
this->evaluateObservable(kBachelorTofPion, bachelor.tofNSigmaPi());
338+
this->evaluateObservable(kBachelorTofKaon, bachelor.tofNSigmaKa());
339+
}
321340

322341
// depending on the charge, we check lambda or antilambda hypothesis
323342
if (cascade.sign() < 0) {
324343
this->evaluateObservable(kPosDauTpc, posDaughter.tpcNSigmaPr());
325344
this->evaluateObservable(kNegDauTpc, negDaughter.tpcNSigmaPi());
345+
if (posDaughter.hasTOF()) {
346+
this->evaluateObservable(kPosDauTof, posDaughter.tofNSigmaPr());
347+
}
348+
if (negDaughter.hasTOF()) {
349+
this->evaluateObservable(kNegDauTof, negDaughter.tofNSigmaPi());
350+
}
326351
} else if (cascade.sign() > 0) {
327352
this->evaluateObservable(kPosDauTpc, posDaughter.tpcNSigmaPi());
328353
this->evaluateObservable(kNegDauTpc, negDaughter.tpcNSigmaPr());
354+
if (posDaughter.hasTOF()) {
355+
this->evaluateObservable(kPosDauTof, posDaughter.tofNSigmaPi());
356+
}
357+
if (negDaughter.hasTOF()) {
358+
this->evaluateObservable(kNegDauTof, negDaughter.tofNSigmaPr());
359+
}
329360
} else {
330361
LOG(warn) << "Encountered Cascade candidate with 0 charge";
331362
}

‎PWGCF/Femto/Core/v0Builder.h‎

Lines changed: 30 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -81,6 +81,10 @@ struct ConfLambdaBits : o2::framework::ConfigurableGroup {
8181
o2::framework::Configurable<std::vector<float>> posDauTpcProton{"posDauTpcProton", {5.f}, "Maximum |nsimga_Proton| TPC for positive daughter tracks"};
8282
o2::framework::Configurable<std::vector<float>> negDauTpcPion{"negDauTpcPion", {5.f}, "Maximum |nsimga_Pion| TPC for negative daughter tracks"};
8383
o2::framework::Configurable<std::vector<float>> negDauTpcProton{"negDauTpcProton", {5.f}, "Maximum |nsimga_Proton| TPC negative for daughter tracks"};
84+
o2::framework::Configurable<std::vector<float>> posDauTofPion{"posDauTofPion", {}, "Maximum |nsigma_Pion| TOF for positive daughter tracks"};
85+
o2::framework::Configurable<std::vector<float>> posDauTofProton{"posDauTofProton", {}, "Maximum |nsigma_Proton| TOF for positive daughter tracks"};
86+
o2::framework::Configurable<std::vector<float>> negDauTofPion{"negDauTofPion", {}, "Maximum |nsigma_Pion| TOF for negative daughter tracks"};
87+
o2::framework::Configurable<std::vector<float>> negDauTofProton{"negDauTofProton", {}, "Maximum |nsigma_Proton| TOF for negative daughter tracks"};
8488
};
8589

8690
// derived selection bits for K0Short
@@ -89,6 +93,8 @@ struct ConfK0shortBits : o2::framework::ConfigurableGroup {
8993
V0_DEFAULT_BITS
9094
o2::framework::Configurable<std::vector<float>> posDauTpcPion{"posDauTpcPion", {5.f}, "Maximum |nsimga_Pion| TPC for positive daughter tracks"};
9195
o2::framework::Configurable<std::vector<float>> negDauTpcPion{"negDauTpcPion", {5.f}, "Maximum |nsimga_Pion| TPC for negative daughter tracks"};
96+
o2::framework::Configurable<std::vector<float>> posDauTofPion{"posDauTofPion", {}, "Maximum |nsigma_Pion| TOF for positive daughter tracks"};
97+
o2::framework::Configurable<std::vector<float>> negDauTofPion{"negDauTofPion", {}, "Maximum |nsigma_Pion| TOF for negative daughter tracks"};
9298
};
9399

94100
#undef V0_DEFAULT_BITS
@@ -154,6 +160,11 @@ enum V0Sels {
154160
kNegDaughTpcPion, ///< TPC Pion PID for negative daughter
155161
kNegDaughTpcProton, ///< TPC Proton PID for negative daughter
156162

163+
kPosDaughTofPion, ///< TOF Pion PID for positive daughter
164+
kPosDaughTofProton, ///< TOF Proton PID for positive daughter
165+
kNegDaughTofPion, ///< TOF Pion PID for negative daughter
166+
kNegDaughTofProton, ///< TOF Proton PID for negative daughter
167+
157168
kV0SelsMax
158169
};
159170

@@ -175,7 +186,11 @@ const std::unordered_map<V0Sels, std::string> v0SelectionNames = {
175186
{kPosDaughTpcPion, "TPC Pion PID for positive daughter"},
176187
{kPosDaughTpcProton, "TPC Proton PID for positive daughter"},
177188
{kNegDaughTpcPion, "TPC Pion PID for negative daughter"},
178-
{kNegDaughTpcProton, "TPC Proton PID for negative daughter"}};
189+
{kNegDaughTpcProton, "TPC Proton PID for negative daughter"},
190+
{kPosDaughTofPion, "TOF Pion PID for positive daughter"},
191+
{kPosDaughTofProton, "TOF Proton PID for positive daughter"},
192+
{kNegDaughTofPion, "TOF Pion PID for negative daughter"},
193+
{kNegDaughTofProton, "TOF Proton PID for negative daughter"}};
179194

180195
// enum for all track filters (loose pre-selection, applied before quality/PID cuts)
181196
enum V0Filters {
@@ -242,11 +257,15 @@ class V0Selection : public baseselection::BaseSelection<float, datatypes::V0Mask
242257
if constexpr (modes::isEqual(v0Type, modes::V0::kLambda)) {
243258
this->addSelection(kPosDaughTpcProton, v0SelectionNames.at(kPosDaughTpcProton), config.posDauTpcProton.value, limits::kAbsUpperLimit, true, true, false);
244259
this->addSelection(kNegDaughTpcPion, v0SelectionNames.at(kNegDaughTpcPion), config.negDauTpcPion.value, limits::kAbsUpperLimit, true, true, false);
260+
this->addSelection(kPosDaughTofProton, v0SelectionNames.at(kPosDaughTofProton), config.posDauTofProton.value, limits::kAbsUpperLimit, true, true, false);
261+
this->addSelection(kNegDaughTofPion, v0SelectionNames.at(kNegDaughTofPion), config.negDauTofPion.value, limits::kAbsUpperLimit, true, true, false);
245262
}
246263

247264
if constexpr (modes::isEqual(v0Type, modes::V0::kAntiLambda)) {
248265
this->addSelection(kPosDaughTpcPion, v0SelectionNames.at(kPosDaughTpcPion), config.posDauTpcPion.value, limits::kAbsUpperLimit, true, true, false);
249266
this->addSelection(kNegDaughTpcProton, v0SelectionNames.at(kNegDaughTpcProton), config.negDauTpcProton.value, limits::kAbsUpperLimit, true, true, false);
267+
this->addSelection(kPosDaughTofPion, v0SelectionNames.at(kPosDaughTofPion), config.posDauTofPion.value, limits::kAbsUpperLimit, true, true, false);
268+
this->addSelection(kNegDaughTofProton, v0SelectionNames.at(kNegDaughTofProton), config.negDauTofProton.value, limits::kAbsUpperLimit, true, true, false);
250269
}
251270
}
252271
if constexpr (modes::isEqual(v0Type, modes::V0::kK0short)) {
@@ -258,6 +277,8 @@ class V0Selection : public baseselection::BaseSelection<float, datatypes::V0Mask
258277

259278
this->addSelection(kPosDaughTpcPion, v0SelectionNames.at(kPosDaughTpcPion), config.posDauTpcPion.value, limits::kAbsUpperLimit, true, true, false);
260279
this->addSelection(kNegDaughTpcPion, v0SelectionNames.at(kNegDaughTpcPion), config.negDauTpcPion.value, limits::kAbsUpperLimit, true, true, false);
280+
this->addSelection(kPosDaughTofPion, v0SelectionNames.at(kPosDaughTofPion), config.posDauTofPion.value, limits::kAbsUpperLimit, true, true, false);
281+
this->addSelection(kNegDaughTofPion, v0SelectionNames.at(kNegDaughTofPion), config.negDauTofPion.value, limits::kAbsUpperLimit, true, true, false);
261282
}
262283

263284
this->addSelection(kDcaDaughMax, v0SelectionNames.at(kDcaDaughMax), config.dcaDauMax.value, limits::kAbsUpperLimit, true, true, false);
@@ -321,6 +342,14 @@ class V0Selection : public baseselection::BaseSelection<float, datatypes::V0Mask
321342
this->evaluateObservable(kPosDaughTpcProton, posDaughter.tpcNSigmaPr());
322343
this->evaluateObservable(kNegDaughTpcPion, negDaughter.tpcNSigmaPi());
323344
this->evaluateObservable(kNegDaughTpcProton, negDaughter.tpcNSigmaPr());
345+
if (posDaughter.hasTOF()) {
346+
this->evaluateObservable(kPosDaughTofPion, posDaughter.tofNSigmaPi());
347+
this->evaluateObservable(kPosDaughTofProton, posDaughter.tofNSigmaPr());
348+
}
349+
if (negDaughter.hasTOF()) {
350+
this->evaluateObservable(kNegDaughTofPion, negDaughter.tofNSigmaPi());
351+
this->evaluateObservable(kNegDaughTofProton, negDaughter.tofNSigmaPr());
352+
}
324353

325354
this->assembleBitmask<SelectionHistName>();
326355
}

0 commit comments

Comments
 (0)