Skip to content

Commit 375831a

Browse files
deependra170598alibuild
deependra170598
and
alibuild
authored
PWGHF: Add more p t bin in selector cut (AliceO2Group#4540)
* Adding more pT bins for Dstar candidates and rebinning Decaylength and Impactparameter Hists * correction to binning option * Correcting default value of impact parameter cut of soft pi * Please consider the following formatting changes * correcting spelling mistake * fix spelling * MegaLinter fixes --------- Co-authored-by: ALICE Action Bot <[email protected]>
1 parent 0accdd3 commit 375831a

File tree

3 files changed

+55
-41
lines changed

3 files changed

+55
-41
lines changed

‎PWGHF/Core/SelectorCuts.h‎

Lines changed: 31 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -356,11 +356,13 @@ static const std::vector<std::string> labelsCutVar = {"m", "DCA", "cos theta*",
356356

357357
namespace hf_cuts_dstar_to_d0_pi
358358
{
359-
static constexpr int nBinsPt = 23;
359+
static constexpr int nBinsPt = 25;
360360
static constexpr int nCutVars = 8;
361361
// default values for the pT bin edges (can be used to configure histogram axis)
362362
// offset by 1 from the bin numbers in cuts array
363363
constexpr double binsPt[nBinsPt + 1] = {
364+
0.,
365+
0.5,
364366
1.0,
365367
1.5,
366368
2.0,
@@ -411,35 +413,39 @@ static const std::vector<std::string> labelsPt = {
411413
"pT bin 19",
412414
"pT bin 20",
413415
"pT bin 21",
414-
"pT bin 22"};
416+
"pT bin 22",
417+
"pT bin 23",
418+
"pT bin 24"};
415419

416420
// column label
417421
static const std::vector<std::string> labelsCutVar = {"ptSoftPiMin", "ptSoftPiMax", "d0SoftPi", "d0SoftPiNormalised", "deltaMInvDstar", "chi2PCA", "d0Prong0Normalised", "d0Prong1Normalised"};
418422

419423
// default values for the cuts
420-
constexpr double cuts[nBinsPt][nCutVars] = {{0.05, 0.3, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
421-
{0.05, 0.3, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
422-
{0.05, 0.4, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
423-
{0.05, 0.4, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
424-
{0.05, 0.6, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
425-
{0.05, 0.6, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
426-
{0.05, 0.6, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
427-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
428-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
429-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
430-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
431-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
432-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
433-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
434-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
435-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
436-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
437-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
438-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
439-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
440-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
441-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0},
442-
{0.05, 100, 0.1, 0.0, 0.2, 300.0, 0.0, 0.0}};
424+
constexpr double cuts[nBinsPt][nCutVars] = {{0.05, 0.2, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
425+
{0.05, 0.2, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
426+
{0.05, 0.3, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
427+
{0.05, 0.3, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
428+
{0.05, 0.4, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
429+
{0.05, 0.4, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
430+
{0.05, 0.6, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
431+
{0.05, 0.6, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
432+
{0.05, 0.6, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
433+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
434+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
435+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
436+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
437+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
438+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
439+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
440+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
441+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
442+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
443+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
444+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
445+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
446+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
447+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0},
448+
{0.05, 100, 0.1, 1000.0, 0.2, 300.0, 0.0, 0.0}};
443449
} // namespace hf_cuts_dstar_to_d0_pi
444450

445451
namespace hf_cuts_lc_to_p_k_pi

‎PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx‎

Lines changed: 23 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -45,6 +45,10 @@ struct HfTaskDstarToD0Pi {
4545
Partition<CandDstarWSelFlag> rowsSelectedCandDstar = aod::hf_sel_candidate_dstar::isSelDstarToD0Pi == selectionFlagDstarToD0Pi;
4646
Partition<CandDstarWSelFlagMcRec> rowsSelectedCandDstarMcRec = aod::hf_sel_candidate_dstar::isRecoD0Flag == selectionFlagHfD0ToPiK;
4747

48+
ConfigurableAxis binningImpactParam{"binningImpactParam", {1000, 0.1, -0.1}, " Bins of Impact Parameter"};
49+
ConfigurableAxis binningDecayLength{"binningDecayLength", {1000, 0.0, 0.7}, "Bins of Decay Length"};
50+
ConfigurableAxis binningNormDecayLength{"binningNormDecayLength", {1000, 0.0, 40.0}, "Bins of Normalised Decay Length"};
51+
4852
HistogramRegistry registry{
4953
"registry",
5054
{{"QA/hPtDstar", "Dstar Candidates; Dstar candidate #it{p}_{T} (GeV/#it{c}); entries", {HistType::kTH1F, {{360, 0., 36.}}}},
@@ -59,28 +63,32 @@ struct HfTaskDstarToD0Pi {
5963
{
6064
auto vecPtBins = (std::vector<double>)ptBins;
6165

66+
AxisSpec axisImpactParam = {binningImpactParam, "impact parameter (cm)"};
67+
AxisSpec axisDecayLength = {binningDecayLength, " decay length (cm)"};
68+
AxisSpec axisNormDecayLength = {binningNormDecayLength, "normalised decay length (cm)"};
69+
6270
registry.add("Yield/hDeltaInvMassDstar2D", "#Delta #it{M}_{inv} D* Candidate; inv. mass ((#pi #pi k) - (#pi k)) (GeV/#it{c}^{2});#it{p}_{T} (GeV/#it{c})", {HistType::kTH2F, {{100, 0.13, 0.16}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}}, true);
6371
registry.add("Yield/hDeltaInvMassDstar1D", "#Delta #it{M}_{inv} D* Candidate; inv. mass ((#pi #pi k) - (#pi k)) (GeV/#it{c}^{2}); entries", {HistType::kTH1F, {{100, 0.13, 0.16}}}, true);
6472
registry.add("Yield/hInvMassDstar", "#Delta #it{M}_{inv} D* Candidate; inv. mass (#pi #pi k) (GeV/#it{c}^{2}); entries", {HistType::kTH1F, {{500, 0., 5.0}}}, true);
6573
registry.add("Yield/hInvMassD0", "#it{M}_{inv}D^{0} candidate;#it{M}_{inv} D^{0} (GeV/#it{c});#it{p}_{T} (GeV/#it{c})", {HistType::kTH2F, {{500, 0., 5.0}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}}, true);
6674
// only QA
6775
registry.add("QA/hEtaDstar", "D* Candidate; D* Candidate #it{#eta};entries", {HistType::kTH2F, {{100, -2., 2.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
68-
registry.add("QA/hCtD0", "D0 Candidate; D0 Candidate's proper life time #it{c}#tau (cm) ; entries ", {HistType::kTH2F, {{120, -20., 100.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
69-
registry.add("QA/hDecayLengthD0", "D0 Candidate; D0 Candidate's decay length (cm); entries", {HistType::kTH2F, {{800, 0., 4.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
70-
registry.add("QA/hDecayLengthXYD0", "D0 Candidate; D0 Candidate's decay length xy (cm); entries", {HistType::kTH2F, {{800, 0., 4.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
71-
registry.add("QA/hDecayLengthNormalisedD0", "D0 Candidates;Do Candidate's normalised decay length (cm); entries", {HistType::kTH2F, {{800, 0., 40.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
72-
registry.add("QA/hDecayLengthXYNormalisedD0", "D0 candidate; D0 Candidate's normalised decay length xy (cm); entries", {HistType::kTH2F, {{800, 0., 40.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
76+
registry.add("QA/hCtD0", "D0 Candidate; D0 Candidate's proper life time #it{c}#tau (cm) ; entries ", {HistType::kTH2F, {{1000, -0.1, 14.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
77+
registry.add("QA/hDecayLengthD0", "D0 Candidate; D0 Candidate's decay length (cm); entries", {HistType::kTH2F, {axisDecayLength, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
78+
registry.add("QA/hDecayLengthXYD0", "D0 Candidate; D0 Candidate's decay length xy (cm); entries", {HistType::kTH2F, {axisDecayLength, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
79+
registry.add("QA/hDecayLengthNormalisedD0", "D0 Candidates;Do Candidate's normalised decay length (cm); entries", {HistType::kTH2F, {axisNormDecayLength, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
80+
registry.add("QA/hDecayLengthXYNormalisedD0", "D0 candidate; D0 Candidate's normalised decay length xy (cm); entries", {HistType::kTH2F, {axisNormDecayLength, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
7381
registry.add("QA/hCPAD0", "D0 Candidates; D0 Candidate's cosine pointing angle; entries", {HistType::kTH2F, {{110, -1., 1.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
7482
registry.add("QA/hCPAxyD0", "D0 candidates; D0 Candidate's cosine of pointing angle xy; entries", {HistType::kTH2F, {{110, -1., 1.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
75-
registry.add("QA/hImpactParameterXYD0", "D0 Candidates; D0 Candidate's reconstructed impact parameter xy (cm); entries", {HistType::kTH2F, {{200, -1., 1.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
76-
registry.add("QA/hDeltaIPMaxNormalisedD0", "D0 Candidate; D0 Candidate's Norm. Delta IP; entries", {HistType::kTH2F, {{200, -20., 20.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
77-
registry.add("QA/hSqSumProngesImpactParameterD0", "D0 Candidates; Sqr Sum of Impact params of D0 Pronges; enteries ", {HistType::kTH2F, {{100, 0., 1.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
78-
registry.add("QA/hDecayLengthErrorD0", "D0 Candidates; D0 Candidate's Decay Lenth Error (cm); entries", {HistType::kTH2F, {{100, 0., 1.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
79-
registry.add("QA/hDecayLengthXYErrorD0", "D0 Candidates; D0 Candidate's Decay Length Error XY (cm); entries", {HistType::kTH2F, {{100, 0., 1.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
80-
registry.add("QA/hImpactParameterError", "D0 Pronges; Impactparam error of different D0 Pronges (cm); entries", {HistType::kTH2F, {{100, 0., 1.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
81-
registry.add("QA/hd0Prong0", "Prong0; DCAxy of Prong0 (cm); entries", {HistType::kTH2F, {{100, -1., 1.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
82-
registry.add("QA/hd0Prong1", "Prong1; DCAxy of Prong1 (cm); entries", {HistType::kTH2F, {{100, -1., 1.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
83-
registry.add("QA/hd0ProngSoftPi", "ProngSoftPi; DCAxy of Prong Soft Pi (cm); entries", {HistType::kTH2F, {{100, -1., 1.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
83+
registry.add("QA/hImpactParameterXYD0", "D0 Candidates; D0 Candidate's reconstructed impact parameter xy (cm); entries", {HistType::kTH2F, {axisImpactParam, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
84+
registry.add("QA/hDeltaIPMaxNormalisedD0", "D0 Candidate; D0 Candidate's Norm. Delta IP; entries", {HistType::kTH2F, {{1000, -20., 20.}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
85+
registry.add("QA/hSqSumProngsImpactParameterD0", "D0 Candidates; Sqr Sum of Impact params of D0 Prongs; entries ", {HistType::kTH2F, {{1000, 0., 0.25}, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
86+
registry.add("QA/hDecayLengthErrorD0", "D0 Candidates; D0 Candidate's Decay Length Error (cm); entries", {HistType::kTH2F, {axisDecayLength, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
87+
registry.add("QA/hDecayLengthXYErrorD0", "D0 Candidates; D0 Candidate's Decay Length Error XY (cm); entries", {HistType::kTH2F, {axisDecayLength, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
88+
registry.add("QA/hImpactParameterError", "D0 Prongs; Impact param error of different D0 Prongs (cm); entries", {HistType::kTH2F, {axisImpactParam, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
89+
registry.add("QA/hd0Prong0", "Prong0; DCAxy of Prong0 (cm); entries", {HistType::kTH2F, {axisImpactParam, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
90+
registry.add("QA/hd0Prong1", "Prong1; DCAxy of Prong1 (cm); entries", {HistType::kTH2F, {axisImpactParam, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
91+
registry.add("QA/hd0ProngSoftPi", "ProngSoftPi; DCAxy of Prong Soft Pi (cm); entries", {HistType::kTH2F, {axisImpactParam, {vecPtBins, "#it{p}_{T} (GeV/#it{c})"}}});
8492
// MC Matching at Reconstruction Level Successful
8593
registry.add("QA/hCPASkimD0RecSig", "MC Matched Skimed D* Candidates at Reconstruction Level; cosine of pointing angle", {HistType::kTH1F, {{110, -1.1, 1.1}}});
8694
registry.add("QA/hEtaSkimD0RecSig", "MC Matched Skimed D* Candidates at Reconstruction Level; #it{#eta} of D0 Prong", {HistType::kTH1F, {{100, -2., 2.}}});
@@ -139,7 +147,7 @@ struct HfTaskDstarToD0Pi {
139147
registry.fill(HIST("QA/hCPAxyD0"), candDstar.cpaXYD0(), candDstar.pt());
140148
registry.fill(HIST("QA/hImpactParameterXYD0"), candDstar.impactParameterXYD0(), candDstar.pt());
141149
registry.fill(HIST("QA/hDeltaIPMaxNormalisedD0"), candDstar.deltaIPNormalisedMaxD0(), candDstar.pt());
142-
registry.fill(HIST("QA/hSqSumProngesImpactParameterD0"), candDstar.impactParameterProngSqSumD0(), candDstar.pt());
150+
registry.fill(HIST("QA/hSqSumProngsImpactParameterD0"), candDstar.impactParameterProngSqSumD0(), candDstar.pt());
143151
registry.fill(HIST("QA/hDecayLengthErrorD0"), candDstar.errorDecayLengthD0(), candDstar.pt());
144152
registry.fill(HIST("QA/hDecayLengthXYErrorD0"), candDstar.errorDecayLengthXYD0(), candDstar.pt());
145153
registry.fill(HIST("QA/hImpactParameterError"), candDstar.errorImpactParameter0(), candDstar.pt());

‎PWGHF/TableProducer/candidateSelectorDstarToD0Pi.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -137,7 +137,7 @@ struct HfCandidateSelectorDstarToD0Pi {
137137
if (candidate.chi2PCAD0() > cutsDstar->get(binPt, "chi2PCA")) {
138138
return false;
139139
}
140-
if (candidate.impactParameterXYD0() > cutsD0->get(binPt, "DCA")) {
140+
if (std::abs(candidate.impactParameterXYD0()) > cutsD0->get(binPt, "DCA")) {
141141
return false;
142142
}
143143
// d0Prong0Normalised,1

0 commit comments

Comments
 (0)