Skip to content

Commit b363810

Browse files
authored
Fix match of decays including resonant channels (#985)
1 parent bc358c4 commit b363810

2 files changed

Lines changed: 4 additions & 4 deletions

File tree

EventFiltering/PWGHF/HFFilter.cxx

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1012,13 +1012,13 @@ struct HfFilter { // Main struct for HF triggers
10121012
int8_t origin = 0;
10131013

10141014
// D± → π± K∓ π±
1015-
auto indexRec = RecoDecay::getMatchedMCRec(particlesMC, arrayDaughters, pdg::Code::kDPlus, array{+kPiPlus, -kKPlus, +kPiPlus}, true, &sign);
1015+
auto indexRec = RecoDecay::getMatchedMCRec(particlesMC, arrayDaughters, pdg::Code::kDPlus, array{+kPiPlus, -kKPlus, +kPiPlus}, true, &sign, 2);
10161016
if (indexRec >= 0) {
10171017
channel = kDplus;
10181018
}
10191019
if (indexRec < 0) {
10201020
// Ds± → K± K∓ π±
1021-
indexRec = RecoDecay::getMatchedMCRec(particlesMC, arrayDaughters, 431, array{+kKPlus, -kKPlus, +kPiPlus}, true, &sign); // TODO: replace hard coded pdg code
1021+
indexRec = RecoDecay::getMatchedMCRec(particlesMC, arrayDaughters, 431, array{+kKPlus, -kKPlus, +kPiPlus}, true, &sign, 2); // TODO: replace hard coded pdg code
10221022
if (indexRec >= 0) {
10231023
channel = kDs;
10241024
}
@@ -1032,7 +1032,7 @@ struct HfFilter { // Main struct for HF triggers
10321032
}
10331033
if (indexRec < 0) {
10341034
// Ξc± → p± K∓ π±
1035-
indexRec = RecoDecay::getMatchedMCRec(particlesMC, arrayDaughters, pdg::Code::kXiCPlus, array{+kProton, -kKPlus, +kPiPlus}, true, &sign);
1035+
indexRec = RecoDecay::getMatchedMCRec(particlesMC, arrayDaughters, pdg::Code::kXiCPlus, array{+kProton, -kKPlus, +kPiPlus}, true, &sign, 2);
10361036
if (indexRec >= 0) {
10371037
channel = kXic;
10381038
}

EventFiltering/PWGHF/HFFilterQC.cxx

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -110,7 +110,7 @@ struct HfFilterQc { // Main struct for HF trigger QC
110110
return;
111111
} else {
112112
int8_t sign = 0;
113-
if (RecoDecay::isMatchedMCGen(particlesMC, particle, pdgCodes[I], std::get<I>(pdgDau), true, &sign)) {
113+
if (RecoDecay::isMatchedMCGen(particlesMC, particle, pdgCodes[I], std::get<I>(pdgDau), true, &sign, 2)) {
114114
nParticles[I]++;
115115
hPtDistr[0]->Fill(static_cast<double>(I), static_cast<double>(particle.pt()));
116116
for (auto iTrig = 0; iTrig < kNtriggersHF; ++iTrig) {

0 commit comments

Comments
 (0)