From 40ff49e61d6300fbaa14eac3ed1d640e6e98a469 Mon Sep 17 00:00:00 2001 From: Fabio Colamaria Date: Wed, 6 Oct 2021 16:50:53 +0200 Subject: [PATCH 1/3] PWGHF: Adding centrality selection to DDbar correlation framework --- PWGHF/Tasks/HFCorrelatorD0D0bar.cxx | 23 +++++++++++++++++++++-- PWGHF/Tasks/HFCorrelatorDplusDminus.cxx | 23 +++++++++++++++++++++-- 2 files changed, 42 insertions(+), 4 deletions(-) diff --git a/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx b/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx index 43fa233babd..5bfa0d8af75 100644 --- a/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx +++ b/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx @@ -18,6 +18,7 @@ #include "Framework/HistogramRegistry.h" #include "PWGHF/DataModel/HFSecondaryVertex.h" #include "PWGHF/DataModel/HFCandidateSelectionTables.h" +#include "Common/DataModel/Centrality.h" using namespace o2; using namespace o2::framework; @@ -77,6 +78,7 @@ struct HfCorrelatorD0D0bar { {"hEta", "D0,D0bar candidates;candidate #it{#eta};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, {"hPhi", "D0,D0bar candidates;candidate #it{#varphi};entries", {HistType::kTH1F, {{32, 0., 2. * o2::constants::math::PI}}}}, {"hY", "D0,D0bar candidates;candidate #it{#y};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, + {"hCentrality", "centrality;centrality percentile;entries", {HistType::kTH1F, {{100, 0., 100.}}}}, {"hDDbarVsEtaCut", "D0,D0bar pairs vs #eta cut;#eta_{max};entries", {HistType::kTH2F, {{(int)(maxEtaCut / incrementEtaCut), 0., maxEtaCut}, {(int)(ptThresholdForMaxEtaCut / incrementPtThreshold), 0., ptThresholdForMaxEtaCut}}}}}}; Configurable selectionFlagD0{"selectionFlagD0", 1, "Selection Flag for D0"}; @@ -86,6 +88,8 @@ struct HfCorrelatorD0D0bar { Configurable> bins{"ptBinsForMassAndEfficiency", std::vector{o2::analysis::hf_cuts_d0_topik::pTBins_v}, "pT bin limits for candidate mass plots and efficiency"}; Configurable> efficiencyDmeson{"efficiencyDmeson", std::vector{efficiencyDmeson_v}, "Efficiency values for D0 meson"}; Configurable flagApplyEfficiency{"efficiencyFlagD", 1, "Flag for applying D-meson efficiency weights"}; + Configurable centMin{"centMin", 0., "minimum centrality accepted"}; + Configurable centMax{"centMax", 100., "maximum centrality accepted"}; void init(o2::framework::InitContext&) { @@ -96,8 +100,14 @@ struct HfCorrelatorD0D0bar { Filter filterSelectCandidates = (aod::hf_selcandidate_d0::isSelD0 >= selectionFlagD0 || aod::hf_selcandidate_d0::isSelD0bar >= selectionFlagD0bar); - void process(aod::Collision const& collision, soa::Filtered> const& candidates) + void process(soa::Join::iterator const& collision, soa::Filtered> const& candidates) { + float centrality = collision.centV0M(); + if (centrality < centMin || centrality > centMax) { + return; + } + registry.fill(HIST("hCentrality"), centrality); + for (auto& candidate1 : candidates) { if (cutYCandMax >= 0. && std::abs(YD0(candidate1)) > cutYCandMax) { continue; @@ -199,6 +209,7 @@ struct HfCorrelatorD0D0barMcRec { {"hEtaMCRec", "D0,D0bar candidates - MC reco;candidate #it{#eta};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, {"hPhiMCRec", "D0,D0bar candidates - MC reco;candidate #it{#varphi};entries", {HistType::kTH1F, {{32, 0., 2. * o2::constants::math::PI}}}}, {"hYMCRec", "D0,D0bar candidates - MC reco;candidate #it{#y};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, + {"hCentralityMCRec", "centrality;centrality percentile;entries", {HistType::kTH1F, {{100, 0., 100.}}}}, {"hDDbarVsEtaCut", "D0,D0bar pairs vs #eta cut;#eta_{max};entries", {HistType::kTH2F, {{(int)(maxEtaCut / incrementEtaCut), 0., maxEtaCut}, {(int)(ptThresholdForMaxEtaCut / incrementPtThreshold), 0., ptThresholdForMaxEtaCut}}}}}}; Configurable selectionFlagD0{"selectionFlagD0", 1, "Selection Flag for D0"}; @@ -208,6 +219,8 @@ struct HfCorrelatorD0D0barMcRec { Configurable> bins{"ptBinsForMassAndEfficiency", std::vector{o2::analysis::hf_cuts_d0_topik::pTBins_v}, "pT bin limits for candidate mass plots and efficiency"}; Configurable> efficiencyDmeson{"efficiencyDmeson", std::vector{efficiencyDmeson_v}, "Efficiency values for D0 meson"}; Configurable flagApplyEfficiency{"efficiencyFlagD", 1, "Flag for applying D-meson efficiency weights"}; + Configurable centMin{"centMin", 0., "minimum centrality accepted"}; + Configurable centMax{"centMax", 100., "maximum centrality accepted"}; void init(o2::framework::InitContext&) { @@ -221,8 +234,14 @@ struct HfCorrelatorD0D0barMcRec { Filter filterSelectCandidates = (aod::hf_selcandidate_d0::isSelD0 >= selectionFlagD0 || aod::hf_selcandidate_d0::isSelD0bar >= selectionFlagD0bar); - void process(aod::Collision const& collision, soa::Filtered> const& candidates) + void process(soa::Join::iterator const& collision, soa::Filtered> const& candidates) { + float centrality = collision.centV0M(); + if (centrality < centMin || centrality > centMax) { + return; + } + registry.fill(HIST("hCentralityMCRec"), centrality); + //MC reco level bool flagD0Signal = false; bool flagD0Reflection = false; diff --git a/PWGHF/Tasks/HFCorrelatorDplusDminus.cxx b/PWGHF/Tasks/HFCorrelatorDplusDminus.cxx index 3a82f84ecdd..a96443ea576 100644 --- a/PWGHF/Tasks/HFCorrelatorDplusDminus.cxx +++ b/PWGHF/Tasks/HFCorrelatorDplusDminus.cxx @@ -18,6 +18,7 @@ #include "Framework/HistogramRegistry.h" #include "PWGHF/DataModel/HFSecondaryVertex.h" #include "PWGHF/DataModel/HFCandidateSelectionTables.h" +#include "Common/DataModel/Centrality.h" using namespace o2; using namespace o2::framework; @@ -76,6 +77,7 @@ struct HfCorrelatorDplusDminus { {"hEta", "Dplus,Dminus candidates;candidate #it{#eta};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, {"hPhi", "Dplus,Dminus candidates;candidate #it{#varphi};entries", {HistType::kTH1F, {{32, 0., 2. * o2::constants::math::PI}}}}, {"hY", "Dplus,Dminus candidates;candidate #it{#y};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, + {"hCentrality", "centrality;centrality percentile;entries", {HistType::kTH1F, {{100, 0., 100.}}}}, {"hDDbarVsEtaCut", "Dplus,Dminus pairs vs #eta cut;#eta_{max};entries", {HistType::kTH2F, {{(int)(maxEtaCut / incrementEtaCut), 0., maxEtaCut}, {(int)(ptThresholdForMaxEtaCut / incrementPtThreshold), 0., ptThresholdForMaxEtaCut}}}}}}; Configurable selectionFlagDplus{"selectionFlagDplus", 1, "Selection Flag for Dplus,Dminus"}; @@ -84,6 +86,8 @@ struct HfCorrelatorDplusDminus { Configurable> bins{"ptBinsForMassAndEfficiency", std::vector{o2::analysis::hf_cuts_dplus_topikpi::pTBins_v}, "pT bin limits for candidate mass plots and efficiency"}; Configurable> efficiencyDmeson{"efficiencyDmeson", std::vector{efficiencyDmeson_v}, "Efficiency values for Dplus meson"}; Configurable flagApplyEfficiency{"efficiencyFlagD", 1, "Flag for applying D-meson efficiency weights"}; + Configurable centMin{"centMin", 0., "minimum centrality accepted"}; + Configurable centMax{"centMax", 100., "maximum centrality accepted"}; void init(o2::framework::InitContext&) { @@ -94,8 +98,14 @@ struct HfCorrelatorDplusDminus { Filter filterSelectCandidates = (aod::hf_selcandidate_dplus::isSelDplusToPiKPi >= selectionFlagDplus); - void process(aod::Collision const& collision, soa::Filtered> const& candidates, aod::BigTracks const& tracks) + void process(soa::Join::iterator const& collision, soa::Filtered> const& candidates, aod::BigTracks const& tracks) { + float centrality = collision.centV0M(); + if (centrality < centMin || centrality > centMax) { + return; + } + registry.fill(HIST("hCentrality"), centrality); + for (auto& candidate1 : candidates) { if (cutYCandMax >= 0. && std::abs(YDPlus(candidate1)) > cutYCandMax) { continue; @@ -196,6 +206,7 @@ struct HfCorrelatorDplusDminusMcRec { {"hEtaMCRec", "Dplus,Dminus candidates - MC reco;candidate #it{#eta};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, {"hPhiMCRec", "Dplus,Dminus candidates - MC reco;candidate #it{#varphi};entries", {HistType::kTH1F, {{32, 0., 2. * o2::constants::math::PI}}}}, {"hYMCRec", "Dplus,Dminus candidates - MC reco;candidate #it{#y};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, + {"hCentralityMCRec", "centrality;centrality percentile;entries", {HistType::kTH1F, {{100, 0., 100.}}}}, {"hDDbarVsEtaCut", "Dplus,Dminus pairs vs #eta cut;#eta_{max};entries", {HistType::kTH2F, {{(int)(maxEtaCut / incrementEtaCut), 0., maxEtaCut}, {(int)(ptThresholdForMaxEtaCut / incrementPtThreshold), 0., ptThresholdForMaxEtaCut}}}}}}; Configurable selectionFlagDplus{"selectionFlagDplus", 1, "Selection Flag for Dplus,Dminus"}; @@ -204,6 +215,8 @@ struct HfCorrelatorDplusDminusMcRec { Configurable> bins{"ptBinsForMassAndEfficiency", std::vector{o2::analysis::hf_cuts_dplus_topikpi::pTBins_v}, "pT bin limits for candidate mass plots and efficiency"}; Configurable> efficiencyDmeson{"efficiencyDmeson", std::vector{efficiencyDmeson_v}, "Efficiency values for D0 meson"}; Configurable flagApplyEfficiency{"efficiencyFlagD", 1, "Flag for applying D-meson efficiency weights"}; + Configurable centMin{"centMin", 0., "minimum centrality accepted"}; + Configurable centMax{"centMax", 100., "maximum centrality accepted"}; void init(o2::framework::InitContext&) { @@ -215,8 +228,14 @@ struct HfCorrelatorDplusDminusMcRec { Filter filterSelectCandidates = (aod::hf_selcandidate_dplus::isSelDplusToPiKPi >= selectionFlagDplus); - void process(aod::Collision const& collision, soa::Filtered> const& candidates, aod::BigTracks const& tracks) + void process(soa::Join::iterator const& collision, soa::Filtered> const& candidates, aod::BigTracks const& tracks) { + float centrality = collision.centV0M(); + if (centrality < centMin || centrality > centMax) { + return; + } + registry.fill(HIST("hCentralityMCRec"), centrality); + //MC reco level bool flagDplusSignal = false; bool flagDminusSignal = false; From 06e9e01565767ead6ce4bceaa602c0553482d608 Mon Sep 17 00:00:00 2001 From: Fabio Colamaria Date: Wed, 6 Oct 2021 17:18:44 +0200 Subject: [PATCH 2/3] Update to previous commit on DDbar centrality selection --- PWGHF/Tasks/HFCorrelatorD0D0bar.cxx | 8 ++++++-- PWGHF/Tasks/HFCorrelatorDplusDminus.cxx | 12 ++++++++---- 2 files changed, 14 insertions(+), 6 deletions(-) diff --git a/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx b/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx index 5bfa0d8af75..094e8ebc4ab 100644 --- a/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx +++ b/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx @@ -78,6 +78,7 @@ struct HfCorrelatorD0D0bar { {"hEta", "D0,D0bar candidates;candidate #it{#eta};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, {"hPhi", "D0,D0bar candidates;candidate #it{#varphi};entries", {HistType::kTH1F, {{32, 0., 2. * o2::constants::math::PI}}}}, {"hY", "D0,D0bar candidates;candidate #it{#y};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, + {"hCentralityPreSelection", "centrality prior to selection;centrality percentile;entries", {HistType::kTH1F, {{101, -1., 100.}}}}, {"hCentrality", "centrality;centrality percentile;entries", {HistType::kTH1F, {{100, 0., 100.}}}}, {"hDDbarVsEtaCut", "D0,D0bar pairs vs #eta cut;#eta_{max};entries", {HistType::kTH2F, {{(int)(maxEtaCut / incrementEtaCut), 0., maxEtaCut}, {(int)(ptThresholdForMaxEtaCut / incrementPtThreshold), 0., ptThresholdForMaxEtaCut}}}}}}; @@ -88,7 +89,7 @@ struct HfCorrelatorD0D0bar { Configurable> bins{"ptBinsForMassAndEfficiency", std::vector{o2::analysis::hf_cuts_d0_topik::pTBins_v}, "pT bin limits for candidate mass plots and efficiency"}; Configurable> efficiencyDmeson{"efficiencyDmeson", std::vector{efficiencyDmeson_v}, "Efficiency values for D0 meson"}; Configurable flagApplyEfficiency{"efficiencyFlagD", 1, "Flag for applying D-meson efficiency weights"}; - Configurable centMin{"centMin", 0., "minimum centrality accepted"}; + Configurable centMin{"centMin", -1., "minimum centrality accepted"}; Configurable centMax{"centMax", 100., "maximum centrality accepted"}; void init(o2::framework::InitContext&) @@ -103,6 +104,7 @@ struct HfCorrelatorD0D0bar { void process(soa::Join::iterator const& collision, soa::Filtered> const& candidates) { float centrality = collision.centV0M(); + registry.fill(HIST("hCentralityPreSelection"), centrality); if (centrality < centMin || centrality > centMax) { return; } @@ -209,6 +211,7 @@ struct HfCorrelatorD0D0barMcRec { {"hEtaMCRec", "D0,D0bar candidates - MC reco;candidate #it{#eta};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, {"hPhiMCRec", "D0,D0bar candidates - MC reco;candidate #it{#varphi};entries", {HistType::kTH1F, {{32, 0., 2. * o2::constants::math::PI}}}}, {"hYMCRec", "D0,D0bar candidates - MC reco;candidate #it{#y};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, + {"hCentralityPreSelectionMCRec", "centrality prior to selection;centrality percentile;entries", {HistType::kTH1F, {{101, -1., 100.}}}}, {"hCentralityMCRec", "centrality;centrality percentile;entries", {HistType::kTH1F, {{100, 0., 100.}}}}, {"hDDbarVsEtaCut", "D0,D0bar pairs vs #eta cut;#eta_{max};entries", {HistType::kTH2F, {{(int)(maxEtaCut / incrementEtaCut), 0., maxEtaCut}, {(int)(ptThresholdForMaxEtaCut / incrementPtThreshold), 0., ptThresholdForMaxEtaCut}}}}}}; @@ -219,7 +222,7 @@ struct HfCorrelatorD0D0barMcRec { Configurable> bins{"ptBinsForMassAndEfficiency", std::vector{o2::analysis::hf_cuts_d0_topik::pTBins_v}, "pT bin limits for candidate mass plots and efficiency"}; Configurable> efficiencyDmeson{"efficiencyDmeson", std::vector{efficiencyDmeson_v}, "Efficiency values for D0 meson"}; Configurable flagApplyEfficiency{"efficiencyFlagD", 1, "Flag for applying D-meson efficiency weights"}; - Configurable centMin{"centMin", 0., "minimum centrality accepted"}; + Configurable centMin{"centMin", -1., "minimum centrality accepted"}; Configurable centMax{"centMax", 100., "maximum centrality accepted"}; void init(o2::framework::InitContext&) @@ -237,6 +240,7 @@ struct HfCorrelatorD0D0barMcRec { void process(soa::Join::iterator const& collision, soa::Filtered> const& candidates) { float centrality = collision.centV0M(); + registry.fill(HIST("hCentralityPreSelectionMCRec"), centrality); if (centrality < centMin || centrality > centMax) { return; } diff --git a/PWGHF/Tasks/HFCorrelatorDplusDminus.cxx b/PWGHF/Tasks/HFCorrelatorDplusDminus.cxx index a96443ea576..d3c99d2d595 100644 --- a/PWGHF/Tasks/HFCorrelatorDplusDminus.cxx +++ b/PWGHF/Tasks/HFCorrelatorDplusDminus.cxx @@ -77,7 +77,8 @@ struct HfCorrelatorDplusDminus { {"hEta", "Dplus,Dminus candidates;candidate #it{#eta};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, {"hPhi", "Dplus,Dminus candidates;candidate #it{#varphi};entries", {HistType::kTH1F, {{32, 0., 2. * o2::constants::math::PI}}}}, {"hY", "Dplus,Dminus candidates;candidate #it{#y};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, - {"hCentrality", "centrality;centrality percentile;entries", {HistType::kTH1F, {{100, 0., 100.}}}}, + {"hCentralityPreSelection", "centrality prior selection;centrality percentile;entries", {HistType::kTH1F, {{101, -1., 100.}}}}, + {"hCentrality", "centrality;centrality percentile;entries", {HistType::kTH1F, {{101, -1., 100.}}}}, {"hDDbarVsEtaCut", "Dplus,Dminus pairs vs #eta cut;#eta_{max};entries", {HistType::kTH2F, {{(int)(maxEtaCut / incrementEtaCut), 0., maxEtaCut}, {(int)(ptThresholdForMaxEtaCut / incrementPtThreshold), 0., ptThresholdForMaxEtaCut}}}}}}; Configurable selectionFlagDplus{"selectionFlagDplus", 1, "Selection Flag for Dplus,Dminus"}; @@ -86,7 +87,7 @@ struct HfCorrelatorDplusDminus { Configurable> bins{"ptBinsForMassAndEfficiency", std::vector{o2::analysis::hf_cuts_dplus_topikpi::pTBins_v}, "pT bin limits for candidate mass plots and efficiency"}; Configurable> efficiencyDmeson{"efficiencyDmeson", std::vector{efficiencyDmeson_v}, "Efficiency values for Dplus meson"}; Configurable flagApplyEfficiency{"efficiencyFlagD", 1, "Flag for applying D-meson efficiency weights"}; - Configurable centMin{"centMin", 0., "minimum centrality accepted"}; + Configurable centMin{"centMin", -1., "minimum centrality accepted"}; Configurable centMax{"centMax", 100., "maximum centrality accepted"}; void init(o2::framework::InitContext&) @@ -101,6 +102,7 @@ struct HfCorrelatorDplusDminus { void process(soa::Join::iterator const& collision, soa::Filtered> const& candidates, aod::BigTracks const& tracks) { float centrality = collision.centV0M(); + registry.fill(HIST("hCentralityPreSelection"), centrality); if (centrality < centMin || centrality > centMax) { return; } @@ -206,7 +208,8 @@ struct HfCorrelatorDplusDminusMcRec { {"hEtaMCRec", "Dplus,Dminus candidates - MC reco;candidate #it{#eta};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, {"hPhiMCRec", "Dplus,Dminus candidates - MC reco;candidate #it{#varphi};entries", {HistType::kTH1F, {{32, 0., 2. * o2::constants::math::PI}}}}, {"hYMCRec", "Dplus,Dminus candidates - MC reco;candidate #it{#y};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, - {"hCentralityMCRec", "centrality;centrality percentile;entries", {HistType::kTH1F, {{100, 0., 100.}}}}, + {"hCentralityPreSelectionMCRec", "centrality prior selection;centrality percentile;entries", {HistType::kTH1F, {{101, -1., 100.}}}}, + {"hCentralityMCRec", "centrality;centrality percentile;entries", {HistType::kTH1F, {{101, -1., 100.}}}}, {"hDDbarVsEtaCut", "Dplus,Dminus pairs vs #eta cut;#eta_{max};entries", {HistType::kTH2F, {{(int)(maxEtaCut / incrementEtaCut), 0., maxEtaCut}, {(int)(ptThresholdForMaxEtaCut / incrementPtThreshold), 0., ptThresholdForMaxEtaCut}}}}}}; Configurable selectionFlagDplus{"selectionFlagDplus", 1, "Selection Flag for Dplus,Dminus"}; @@ -215,7 +218,7 @@ struct HfCorrelatorDplusDminusMcRec { Configurable> bins{"ptBinsForMassAndEfficiency", std::vector{o2::analysis::hf_cuts_dplus_topikpi::pTBins_v}, "pT bin limits for candidate mass plots and efficiency"}; Configurable> efficiencyDmeson{"efficiencyDmeson", std::vector{efficiencyDmeson_v}, "Efficiency values for D0 meson"}; Configurable flagApplyEfficiency{"efficiencyFlagD", 1, "Flag for applying D-meson efficiency weights"}; - Configurable centMin{"centMin", 0., "minimum centrality accepted"}; + Configurable centMin{"centMin", -1., "minimum centrality accepted"}; Configurable centMax{"centMax", 100., "maximum centrality accepted"}; void init(o2::framework::InitContext&) @@ -231,6 +234,7 @@ struct HfCorrelatorDplusDminusMcRec { void process(soa::Join::iterator const& collision, soa::Filtered> const& candidates, aod::BigTracks const& tracks) { float centrality = collision.centV0M(); + registry.fill(HIST("hCentralityPreSelectionMCRec"), centrality); if (centrality < centMin || centrality > centMax) { return; } From 2d6db688b668f5f6075cf99086afe7dbfb3d67d1 Mon Sep 17 00:00:00 2001 From: Fabio Colamaria Date: Wed, 6 Oct 2021 17:23:14 +0200 Subject: [PATCH 3/3] Further update to a histogram axis --- PWGHF/Tasks/HFCorrelatorD0D0bar.cxx | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx b/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx index 094e8ebc4ab..2c1761bcb60 100644 --- a/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx +++ b/PWGHF/Tasks/HFCorrelatorD0D0bar.cxx @@ -79,7 +79,7 @@ struct HfCorrelatorD0D0bar { {"hPhi", "D0,D0bar candidates;candidate #it{#varphi};entries", {HistType::kTH1F, {{32, 0., 2. * o2::constants::math::PI}}}}, {"hY", "D0,D0bar candidates;candidate #it{#y};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, {"hCentralityPreSelection", "centrality prior to selection;centrality percentile;entries", {HistType::kTH1F, {{101, -1., 100.}}}}, - {"hCentrality", "centrality;centrality percentile;entries", {HistType::kTH1F, {{100, 0., 100.}}}}, + {"hCentrality", "centrality;centrality percentile;entries", {HistType::kTH1F, {{101, -1., 100.}}}}, {"hDDbarVsEtaCut", "D0,D0bar pairs vs #eta cut;#eta_{max};entries", {HistType::kTH2F, {{(int)(maxEtaCut / incrementEtaCut), 0., maxEtaCut}, {(int)(ptThresholdForMaxEtaCut / incrementPtThreshold), 0., ptThresholdForMaxEtaCut}}}}}}; Configurable selectionFlagD0{"selectionFlagD0", 1, "Selection Flag for D0"}; @@ -212,7 +212,7 @@ struct HfCorrelatorD0D0barMcRec { {"hPhiMCRec", "D0,D0bar candidates - MC reco;candidate #it{#varphi};entries", {HistType::kTH1F, {{32, 0., 2. * o2::constants::math::PI}}}}, {"hYMCRec", "D0,D0bar candidates - MC reco;candidate #it{#y};entries", {HistType::kTH1F, {{100, -5., 5.}}}}, {"hCentralityPreSelectionMCRec", "centrality prior to selection;centrality percentile;entries", {HistType::kTH1F, {{101, -1., 100.}}}}, - {"hCentralityMCRec", "centrality;centrality percentile;entries", {HistType::kTH1F, {{100, 0., 100.}}}}, + {"hCentralityMCRec", "centrality;centrality percentile;entries", {HistType::kTH1F, {{101, -1., 100.}}}}, {"hDDbarVsEtaCut", "D0,D0bar pairs vs #eta cut;#eta_{max};entries", {HistType::kTH2F, {{(int)(maxEtaCut / incrementEtaCut), 0., maxEtaCut}, {(int)(ptThresholdForMaxEtaCut / incrementPtThreshold), 0., ptThresholdForMaxEtaCut}}}}}}; Configurable selectionFlagD0{"selectionFlagD0", 1, "Selection Flag for D0"};