|
52 | 52 | #include <THnSparse.h> |
53 | 53 | #include <TString.h> |
54 | 54 |
|
| 55 | +#include <array> |
55 | 56 | #include <chrono> |
56 | 57 | #include <cstddef> |
57 | 58 | #include <cstdint> |
@@ -87,20 +88,37 @@ const TString stringMcGenLcFd = "MC gen, non-prompt #Lambda_c;"; |
87 | 88 |
|
88 | 89 | // definition of vectors for standard ptbin and invariant mass configurables |
89 | 90 | const int nPtBinsCorrelations = 8; |
90 | | -const double pTBinsCorrelations[nPtBinsCorrelations + 1] = {0., 2., 4., 6., 8., 12., 16., 24., 99.}; |
91 | | -const auto vecBinsPtCorrelations = std::vector<double>{pTBinsCorrelations, pTBinsCorrelations + nPtBinsCorrelations + 1}; |
92 | | -const double signalRegionInnerDefault[nPtBinsCorrelations] = {2.269, 2.269, 2.269, 2.269, 2.269, 2.269, 2.269, 2.269}; |
93 | | -const double signalRegionOuterDefault[nPtBinsCorrelations] = {2.309, 2.309, 2.309, 2.309, 2.309, 2.309, 2.309, 2.309}; |
94 | | -const double sidebandLeftOuterDefault[nPtBinsCorrelations] = {2.209, 2.209, 2.209, 2.209, 2.209, 2.209, 2.209, 2.209}; |
95 | | -const double sidebandLeftInnerDefault[nPtBinsCorrelations] = {2.249, 2.249, 2.249, 2.249, 2.249, 2.249, 2.249, 2.249}; |
96 | | -const double sidebandRightInnerDefault[nPtBinsCorrelations] = {2.329, 2.329, 2.329, 2.329, 2.329, 2.329, 2.329, 2.329}; |
97 | | -const double sidebandRightOuterDefault[nPtBinsCorrelations] = {2.369, 2.369, 2.369, 2.369, 2.369, 2.369, 2.369, 2.369}; |
98 | | -const auto vecSignalRegionInner = std::vector<double>{signalRegionInnerDefault, signalRegionInnerDefault + nPtBinsCorrelations}; |
99 | | -const auto vecSignalRegionOuter = std::vector<double>{signalRegionOuterDefault, signalRegionOuterDefault + nPtBinsCorrelations}; |
100 | | -const auto vecSidebandLeftInner = std::vector<double>{sidebandLeftInnerDefault, sidebandLeftInnerDefault + nPtBinsCorrelations}; |
101 | | -const auto vecSidebandLeftOuter = std::vector<double>{sidebandLeftOuterDefault, sidebandLeftOuterDefault + nPtBinsCorrelations}; |
102 | | -const auto vecSidebandRightInner = std::vector<double>{sidebandRightInnerDefault, sidebandRightInnerDefault + nPtBinsCorrelations}; |
103 | | -const auto vecSidebandRightOuter = std::vector<double>{sidebandRightOuterDefault, sidebandRightOuterDefault + nPtBinsCorrelations}; |
| 91 | + |
| 92 | +const std::array<double, nPtBinsCorrelations + 1> pTBinsCorrelations = { |
| 93 | + 0., 2., 4., 6., 8., 12., 16., 24., 99.}; |
| 94 | +const auto vecBinsPtCorrelations = |
| 95 | + std::vector<double>{pTBinsCorrelations.begin(), pTBinsCorrelations.end()}; |
| 96 | + |
| 97 | +const std::array<double, nPtBinsCorrelations> signalRegionInnerDefault = { |
| 98 | + 2.269, 2.269, 2.269, 2.269, 2.269, 2.269, 2.269, 2.269}; |
| 99 | +const std::array<double, nPtBinsCorrelations> signalRegionOuterDefault = { |
| 100 | + 2.309, 2.309, 2.309, 2.309, 2.309, 2.309, 2.309, 2.309}; |
| 101 | +const std::array<double, nPtBinsCorrelations> sidebandLeftOuterDefault = { |
| 102 | + 2.209, 2.209, 2.209, 2.209, 2.209, 2.209, 2.209, 2.209}; |
| 103 | +const std::array<double, nPtBinsCorrelations> sidebandLeftInnerDefault = { |
| 104 | + 2.249, 2.249, 2.249, 2.249, 2.249, 2.249, 2.249, 2.249}; |
| 105 | +const std::array<double, nPtBinsCorrelations> sidebandRightInnerDefault = { |
| 106 | + 2.329, 2.329, 2.329, 2.329, 2.329, 2.329, 2.329, 2.329}; |
| 107 | +const std::array<double, nPtBinsCorrelations> sidebandRightOuterDefault = { |
| 108 | + 2.369, 2.369, 2.369, 2.369, 2.369, 2.369, 2.369, 2.369}; |
| 109 | + |
| 110 | +const auto vecSignalRegionInner = |
| 111 | + std::vector<double>{signalRegionInnerDefault.begin(), signalRegionInnerDefault.end()}; |
| 112 | +const auto vecSignalRegionOuter = |
| 113 | + std::vector<double>{signalRegionOuterDefault.begin(), signalRegionOuterDefault.end()}; |
| 114 | +const auto vecSidebandLeftInner = |
| 115 | + std::vector<double>{sidebandLeftInnerDefault.begin(), sidebandLeftInnerDefault.end()}; |
| 116 | +const auto vecSidebandLeftOuter = |
| 117 | + std::vector<double>{sidebandLeftOuterDefault.begin(), sidebandLeftOuterDefault.end()}; |
| 118 | +const auto vecSidebandRightInner = |
| 119 | + std::vector<double>{sidebandRightInnerDefault.begin(), sidebandRightInnerDefault.end()}; |
| 120 | +const auto vecSidebandRightOuter = |
| 121 | + std::vector<double>{sidebandRightOuterDefault.begin(), sidebandRightOuterDefault.end()}; |
104 | 122 |
|
105 | 123 | /// Lc-Hadron correlation pair filling task, from pair tables - for real data and data-like analysis (i.e. reco-level w/o matching request via Mc truth) |
106 | 124 | struct HfTaskCorrelationLcHadrons { |
@@ -405,7 +423,13 @@ struct HfTaskCorrelationLcHadrons { |
405 | 423 |
|
406 | 424 | if (activateQA) { |
407 | 425 | const int regionLimits = 6; |
408 | | - std::string labels[regionLimits] = {"SigReg Left", "SigReg Right", "Left SB Low", "Left SB Up", "Right SB Low", "Right SB Up"}; |
| 426 | + const std::array<std::string, regionLimits> labels = { |
| 427 | + "SigReg Left", |
| 428 | + "SigReg Right", |
| 429 | + "Left SB Low", |
| 430 | + "Left SB Up", |
| 431 | + "Right SB Low", |
| 432 | + "Right SB Up"}; |
409 | 433 | static const AxisSpec axisSidebandLimits = {regionLimits, 0.5, 6.5, ""}; |
410 | 434 | auto hSigSidebandLimits = registry.add<TH2>("Inputs/hSigSidebandLimits", "Signal and Sideband Limits;;#it{p}_{T} (GeV/#it{c})", {HistType::kTH2F, {axisSidebandLimits, {(std::vector<double>)binsPtCorrelations, "#it{p}_{T} (GeV/#it{c})"}}}); |
411 | 435 | for (int iLim = 0; iLim < regionLimits; iLim++) { |
@@ -521,7 +545,7 @@ struct HfTaskCorrelationLcHadrons { |
521 | 545 | registry.fill(HIST("hCorrel2DVsPtSignMass"), deltaPhi, deltaEta, ptLc, ptHadron, massLc, signPair, poolBin, efficiencyWeight); |
522 | 546 | } |
523 | 547 | // check if correlation entry belongs to signal region, sidebands or is outside both, and fill correlation plots |
524 | | - if (storeMass) { |
| 548 | + if (storeMass != 0) { |
525 | 549 | registry.fill(HIST("hCorrel2DVsPtGlobalRegion"), deltaPhi, deltaEta, ptLc, ptHadron, poolBin, massLc, efficiencyWeight); |
526 | 550 | continue; |
527 | 551 | } |
|
0 commit comments