[PWGHF] Refactor of resonances workflow, pt 2

PWGHF/Core/DecayChannels.h

@@ -219,6 +219,40 @@ enum DecayChannelToJpsiResonant : int8_t {
  NChannelsToJpsiResonant = BsToJpsiPhi // last channel
 };
 } // namespace hf_cand_beauty
+
+namespace hf_cand_reso
+{
+/// @brief resonance candidates: main channels
+enum DecayChannelMain : int8_t {
+ // D1(2420)0
+ D1zeroToDstarPi = 1, // D*+ π-
+ // D2*(2460)0
+ D2starzeroToDplusPi = 2, // D+ π−
+ D2starzeroToDstarPi = 3, // D*+ π-
+ // D2*(2460)+
+ D2starplusToD0Pi = 4, // D0 π+
+ // Ds1(2536)+
+ Ds1ToDstarK0s = 5, // D*+ K0s
+ // Ds2*(2573)+
+ Ds2starToD0Kplus = 6, // D0 K+
+ Ds2starToDplusK0s = 7, // D+ K0s
+ Ds2starToDstarK0s = 8, // D*+ K0s
+ // Ds1*(2700)+
+ Ds1star2700ToDstarK0s = 9, // D*+ K0s
+ // Ds1*(2860)+
+ Ds1star2860ToDstarK0s = 10, // D*+ K0s
+ // Ds3*(2860)+
+ Ds3star2860ToDstarK0s = 11, // D*+ K0s
+ // Xic(3055)0
+ Xic3055zeroToD0Lambda = 12, // D0 Λ
+ // Xic(3055)+
+ Xic3055plusToDplusLambda = 13, // D+ Λ
+ // Xic(3080)0
+ Xic3080zeroToD0Lambda = 14, // D0 Λ
+ // Xic(3080)+
+ Xic3080plusToDplusLambda = 15 // D+ Λ
+};
+} // namespace hf_cand_reso
 } // namespace o2::hf_decay
 
 #endif // PWGHF_CORE_DECAYCHANNELS_H_

PWGHF/D2H/Core/SelectorCutsRedDataFormat.h

@@ -29,7 +29,7 @@ namespace hf_cuts_d_daughter
 static constexpr int NBinsPt = 7;
 static constexpr int NCutVars = 6;
 constexpr double BinsPt[NBinsPt + 1] = {
-1.,
+0.,
  2.,
  4.,
  6.,
@@ -87,5 +87,15 @@ static const std::vector<std::string> labelsPt{};
 // column labels
 static const std::vector<std::string> labelsCutVar = {"invMassLow", "invMassHigh", "cpaMin", "dcaMax", "radiusMin"};
 } // namespace hf_cuts_v0_daughter
+
+namespace hf_cuts_track_daughter
+{
+static constexpr int NCutVars = 7;
+// default values for the cuts
+constexpr double Cuts[1][NCutVars] = {{0.1, 3, 40, 4, 3, -1, -1}}; // nSigmaTpc, nSigmaTof, nSigmaCombined
+// row labels
+static const std::vector<std::string> labelsCutVar = {"ptMin", "itsNClsMin", "tpcNCrossedRowsMin", "tpcChi2Max", "nSigmaTpc", "nSigmaTof", "nSigmaComb"};
+} // namespace hf_cuts_track_daughter
+
 } // namespace o2::analysis
 #endif // PWGHF_D2H_CORE_SELECTORCUTSREDDATAFORMAT_H_

PWGHF/D2H/Tasks/CMakeLists.txt

@@ -1,4 +1,4 @@
 # Copyright 2019-2020 CERN and copyright holders of ALICE O2.
 # See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
 # All rights not expressly granted are reserved.
 #
@@ -54,8 +54,13 @@
  PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
  COMPONENT_NAME Analysis)
 
-o2physics_add_dpl_workflow(task-charm-reso-reduced
- SOURCES taskCharmResoReduced.cxx
+o2physics_add_dpl_workflow(task-charm-reso-to-d-v0-reduced
+ SOURCES taskCharmResoToDV0Reduced.cxx
+ PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
+ COMPONENT_NAME Analysis)
+
+o2physics_add_dpl_workflow(task-charm-reso-to-d-trk-reduced
+ SOURCES taskCharmResoToDTrkReduced.cxx
  PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
  COMPONENT_NAME Analysis)
 
@@ -109,7 +114,7 @@
  PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
  COMPONENT_NAME Analysis)
 
 o2physics_add_dpl_workflow(task-omegac0-to-omega-pi
  SOURCES taskOmegac0ToOmegapi.cxx
  PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
  COMPONENT_NAME Analysis)

PWGHF/Utils/utilsMcMatching.h

@@ -176,6 +176,45 @@ inline std::unordered_map<DecayChannelMain, const std::vector<int>> getDecayChan
 }
 } // namespace hf_cand_3prong
 
+namespace hf_cand_reso
+{
+const std::unordered_map<int, int> particlesToDstarK0s = {
+ {DecayChannelMain::Ds1ToDstarK0s, constants::physics::Pdg::kDS1},
+ {DecayChannelMain::Ds2starToDstarK0s, constants::physics::Pdg::kDS2Star},
+ {DecayChannelMain::Ds1star2700ToDstarK0s, constants::physics::Pdg::kDS1Star2700},
+ {DecayChannelMain::Ds1star2860ToDstarK0s, constants::physics::Pdg::kDS1Star2860},
+ {DecayChannelMain::Ds3star2860ToDstarK0s, constants::physics::Pdg::kDS3Star2860}};
+const std::unordered_map<int, int> particlesToDplusK0s = {
+ {DecayChannelMain::Ds2starToDplusK0s, constants::physics::Pdg::kDS2Star}};
+const std::unordered_map<int, int> particlesToDplusLambda = {
+ {DecayChannelMain::Xic3055plusToDplusLambda, constants::physics::Pdg::kXiC3055Plus},
+ {DecayChannelMain::Xic3080plusToDplusLambda, constants::physics::Pdg::kXiC3080Plus}};
+const std::unordered_map<int, int> particlesToD0Lambda = {
+ {DecayChannelMain::Xic3055zeroToD0Lambda, constants::physics::Pdg::kXiC3055_0},
+ {DecayChannelMain::Xic3080zeroToD0Lambda, constants::physics::Pdg::kXiC3080_0}};
+const std::unordered_map<int, int> particlesToDstarPi = {
+ {DecayChannelMain::D1zeroToDstarPi, constants::physics::Pdg::kD10},
+ {DecayChannelMain::D2starzeroToDstarPi, constants::physics::Pdg::kD2Star0}};
+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}};
+const std::unordered_map<int, int> particlesToD0Kplus = {
+ {DecayChannelMain::Ds2starToD0Kplus, constants::physics::Pdg::kDS2Star}};
+
+enum PartialMatchMc : uint8_t {
+ D0Matched = 0,
+ DstarMatched,
+ DplusMatched,
+ K0Matched,
+ LambdaMatched,
+ PionMatched,
+ KaonMatched,
+ ProtonMatched,
+ ResoPartlyMatched
+};
+} // namespace hf_cand_reso
+
 /// Compare an array of PDG codes with an expected array
 /// \tparam N size of the arrays to be compared
 /// \param arrPdgTested array of PDG codes to be tested