Skip to content

PWGHF: add info about K-pi pairs from Lc->pKpi decays, and add direct/reso. decay channel info #6590

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
fgrosa merged 7 commits into from
Jun 21, 2024
Merged

PWGHF: add info about K-pi pairs from Lc->pKpi decays, and add direct/reso. decay channel info #6590

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
6125b7d
Add KPi mass selection for Lc->pKpi.
Jun 19, 2024
fc48928
Add KPi inv. mass to Lc sparses.
Jun 19, 2024
0ca0ebc
Fix pion mass in KPi mass calculation.
Jun 19, 2024
b6d7aed
Add reso. channel info in LC sparses.
Jun 19, 2024
3f6d62c
Add reso. channel info in gen. Lc sparses.
Jun 19, 2024
cc441dd
Add helper function for KPi pair mass.
Jun 20, 2024
0f18d4e
Take K*0(892) mass from PhysicsConstants.h in O2.
Jun 20, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 12 additions & 0 deletions PWGHF/Core/HfHelper.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -358,6 +358,18 @@ class HfHelper
return candidate.m(std::array{o2::constants::physics::MassPiPlus, o2::constants::physics::MassKPlus, o2::constants::physics::MassProton});
}

template <typename T>
auto invMassKPiPairLcToPKPi(const T& candidate)
{
return RecoDecay::m(std::array{candidate.pVectorProng1(), candidate.pVectorProng2()}, std::array{o2::constants::physics::MassKPlus, o2::constants::physics::MassPiPlus});
}

template <typename T>
auto invMassKPiPairLcToPiKP(const T& candidate)
{
return RecoDecay::m(std::array{candidate.pVectorProng1(), candidate.pVectorProng0()}, std::array{o2::constants::physics::MassKPlus, o2::constants::physics::MassPiPlus});
}

// Ξc± → p± K∓ π±

template <typename T>
Expand Down
24 changes: 12 additions & 12 deletions PWGHF/Core/SelectorCuts.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -462,7 +462,7 @@ constexpr double cuts[nBinsPt][nCutVars] = {{0.05, 0.2, 0.1, 1000.0, 0.2, 300.0,
namespace hf_cuts_lc_to_p_k_pi
{
static constexpr int nBinsPt = 10;
static constexpr int nCutVars = 7;
static constexpr int nCutVars = 8;
// default values for the pT bin edges (can be used to configure histogram axis)
// offset by 1 from the bin numbers in cuts array
constexpr double binsPt[nBinsPt + 1] = {
Expand All @@ -480,16 +480,16 @@ constexpr double binsPt[nBinsPt + 1] = {
auto vecBinsPt = std::vector<double>{binsPt, binsPt + nBinsPt + 1};

// default values for the cuts
constexpr double cuts[nBinsPt][nCutVars] = {{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0.}, /* 0 < pT < 1 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0.}, /* 1 < pT < 2 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0.}, /* 2 < pT < 3 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0.}, /* 3 < pT < 4 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0.}, /* 4 < pT < 5 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0.}, /* 5 < pT < 6 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0.}, /* 6 < pT < 8 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0.}, /* 8 < pT < 12 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0.}, /* 12 < pT < 24 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0.}}; /* 24 < pT < 36 */
constexpr double cuts[nBinsPt][nCutVars] = {{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0., -1.}, /* 0 < pT < 1 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0., -1.}, /* 1 < pT < 2 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0., -1.}, /* 2 < pT < 3 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0., -1.}, /* 3 < pT < 4 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0., -1.}, /* 4 < pT < 5 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0., -1.}, /* 5 < pT < 6 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0., -1.}, /* 6 < pT < 8 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0., -1.}, /* 8 < pT < 12 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0., -1.}, /* 12 < pT < 24 */
{0.400, 0.4, 0.4, 0.4, 0., 0.005, 0., -1.}}; /* 24 < pT < 36 */

// row labels
static const std::vector<std::string> labelsPt = {
Expand All @@ -505,7 +505,7 @@ static const std::vector<std::string> labelsPt = {
"pT bin 9"};

// column labels
static const std::vector<std::string> labelsCutVar = {"m", "pT p", "pT K", "pT Pi", "Chi2PCA", "decay length", "cos pointing angle"};
static const std::vector<std::string> labelsCutVar = {"m", "pT p", "pT K", "pT Pi", "Chi2PCA", "decay length", "cos pointing angle", "mass (Kpi)"};
} // namespace hf_cuts_lc_to_p_k_pi

namespace hf_cuts_lc_to_k0s_p
Expand Down
Loading