Skip to content

Commit 4d4d2ac

Browse files
mfagginmfaggin
committed
Replace legacy decay channels for SigmaC.
1 parent baa254a commit 4d4d2ac

File tree

3 files changed

+32
-18
lines changed

3 files changed

+32
-18
lines changed

‎PWGHF/Core/DecayChannels.h‎

Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -259,6 +259,22 @@ enum DecayChannelMain : HfDecayChannel {
259259
NChannelsMain = DstarToD0Pi // last channel
260260
};
261261
} // namespace hf_cand_reso
262+
263+
/// @brief Sigmac candidates: main channels
264+
namespace hf_cand_sigmac
265+
{
266+
enum DecayChannelMain : HfDecayChannel {
267+
// Σc0(2455) → Λc+ π- → p K- π+ π-
268+
Sc0ToPKPiPi = 0,
269+
// Σc++(2455) → Λc+ π+ → p K- π+ π+
270+
ScplusplusToPKPiPi,
271+
// Σc0(2520) → Λc+ π- → p K- π+ π-
272+
ScStar0ToPKPiPi,
273+
// Σc++(2520) → Λc+ π+ → p K- π+ π+
274+
ScStarPlusPlusToPKPiPi
275+
};
276+
} // namespace hf_cand_sigmac
277+
262278
} // namespace o2::hf_decay
263279

264280
#endif // PWGHF_CORE_DECAYCHANNELS_H_

‎PWGHF/D2H/Tasks/taskSigmac.cxx‎

Lines changed: 8 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -16,7 +16,6 @@
1616
/// \author Mattia Faggin <mfaggin@cern.ch>, University and INFN PADOVA
1717

1818
#include "PWGHF/Core/DecayChannels.h"
19-
#include "PWGHF/Core/DecayChannelsLegacy.h"
2019
#include "PWGHF/Core/HfHelper.h"
2120
#include "PWGHF/D2H/Utils/utilsSigmac.h"
2221
#include "PWGHF/DataModel/AliasTables.h"
@@ -603,10 +602,10 @@ struct HfTaskSigmac {
603602
for (const auto& particle : mcParticlesSc) {
604603

605604
/// reject immediately particles different from Σc0,++
606-
bool const isSc0Gen = (std::abs(particle.flagMcMatchGen()) == BIT(aod::hf_cand_sigmac::DecayType::Sc0ToPKPiPi));
607-
bool const isScStar0Gen = (std::abs(particle.flagMcMatchGen()) == BIT(aod::hf_cand_sigmac::DecayType::ScStar0ToPKPiPi));
608-
bool const isScPlusPlusGen = (std::abs(particle.flagMcMatchGen()) == BIT(aod::hf_cand_sigmac::DecayType::ScplusplusToPKPiPi));
609-
bool const isScStarPlusPlusGen = (std::abs(particle.flagMcMatchGen()) == BIT(aod::hf_cand_sigmac::DecayType::ScStarPlusPlusToPKPiPi));
605+
bool const isSc0Gen = (std::abs(particle.flagMcMatchGen()) == BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::Sc0ToPKPiPi));
606+
bool const isScStar0Gen = (std::abs(particle.flagMcMatchGen()) == BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStar0ToPKPiPi));
607+
bool const isScPlusPlusGen = (std::abs(particle.flagMcMatchGen()) == BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScplusplusToPKPiPi));
608+
bool const isScStarPlusPlusGen = (std::abs(particle.flagMcMatchGen()) == BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStarPlusPlusToPKPiPi));
610609
if (!isSc0Gen && !isScPlusPlusGen && !isScStar0Gen && !isScStarPlusPlusGen) {
611610
continue;
612611
}
@@ -799,10 +798,10 @@ struct HfTaskSigmac {
799798

800799
// candidateLc.flagMcDecayChanRec();
801800

802-
bool const isTrueSc0Reco = std::abs(candSc.flagMcMatchRec()) == BIT(aod::hf_cand_sigmac::DecayType::Sc0ToPKPiPi);
803-
bool const isTrueScStar0Reco = std::abs(candSc.flagMcMatchRec()) == BIT(aod::hf_cand_sigmac::DecayType::ScStar0ToPKPiPi);
804-
bool const isTrueScPlusPlusReco = std::abs(candSc.flagMcMatchRec()) == BIT(aod::hf_cand_sigmac::DecayType::ScplusplusToPKPiPi);
805-
bool const isTrueScStarPlusPlusReco = std::abs(candSc.flagMcMatchRec()) == BIT(aod::hf_cand_sigmac::DecayType::ScStarPlusPlusToPKPiPi);
801+
bool const isTrueSc0Reco = std::abs(candSc.flagMcMatchRec()) == BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::Sc0ToPKPiPi);
802+
bool const isTrueScStar0Reco = std::abs(candSc.flagMcMatchRec()) == BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStar0ToPKPiPi);
803+
bool const isTrueScPlusPlusReco = std::abs(candSc.flagMcMatchRec()) == BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScplusplusToPKPiPi);
804+
bool const isTrueScStarPlusPlusReco = std::abs(candSc.flagMcMatchRec()) == BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStarPlusPlusToPKPiPi);
806805
int sigmacSpecies = -1;
807806
if ((isTrueSc0Reco || isTrueScStar0Reco) && (chargeSc == o2::aod::hf_cand_sigmac::ChargeNull)) {
808807
/// Reconstructed Σc0 signal

‎PWGHF/TableProducer/candidateCreatorSigmac0plusplus.cxx‎

Lines changed: 8 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -17,7 +17,6 @@
1717

1818
#include "PWGHF/Core/CentralityEstimation.h"
1919
#include "PWGHF/Core/DecayChannels.h"
20-
#include "PWGHF/Core/DecayChannelsLegacy.h"
2120
#include "PWGHF/Core/HfHelper.h"
2221
#include "PWGHF/Core/SelectorCuts.h"
2322
#include "PWGHF/DataModel/AliasTables.h"
@@ -519,7 +518,7 @@ struct HfCandidateSigmac0plusplusMc {
519518
/// look for Σc0(2455)
520519
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughters, Pdg::kSigmaC0, std::array{+kProton, -kKPlus, +kPiPlus, -kPiPlus}, true, &sign, 3);
521520
if (indexRec > -1) { /// due to (*) no need to check anything for LambdaC
522-
flag = sign * BIT(aod::hf_cand_sigmac::DecayType::Sc0ToPKPiPi);
521+
flag = sign * BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::Sc0ToPKPiPi);
523522
}
524523
auto particle = mcParticles.rawIteratorAt(indexRec);
525524
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaC0);
@@ -528,7 +527,7 @@ struct HfCandidateSigmac0plusplusMc {
528527
if (flag == 0) {
529528
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughters, Pdg::kSigmaCStar0, std::array{+kProton, -kKPlus, +kPiPlus, -kPiPlus}, true, &sign, 3);
530529
if (indexRec > -1) { /// due to (*) no need to check anything for LambdaC
531-
flag = sign * BIT(aod::hf_cand_sigmac::DecayType::ScStar0ToPKPiPi);
530+
flag = sign * BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStar0ToPKPiPi);
532531
}
533532
if (particleAntiparticle < 0) {
534533
auto particle = mcParticles.rawIteratorAt(indexRec);
@@ -546,7 +545,7 @@ struct HfCandidateSigmac0plusplusMc {
546545
/// look for Σc++(2455)
547546
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughters, Pdg::kSigmaCPlusPlus, std::array{+kProton, -kKPlus, +kPiPlus, +kPiPlus}, true, &sign, 3);
548547
if (indexRec > -1) { /// due to (*) no need to check anything for LambdaC
549-
flag = sign * BIT(aod::hf_cand_sigmac::DecayType::ScplusplusToPKPiPi);
548+
flag = sign * BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScplusplusToPKPiPi);
550549
}
551550
auto particle = mcParticles.rawIteratorAt(indexRec);
552551
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaCPlusPlus);
@@ -555,7 +554,7 @@ struct HfCandidateSigmac0plusplusMc {
555554
if (flag == 0) {
556555
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughters, Pdg::kSigmaCStarPlusPlus, std::array{+kProton, -kKPlus, +kPiPlus, +kPiPlus}, true, &sign, 3);
557556
if (indexRec > -1) { /// due to (*) no need to check anything for LambdaC
558-
flag = sign * BIT(aod::hf_cand_sigmac::DecayType::ScStarPlusPlusToPKPiPi);
557+
flag = sign * BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStarPlusPlusToPKPiPi);
559558
}
560559
if (particleAntiparticle < 0) {
561560
auto particle = mcParticles.rawIteratorAt(indexRec);
@@ -614,7 +613,7 @@ struct HfCandidateSigmac0plusplusMc {
614613
}
615614
if (std::abs(daughter.flagMcMatchGen()) == o2::hf_decay::hf_cand_3prong::DecayChannelMain::LcToPKPi) {
616615
/// Λc+ daughter decaying in pK-π+ found!
617-
flag = sign * BIT(aod::hf_cand_sigmac::DecayType::Sc0ToPKPiPi);
616+
flag = sign * BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::Sc0ToPKPiPi);
618617
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaC0);
619618
break;
620619
}
@@ -628,7 +627,7 @@ struct HfCandidateSigmac0plusplusMc {
628627
}
629628
if (std::abs(daughter.flagMcMatchGen()) == o2::hf_decay::hf_cand_3prong::DecayChannelMain::LcToPKPi) {
630629
/// Λc+ daughter decaying in pK-π+ found!
631-
flag = sign * BIT(aod::hf_cand_sigmac::DecayType::ScplusplusToPKPiPi);
630+
flag = sign * BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScplusplusToPKPiPi);
632631
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaCPlusPlus);
633632
break;
634633
}
@@ -646,7 +645,7 @@ struct HfCandidateSigmac0plusplusMc {
646645
}
647646
if (std::abs(daughter.flagMcMatchGen()) == o2::hf_decay::hf_cand_3prong::DecayChannelMain::LcToPKPi) {
648647
/// Λc+ daughter decaying in pK-π+ found!
649-
flag = sign * BIT(aod::hf_cand_sigmac::DecayType::ScStar0ToPKPiPi);
648+
flag = sign * BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStar0ToPKPiPi);
650649
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaCStar0);
651650
break;
652651
}
@@ -660,7 +659,7 @@ struct HfCandidateSigmac0plusplusMc {
660659
}
661660
if (std::abs(daughter.flagMcMatchGen()) == o2::hf_decay::hf_cand_3prong::DecayChannelMain::LcToPKPi) {
662661
/// Λc+ daughter decaying in pK-π+ found!
663-
flag = sign * BIT(aod::hf_cand_sigmac::DecayType::ScStarPlusPlusToPKPiPi);
662+
flag = sign * BIT(o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStarPlusPlusToPKPiPi);
664663
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaCStarPlusPlus);
665664
break;
666665
}

0 commit comments

Comments
 (0)