Skip to content

Commit ce38558

Browse files
pstahlhupstahlhu
committed
Introduce separate eta cuts and change to TH1D for counter histograms
1 parent 18681b8 commit ce38558

File tree

3 files changed

+46
-30
lines changed

3 files changed

+46
-30
lines changed

‎PWGHF/Core/SelectorCuts.h‎

Lines changed: 16 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -788,7 +788,7 @@ static const std::vector<std::string> labelsCutVar = {"m", "pT p", "pT K", "pT P
788788
namespace hf_cuts_xic_to_xi_pi_pi
789789
{
790790
static constexpr int NBinsPt = 13;
791-
static constexpr int NCutVars = 12;
791+
static constexpr int NCutVars = 13;
792792
// default values for the pT bin edges (can be used to configure histogram axis)
793793
// offset by 1 from the bin numbers in cuts array
794794
constexpr double BinsPt[NBinsPt + 1] = {
@@ -808,20 +808,20 @@ constexpr double BinsPt[NBinsPt + 1] = {
808808
20.};
809809
auto vecBinsPt = std::vector<double>{BinsPt, BinsPt + NBinsPt + 1};
810810

811-
// default values for the cuts m Y Eta EtaDau pT Pi0 Pi1 Sum chi2SV dL dLXY invMass Xi-Pi pairs
812-
constexpr double Cuts[NBinsPt][NCutVars] = {{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 0 < pT < 1 */
813-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 1 < pT < 2 */
814-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 2 < pT < 3 */
815-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 3 < pT < 4 */
816-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 4 < pT < 5 */
817-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 5 < pT < 6 */
818-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 6 < pT < 7 */
819-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 7 < pT < 8 */
820-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 8 < pT < 9 */
821-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 9 < pT < 10 */
822-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 10 < pT < 11 */
823-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 11 < pT < 12 */
824-
{0.4, 0.8, 0.8, 0.8, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}}; /* 12 < pT < 20 */
811+
// default values for the cuts m Y Eta EtaPi EtaXi pT Pi0 Pi1 Sum chi2SV dL dLXY invMass Xi-Pi pairs
812+
constexpr double Cuts[NBinsPt][NCutVars] = {{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 0 < pT < 1 */
813+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 1 < pT < 2 */
814+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 2 < pT < 3 */
815+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 3 < pT < 4 */
816+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 4 < pT < 5 */
817+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 5 < pT < 6 */
818+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 6 < pT < 7 */
819+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 7 < pT < 8 */
820+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 8 < pT < 9 */
821+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 9 < pT < 10 */
822+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 10 < pT < 11 */
823+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}, /* 11 < pT < 12 */
824+
{0.4, 0.8, 0.8, 0.8, 1.0, 0.1, 0.1, 0.2, 100, 0.0, 0.0, 2.4, 2.4}}; /* 12 < pT < 20 */
825825

826826
// row labels
827827
static const std::vector<std::string> labelsPt = {
@@ -840,7 +840,7 @@ static const std::vector<std::string> labelsPt = {
840840
"pT bin 12"};
841841

842842
// column labels
843-
static const std::vector<std::string> labelsCutVar = {"m", "y", "eta", "eta Daughters", "pT Pi0", "pT Pi1", "pT Pi0 + Pi1", "chi2SV", "min decay length", "min decay length XY", "max inv mass Xi-Pi0", "max inv mass Xi-Pi1"};
843+
static const std::vector<std::string> labelsCutVar = {"m", "y", "eta", "eta Pi from XicPlus", "eta Xi Daughters", "pT Pi0", "pT Pi1", "pT Pi0 + Pi1", "chi2SV", "min decay length", "min decay length XY", "max inv mass Xi-Pi0", "max inv mass Xi-Pi1"};
844844
} // namespace hf_cuts_xic_to_xi_pi_pi
845845

846846
namespace hf_cuts_xicc_to_p_k_pi_pi

‎PWGHF/TableProducer/candidateCreatorXicToXiPiPi.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -146,7 +146,7 @@ struct HfCandidateCreatorXicToXiPiPi {
146146
registry.add("hVertexerType", "Use DCAFitter or KFParticle;;entries", {HistType::kTH1F, {{2, -0.5, 1.5}}});
147147
registry.get<TH1>(HIST("hVertexerType"))->GetXaxis()->SetBinLabel(1 + aod::hf_cand::VertexerType::DCAFitter, "DCAFitter");
148148
registry.get<TH1>(HIST("hVertexerType"))->GetXaxis()->SetBinLabel(1 + aod::hf_cand::VertexerType::KfParticle, "KFParticle");
149-
registry.add("hCandCounter", "hCandCounter", {HistType::kTH1F, {{4, -0.5, 3.5}}});
149+
registry.add("hCandCounter", "hCandCounter", {HistType::kTH1D, {{4, -0.5, 3.5}}});
150150
registry.get<TH1>(HIST("hCandCounter"))->GetXaxis()->SetBinLabel(1 + TotalSkimmedTriplets, "total");
151151
registry.get<TH1>(HIST("hCandCounter"))->GetXaxis()->SetBinLabel(1 + SelEvent, "Event selected");
152152
registry.get<TH1>(HIST("hCandCounter"))->GetXaxis()->SetBinLabel(1 + CascPreSel, "Cascade preselection");

‎PWGHF/TableProducer/candidateSelectorXicToXiPiPi.cxx‎

Lines changed: 29 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -101,7 +101,8 @@ struct HfCandidateSelectorXicToXiPiPi {
101101
Mass,
102102
Rapidity,
103103
Eta,
104-
EtaDaughters,
104+
EtaPionFromXicPlus,
105+
EtaXiDaughters,
105106
PtPionFromXicPlus,
106107
Chi2SV,
107108
MinDecayLength,
@@ -146,25 +147,26 @@ struct HfCandidateSelectorXicToXiPiPi {
146147
labels[Mass] = "#Delta M";
147148
labels[Rapidity] = "y";
148149
labels[Eta] = "#eta";
149-
labels[EtaDaughters] = "#eta final state daughters";
150-
labels[PtPionFromXicPlus] = "#it{p}_{T} (#pi #leftarrow #Xi_{c}^{+})";
150+
labels[EtaPionFromXicPlus] = "#eta (#pi #leftarrow #Xi_{c}^{#plus})";
151+
labels[EtaXiDaughters] = "#eta (#Xi daughters)";
152+
labels[PtPionFromXicPlus] = "#it{p}_{T} (#pi #leftarrow #Xi_{c}^{#plus})";
151153
labels[Chi2SV] = "#chi^{2}_{SV}";
152154
labels[MinDecayLength] = "Decay length";
153155
labels[MaxInvMassXiPiPairs] = "M_{#Xi #pi}";
154156
labels[TpcTrackQualityXiDaughters] = "TPC track quality selection on #Xi daughters";
155-
labels[TpcTrackQualityPiFromCharm] = "TPC track quality selection on #pi #leftarrow #Xi_{c}^{+}";
156-
labels[ItsTrackQualityPiFromCharm] = "ITS track quality selection on #pi #leftarrow #Xi_{c}^{+}";
157+
labels[TpcTrackQualityPiFromCharm] = "TPC track quality selection on #pi #leftarrow #Xi_{c}^{#plus}";
158+
labels[ItsTrackQualityPiFromCharm] = "ITS track quality selection on #pi #leftarrow #Xi_{c}^{#plus}";
157159
labels[PidSelected] = "PID selection";
158160
labels[BdtSelected] = "BDT selection";
159161

160162
if (doprocessData) {
161-
registry.add("hSelCandidates", ";;entries", {HistType::kTH1F, {{NSelectionCriteria, -0.5f, +NSelectionCriteria - 0.5f}}});
163+
registry.add("hSelCandidates", ";;entries", {HistType::kTH1D, {{NSelectionCriteria, -0.5f, +NSelectionCriteria - 0.5f}}});
162164
for (int iBin = 0; iBin < NSelectionCriteria; ++iBin) {
163165
registry.get<TH1>(HIST("hSelCandidates"))->GetXaxis()->SetBinLabel(iBin + 1, labels[iBin].data());
164166
}
165167
} else if (doprocessMc) {
166-
registry.add("hSelCandidatesRecSig", ";;entries", {HistType::kTH1F, {{NSelectionCriteria, -0.5f, +NSelectionCriteria - 0.5f}}});
167-
registry.add("hSelCandidatesRecBkg", ";;entries", {HistType::kTH1F, {{NSelectionCriteria, -0.5f, +NSelectionCriteria - 0.5f}}});
168+
registry.add("hSelCandidatesRecSig", ";;entries", {HistType::kTH1D, {{NSelectionCriteria, -0.5f, +NSelectionCriteria - 0.5f}}});
169+
registry.add("hSelCandidatesRecBkg", ";;entries", {HistType::kTH1D, {{NSelectionCriteria, -0.5f, +NSelectionCriteria - 0.5f}}});
168170
for (int iBin = 0; iBin < NSelectionCriteria; ++iBin) {
169171
registry.get<TH1>(HIST("hSelCandidatesRecSig"))->GetXaxis()->SetBinLabel(iBin + 1, labels[iBin].data());
170172
registry.get<TH1>(HIST("hSelCandidatesRecBkg"))->GetXaxis()->SetBinLabel(iBin + 1, labels[iBin].data());
@@ -443,18 +445,32 @@ struct HfCandidateSelectorXicToXiPiPi {
443445
registry.fill(HIST("hSelCandidates"), Eta);
444446
}
445447

446-
// cut on pseudorapidity of final state daughters
447-
if (std::abs(trackPi0.eta()) > cuts->get(pTBin, "eta Daughters") || std::abs(trackPi1.eta()) > cuts->get(pTBin, "eta Daughters") || std::abs(trackPiFromXi.eta()) > cuts->get(pTBin, "eta Daughters") || std::abs(trackV0PosDau.eta()) > cuts->get(pTBin, "eta Daughters") || std::abs(trackV0NegDau.eta()) > cuts->get(pTBin, "eta Daughters")) {
448+
// cut on pseudorapidity of pions from XicPlus
449+
if (std::abs(trackPi0.eta()) > cuts->get(pTBin, "eta Pi from XicPlus") || std::abs(trackPi1.eta()) > cuts->get(pTBin, "eta Pi from XicPlus")) {
448450
return false;
449451
}
450452
if constexpr (isMc) {
451453
if (isMatchedSignal) {
452-
registry.fill(HIST("hSelCandidatesRecSig"), EtaDaughters);
454+
registry.fill(HIST("hSelCandidatesRecSig"), EtaPionFromXicPlus);
453455
} else {
454-
registry.fill(HIST("hSelCandidatesRecBkg"), EtaDaughters);
456+
registry.fill(HIST("hSelCandidatesRecBkg"), EtaPionFromXicPlus);
455457
}
456458
} else {
457-
registry.fill(HIST("hSelCandidates"), EtaDaughters);
459+
registry.fill(HIST("hSelCandidates"), EtaPionFromXicPlus);
460+
}
461+
462+
// cut on pseudorapidity of Xi daughters
463+
if (std::abs(trackPiFromXi.eta()) > cuts->get(pTBin, "eta Xi Daughters") || std::abs(trackV0PosDau.eta()) > cuts->get(pTBin, "eta Xi Daughters") || std::abs(trackV0NegDau.eta()) > cuts->get(pTBin, "eta Xi Daughters")) {
464+
return false;
465+
}
466+
if constexpr (isMc) {
467+
if (isMatchedSignal) {
468+
registry.fill(HIST("hSelCandidatesRecSig"), EtaXiDaughters);
469+
} else {
470+
registry.fill(HIST("hSelCandidatesRecBkg"), EtaXiDaughters);
471+
}
472+
} else {
473+
registry.fill(HIST("hSelCandidates"), EtaXiDaughters);
458474
}
459475

460476
// cut on pion pT

0 commit comments

Comments
 (0)