Skip to content

Commit d222cd2

Browse files
committed
Use old flag scheme for 2prongs
1 parent f5abc56 commit d222cd2

File tree

4 files changed

+107
-107
lines changed

4 files changed

+107
-107
lines changed

PWGHF/TableProducer/candidateCreator2Prong.cxx

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -61,8 +61,8 @@ using namespace o2::analysis;
6161
using namespace o2::hf_evsel;
6262
using namespace o2::hf_trkcandsel;
6363
using namespace o2::aod::hf_cand_2prong;
64+
using namespace o2::hf_decay;
6465
using namespace o2::hf_decay::hf_cand_2prong;
65-
using namespace o2::hf_corrbkg;
6666
using namespace o2::hf_centrality;
6767
using namespace o2::hf_occupancy;
6868
using namespace o2::constants::physics;
@@ -779,7 +779,7 @@ struct HfCandidateCreator2ProngExpressions {
779779
if (matchCorrBkgs) {
780780
indexRec = -1; // Index of the matched reconstructed candidate
781781
int depth = 2;
782-
for (const auto& [chn, finalState] : hf_chns_2prong::finalStates2Prongs) {
782+
for (const auto& [chn, finalState] : hf_cand_2prong::finalStates2Prongs) {
783783
std::array<int, 2> finalStateParts2Prong = std::array{finalState[0], finalState[1]};
784784
if (finalState.size() == 3) { // Partly Reco 2-prong decays
785785
if (matchKinkedDecayTopology && matchInteractionsWithMaterial) {
@@ -821,7 +821,7 @@ struct HfCandidateCreator2ProngExpressions {
821821
continue;
822822
}
823823
if (indexRec > -1) {
824-
flag = sign * chn;
824+
flag = sign * (1 << chn);
825825

826826
// Flag the resonant decay channel
827827
int resoMaxDepth = 1;

PWGHF/TableProducer/candidateCreator3Prong.cxx

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -59,14 +59,14 @@ using namespace o2;
5959
using namespace o2::hf_evsel;
6060
using namespace o2::hf_trkcandsel;
6161
using namespace o2::aod::hf_cand_3prong;
62+
using namespace o2::hf_decay;
6263
using namespace o2::hf_decay::hf_cand_3prong;
6364
using namespace o2::hf_centrality;
6465
using namespace o2::hf_occupancy;
6566
using namespace o2::constants::physics;
6667
using namespace o2::framework;
6768
using namespace o2::framework::expressions;
6869
using namespace o2::aod::pid_tpc_tof_utils;
69-
using namespace o2::hf_corrbkg;
7070

7171
/// Reconstruction of heavy-flavour 3-prong decay candidates
7272
struct HfCandidateCreator3Prong {

PWGHF/Utils/utilsMcGen.h

Lines changed: 5 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -35,8 +35,6 @@
3535
namespace hf_mc_gen
3636
{
3737

38-
using namespace o2::hf_corrbkg;
39-
4038
template <bool matchCorrBkgs = false, typename T, typename U, typename V>
4139
void fillMcMatchGen2Prong(T const& mcParticles, U const& mcParticlesPerMcColl, V& rowMcMatchGen, bool rejectBackground)
4240
{
@@ -59,7 +57,7 @@ void fillMcMatchGen2Prong(T const& mcParticles, U const& mcParticlesPerMcColl, V
5957
if (matchCorrBkgs) {
6058
bool matched = false;
6159

62-
for (const auto& [chn, finalState] : o2::hf_corrbkg::hf_chns_2prong::finalStates2Prongs) {
60+
for (const auto& [chn, finalState] : o2::hf_decay::hf_cand_2prong::finalStates2Prongs) {
6361
if (finalState.size() == 3) { // Partly Reco 3-prong decays
6462
std::array<int, 3> finalStateParts = std::array{finalState[0], finalState[1], finalState[2]};
6563
if (particle.pdgCode() < 0) {
@@ -78,7 +76,7 @@ void fillMcMatchGen2Prong(T const& mcParticles, U const& mcParticlesPerMcColl, V
7876
continue;
7977
}
8078
if (matched) {
81-
flag = sign * chn;
79+
flag = sign * (1 << chn);
8280

8381
// Flag the resonant decay channel
8482
int resoMaxDepth = 1;
@@ -90,7 +88,7 @@ void fillMcMatchGen2Prong(T const& mcParticles, U const& mcParticlesPerMcColl, V
9088
auto daughI = mcParticles.rawIteratorAt(arrResoDaughIndex[iProng]);
9189
arrPDGDaugh[iProng] = daughI.pdgCode();
9290
}
93-
flagResonantDecay(Pdg::kD0, &channel, arrPDGDaugh);
91+
o2::hf_decay::flagResonantDecay(Pdg::kD0, &channel, arrPDGDaugh);
9492
}
9593
break;
9694
}
@@ -159,7 +157,7 @@ void fillMcMatchGen3Prong(T const& mcParticles, U const& mcParticlesPerMcColl, V
159157
if (std::abs(particle.pdgCode()) != motherPdgCode) {
160158
continue; // Skip if the particle PDG code does not match the mother PDG code
161159
}
162-
auto finalStates = getDecayChannel3Prong(motherPdgCode);
160+
auto finalStates = o2::hf_decay::hf_cand_3prong::getDecayChannel3Prong(motherPdgCode);
163161
int maxDepth = 2;
164162
bool matched = false;
165163
if (motherPdgCode == Pdg::kDStar) {
@@ -224,7 +222,7 @@ void fillMcMatchGen3Prong(T const& mcParticles, U const& mcParticlesPerMcColl, V
224222
auto daughI = mcParticles.rawIteratorAt(arrResoDaughIndex[iProng]);
225223
arrPDGDaugh[iProng] = daughI.pdgCode();
226224
}
227-
flagResonantDecay<true>(motherPdgCode, &channel, arrPDGDaugh);
225+
o2::hf_decay::flagResonantDecay<true>(motherPdgCode, &channel, arrPDGDaugh);
228226
}
229227
break; // Exit loop if a match is found
230228
}

0 commit comments

Comments
 (0)