@@ -58,10 +58,13 @@ enum V0Type : uint8_t {
5858};
5959
6060enum DecayTypeMc : uint8_t {
61- Ds1ToDStarK0ToD0PiK0s = 0 ,
61+ Ds1ToDStarK0ToD0PiK0s = 1 ,
6262 Ds2StarToDplusK0sToPiKaPiPiPi,
6363 Ds1ToDStarK0ToDPlusPi0K0s,
64- Ds1ToDStarK0ToDPlusGammaK0s
64+ Ds1ToDStarK0ToD0PiK0sPart,
65+ Ds1ToDStarK0ToD0NoPiK0sPart,
66+ Ds1ToDStarK0ToD0PiK0sOneMu,
67+ Ds2StarToDplusK0sOneMu
6568};
6669
6770const int nBinsPt = 7 ;
@@ -676,13 +679,14 @@ struct HfCandidateCreatorCharmResoReducedExpressions {
676679 }
677680 rowResoMcRec (rowDV0McRec.flagMcMatchRec (), rowDV0McRec.debugMcRec (), rowDV0McRec.origin (), rowDV0McRec.ptMother ());
678681 filledMcInfo = true ;
679- if (TESTBIT (rowDV0McRec.flagMcMatchRec (), DecayTypeMc::Ds1ToDStarK0ToD0PiK0s) || TESTBIT (rowDV0McRec.flagMcMatchRec (), DecayTypeMc::Ds2StarToDplusK0sToPiKaPiPiPi)) {
682+ if (std::abs (rowDV0McRec.flagMcMatchRec ()) == DecayTypeMc::Ds1ToDStarK0ToD0PiK0s || std::abs (rowDV0McRec.flagMcMatchRec ()) == DecayTypeMc::Ds2StarToDplusK0sToPiKaPiPiPi ||
683+ std::abs (rowDV0McRec.flagMcMatchRec ()) == DecayTypeMc::Ds1ToDStarK0ToD0PiK0sOneMu || std::abs (rowDV0McRec.flagMcMatchRec ()) == DecayTypeMc::Ds2StarToDplusK0sOneMu) {
680684 registry.fill (HIST (" hMassMcMatched" ), candReso.invMass (), candReso.pt ());
681685 registry.fill (HIST (" hMassMcMatchedVsBach0Mass" ), candReso.invMass (), candReso.invMassProng0 () - candReso.invMassD0 ());
682686 registry.fill (HIST (" hMassMcMatchedVsBach1Mass" ), candReso.invMass (), candReso.invMassProng1 ());
683687 registry.fill (HIST (" hMassMcMatchedVsD0Mass" ), candReso.invMass (), candReso.invMassD0 ());
684688
685- } else if (TESTBIT (rowDV0McRec.flagMcMatchRec (), DecayTypeMc::Ds1ToDStarK0ToDPlusGammaK0s) || TESTBIT (rowDV0McRec.flagMcMatchRec (), DecayTypeMc::Ds1ToDStarK0ToDPlusPi0K0s) ) {
689+ } else if (std::abs (rowDV0McRec.flagMcMatchRec ()) == DecayTypeMc::Ds1ToDStarK0ToD0NoPiK0sPart || std::abs (rowDV0McRec.flagMcMatchRec ()) == DecayTypeMc::Ds1ToDStarK0ToDPlusPi0K0s) {
686690 registry.fill (HIST (" hMassMcMatchedIncomplete" ), candReso.invMass (), candReso.pt ());
687691 } else {
688692 registry.fill (HIST (" hMassMcUnmatched" ), candReso.invMass (), candReso.pt ());
0 commit comments