Skip to content

PWGHF: add pt-dependent Lc-window cut and fix ConfigAxes. #3826

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 1 commit into from
Nov 14, 2023
Merged

PWGHF: add pt-dependent Lc-window cut and fix ConfigAxes. #3826

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
53 changes: 51 additions & 2 deletions PWGHF/Core/SelectorCuts.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -15,8 +15,8 @@
#ifndef PWGHF_CORE_SELECTORCUTS_H_
#define PWGHF_CORE_SELECTORCUTS_H_

#include <string> // std::string
#include <vector> // std::vector
#include <string> // std::string
#include <vector> // std::vector

namespace o2::analysis
{
Expand Down Expand Up @@ -927,6 +927,55 @@ static const std::vector<std::string> labelsPt = {
static const std::vector<std::string> labelsCutVar = {"m", "CPA", "d0 Jpsi", "d0 Gamma", "pT Jpsi", "pT Gamma", "chi2PCA"};
} // namespace hf_cuts_chic_to_jpsi_gamma

namespace hf_cuts_sigmac_to_p_k_pi
{
static constexpr int nBinsPt = 10;
static constexpr int nCutVars = 2;
// 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] = {
0.,
1.,
2.,
3.,
4.,
5.,
6.,
8.,
12.,
24.,
36.};
auto vecBinsPt = std::vector<double>{binsPt, binsPt + nBinsPt + 1};

// default values for the cuts
constexpr double cuts[nBinsPt][nCutVars] = {{0.03, 0.03}, /* 0 < pT < 1 */
Comment on lines +936 to +951
Copy link
Collaborator

@fgrosa fgrosa Nov 14, 2023

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I know that the cut values defined for all the other channels are double, however we don't need double precision for them. It is ok for the moment to use double for consistency, but I would keep it in mind and modify all of them to float in the future.

{0.03, 0.03}, /* 1 < pT < 2 */
{0.03, 0.03}, /* 2 < pT < 3 */
{0.03, 0.03}, /* 3 < pT < 4 */
{0.03, 0.03}, /* 4 < pT < 5 */
{0.03, 0.03}, /* 5 < pT < 6 */
{0.03, 0.03}, /* 6 < pT < 8 */
{0.03, 0.03}, /* 8 < pT < 12 */
{0.03, 0.03}, /* 12 < pT < 24 */
{0.03, 0.03}}; /* 24 < pT < 36 */

// row labels
static const std::vector<std::string> labelsPt = {
"pT bin 0 Lc",
"pT bin 1 Lc",
"pT bin 2 Lc",
"pT bin 3 Lc",
"pT bin 4 Lc",
"pT bin 5 Lc",
"pT bin 6 Lc",
"pT bin 7 Lc",
"pT bin 8 Lc",
"pT bin 9 Lc"};

// column labels
static const std::vector<std::string> labelsCutVar = {"max pKpi mass Lc", "max piKp mass Lc"};
} // namespace hf_cuts_sigmac_to_p_k_pi

} // namespace o2::analysis

#endif // PWGHF_CORE_SELECTORCUTS_H_
Loading