Skip to content

Commit

Permalink
optimized binning (AliceO2Group#3507)
Browse files Browse the repository at this point in the history
  • Loading branch information
@alcaliva
alcaliva authored Sep 26, 2023
1 parent 124dd23 commit 50be35f
Showing 1 changed file with 55 additions and 52 deletions.
107 changes: 55 additions & 52 deletions PWGLF/Tasks/QC/tpc_dedx_qa.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -19,6 +19,7 @@
#include "Common/DataModel/PIDResponse.h"
#include "Common/DataModel/TrackSelectionTables.h"
#include "Common/DataModel/Centrality.h"
#include "Common/DataModel/Multiplicity.h"
#include "Framework/ASoAHelpers.h"
#include "Framework/AnalysisDataModel.h"
#include "Framework/AnalysisTask.h"
Expand Down Expand Up @@ -115,191 +116,193 @@ struct tpc_dedx_qa {
Configurable<float> maxDCAxy{"maxDCAxy", 0.1f, "maxDCAxy"};
Configurable<float> maxDCAz{"maxDCAz", 0.1f, "maxDCAz"};
Configurable<bool> eventSelection{"eventSelection", true, "event selection"};
ConfigurableAxis pBins_Pi{"pBins_Pi", {VARIABLE_WIDTH, -10.0, -5.0, -4.5, -4.0, -3.8, -3.6, -3.4, -3.2, -3.0, -2.8, -2.6, -2.4, -2.2, -2.0, -1.9, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.9, -0.85, -0.8, -0.76, -0.72, -0.68, -0.64, -0.6, -0.58, -0.56, -0.54, -0.52, -0.5, -0.48, -0.46, -0.44, -0.42, -0.4, -0.38, -0.36, -0.34, -0.32, -0.3, -0.28, -0.26, -0.24, -0.22, -0.2, 0.0, 0.2, 0.22, 0.24, 0.26, 0.28, 0.3, 0.32, 0.34, 0.36, 0.38, 0.4, 0.42, 0.44, 0.46, 0.48, 0.5, 0.52, 0.54, 0.56, 0.58, 0.6, 0.64, 0.68, 0.72, 0.76, 0.8, 0.85, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.5, 5.0, 10}, "p bins for pions"};
ConfigurableAxis pBins_Ka{"pBins_Ka", {VARIABLE_WIDTH, -2.0, -1.9, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.9, -0.85, -0.8, -0.78, -0.76, -0.74, -0.72, -0.7, -0.68, -0.66, -0.64, -0.62, -0.6, -0.58, -0.56, -0.54, -0.52, -0.5, -0.49, -0.48, -0.47, -0.46, -0.45, -0.44, -0.43, -0.42, -0.41, -0.4, -0.39, -0.38, -0.37, -0.36, -0.35, -0.34, -0.33, -0.32, -0.31, -0.3, 0.0, 0.3, 0.31, 0.32, 0.33, 0.34, 0.35, 0.36, 0.37, 0.38, 0.39, 0.4, 0.41, 0.42, 0.43, 0.44, 0.45, 0.46, 0.47, 0.48, 0.49, 0.5, 0.52, 0.54, 0.56, 0.58, 0.6, 0.62, 0.64, 0.66, 0.68, 0.70, 0.72, 0.74, 0.76, 0.78, 0.8, 0.85, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0}, "p bins for kaons"};
ConfigurableAxis pBins_Pr{"pBins_Pr", {VARIABLE_WIDTH, -10.0, -8.0, -7.0, -6.5, -6.0, -5.5, -5.0, -4.5, -4.0, -3.8, -3.6, -3.4, -3.2, -3.0, -2.8, -2.6, -2.4, -2.2, -2.0, -1.9, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.9, -0.85, -0.8, -0.78, -0.76, -0.74, -0.72, -0.70, -0.68, -0.66, -0.64, -0.62, -0.6, -0.58, -0.56, -0.54, -0.52, -0.5, -0.49, -0.48, -0.47, -0.46, -0.45, -0.44, -0.43, -0.42, -0.41, -0.4, -0.39, -0.38, -0.37, -0.36, -0.35, -0.34, -0.33, -0.32, -0.31, -0.3, -0.29, -0.28, -0.27, -0.26, -0.25, -0.24, -0.23, -0.22, -0.21, -0.2, 0.0, 0.2, 0.21, 0.22, 0.23, 0.24, 0.25, 0.26, 0.27, 0.28, 0.29, 0.3, 0.31, 0.32, 0.33, 0.34, 0.35, 0.36, 0.37, 0.38, 0.39, 0.4, 0.41, 0.42, 0.43, 0.44, 0.45, 0.46, 0.47, 0.48, 0.49, 0.5, 0.52, 0.54, 0.56, 0.58, 0.6, 0.62, 0.64, 0.66, 0.68, 0.70, 0.72, 0.74, 0.76, 0.78, 0.8, 0.85, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5, 7.0, 8.0, 10.0}, "p bins for protons"};
ConfigurableAxis pBins_De{"pBins_De", {VARIABLE_WIDTH, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1, -0.95, -0.9, -0.88, -0.86, -0.84, -0.82, -0.8, -0.78, -0.76, -0.74, -0.72, -0.7, -0.68, -0.66, -0.64, -0.62, -0.6, -0.58, -0.56, -0.54, -0.52, -0.5, 0.0, 0.5, 0.52, 0.54, 0.56, 0.58, 0.6, 0.62, 0.64, 0.66, 0.68, 0.70, 0.72, 0.74, 0.76, 0.78, 0.8, 0.82, 0.84, 0.86, 0.88, 0.9, 0.95, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8}, "p bins for deuterons"};
ConfigurableAxis pBins_He{"pBins_He", {VARIABLE_WIDTH, -10, -8, -7, -6, -5.5, -5, -4.9, -4.8, -4.7, -4.6, -4.5, -4.4, -4.3, -4.2, -4.1, -4, -3.9, -3.8, -3.7, -3.6, -3.5, -3.4, -3.3, -3.2, -3.1, -3, -2.9, -2.8, -2.7, -2.6, -2.5, -2.4, -2.3, -2.2, -2.1, -2, -1.9, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1, -0.98, -0.96, -0.94, -0.92, -0.9, -0.88, -0.86, -0.84, -0.82, -0.8, 0.0, 0.8, 0.82, 0.84, 0.86, 0.88, 0.9, 0.92, 0.94, 0.96, 0.98, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.1, 3.2, 3.3, 3.4, 3.5, 3.6, 3.7, 3.8, 3.9, 4.0, 4.1, 4.2, 4.3, 4.4, 4.5, 4.6, 4.7, 4.8, 4.9, 5.0, 5.5, 6.0, 7.0, 8.0, 10.0}, "p bins for helium-3"};

// Reduce binning
ConfigurableAxis pBins_Pi{"pBins_Pi", {VARIABLE_WIDTH, -10.0, -5.0, -4.5, -4.0, -3.8, -3.6, -3.4, -3.2, -3.0, -2.8, -2.6, -2.4, -2.2, -2.0, -1.9, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.9, -0.8, -0.75, -0.7, -0.65, -0.6, -0.55, -0.5, -0.45, -0.4, -0.35, -0.3, -0.28, -0.26, -0.24, -0.22, -0.2, 0.0, 0.2, 0.22, 0.24, 0.26, 0.28, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.5, 5.0, 10}, "p bins for pions"};
ConfigurableAxis pBins_Ka{"pBins_Ka", {VARIABLE_WIDTH, -2.0, -1.9, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.9, -0.8, -0.75, -0.7, -0.65, -0.6, -0.55, -0.5, -0.45, -0.4, -0.35, -0.3, -0.28, -0.26, -0.24, -0.22, -0.2, 0.0, 0.2, 0.22, 0.24, 0.26, 0.28, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0}, "p bins for kaons"};
ConfigurableAxis pBins_Pr{"pBins_Pr", {VARIABLE_WIDTH, -5.0, -4.5, -4.0, -3.5, -3.0, -2.8, -2.6, -2.4, -2.2, -2.0, -1.9, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.9, -0.8, -0.75, -0.7, -0.65, -0.6, -0.55, -0.5, -0.48, -0.46, -0.44, -0.42, -0.4, -0.38, -0.36, -0.34, -0.32, -0.3, 0.0, 0.3, 0.32, 0.34, 0.36, 0.38, 0.4, 0.42, 0.44, 0.46, 0.48, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.5, 4.0, 4.5, 5.0}, "p bins for protons"};
ConfigurableAxis pBins_De{"pBins_De", {VARIABLE_WIDTH, -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.9, -0.8, -0.75, -0.7, -0.65, -0.6, -0.55, -0.5, -0.48, -0.46, -0.44, -0.42, -0.4, -0.38, -0.36, -0.34, -0.32, -0.3, 0.0, 0.3, 0.32, 0.34, 0.36, 0.38, 0.4, 0.42, 0.44, 0.46, 0.48, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5}, "p bins for deuterons"};
ConfigurableAxis pBins_He{"pBins_He", {VARIABLE_WIDTH, -5.0, -4.5, -4.0, -3.5, -3.0, -2.8, -2.6, -2.4, -2.2, -2.0, -1.9, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, 0.0, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.5, 4.0, 4.5, 5.0}, "p bins for helium-3"};

void init(InitContext const&)
{

// Variable binning
AxisSpec pAxis_Pi{pBins_Pi, "z#upoint p (GeV/c)"};
AxisSpec pAxis_Ka{pBins_Ka, "z#upoint p (GeV/c)"};
AxisSpec pAxis_Pr{pBins_Pr, "z#upoint p (GeV/c)"};
AxisSpec pAxis_De{pBins_De, "z#upoint p (GeV/c)"};
AxisSpec pAxis_He{pBins_He, "z#upoint p (GeV/c)"};
AxisSpec pAxis_Pi{pBins_Pi, "#it{p}/Z (GeV/c)"};
AxisSpec pAxis_Ka{pBins_Ka, "#it{p}/Z (GeV/c)"};
AxisSpec pAxis_Pr{pBins_Pr, "#it{p}/Z (GeV/c)"};
AxisSpec pAxis_De{pBins_De, "#it{p}/Z (GeV/c)"};
AxisSpec pAxis_He{pBins_He, "#it{p}/Z (GeV/c)"};

// Charged Particles
registryCh.add(
"dEdx_vs_Momentum", "dE/dx", HistType::kTH2F,
{{200, -10.0, 10.0, "z#upoint p (GeV/c)"}, {1500, 0, 1500, "dE/dx (a. u.)"}});
{{100, -10.0, 10.0, "#it{p}/Z (GeV/c)"}, {100, 0.0, 600.0, "dE/dx (a. u.)"}});

// Pions
registryPi.add(
"dEdx_vs_Momentum_Pi_m0906", "dE/dx", HistType::kTH3F,
{pAxis_Pi, {200, 0, 200, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pi, {150, 0.0, 150.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPi.add(
"dEdx_vs_Momentum_Pi_m0604", "dE/dx", HistType::kTH3F,
{pAxis_Pi, {200, 0, 200, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pi, {150, 0.0, 150.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPi.add(
"dEdx_vs_Momentum_Pi_m0402", "dE/dx", HistType::kTH3F,
{pAxis_Pi, {200, 0, 200, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pi, {150, 0.0, 150.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPi.add(
"dEdx_vs_Momentum_Pi_m0200", "dE/dx", HistType::kTH3F,
{pAxis_Pi, {200, 0, 200, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pi, {150, 0.0, 150.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPi.add(
"dEdx_vs_Momentum_Pi_p0002", "dE/dx", HistType::kTH3F,
{pAxis_Pi, {200, 0, 200, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pi, {150, 0.0, 150.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPi.add(
"dEdx_vs_Momentum_Pi_p0204", "dE/dx", HistType::kTH3F,
{pAxis_Pi, {200, 0, 200, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pi, {150, 0.0, 150.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPi.add(
"dEdx_vs_Momentum_Pi_p0406", "dE/dx", HistType::kTH3F,
{pAxis_Pi, {200, 0, 200, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pi, {150, 0.0, 150.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPi.add(
"dEdx_vs_Momentum_Pi_p0609", "dE/dx", HistType::kTH3F,
{pAxis_Pi, {200, 0, 200, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pi, {150, 0.0, 150.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

// Kaons
registryKa.add(
"dEdx_vs_Momentum_Ka_m0906", "dE/dx", HistType::kTH3F,
{pAxis_Ka, {600, 0, 600, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Ka, {200, 0.0, 400.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryKa.add(
"dEdx_vs_Momentum_Ka_m0604", "dE/dx", HistType::kTH3F,
{pAxis_Ka, {600, 0, 600, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Ka, {200, 0.0, 400.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryKa.add(
"dEdx_vs_Momentum_Ka_m0402", "dE/dx", HistType::kTH3F,
{pAxis_Ka, {600, 0, 600, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Ka, {200, 0.0, 400.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryKa.add(
"dEdx_vs_Momentum_Ka_m0200", "dE/dx", HistType::kTH3F,
{pAxis_Ka, {600, 0, 600, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Ka, {200, 0.0, 400.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryKa.add(
"dEdx_vs_Momentum_Ka_p0002", "dE/dx", HistType::kTH3F,
{pAxis_Ka, {600, 0, 600, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Ka, {200, 0.0, 400.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryKa.add(
"dEdx_vs_Momentum_Ka_p0204", "dE/dx", HistType::kTH3F,
{pAxis_Ka, {600, 0, 600, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Ka, {200, 0.0, 400.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryKa.add(
"dEdx_vs_Momentum_Ka_p0406", "dE/dx", HistType::kTH3F,
{pAxis_Ka, {600, 0, 600, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Ka, {200, 0.0, 400.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryKa.add(
"dEdx_vs_Momentum_Ka_p0609", "dE/dx", HistType::kTH3F,
{pAxis_Ka, {600, 0, 600, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Ka, {200, 0.0, 400.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

// Protons
registryPr.add(
"dEdx_vs_Momentum_Pr_m0906", "dE/dx", HistType::kTH3F,
{pAxis_Pr, {1000, 0, 1000, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pr, {200, 0.0, 500.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPr.add(
"dEdx_vs_Momentum_Pr_m0604", "dE/dx", HistType::kTH3F,
{pAxis_Pr, {1000, 0, 1000, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pr, {200, 0.0, 500.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPr.add(
"dEdx_vs_Momentum_Pr_m0402", "dE/dx", HistType::kTH3F,
{pAxis_Pr, {1000, 0, 1000, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pr, {200, 0.0, 500.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPr.add(
"dEdx_vs_Momentum_Pr_m0200", "dE/dx", HistType::kTH3F,
{pAxis_Pr, {1000, 0, 1000, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pr, {200, 0.0, 500.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPr.add(
"dEdx_vs_Momentum_Pr_p0002", "dE/dx", HistType::kTH3F,
{pAxis_Pr, {1000, 0, 1000, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pr, {200, 0.0, 500.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPr.add(
"dEdx_vs_Momentum_Pr_p0204", "dE/dx", HistType::kTH3F,
{pAxis_Pr, {1000, 0, 1000, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pr, {200, 0.0, 500.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPr.add(
"dEdx_vs_Momentum_Pr_p0406", "dE/dx", HistType::kTH3F,
{pAxis_Pr, {1000, 0, 1000, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pr, {200, 0.0, 500.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryPr.add(
"dEdx_vs_Momentum_Pr_p0609", "dE/dx", HistType::kTH3F,
{pAxis_Pr, {1000, 0, 1000, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_Pr, {200, 0.0, 500.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

// Deuterons
registryDe.add(
"dEdx_vs_Momentum_De_m0906", "dE/dx", HistType::kTH3F,
{pAxis_De, {1500, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_De, {300, 0.0, 600.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryDe.add(
"dEdx_vs_Momentum_De_m0604", "dE/dx", HistType::kTH3F,
{pAxis_De, {1500, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_De, {300, 0.0, 600.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryDe.add(
"dEdx_vs_Momentum_De_m0402", "dE/dx", HistType::kTH3F,
{pAxis_De, {1500, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_De, {300, 0.0, 600.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryDe.add(
"dEdx_vs_Momentum_De_m0200", "dE/dx", HistType::kTH3F,
{pAxis_De, {1500, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_De, {300, 0.0, 600.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryDe.add(
"dEdx_vs_Momentum_De_p0002", "dE/dx", HistType::kTH3F,
{pAxis_De, {1500, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_De, {300, 0.0, 600.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryDe.add(
"dEdx_vs_Momentum_De_p0204", "dE/dx", HistType::kTH3F,
{pAxis_De, {1500, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_De, {300, 0.0, 600.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryDe.add(
"dEdx_vs_Momentum_De_p0406", "dE/dx", HistType::kTH3F,
{pAxis_De, {1500, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_De, {300, 0.0, 600.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryDe.add(
"dEdx_vs_Momentum_De_p0609", "dE/dx", HistType::kTH3F,
{pAxis_De, {1500, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_De, {300, 0.0, 600.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

// Helium-3
registryHe.add(
"dEdx_vs_Momentum_He_m0906", "dE/dx", HistType::kTH3F,
{pAxis_He, {1000, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_He, {400, 0.0, 800.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryHe.add(
"dEdx_vs_Momentum_He_m0604", "dE/dx", HistType::kTH3F,
{pAxis_He, {1000, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_He, {400, 0.0, 800.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryHe.add(
"dEdx_vs_Momentum_He_m0402", "dE/dx", HistType::kTH3F,
{pAxis_He, {1000, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_He, {400, 0.0, 800.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryHe.add(
"dEdx_vs_Momentum_He_m0200", "dE/dx", HistType::kTH3F,
{pAxis_He, {1000, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_He, {400, 0.0, 800.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryHe.add(
"dEdx_vs_Momentum_He_p0002", "dE/dx", HistType::kTH3F,
{pAxis_He, {1000, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_He, {400, 0.0, 800.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryHe.add(
"dEdx_vs_Momentum_He_p0204", "dE/dx", HistType::kTH3F,
{pAxis_He, {1000, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_He, {400, 0.0, 800.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryHe.add(
"dEdx_vs_Momentum_He_p0406", "dE/dx", HistType::kTH3F,
{pAxis_He, {1000, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_He, {400, 0.0, 800.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

registryHe.add(
"dEdx_vs_Momentum_He_p0609", "dE/dx", HistType::kTH3F,
{pAxis_He, {1000, 0, 1500, "dE/dx (a. u.)"}, {10, 0, 100, "centrality"}});
{pAxis_He, {400, 0.0, 800.0, "dE/dx (a. u.)"}, {10, 0.0, 100.0, "centrality"}});

// Event Counter
registryCh.add("histRecVtxZData", "collision z position", HistType::kTH1F, {{100, -20.0, +20.0, "z_{vtx} (cm)"}});
Expand Down Expand Up @@ -540,7 +543,7 @@ struct tpc_dedx_qa {
}

// Helium-3
if (trk.tpcSignal() > 180 && trk.tpcInnerParam() > 0.8 && trk.tpcNSigmaDe() > 2.0) {
if (trk.tpcSignal() > 180.0 && trk.tpcInnerParam() > 0.8 && trk.tpcNSigmaDe() > 2.0) {
if (trk.eta() > -0.9 && trk.eta() < -0.6)
registryHe.fill(HIST("dEdx_vs_Momentum_He_m0906"), signedP, trk.tpcSignal(), centrality);
if (trk.eta() > -0.6 && trk.eta() < -0.4)
Expand Down

0 comments on commit 50be35f

Please sign in to comment.