Matching of D*->D0pi in resonances data Creator for backgoround studies

Author: Luca610

PWGHF/Core/DecayChannels.h

@@ -250,7 +250,10 @@ enum DecayChannelMain : int8_t {
   // Xic(3080)0
   Xic3080zeroToD0Lambda = 14, // D0 Λ
   // Xic(3080)+
-  Xic3080plusToDplusLambda = 15 // D+ Λ
+  Xic3080plusToDplusLambda = 15, // D+ Λ
+  // D*+
+  DstarToD0Pi = 16, // D0 π+
+  NChannelsMain = DstarToD0Pi // last channel
 };
 } // namespace hf_cand_reso
 } // namespace o2::hf_decay

PWGHF/D2H/TableProducer/dataCreatorCharmResoReduced.cxx

@@ -350,7 +350,8 @@ struct HfDataCreatorCharmResoReduced {
         doprocessDstarV0MC || doprocessDstarTrackMC || doprocessDstarV0AndTrackMC || doprocessDstarV0MCWithMl || doprocessDstarTrackMCWithMl || doprocessDstarV0AndTrackMCWithMl ||
         doprocessDplusV0MC || doprocessDplusTrackMC || doprocessDplusV0AndTrackMC || doprocessDplusV0MCWithMl || doprocessDplusTrackMCWithMl || doprocessDplusV0AndTrackMCWithMl) {
       // MC Rec
-      registry.add("hMCRecCounter", "Number of Reconstructed MC Matched candidates per channel", {HistType::kTH1D, {{31, -15.5, 15.5}}});
+      int nChannels = hf_decay::hf_cand_reso::DecayChannelMain::NChannelsMain;
+      registry.add("hMCRecCounter", "Number of Reconstructed MC Matched candidates per channel", {HistType::kTH1D, {{2*nChannels + 1, -(nChannels + 0.5), nChannels + 0.5}}});
       registry.add("hMCRecDebug", "Debug of MC Reco", {HistType::kTH1D, {{551, -0.5, 550.5}}});
       registry.add("hMCRecOrigin", "Origin of Matched particles", {HistType::kTH1D, {{3, -0.5, 2.5}}});
       registry.add("hMCRecMassGen", "Generated inv. mass of resoncances", {HistType::kTH1D, {{2000, 1.8, 3.8}}});
@@ -931,6 +932,9 @@ struct HfDataCreatorCharmResoReduced {
         origin = candCharmBach.originMcRec();
       }
       flagTrack = getMatchingFlagTrack(bachelorTrack);
+      if (std::abs(flagTrack) > 0) {
+        SETBIT(debugMcRec, flagTrack);
+      }
       if (hf_decay::hf_cand_2prong::daughtersD0Main.contains(static_cast<hf_decay::hf_cand_2prong::DecayChannelMain>(std::abs(flagCharmBach))) && flagTrack == hf_decay::hf_cand_reso::PartialMatchMc::PionMatched) {
         auto arrDaughtersReso = std::array{vecDaughtersReso[0], vecDaughtersReso[1], bachelorTrack};
         auto pdgCodesDzeroDaughters = hf_decay::hf_cand_2prong::daughtersD0Main.at(static_cast<hf_decay::hf_cand_2prong::DecayChannelMain>(std::abs(flagCharmBach)));

PWGHF/Utils/utilsMcMatching.h

@@ -198,7 +198,8 @@ const std::unordered_map<int, int> particlesToDstarPi = {
 const std::unordered_map<int, int> particlesToDplusPi = {
   {DecayChannelMain::D2starzeroToDplusPi, constants::physics::Pdg::kD2Star0}};
 const std::unordered_map<int, int> particlesToD0Pi = {
-  {DecayChannelMain::D2starplusToD0Pi, constants::physics::Pdg::kD2StarPlus}};
+  {DecayChannelMain::D2starplusToD0Pi, constants::physics::Pdg::kD2StarPlus},
+  {DecayChannelMain::DstarToD0Pi, constants::physics::Pdg::kDStar}};
 const std::unordered_map<int, int> particlesToD0Kplus = {
   {DecayChannelMain::Ds2starToD0Kplus, constants::physics::Pdg::kDS2Star}};