Skip to content

Commit 67f45b1

Browse files
authored
[PWGCF] FemtoUniverse: Add eta, pT cuts on MCTruth tracks (#13457)
1 parent 4404428 commit 67f45b1

File tree

1 file changed

+6
-1
lines changed

1 file changed

+6
-1
lines changed

PWGCF/FemtoUniverse/TableProducer/femtoUniverseProducerTask.cxx

Lines changed: 6 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -240,6 +240,7 @@ struct FemtoUniverseProducerTask {
240240
Configurable<bool> confDcaXYCustom1Cut{"confDcaXYCustom1Cut", true, "Enable Custom |DCAxy| < [1] + [2]/pt cut."};
241241
Configurable<float> confDcaXYCustom11FilterCut{"confDcaXYCustom11FilterCut", 0.004, "Value for [1] custom DCAxy cut -> |DCAxy| < [1] + [2]/pT"};
242242
Configurable<float> confDcaXYCustom12FilterCut{"confDcaXYCustom12FilterCut", 0.013, "Value for [2] custom DCAxy cut -> |DCAxy| < [1] + [2]/pT"};
243+
Configurable<bool> confIsApplyTrkCutMCTruth{"confIsApplyTrkCutMCTruth", false, "Apply eta, pT selection cut on MCTruth tracks "};
243244
} ConfFilterCuts;
244245

245246
Filter globalCutFilter = requireGlobalTrackInFilter();
@@ -1361,7 +1362,11 @@ struct FemtoUniverseProducerTask {
13611362
void fillTracksMCTruth(MCParticlesType const& mcParticles)
13621363
{
13631364
for (const auto& mc : mcParticles) { // Loop over all MC Truth particles
1364-
1365+
if (ConfFilterCuts.confIsApplyTrkCutMCTruth) {
1366+
if (std::abs(mc.eta()) > ConfFilterCuts.confEtaFilterCut || mc.pt() < ConfFilterCuts.confPtLowFilterCut || mc.pt() > ConfFilterCuts.confPtHighFilterCut) {
1367+
continue;
1368+
}
1369+
}
13651370
std::vector<int> childIDs = {0, 0};
13661371
outputParts(outputCollision.lastIndex(),
13671372
mc.pt(),

0 commit comments

Comments
 (0)