Skip to content

Commit e8aa7e0

Browse files
deependra170598deependra170598
authored
PWGHF: Change in number of pT Bins and selection criteria for d0SoftPiNormalised (#4336)
* Increasing pT Bins for Dstar with labels & changing inequality condition for d0SoftPiNormalised * Restoring d0SoftPiNormalised inequality condition
1 parent c0478bb commit e8aa7e0

2 files changed

Lines changed: 54 additions & 15 deletions

File tree

‎PWGHF/Core/SelectorCuts.h‎

Lines changed: 52 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -356,7 +356,7 @@ 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 = 10;
359+
static constexpr int nBinsPt = 23;
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
@@ -371,7 +371,20 @@ constexpr double binsPt[nBinsPt + 1] = {
371371
4.5,
372372
5.0,
373373
5.5,
374-
6.0};
374+
6.0,
375+
6.5,
376+
7.0,
377+
7.5,
378+
8.0,
379+
9.0,
380+
10.0,
381+
12.0,
382+
16.0,
383+
20.0,
384+
24.0,
385+
36.0,
386+
50.0,
387+
100.0};
375388
auto vecBinsPt = std::vector<double>{binsPt, binsPt + nBinsPt + 1};
376389

377390
// row labels
@@ -385,22 +398,48 @@ static const std::vector<std::string> labelsPt = {
385398
"pT bin 6",
386399
"pT bin 7",
387400
"pT bin 8",
388-
"pT bin 9"};
401+
"pT bin 9",
402+
"pT bin 10",
403+
"pT bin 11",
404+
"pT bin 12",
405+
"pT bin 13",
406+
"pT bin 14",
407+
"pT bin 15",
408+
"pT bin 16",
409+
"pT bin 17",
410+
"pT bin 18",
411+
"pT bin 19",
412+
"pT bin 20",
413+
"pT bin 21",
414+
"pT bin 22"};
389415

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

393419
// default values for the cuts
394-
constexpr double cuts[nBinsPt][nCutVars] = {{0.05, 0.3, 0.1, 0.5, 0.2, 300.0, 0.5, 0.5},
395-
{0.05, 0.3, 0.1, 0.5, 0.2, 300.0, 0.5, 0.5},
396-
{0.05, 0.4, 0.1, 0.5, 0.2, 300.0, 0.5, 0.5},
397-
{0.05, 0.4, 0.1, 0.5, 0.2, 300.0, 0.5, 0.5},
398-
{0.05, 0.6, 0.1, 0.5, 0.2, 300.0, 0.5, 0.5},
399-
{0.05, 0.6, 0.1, 0.5, 0.2, 300.0, 0.5, 0.5},
400-
{0.05, 0.6, 0.1, 0.5, 0.2, 300.0, 0.5, 0.5},
401-
{0.05, 100, 0.1, 0.5, 0.2, 300.0, 0.5, 0.5},
402-
{0.05, 100, 0.1, 0.5, 0.2, 300.0, 0.5, 0.5},
403-
{0.05, 100, 0.1, 0.5, 0.2, 300.0, 0.5, 0.5}};
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}};
404443
} // namespace hf_cuts_dstar_to_d0_pi
405444

406445
namespace hf_cuts_lc_to_p_k_pi

‎PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -59,8 +59,8 @@ struct HfTaskDstarToD0Pi {
5959
{
6060
auto vecPtBins = (std::vector<double>)ptBins;
6161

62-
registry.add("Yield/hDeltaInvMassDstar2D", "#Delta #it{M}_{inv} D* Candidate; inv. mass (#pi #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);
63-
registry.add("Yield/hDeltaInvMassDstar1D", "#Delta #it{M}_{inv} D* Candidate; inv. mass (#pi #pi k) (GeV/#it{c}^{2}); entries", {HistType::kTH1F, {{100, 0.13, 0.16}}}, true);
62+
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);
63+
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);
6464
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);
6565
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);
6666
// only QA

0 commit comments

Comments
 (0)