@@ -476,7 +476,7 @@ struct HfDataCreatorCharmHadPiReduced {
476476 // B+ → D0(bar) π+ → (K+ π-) π+
477477 auto indexRec = RecoDecay::getMatchedMCRec (particlesMc, std::array{vecDaughtersB[0 ], vecDaughtersB[1 ], vecDaughtersB[2 ]}, Pdg::kBPlus , std::array{+kPiPlus , +kKPlus , -kPiPlus }, true , &sign, 2 );
478478 if (indexRec > -1 ) {
479- // Printf("Checking D0bar → K+ π-") ;
479+ // D0(bar) → K+ π-;
480480 indexRec = RecoDecay::getMatchedMCRec (particlesMc, std::array{vecDaughtersB[0 ], vecDaughtersB[1 ]}, Pdg::kD0 , std::array{+kPiPlus , -kKPlus }, true , &sign, 1 );
481481 if (indexRec > -1 ) {
482482 flag = sign * BIT (hf_cand_bplus::DecayType::BplusToD0Pi);
@@ -492,8 +492,6 @@ struct HfDataCreatorCharmHadPiReduced {
492492 checkWrongCollision (particleMother, collision, indexCollisionMaxNumContrib, flagWrongCollision);
493493 }
494494 }
495- rowHfD0PiMcRecReduced (indexHfCandCharm, selectedTracksPion[vecDaughtersB.back ().globalIndex ()], flag, flagWrongCollision, debug, motherPt);
496-
497495 // additional checks for correlated backgrounds
498496 if (checkDecayTypeMc) {
499497 // Partly reconstructed decays, i.e. the 3 prongs have a common b-hadron ancestor
@@ -526,6 +524,7 @@ struct HfDataCreatorCharmHadPiReduced {
526524 }
527525 rowHfD0PiMcCheckReduced (pdgCodeBeautyMother, pdgCodeProng0, pdgCodeProng1, pdgCodeProng2);
528526 }
527+ rowHfD0PiMcRecReduced (indexHfCandCharm, selectedTracksPion[vecDaughtersB.back ().globalIndex ()], flag, flagWrongCollision, debug, motherPt);
529528 }
530529 }
531530
0 commit comments