Skip to content

Commit f867878

Browse files
vkuceravkucera
committed
tmp
1 parent 846cd4b commit f867878

File tree

7 files changed

+46
-38
lines changed

7 files changed

+46
-38
lines changed

‎PWGHF/D2H/TableProducer/dataCreatorCharmHadPiReduced.cxx‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -598,7 +598,7 @@ struct HfDataCreatorCharmHadPiReduced {
598598
}
599599
// Ds- → Phi π- → K- K+ π- and Ds- → K0* K- → K- K+ π-
600600
if ((arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[0] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[1]) || (arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[1] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[0])) {
601-
flag = sign * BIT(hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi);
601+
flag = sign * DecayChannelMain::BsToDsPi;
602602
} else if ((arrPDGDaughDs[0] == arrPDGResonantDKstarK[0] && arrPDGDaughDs[1] == arrPDGResonantDKstarK[1]) || (arrPDGDaughDs[0] == arrPDGResonantDKstarK[1] && arrPDGDaughDs[1] == arrPDGResonantDKstarK[0])) {
603603
flag = sign * BIT(hf_cand_bs::DecayTypeMc::BsToDsPiToK0starKPiToKKPiPi);
604604
}
@@ -1407,7 +1407,7 @@ struct HfDataCreatorCharmHadPiReduced {
14071407
}
14081408
// Ds- → Phi π- → K- K+ π- and Ds- → K0* K- → K- K+ π-
14091409
if ((arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[0] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[1]) || (arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[1] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[0])) {
1410-
flag = sign * BIT(hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi);
1410+
flag = sign * DecayChannelMain::BsToDsPi;
14111411
} else if ((arrPDGDaughDs[0] == arrPDGResonantDKstarK[0] && arrPDGDaughDs[1] == arrPDGResonantDKstarK[1]) || (arrPDGDaughDs[0] == arrPDGResonantDKstarK[1] && arrPDGDaughDs[1] == arrPDGResonantDKstarK[0])) {
14121412
flag = sign * BIT(hf_cand_bs::DecayTypeMc::BsToDsPiToK0starKPiToKKPiPi);
14131413
}
@@ -1444,7 +1444,7 @@ struct HfDataCreatorCharmHadPiReduced {
14441444
}
14451445

14461446
// save information for Bs task
1447-
if (!TESTBIT(std::abs(flag), hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi) && !TESTBIT(std::abs(flag), hf_cand_bs::DecayTypeMc::BsToDsPiToK0starKPiToKKPiPi) &&
1447+
if (!TESTBIT(std::abs(flag), DecayChannelMain::BsToDsPi) && !TESTBIT(std::abs(flag), hf_cand_bs::DecayTypeMc::BsToDsPiToK0starKPiToKKPiPi) &&
14481448
!TESTBIT(std::abs(flag), hf_cand_bs::DecayTypeMc::B0ToDsPiToPhiPiPiToKKPiPi) && !TESTBIT(std::abs(flag), hf_cand_bs::DecayTypeMc::B0ToDsPiToK0starKPiToKKPiPi)) {
14491449
continue;
14501450
}

‎PWGHF/D2H/Tasks/taskBs.cxx‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -218,7 +218,7 @@ struct HfTaskBs {
218218
auto invMassCandBs = hfHelper.invMassBsToDsPi(candidate);
219219
int flagMcMatchRecBs = std::abs(candidate.flagMcMatchRec());
220220

221-
if (TESTBIT(flagMcMatchRecBs, hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi)) {
221+
if (TESTBIT(flagMcMatchRecBs, DecayChannelMain::BsToDsPi)) {
222222
auto indexMother = RecoDecay::getMother(mcParticles, candidate.prong1_as<aod::TracksWMc>().mcParticle_as<soa::Join<aod::McParticles, aod::HfCandBsMcGen>>(), o2::constants::physics::Pdg::kBS, true);
223223
auto particleMother = mcParticles.rawIteratorAt(indexMother);
224224

@@ -277,7 +277,7 @@ struct HfTaskBs {
277277

278278
// MC gen. level
279279
for (const auto& particle : mcParticles) {
280-
if (TESTBIT(std::abs(particle.flagMcMatchGen()), hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi)) {
280+
if (TESTBIT(std::abs(particle.flagMcMatchGen()), DecayChannelMain::BsToDsPi)) {
281281

282282
auto ptParticle = particle.pt();
283283
auto yParticle = RecoDecay::y(particle.pVector(), o2::constants::physics::MassBS);

‎PWGHF/D2H/Tasks/taskBsReduced.cxx‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -368,7 +368,7 @@ struct HfTaskBsReduced {
368368
if constexpr (doMc) {
369369
flagMcMatchRec = candidate.flagMcMatchRec();
370370
flagWrongCollision = candidate.flagWrongCollision();
371-
isSignal = TESTBIT(std::abs(flagMcMatchRec), hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi);
371+
isSignal = TESTBIT(std::abs(flagMcMatchRec), DecayChannelMain::BsToDsPi);
372372
}
373373

374374
if (fillHistograms) {
@@ -572,7 +572,7 @@ struct HfTaskBsReduced {
572572
void fillCandMcGen(aod::HfMcGenRedBss::iterator const& particle)
573573
{
574574
// keep only generated Bs with the analysis decay channel
575-
if (!TESTBIT(std::abs(particle.flagMcMatchGen()), hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi)) {
575+
if (!TESTBIT(std::abs(particle.flagMcMatchGen()), DecayChannelMain::BsToDsPi)) {
576576
return;
577577
}
578578
auto ptParticle = particle.ptTrack();

‎PWGHF/DataModel/CandidateReconstructionTables.h‎

Lines changed: 14 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -876,7 +876,7 @@ namespace hf_cand_bplus
876876
DECLARE_SOA_INDEX_COLUMN_FULL(Prong0, prong0, int, HfCand2Prong, "_0"); // D0 index
877877
// MC matching result:
878878
DECLARE_SOA_COLUMN(FlagMcMatchRec, flagMcMatchRec, int8_t); // reconstruction level
879-
DECLARE_SOA_COLUMN(ChannelMcMatchRec, channelMcMatchRec, int8_t); // reconstruction level
879+
DECLARE_SOA_COLUMN(FlagMcDecayChanRec, flagMcDecayChanRec, int8_t); // reconstruction level
880880
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); // reconstruction level
881881
DECLARE_SOA_COLUMN(FlagMcMatchGen, flagMcMatchGen, int8_t); // generator level
882882
DECLARE_SOA_COLUMN(OriginMcRec, originMcRec, int8_t); // particle origin, reconstruction level
@@ -2107,11 +2107,12 @@ namespace hf_cand_b0
21072107
DECLARE_SOA_INDEX_COLUMN_FULL(Prong0, prong0, int, HfCand3Prong, "_0"); // D index
21082108
// MC matching result:
21092109
DECLARE_SOA_COLUMN(FlagMcMatchRec, flagMcMatchRec, int8_t); // reconstruction level
2110-
DECLARE_SOA_COLUMN(ChannelMcMatchRec, channelMcMatchRec, int8_t); // reconstruction level
2111-
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); // reconstruction level
2110+
DECLARE_SOA_COLUMN(FlagMcDecayChanRec, flagMcDecayChanRec, int8_t); // reconstruction level
21122111
DECLARE_SOA_COLUMN(FlagMcMatchGen, flagMcMatchGen, int8_t); // generator level
2112+
DECLARE_SOA_COLUMN(FlagMcDecayChanGen, flagMcDecayChanGen, int8_t); // generator level
21132113
DECLARE_SOA_COLUMN(OriginMcRec, originMcRec, int8_t); // particle origin, reconstruction level
21142114
DECLARE_SOA_COLUMN(OriginMcGen, originMcGen, int8_t); // particle origin, generator level
2115+
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); // reconstruction level
21152116
DECLARE_SOA_COLUMN(DebugMcRec, debugMcRec, int8_t); // debug flag for mis-association reconstruction level
21162117

21172118
// mapping of decay types
@@ -2171,24 +2172,27 @@ using HfCandB0 = soa::Join<HfCandB0Ext, HfCandB0Prongs>;
21712172
DECLARE_SOA_TABLE(HfCandB0McRec, "AOD", "HFCANDB0MCREC",
21722173
hf_cand_b0::FlagMcMatchRec,
21732174
hf_cand_b0::OriginMcRec,
2175+
hf_cand_b0::FlagMcDecayChanRec,
21742176
hf_cand_b0::DebugMcRec);
21752177

21762178
// table with results of generator level MC matching
21772179
DECLARE_SOA_TABLE(HfCandB0McGen, "AOD", "HFCANDB0MCGEN",
21782180
hf_cand_b0::FlagMcMatchGen,
2179-
hf_cand_b0::OriginMcGen);
2181+
hf_cand_b0::OriginMcGen,
2182+
hf_cand_b0::FlagMcDecayChanGen);
21802183

21812184
// specific Bs candidate properties
21822185
namespace hf_cand_bs
21832186
{
21842187
DECLARE_SOA_INDEX_COLUMN_FULL(Prong0, prong0, int, HfCand3Prong, "_0"); // Ds index
21852188
// MC matching result:
21862189
DECLARE_SOA_COLUMN(FlagMcMatchRec, flagMcMatchRec, int8_t); // reconstruction level
2187-
DECLARE_SOA_COLUMN(ChannelMcMatchRec, channelMcMatchRec, int8_t); // reconstruction level
2188-
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); // reconstruction level
2190+
DECLARE_SOA_COLUMN(FlagMcDecayChanRec, flagMcDecayChanRec, int8_t); // reconstruction level
21892191
DECLARE_SOA_COLUMN(FlagMcMatchGen, flagMcMatchGen, int8_t); // generator level
2192+
DECLARE_SOA_COLUMN(FlagMcDecayChanGen, flagMcDecayChanGen, int8_t); // reconstruction level
21902193
DECLARE_SOA_COLUMN(OriginMcRec, originMcRec, int8_t); // particle origin, reconstruction level
21912194
DECLARE_SOA_COLUMN(OriginMcGen, originMcGen, int8_t); // particle origin, generator level
2195+
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); // reconstruction level
21922196
DECLARE_SOA_COLUMN(DebugMcRec, debugMcRec, int8_t); // debug flag for mis-association reconstruction level
21932197
DECLARE_SOA_DYNAMIC_COLUMN(ImpactParameterProduct, impactParameterProduct, // Impact parameter product for Bs -> J/Psi phi
21942198
[](float pxJpsiDauPos, float pyJpsiDauPos, float pzJpsiDauPos, float pxJpsiDauNeg, float pyJpsiDauNeg, float pzJpsiDauNeg, float pxLfTrack0, float pyLfTrack0, float pzLfTrack0, float pxLfTrack1, float pyLfTrack1, float pzLfTrack1, float xVtxP, float yVtxP, float zVtxP, float xVtxS, float yVtxS, float zVtxS) -> float {
@@ -2266,11 +2270,13 @@ using HfCandBs = soa::Join<HfCandBsExt, HfCandBsProngs>;
22662270

22672271
// table with results of reconstruction level MC matching
22682272
DECLARE_SOA_TABLE(HfCandBsMcRec, "AOD", "HFCANDBSMCREC",
2269-
hf_cand_bs::FlagMcMatchRec);
2273+
hf_cand_bs::FlagMcMatchRec,
2274+
hf_cand_bs::FlagMcDecayChanRec);
22702275

22712276
// table with results of generator level MC matching
22722277
DECLARE_SOA_TABLE(HfCandBsMcGen, "AOD", "HFCANDBSMCGEN",
2273-
hf_cand_bs::FlagMcMatchGen);
2278+
hf_cand_bs::FlagMcMatchGen,
2279+
hf_cand_bs::FlagMcDecayChanGen);
22742280

22752281
namespace hf_cand_4prong
22762282
{

‎PWGHF/TableProducer/candidateCreatorB0.cxx‎

Lines changed: 9 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -406,13 +406,14 @@ struct HfCandidateCreatorB0Expressions {
406406

407407
int indexRec = -1;
408408
int8_t sign = 0;
409-
int8_t flag = 0;
409+
int8_t flagChannelMain = 0;
410+
int8_t flagChannelReso = 0;
410411
int8_t origin = 0;
411412
int8_t debug = 0;
412413

413414
// Match reconstructed candidates.
414415
for (const auto& candidate : candsB0) {
415-
flag = 0;
416+
flagChannelMain = 0;
416417
origin = 0;
417418
debug = 0;
418419
auto candD = candidate.prong0();
@@ -430,28 +431,28 @@ struct HfCandidateCreatorB0Expressions {
430431
// D- → π- K+ π-
431432
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughtersD, Pdg::kDMinus, std::array{-kPiPlus, +kKPlus, -kPiPlus}, true, &sign, 2);
432433
if (indexRec > -1) {
433-
flag = sign * o2::hf_decay::hf_cand_beauty::DecayChannelMain::B0ToDminusPi;
434+
flagChannelMain = sign * o2::hf_decay::hf_cand_beauty::DecayChannelMain::B0ToDminusPi;
434435
} else {
435436
debug = 1;
436437
LOGF(debug, "WARNING: B0 in decays in the expected final state but the condition on the intermediate state is not fulfilled");
437438
}
438439
}
439440

440441
// B0 → Ds- π+ → (K- K+ π-) π+
441-
if (!flag) {
442+
if (!flagChannelMain) {
442443
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughtersB0, Pdg::kB0, std::array{-kKPlus, +kKPlus, -kPiPlus, +kPiPlus}, true, &sign, 3);
443444
if (indexRec > -1) {
444445
// Ds- → K- K+ π-
445446
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughtersD, -Pdg::kDS, std::array{-kKPlus, +kKPlus, -kPiPlus}, true, &sign, 2);
446447
if (indexRec > -1) {
447-
flag = sign * o2::hf_decay::hf_cand_beauty::DecayChannelMain::B0ToDsPi;
448+
flagChannelMain = sign * o2::hf_decay::hf_cand_beauty::DecayChannelMain::B0ToDsPi;
448449
}
449450
}
450451
}
451452

452453
// Partly reconstructed decays, i.e. the 4 prongs have a common b-hadron ancestor
453454
// convention: final state particles are prong0,1,2,3
454-
if (!flag) {
455+
if (!flagChannelMain) {
455456
auto particleProng0 = arrayDaughtersB0[0].mcParticle();
456457
auto particleProng1 = arrayDaughtersB0[1].mcParticle();
457458
auto particleProng2 = arrayDaughtersB0[2].mcParticle();
@@ -467,13 +468,13 @@ struct HfCandidateCreatorB0Expressions {
467468

468469
// look for common b-hadron ancestor
469470
if (index0Mother > -1 && index0Mother == index1Mother && index1Mother == index2Mother && index2Mother == index3Mother) {
470-
flag = hf_cand_b0::DecayTypeMc::PartlyRecoDecay; // FIXME
471+
flagChannelMain = hf_cand_b0::DecayTypeMc::PartlyRecoDecay; // FIXME
471472
break;
472473
}
473474
}
474475
}
475476

476-
rowMcMatchRec(flag, origin, debug);
477+
rowMcMatchRec(flagChannelMain, origin, debug);
477478
} // rec
478479

479480
hf_mc_gen::fillMcMatchGenB0(mcParticles, rowMcMatchGen); // gen

‎PWGHF/TableProducer/candidateCreatorBs.cxx‎

Lines changed: 14 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -391,14 +391,15 @@ struct HfCandidateCreatorBsExpressions {
391391
{
392392
int indexRec = -1;
393393
int8_t sign = 0;
394-
int8_t flag = 0;
394+
int8_t flagChannelMain = 0;
395+
int8_t flagChannelReso = 0;
395396
std::vector<int> arrDaughDsIndex;
396397
std::array<int, 2> arrPDGDaughDs;
397398
std::array<int, 2> arrPDGResonantDsPhiPi = {Pdg::kPhi, kPiPlus}; // Ds± → Phi π±
398399

399400
// Match reconstructed candidates.
400401
for (const auto& candidate : candsBs) {
401-
flag = 0;
402+
flagChannelMain = 0;
402403
arrDaughDsIndex.clear();
403404
auto candDs = candidate.prong0();
404405
auto arrayDaughtersBs = std::array{candDs.prong0_as<aod::TracksWMc>(),
@@ -422,13 +423,13 @@ struct HfCandidateCreatorBsExpressions {
422423
arrPDGDaughDs[iProng] = std::abs(daughI.pdgCode());
423424
}
424425
if ((arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[0] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[1]) || (arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[1] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[0])) {
425-
flag = sign * BIT(hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi);
426+
flagChannelMain = sign * DecayChannelMain::BsToDsPi;
426427
}
427428
}
428429
}
429430
}
430431

431-
if (!flag) {
432+
if (!flagChannelMain) {
432433
// Checking B0(bar) → Ds± π∓ → (K- K+ π±) π∓
433434
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughtersBs, Pdg::kB0, std::array{-kKPlus, +kKPlus, +kPiPlus, -kPiPlus}, true, &sign, 3);
434435
if (indexRec > -1) {
@@ -442,7 +443,7 @@ struct HfCandidateCreatorBsExpressions {
442443
arrPDGDaughDs[iProng] = std::abs(daughI.pdgCode());
443444
}
444445
if ((arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[0] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[1]) || (arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[1] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[0])) {
445-
flag = sign * BIT(hf_cand_bs::DecayTypeMc::B0ToDsPiToPhiPiPiToKKPiPi);
446+
flagChannelMain = sign * BIT(hf_cand_bs::DecayTypeMc::B0ToDsPiToPhiPiPiToKKPiPi);
446447
}
447448
}
448449
}
@@ -451,7 +452,7 @@ struct HfCandidateCreatorBsExpressions {
451452

452453
// Partly reconstructed decays, i.e. the 4 prongs have a common b-hadron ancestor
453454
// convention: final state particles are prong0,1,2,3
454-
if (!flag) {
455+
if (!flagChannelMain) {
455456
auto particleProng0 = arrayDaughtersBs[0].mcParticle();
456457
auto particleProng1 = arrayDaughtersBs[1].mcParticle();
457458
auto particleProng2 = arrayDaughtersBs[2].mcParticle();
@@ -468,19 +469,19 @@ struct HfCandidateCreatorBsExpressions {
468469
// look for common b-hadron ancestor
469470
if (index0Mother > -1 && index1Mother > -1 && index2Mother > -1 && index3Mother > -1) {
470471
if (index0Mother == index1Mother && index1Mother == index2Mother && index2Mother == index3Mother) {
471-
flag = BIT(hf_cand_bs::DecayTypeMc::PartlyRecoDecay); // FIXME
472+
flagChannelMain = BIT(hf_cand_bs::DecayTypeMc::PartlyRecoDecay); // FIXME
472473
break;
473474
}
474475
}
475476
}
476477
}
477478

478-
rowMcMatchRec(flag);
479+
rowMcMatchRec(flagChannelMain);
479480
} // rec
480481

481482
// Match generated particles.
482483
for (const auto& particle : mcParticles) {
483-
flag = 0;
484+
flagChannelMain = 0;
484485
arrDaughDsIndex.clear();
485486

486487
// Checking Bs0(bar) → Ds∓ π± → (K- K+ π∓) π±
@@ -495,13 +496,13 @@ struct HfCandidateCreatorBsExpressions {
495496
arrPDGDaughDs[jProng] = std::abs(daughJ.pdgCode());
496497
}
497498
if ((arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[0] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[1]) || (arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[1] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[0])) {
498-
flag = sign * BIT(hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi);
499+
flag = sign * DecayChannelMain::BsToDsPi;
499500
}
500501
}
501502
}
502503
}
503504

504-
if (!flag) {
505+
if (!flagChannelMain) {
505506
// Checking B0(bar) → Ds± π∓ → (K- K+ π±) π∓
506507
if (RecoDecay::isMatchedMCGen(mcParticles, particle, Pdg::kB0, std::array{+Pdg::kDS, -kPiPlus}, true)) {
507508
// Checking Ds± → K- K+ π±
@@ -514,14 +515,14 @@ struct HfCandidateCreatorBsExpressions {
514515
arrPDGDaughDs[jProng] = std::abs(daughJ.pdgCode());
515516
}
516517
if ((arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[0] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[1]) || (arrPDGDaughDs[0] == arrPDGResonantDsPhiPi[1] && arrPDGDaughDs[1] == arrPDGResonantDsPhiPi[0])) {
517-
flag = sign * BIT(hf_cand_bs::DecayTypeMc::B0ToDsPiToPhiPiPiToKKPiPi);
518+
flagChannelMain = sign * BIT(hf_cand_bs::DecayTypeMc::B0ToDsPiToPhiPiPiToKKPiPi);
518519
}
519520
}
520521
}
521522
}
522523
}
523524

524-
rowMcMatchGen(flag);
525+
rowMcMatchGen(flagChannelMain);
525526
} // gen
526527
} // processMc
527528
PROCESS_SWITCH(HfCandidateCreatorBsExpressions, processMc, "Process MC", false);

‎PWGHF/TableProducer/treeCreatorBsToDsPi.cxx‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -184,8 +184,8 @@ struct HfTreeCreatorBsToDsPi {
184184

185185
Filter filterSelectCandidates = aod::hf_sel_candidate_bs::isSelBsToDsPi >= selectionFlagBs;
186186

187-
Partition<SelectedCandidatesMc> recSig = nabs(aod::hf_cand_bs::flagMcMatchRec) == (int8_t)BIT(aod::hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi);
188-
Partition<SelectedCandidatesMc> recBg = nabs(aod::hf_cand_bs::flagMcMatchRec) != (int8_t)BIT(aod::hf_cand_bs::DecayTypeMc::BsToDsPiToPhiPiPiToKKPiPi);
187+
Partition<SelectedCandidatesMc> recSig = nabs(aod::hf_cand_bs::flagMcMatchRec) == (int8_t)DecayChannelMain::BsToDsPi;
188+
Partition<SelectedCandidatesMc> recBg = nabs(aod::hf_cand_bs::flagMcMatchRec) != (int8_t)DecayChannelMain::BsToDsPi;
189189

190190
void init(InitContext const&)
191191
{

0 commit comments

Comments
 (0)