Skip to content

Commit 273539c

Browse files
committed
tmp
1 parent b828745 commit 273539c

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
@@ -2183,11 +2183,12 @@ namespace hf_cand_b0
21832183
DECLARE_SOA_INDEX_COLUMN_FULL(Prong0, prong0, int, HfCand3Prong, "_0"); // D index
21842184
// MC matching result:
21852185
DECLARE_SOA_COLUMN(FlagMcMatchRec, flagMcMatchRec, int8_t); // reconstruction level
2186-
DECLARE_SOA_COLUMN(ChannelMcMatchRec, channelMcMatchRec, int8_t); // reconstruction level
2187-
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); // reconstruction level
2186+
DECLARE_SOA_COLUMN(FlagMcDecayChanRec, flagMcDecayChanRec, int8_t); // reconstruction level
21882187
DECLARE_SOA_COLUMN(FlagMcMatchGen, flagMcMatchGen, int8_t); // generator level
2188+
DECLARE_SOA_COLUMN(FlagMcDecayChanGen, flagMcDecayChanGen, int8_t); // generator level
21892189
DECLARE_SOA_COLUMN(OriginMcRec, originMcRec, int8_t); // particle origin, reconstruction level
21902190
DECLARE_SOA_COLUMN(OriginMcGen, originMcGen, int8_t); // particle origin, generator level
2191+
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); // reconstruction level
21912192
DECLARE_SOA_COLUMN(DebugMcRec, debugMcRec, int8_t); // debug flag for mis-association reconstruction level
21922193

21932194
// mapping of decay types
@@ -2247,24 +2248,27 @@ using HfCandB0 = soa::Join<HfCandB0Ext, HfCandB0Prongs>;
22472248
DECLARE_SOA_TABLE(HfCandB0McRec, "AOD", "HFCANDB0MCREC",
22482249
hf_cand_b0::FlagMcMatchRec,
22492250
hf_cand_b0::OriginMcRec,
2251+
hf_cand_b0::FlagMcDecayChanRec,
22502252
hf_cand_b0::DebugMcRec);
22512253

22522254
// table with results of generator level MC matching
22532255
DECLARE_SOA_TABLE(HfCandB0McGen, "AOD", "HFCANDB0MCGEN",
22542256
hf_cand_b0::FlagMcMatchGen,
2255-
hf_cand_b0::OriginMcGen);
2257+
hf_cand_b0::OriginMcGen,
2258+
hf_cand_b0::FlagMcDecayChanGen);
22562259

22572260
// specific Bs candidate properties
22582261
namespace hf_cand_bs
22592262
{
22602263
DECLARE_SOA_INDEX_COLUMN_FULL(Prong0, prong0, int, HfCand3Prong, "_0"); // Ds index
22612264
// MC matching result:
22622265
DECLARE_SOA_COLUMN(FlagMcMatchRec, flagMcMatchRec, int8_t); // reconstruction level
2263-
DECLARE_SOA_COLUMN(ChannelMcMatchRec, channelMcMatchRec, int8_t); // reconstruction level
2264-
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); // reconstruction level
2266+
DECLARE_SOA_COLUMN(FlagMcDecayChanRec, flagMcDecayChanRec, int8_t); // reconstruction level
22652267
DECLARE_SOA_COLUMN(FlagMcMatchGen, flagMcMatchGen, int8_t); // generator level
2268+
DECLARE_SOA_COLUMN(FlagMcDecayChanGen, flagMcDecayChanGen, int8_t); // reconstruction level
22662269
DECLARE_SOA_COLUMN(OriginMcRec, originMcRec, int8_t); // particle origin, reconstruction level
22672270
DECLARE_SOA_COLUMN(OriginMcGen, originMcGen, int8_t); // particle origin, generator level
2271+
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); // reconstruction level
22682272
DECLARE_SOA_COLUMN(DebugMcRec, debugMcRec, int8_t); // debug flag for mis-association reconstruction level
22692273
DECLARE_SOA_DYNAMIC_COLUMN(ImpactParameterProduct, impactParameterProduct, // Impact parameter product for Bs -> J/Psi phi
22702274
[](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 {
@@ -2342,11 +2346,13 @@ using HfCandBs = soa::Join<HfCandBsExt, HfCandBsProngs>;
23422346

23432347
// table with results of reconstruction level MC matching
23442348
DECLARE_SOA_TABLE(HfCandBsMcRec, "AOD", "HFCANDBSMCREC",
2345-
hf_cand_bs::FlagMcMatchRec);
2349+
hf_cand_bs::FlagMcMatchRec,
2350+
hf_cand_bs::FlagMcDecayChanRec);
23462351

23472352
// table with results of generator level MC matching
23482353
DECLARE_SOA_TABLE(HfCandBsMcGen, "AOD", "HFCANDBSMCGEN",
2349-
hf_cand_bs::FlagMcMatchGen);
2354+
hf_cand_bs::FlagMcMatchGen,
2355+
hf_cand_bs::FlagMcDecayChanGen);
23502356

23512357
namespace hf_cand_4prong
23522358
{

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)