diff --git a/PWGHF/Utils/utilsMcMatching.h b/PWGHF/Utils/utilsMcMatching.h index 02522a40b4b..22d035c6059 100644 --- a/PWGHF/Utils/utilsMcMatching.h +++ b/PWGHF/Utils/utilsMcMatching.h @@ -33,122 +33,122 @@ namespace hf_chns_2prong using namespace o2::hf_decay::hf_cand_2prong; std::unordered_map> finalStates2Prongs = -{ - {DecayChannelMain::D0ToPiK, {+kKMinus, +kPiPlus}}, - {DecayChannelMain::D0ToPiKPi0, {+kKMinus, +kPiPlus, +kPi0}}, - {DecayChannelMain::D0ToPiPi, {+kPiMinus, +kPiPlus}}, - {DecayChannelMain::D0ToPiPiPi0, {+kPiMinus, +kPiPlus, +kPi0}}, - {DecayChannelMain::D0ToKK, {+kKMinus, +kKPlus}}, + { + {DecayChannelMain::D0ToPiK, {+kKMinus, +kPiPlus}}, + {DecayChannelMain::D0ToPiKPi0, {+kKMinus, +kPiPlus, +kPi0}}, + {DecayChannelMain::D0ToPiPi, {+kPiMinus, +kPiPlus}}, + {DecayChannelMain::D0ToPiPiPi0, {+kPiMinus, +kPiPlus, +kPi0}}, + {DecayChannelMain::D0ToKK, {+kKMinus, +kKPlus}}, }; std::unordered_map> resoStatesD0 = -{ - {DecayChannelResonant::D0ToRhoplusPi, {kRho770Plus, +kPiMinus}}, - {DecayChannelResonant::D0ToRhoplusK, {kRho770Plus, +kKMinus}}, - {DecayChannelResonant::D0ToKstar0Pi0, {-kK0Star892, +kPi0}}, - {DecayChannelResonant::D0ToKstarPi, {-kKPlusStar892, +kPiPlus}}, + { + {DecayChannelResonant::D0ToRhoplusPi, {kRho770Plus, +kPiMinus}}, + {DecayChannelResonant::D0ToRhoplusK, {kRho770Plus, +kKMinus}}, + {DecayChannelResonant::D0ToKstar0Pi0, {-kK0Star892, +kPi0}}, + {DecayChannelResonant::D0ToKstarPi, {-kKPlusStar892, +kPiPlus}}, }; -} // namespace 2prong - +} // namespace hf_chns_2prong + namespace hf_chns_3prong { using namespace o2::hf_decay::hf_cand_3prong; std::unordered_map> resoStatesDPlus = -{ - {DecayChannelResonant::DplusToPhiPi, {+kPhi, +kPiPlus}}, - {DecayChannelResonant::DplusToKstar0K, {-kK0Star892, +kKPlus}}, - {DecayChannelResonant::DplusToKstar1430_0K, {+10311, +kKPlus}}, - {DecayChannelResonant::DplusToRho0Pi, {+kRho770_0, +kPiPlus}}, - {DecayChannelResonant::DplusToF2_1270Pi, {+225, +kPiPlus}}, + { + {DecayChannelResonant::DplusToPhiPi, {+kPhi, +kPiPlus}}, + {DecayChannelResonant::DplusToKstar0K, {-kK0Star892, +kKPlus}}, + {DecayChannelResonant::DplusToKstar1430_0K, {+10311, +kKPlus}}, + {DecayChannelResonant::DplusToRho0Pi, {+kRho770_0, +kPiPlus}}, + {DecayChannelResonant::DplusToF2_1270Pi, {+225, +kPiPlus}}, }; std::unordered_map> finalStatesDPlus = -{ - {DecayChannelMain::DplusToPiKPi, {+kKMinus, +kKPlus, +kPiPlus}}, - {DecayChannelMain::DplusToPiKK, {+kKMinus, +kPiPlus, +kPiPlus}}, - {DecayChannelMain::DplusToPiKPiPi0, {+kKMinus, +kPiPlus, +kPiPlus, +kPi0}}, - {DecayChannelMain::DplusToPiPiPi, {+kPiMinus, +kPiPlus, +kPiPlus}}, + { + {DecayChannelMain::DplusToPiKPi, {+kKMinus, +kKPlus, +kPiPlus}}, + {DecayChannelMain::DplusToPiKK, {+kKMinus, +kPiPlus, +kPiPlus}}, + {DecayChannelMain::DplusToPiKPiPi0, {+kKMinus, +kPiPlus, +kPiPlus, +kPi0}}, + {DecayChannelMain::DplusToPiPiPi, {+kPiMinus, +kPiPlus, +kPiPlus}}, }; // Ds± → K± K∓ π± std::unordered_map> resoStatesDs = -{ - {DecayChannelResonant::DsToPhiPi, {+kPhi, +kPiPlus}}, - {DecayChannelResonant::DsToPhiRhoplus, {+kPhi, +kRho770Plus}}, - {DecayChannelResonant::DsToKstar0K, {-kK0Star892, +kKPlus}}, - {DecayChannelResonant::DsToKstar0Pi, {+kK0Star892, +kPiPlus}}, - {DecayChannelResonant::DsToRho0Pi, {+kRho770_0, +kPiPlus}}, - {DecayChannelResonant::DsToRho0K, {+kRho770_0, +kKPlus}}, - {DecayChannelResonant::DsToF2_1270Pi, {225, +kPiPlus}}, - {DecayChannelResonant::DsToF0_1370K, {10221, +kKPlus}}, - {DecayChannelResonant::DsToEtaPi, {221, +kPiPlus}}, + { + {DecayChannelResonant::DsToPhiPi, {+kPhi, +kPiPlus}}, + {DecayChannelResonant::DsToPhiRhoplus, {+kPhi, +kRho770Plus}}, + {DecayChannelResonant::DsToKstar0K, {-kK0Star892, +kKPlus}}, + {DecayChannelResonant::DsToKstar0Pi, {+kK0Star892, +kPiPlus}}, + {DecayChannelResonant::DsToRho0Pi, {+kRho770_0, +kPiPlus}}, + {DecayChannelResonant::DsToRho0K, {+kRho770_0, +kKPlus}}, + {DecayChannelResonant::DsToF2_1270Pi, {225, +kPiPlus}}, + {DecayChannelResonant::DsToF0_1370K, {10221, +kKPlus}}, + {DecayChannelResonant::DsToEtaPi, {221, +kPiPlus}}, }; std::unordered_map> finalStatesDs = -{ - {DecayChannelMain::DsToPiKK, {+kKMinus, +kKPlus, +kPiPlus}}, - {DecayChannelMain::DsToPiKKPi0, {+kKMinus, +kKPlus, +kPiPlus, +kPi0}}, - {DecayChannelMain::DsToPiPiK, {+kKPlus, +kPiPlus, +kPiMinus}}, - {DecayChannelMain::DsToPiPiPi, {+kPiMinus, +kPiPlus, +kPiPlus}}, - {DecayChannelMain::DsToPiPiPiPi0, {+kPiMinus, +kPiPlus, +kPiPlus, +kPi0}}, + { + {DecayChannelMain::DsToPiKK, {+kKMinus, +kKPlus, +kPiPlus}}, + {DecayChannelMain::DsToPiKKPi0, {+kKMinus, +kKPlus, +kPiPlus, +kPi0}}, + {DecayChannelMain::DsToPiPiK, {+kKPlus, +kPiPlus, +kPiMinus}}, + {DecayChannelMain::DsToPiPiPi, {+kPiMinus, +kPiPlus, +kPiPlus}}, + {DecayChannelMain::DsToPiPiPiPi0, {+kPiMinus, +kPiPlus, +kPiPlus, +kPi0}}, }; // Dstar → K± K∓ π± std::unordered_map> resoStatesDstar = -{ - {DecayChannelResonant::DstarToD0ToRhoplusPi, {kRho770Plus, +kPiMinus}}, - {DecayChannelResonant::DstarToD0ToRhoplusK, {kRho770Plus, +kKMinus}}, - {DecayChannelResonant::DstarToD0ToKstar0Pi0, {-kK0Star892, +kPi0}}, - {DecayChannelResonant::DstarToD0ToKstarPi, {-kKPlusStar892, +kPiPlus}}, - {DecayChannelResonant::DstarToDplusToPhiPi, {+kPhi, +kPiPlus}}, - {DecayChannelResonant::DstarToDplusToKstar0K, {-kK0Star892, +kKPlus}}, - {DecayChannelResonant::DstarToDplusToKstar1430_0K, {+10311, +kKPlus}}, - {DecayChannelResonant::DstarToDplusToRho0Pi, {+kRho770_0, +kPiPlus}}, - {DecayChannelResonant::DstarToDplusToF2_1270Pi, {+225, +kPiPlus}}, + { + {DecayChannelResonant::DstarToD0ToRhoplusPi, {kRho770Plus, +kPiMinus}}, + {DecayChannelResonant::DstarToD0ToRhoplusK, {kRho770Plus, +kKMinus}}, + {DecayChannelResonant::DstarToD0ToKstar0Pi0, {-kK0Star892, +kPi0}}, + {DecayChannelResonant::DstarToD0ToKstarPi, {-kKPlusStar892, +kPiPlus}}, + {DecayChannelResonant::DstarToDplusToPhiPi, {+kPhi, +kPiPlus}}, + {DecayChannelResonant::DstarToDplusToKstar0K, {-kK0Star892, +kKPlus}}, + {DecayChannelResonant::DstarToDplusToKstar1430_0K, {+10311, +kKPlus}}, + {DecayChannelResonant::DstarToDplusToRho0Pi, {+kRho770_0, +kPiPlus}}, + {DecayChannelResonant::DstarToDplusToF2_1270Pi, {+225, +kPiPlus}}, }; std::unordered_map> finalStatesDstar = -{ - {DecayChannelMain::DstarToPiKPi, {+kKMinus, +kPiPlus, +kPiPlus}}, - {DecayChannelMain::DstarToPiKPiPi0, {+kKMinus, +kPiPlus, +kPiPlus, +kPi0}}, - {DecayChannelMain::DstarToPiKPiPi0Pi0, {+kKMinus, +kPiPlus, +kPiPlus, +kPi0, +kPi0}}, - {DecayChannelMain::DstarToPiKK, {+kKMinus, +kKPlus, +kPiPlus}}, - {DecayChannelMain::DstarToPiKKPi0, {+kKMinus, +kKPlus, +kPiPlus, +kPi0}}, - {DecayChannelMain::DstarToPiPiPi, {+kPiMinus, +kPiPlus, +kPiPlus}}, - {DecayChannelMain::DstarToPiPiPiPi0, {+kPiMinus, +kPiPlus, +kPiPlus, +kPi0}}, + { + {DecayChannelMain::DstarToPiKPi, {+kKMinus, +kPiPlus, +kPiPlus}}, + {DecayChannelMain::DstarToPiKPiPi0, {+kKMinus, +kPiPlus, +kPiPlus, +kPi0}}, + {DecayChannelMain::DstarToPiKPiPi0Pi0, {+kKMinus, +kPiPlus, +kPiPlus, +kPi0, +kPi0}}, + {DecayChannelMain::DstarToPiKK, {+kKMinus, +kKPlus, +kPiPlus}}, + {DecayChannelMain::DstarToPiKKPi0, {+kKMinus, +kKPlus, +kPiPlus, +kPi0}}, + {DecayChannelMain::DstarToPiPiPi, {+kPiMinus, +kPiPlus, +kPiPlus}}, + {DecayChannelMain::DstarToPiPiPiPi0, {+kPiMinus, +kPiPlus, +kPiPlus, +kPi0}}, }; // Lc → p K∓ π± std::unordered_map> resoStatesLambdaC = -{ - {DecayChannelResonant::LcToPKstar0, {+kK0Star892, +kProton}}, - {DecayChannelResonant::LcToDeltaplusplusK, {+2224, +kKMinus}}, - {DecayChannelResonant::LcToL1520Pi, {+102134, +kPiPlus}}, + { + {DecayChannelResonant::LcToPKstar0, {+kK0Star892, +kProton}}, + {DecayChannelResonant::LcToDeltaplusplusK, {+2224, +kKMinus}}, + {DecayChannelResonant::LcToL1520Pi, {+102134, +kPiPlus}}, }; std::unordered_map> finalStatesLc = -{ - {DecayChannelMain::LcToPKPi, {+kProton, +kKMinus, +kPiPlus}}, - {DecayChannelMain::LcToPKPiPi0, {+kProton, +kKMinus, +kPiPlus, +kPi0}}, - {DecayChannelMain::LcToPPiPi, {+kProton, +kPiMinus, +kPiPlus}}, - {DecayChannelMain::LcToPKK, {+kProton, +kKMinus, +kKPlus}}}; + { + {DecayChannelMain::LcToPKPi, {+kProton, +kKMinus, +kPiPlus}}, + {DecayChannelMain::LcToPKPiPi0, {+kProton, +kKMinus, +kPiPlus, +kPi0}}, + {DecayChannelMain::LcToPPiPi, {+kProton, +kPiMinus, +kPiPlus}}, + {DecayChannelMain::LcToPKK, {+kProton, +kKMinus, +kKPlus}}}; // Xic → p K∓ π± std::unordered_map> resoStatesXiC = -{ - {DecayChannelResonant::XicToPKstar0, {-kK0Star892, +kProton}}, - {DecayChannelResonant::XicToPPhi, {+kProton, +kPhi}}, + { + {DecayChannelResonant::XicToPKstar0, {-kK0Star892, +kProton}}, + {DecayChannelResonant::XicToPPhi, {+kProton, +kPhi}}, }; std::unordered_map> finalStatesXic = -{ - {DecayChannelMain::XicToPKPi, {+kProton, +kKMinus, +kPiPlus}}, - {DecayChannelMain::XicToPKK, {+kProton, +kKMinus, +kKPlus}}, - {DecayChannelMain::XicToSPiPi, {+kSigmaPlus, +kPiMinus, +kPiPlus}}, + { + {DecayChannelMain::XicToPKPi, {+kProton, +kKMinus, +kPiPlus}}, + {DecayChannelMain::XicToPKK, {+kProton, +kKMinus, +kKPlus}}, + {DecayChannelMain::XicToSPiPi, {+kSigmaPlus, +kPiMinus, +kPiPlus}}, }; -} +} // namespace hf_chns_3prong /// Returns a map of the possible final states for a specific 3-prong particle specie /// \param pdgMother PDG code of the mother particle