Skip to content

Commit 63f0c62

Browse files
author
Marcello Di Costanzo
committed
Fixes to templetization
1 parent 6e87a8d commit 63f0c62

File tree

7 files changed

+74
-63
lines changed

7 files changed

+74
-63
lines changed

‎ALICE3/DataModel/A3DecayFinderTables.h‎

Lines changed: 8 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -393,17 +393,21 @@ DECLARE_SOA_TABLE(Alice3Sel3Ps, "AOD", "ALICE3SEL3P", //!
393393

394394
namespace a3_mc_truth
395395
{
396-
DECLARE_SOA_COLUMN(OriginMcRec, originMcRec, int);
397-
DECLARE_SOA_COLUMN(FlagMcRec, flagMcRec, int);
398-
DECLARE_SOA_COLUMN(OriginMcGen, originMcGen, int);
399-
DECLARE_SOA_COLUMN(FlagMcGen, flagMcGen, int);
396+
DECLARE_SOA_COLUMN(OriginMcRec, originMcRec, int); //!
397+
DECLARE_SOA_COLUMN(BHadMotherPtRec, bHadMotherPtRec, float); //!
398+
DECLARE_SOA_COLUMN(FlagMcRec, flagMcRec, int); //!
399+
DECLARE_SOA_COLUMN(OriginMcGen, originMcGen, int); //!
400+
DECLARE_SOA_COLUMN(BHadMotherPtGen, bHadMotherPtGen, float); //!
401+
DECLARE_SOA_COLUMN(FlagMcGen, flagMcGen, int); //!
400402
} // namespace a3_mc_truth
401403
DECLARE_SOA_TABLE(Alice3McRecFlags, "AOD", "ALICE3MCRECFLAG", //!
402404
a3_mc_truth::OriginMcRec,
405+
a3_mc_truth::BHadMotherPtRec,
403406
a3_mc_truth::FlagMcRec);
404407

405408
DECLARE_SOA_TABLE(Alice3McGenFlags, "AOD", "ALICE3MCGENFLAG", //!
406409
a3_mc_truth::OriginMcGen,
410+
a3_mc_truth::BHadMotherPtGen,
407411
a3_mc_truth::FlagMcGen);
408412

409413
DECLARE_SOA_TABLE(Alice3PidLcs, "AOD", "ALICE3PIDLC", //!

‎ALICE3/ML/alice3-mlresponse3prong.h‎

Lines changed: 20 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -39,9 +39,9 @@
3939
// Specific case of CHECK_AND_FILL_ML_ALICE3_FULL(OBJECT, FEATURE, GETTER)
4040
// where OBJECT is named candidate and FEATURE = GETTER
4141
#define CHECK_AND_FILL_ML_ALICE3(GETTER) \
42-
case static_cast<uint8_t>(InputFeatures3Prong::GETTER): { \
43-
inputFeatures.emplace_back(candidate.GETTER()); \
44-
break; \
42+
case static_cast<uint8_t>(InputFeatures3Prong::GETTER): { \
43+
inputFeatures.emplace_back(candidate.GETTER()); \
44+
break; \
4545
}
4646

4747
namespace o2::analysis
@@ -116,7 +116,7 @@ class Alice3MlResponse3Prong : public MlResponse<TypeOutputScore>
116116
/// \param prong2 is the candidate's prong2
117117
/// \return inputFeatures vector
118118
template <typename T1>
119-
std::vector<float> getInputFeatures(T1 const& candidate, int pdg)
119+
std::vector<float> getInputFeatures(T1 const& candidate)
120120
{
121121
std::vector<float> inputFeatures;
122122

@@ -137,27 +137,22 @@ class Alice3MlResponse3Prong : public MlResponse<TypeOutputScore>
137137
CHECK_AND_FILL_ML_ALICE3(cpa);
138138
CHECK_AND_FILL_ML_ALICE3(cpaXY);
139139
CHECK_AND_FILL_ML_ALICE3(chi2PCA);
140-
switch (pdg) {
141-
case o2::constants::physics::Pdg::kLambdaCPlus: {
142-
// RICH PID variables
143-
CHECK_AND_FILL_ML_ALICE3(nSigRichPr0);
144-
CHECK_AND_FILL_ML_ALICE3(nSigRichKa1);
145-
CHECK_AND_FILL_ML_ALICE3(nSigRichPi2);
146-
// INNER TOF PID variables
147-
CHECK_AND_FILL_ML_ALICE3(nSigInnTofPr0);
148-
CHECK_AND_FILL_ML_ALICE3(nSigInnTofKa1);
149-
CHECK_AND_FILL_ML_ALICE3(nSigInnTofPi2);
150-
// OUTER TOF PID variables
151-
CHECK_AND_FILL_ML_ALICE3(nSigOutTofPr0);
152-
CHECK_AND_FILL_ML_ALICE3(nSigOutTofKa1);
153-
CHECK_AND_FILL_ML_ALICE3(nSigOutTofPi2);
154-
// TRACKER PID variables
155-
CHECK_AND_FILL_ML_ALICE3(nSigTrkPr0);
156-
CHECK_AND_FILL_ML_ALICE3(nSigTrkKa1);
157-
CHECK_AND_FILL_ML_ALICE3(nSigTrkPi2);
158-
break;
159-
}
160-
}
140+
// TRACKER PID variables
141+
CHECK_AND_FILL_ML_ALICE3(nSigTrkPr0);
142+
CHECK_AND_FILL_ML_ALICE3(nSigTrkKa1);
143+
CHECK_AND_FILL_ML_ALICE3(nSigTrkPi2);
144+
// RICH PID variables
145+
CHECK_AND_FILL_ML_ALICE3(nSigRichPr0);
146+
CHECK_AND_FILL_ML_ALICE3(nSigRichKa1);
147+
CHECK_AND_FILL_ML_ALICE3(nSigRichPi2);
148+
// INNER TOF PID variables
149+
CHECK_AND_FILL_ML_ALICE3(nSigInnTofPr0);
150+
CHECK_AND_FILL_ML_ALICE3(nSigInnTofKa1);
151+
CHECK_AND_FILL_ML_ALICE3(nSigInnTofPi2);
152+
// OUTER TOF PID variables
153+
CHECK_AND_FILL_ML_ALICE3(nSigOutTofPr0);
154+
CHECK_AND_FILL_ML_ALICE3(nSigOutTofKa1);
155+
CHECK_AND_FILL_ML_ALICE3(nSigOutTofPi2);
161156
}
162157
}
163158
return inputFeatures;

‎ALICE3/TableProducer/alice3-decayfinder.cxx‎

Lines changed: 26 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -48,8 +48,10 @@
4848
#include "ALICE3/DataModel/OTFRICH.h"
4949
#include "ALICE3/DataModel/RICH.h"
5050
#include "ALICE3/DataModel/A3DecayFinderTables.h"
51+
#include "PWGHF/Utils/utilsAnalysis.h"
5152

5253
using namespace o2;
54+
using namespace o2::analysis;
5355
using namespace o2::framework;
5456
using namespace o2::constants::physics;
5557
using namespace o2::framework::expressions;
@@ -139,10 +141,11 @@ struct alice3decayFinder {
139141
o2::vertexing::DCAFitterN<3> fitter3;
140142

141143
double bz{0.};
142-
const float toMicrometers = 10000.; // from cm to µm
143-
std::array<int, 3> daugsPdgCodes3Prong = {-1, -1, -1};
144-
std::array<float, 3> daughtersMasses3Prong = {-1.f, -1.f, -1.f};
145-
int motherPdgCode = -1;
144+
const float toMicrometers{10000.}; // from cm to µm
145+
std::array<int, 3> daugsPdgCodes3Prong{{-1, -1, -1}};
146+
std::array<float, 3> daughtersMasses3Prong{{-1.f, -1.f, -1.f}};
147+
int motherPdgCode{-1};
148+
int charmHadFlag{0};
146149

147150
HistogramRegistry histos{"histos", {}, OutputObjHandlingPolicy::AnalysisObject};
148151

@@ -241,6 +244,7 @@ struct alice3decayFinder {
241244
float chi2PCA; // normalized 3D decay length
242245
int flagMc; // 0 = bkg, pdg code for signal
243246
int origin; // 1 = prompt, 2 = non-prompt
247+
float ptBMotherRec; // pT of the B hadron mother (reconstructed)
244248
} cand3prong;
245249

246250
template <typename TTrackType>
@@ -417,7 +421,7 @@ struct alice3decayFinder {
417421
if (indexRec < 0) {
418422
cand3prong.flagMc = 0; // bkg
419423
} else {
420-
cand3prong.flagMc = motherPart.pdgCode(); // Particle
424+
cand3prong.flagMc = motherPart.pdgCode() > 0 ? charmHadFlag : -charmHadFlag; // Particle
421425
}
422426

423427
cand3prong.origin = 0;
@@ -428,6 +432,11 @@ struct alice3decayFinder {
428432
int origin = RecoDecay::getCharmHadronOrigin(mcParticles, motherParticle, false, &idxBhadMothers);
429433
LOG(info) << "Origin: " << origin;
430434
cand3prong.origin = origin;
435+
cand3prong.ptBMotherRec = -1.f;
436+
if (origin == RecoDecay::OriginType::NonPrompt) {
437+
auto bHadMother = mcParticles.rawIteratorAt(idxBhadMothers[0]);
438+
cand3prong.ptBMotherRec = bHadMother.pt();
439+
}
431440
}
432441
return true;
433442
}
@@ -590,6 +599,7 @@ struct alice3decayFinder {
590599
daughtersMasses3Prong = {o2::constants::physics::MassProton,
591600
o2::constants::physics::MassKaonCharged,
592601
o2::constants::physics::MassPionCharged};
602+
charmHadFlag = CharmHadAlice3::Lc;
593603
}
594604
}
595605

@@ -637,18 +647,23 @@ struct alice3decayFinder {
637647
LOG(info) << "Processing generated MC particles: total number = " << mcParticles.size();
638648
for (auto const& mcParticle : mcParticles) {
639649
if (std::abs(mcParticle.pdgCode()) != motherPdgCode) {
640-
mcGenFlags(-1, -1);
650+
mcGenFlags(-1, -1, -1);
641651
continue;
642652
}
643653
std::vector<int> idxBhadMothers{};
644654
int origin = RecoDecay::getCharmHadronOrigin(mcParticles, mcParticle, false, &idxBhadMothers);
645-
mcGenFlags(origin, mcParticle.pdgCode());
655+
float ptBMotherGen{-1.f};
656+
if (origin == RecoDecay::OriginType::NonPrompt) {
657+
auto bHadMother = mcParticles.rawIteratorAt(idxBhadMothers[0]);
658+
ptBMotherGen = bHadMother.pt();
659+
}
660+
mcGenFlags(origin, ptBMotherGen, mcParticle.pdgCode() ? charmHadFlag : -charmHadFlag);
646661
if (mcParticle.pdgCode() > 0) {
647-
LOG(info) << "[P] Origin: " << origin << " PDG: " << mcParticle.pdgCode();
662+
LOG(info) << "[P] Origin: " << origin << " PDG: " << charmHadFlag << " " << mcParticle.pdgCode();
648663
histos.fill(HIST("h2dGen3Prong"), mcParticle.pt(), mcParticle.eta());
649664
} else {
650665
histos.fill(HIST("h2dGen3ProngBar"), mcParticle.pt(), mcParticle.eta());
651-
LOG(info) << "[AP] Origin: " << origin << " PDG: " << mcParticle.pdgCode();
666+
LOG(info) << "[AP] Origin: " << origin << " PDG: " << -charmHadFlag << " " << mcParticle.pdgCode();
652667
}
653668
}
654669
}
@@ -1029,7 +1044,7 @@ struct alice3decayFinder {
10291044
false, // is swapped hypothesis
10301045
candPx, candPy, candPz
10311046
);
1032-
mcRecFlags(cand3prong.origin, cand3prong.flagMc); // placeholder for prompt/non-prompt
1047+
mcRecFlags(cand3prong.origin, cand3prong.ptBMotherRec, cand3prong.flagMc); // placeholder for prompt/non-prompt
10331048
fillPidTable(prong0, prong1, prong2);
10341049

10351050
if constexpr (FillSwapHypo) {
@@ -1057,7 +1072,7 @@ struct alice3decayFinder {
10571072
true, // is swapped hypothesis
10581073
candPx, candPy, candPz
10591074
);
1060-
mcRecFlags(cand3prong.origin, cand3prong.flagMc); // placeholder for prompt/non-prompt
1075+
mcRecFlags(cand3prong.origin, cand3prong.ptBMotherRec, cand3prong.flagMc); // placeholder for prompt/non-prompt
10611076
fillPidTable(prong2, prong1, prong0);
10621077
}
10631078
}

‎ALICE3/TableProducer/alice3-selector3prong.cxx‎

Lines changed: 11 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -156,9 +156,6 @@ struct Alice3Selector3Prong {
156156
bool selectionTopol(const T& cand, float candPt)
157157
{
158158
int const ptBin = findBin(binsPt, candPt);
159-
if (ptBin == -1) {
160-
return false;
161-
}
162159

163160
// check that the cand pT is within the analysis range
164161
if (candPt < ptCandMin || candPt >= ptCandMax) {
@@ -209,11 +206,11 @@ struct Alice3Selector3Prong {
209206
return true;
210207
}
211208

212-
template <typename TCandidate>
209+
template <CharmHadAlice3 CharmHad, typename TCandidate>
213210
bool selectionCandidateMass(int const ptBin, const TCandidate& cand)
214211
{
215212
float massCand{0.f};
216-
if (std::abs(cand.flagMcRec()) == o2::constants::physics::Pdg::kLambdaCPlus) {
213+
if constexpr (CharmHad == CharmHadAlice3::Lc) {
217214
if (cand.isSwapped()) {
218215
massCand = hfHelper.invMassLcToPiKP(cand);
219216
} else {
@@ -244,7 +241,7 @@ struct Alice3Selector3Prong {
244241
/// \param pidTrack0 is the PID status of proton cand track
245242
/// \param pidTrack1 is the PID status of kaon cand track
246243
/// \param pidTrack2 is the PID status of pion cand track
247-
template<typename TCand>
244+
template<CharmHadAlice3 CharmHad, typename TCand>
248245
void configurePidMask(const TCand& cand, uint32_t& pidMask)
249246
{
250247

@@ -257,7 +254,7 @@ struct Alice3Selector3Prong {
257254

258255
// prong 0
259256
float ptProng0{cand.ptProng0()};
260-
if (std::abs(cand.flagMcRec()) == o2::constants::physics::Pdg::kLambdaCPlus) {
257+
if constexpr (CharmHad == CharmHadAlice3::Lc) {
261258
isSelPid(cand.nSigTrkPr0(), PidSels::TrkProng0, ptProng0, nSigmaTrkMax, ptPidTrkMin, ptPidTrkMax);
262259
isSelPid(cand.nSigRichPr0(), PidSels::RichProng0, ptProng0, nSigmaRichMax, ptPidRichMin, ptPidRichMax);
263260
isSelPid(cand.nSigInnTofPr0(), PidSels::InnTofProng0, ptProng0, nSigmaInnTofMax, ptPidInnTofMin, ptPidInnTofMax);
@@ -266,7 +263,7 @@ struct Alice3Selector3Prong {
266263

267264
// prong 1
268265
float ptProng1{cand.ptProng1()};
269-
if (std::abs(cand.flagMcRec()) == o2::constants::physics::Pdg::kLambdaCPlus) {
266+
if constexpr (CharmHad == CharmHadAlice3::Lc) {
270267
isSelPid(cand.nSigTrkKa1(), PidSels::TrkProng1, ptProng1, nSigmaTrkMax, ptPidTrkMin, ptPidTrkMax);
271268
isSelPid(cand.nSigRichKa1(), PidSels::RichProng1, ptProng1, nSigmaRichMax, ptPidRichMin, ptPidRichMax);
272269
isSelPid(cand.nSigInnTofKa1(), PidSels::InnTofProng1, ptProng1, nSigmaInnTofMax, ptPidInnTofMin, ptPidInnTofMax);
@@ -275,7 +272,7 @@ struct Alice3Selector3Prong {
275272

276273
// prong 2
277274
float ptProng2{cand.ptProng2()};
278-
if (std::abs(cand.flagMcRec()) == o2::constants::physics::Pdg::kLambdaCPlus) {
275+
if constexpr (CharmHad == CharmHadAlice3::Lc) {
279276
isSelPid(cand.nSigTrkPi2(), PidSels::TrkProng2, ptProng2, nSigmaTrkMax, ptPidTrkMin, ptPidTrkMax);
280277
isSelPid(cand.nSigRichPi2(), PidSels::RichProng2, ptProng2, nSigmaRichMax, ptPidRichMin, ptPidRichMax);
281278
isSelPid(cand.nSigInnTofPi2(), PidSels::InnTofProng2, ptProng2, nSigmaInnTofMax, ptPidInnTofMin, ptPidInnTofMax);
@@ -288,7 +285,7 @@ struct Alice3Selector3Prong {
288285
/// \brief function to apply Lc selections
289286
/// \param cands Lc cand table
290287
/// \param tracks track table
291-
template <typename CandType>
288+
template <CharmHadAlice3 CharmHad, typename CandType>
292289
void runSelect3Prong(CandType const& cands)
293290
{
294291
bool isSel = false;
@@ -312,7 +309,7 @@ struct Alice3Selector3Prong {
312309
}
313310

314311
// Here all cands pass the cut on the mass selection
315-
if (!selectionCandidateMass(ptBin, cand)) {
312+
if (!selectionCandidateMass<CharmHad>(ptBin, cand)) {
316313
candSelFlags(isSel, pidMask);
317314
if (applyMl) {
318315
candMlScores(outputMl[0], outputMl[1], outputMl[2]);
@@ -336,7 +333,7 @@ struct Alice3Selector3Prong {
336333
}
337334

338335
// PID selection
339-
configurePidMask(cand, pidMask);
336+
configurePidMask<CharmHad>(cand, pidMask);
340337
if (pidMask == 0) {
341338
candSelFlags(isSel, pidMask);
342339
if (applyMl) {
@@ -352,7 +349,7 @@ struct Alice3Selector3Prong {
352349
// ML selections
353350
if (applyMl) {
354351

355-
std::vector<float> inputFeaturesMassHypo0 = mlResponse.getInputFeatures(cand, std::abs(cand.flagMcRec()));
352+
std::vector<float> inputFeaturesMassHypo0 = mlResponse.getInputFeatures(cand);
356353
isSelectedMl = mlResponse.isSelectedMl(inputFeaturesMassHypo0, ptCand, outputMl);
357354
candMlScores(outputMl[0], outputMl[1], outputMl[2]);
358355
if (!isSelectedMl) {
@@ -376,7 +373,7 @@ struct Alice3Selector3Prong {
376373
void processLc(CandsLcWMcTruth const& cands)
377374
{
378375
LOG(info) << "Starting Lc 3-prong cand selection: " << cands.size() << " cands to be processed.";
379-
runSelect3Prong(cands);
376+
runSelect3Prong<CharmHadAlice3::Lc>(cands);
380377
}
381378
PROCESS_SWITCH(Alice3Selector3Prong, processLc, "Process 3 prong selection for Lc", true);
382379
};

‎ALICE3/TableProducer/alice3-treecreator3prong.cxx‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -288,11 +288,11 @@ struct Alice3TreeCreator3Prong {
288288
using CandsMcGen = soa::Join<aod::McParticles, aod::Alice3McGenFlags>;
289289

290290
Filter filterSelectCandidates = aod::a3_hf_sel_3prong::isSel == true;
291-
Filter filterSelectGenCands = nabs(aod::a3_mc_truth::flagMcGen) == static_cast<int>(o2::constants::physics::Pdg::kLambdaCPlus);
291+
Filter filterSelectGenCands = nabs(aod::a3_mc_truth::flagMcGen) == static_cast<int>(CharmHadAlice3::Lc);
292292

293-
Partition<CandsLcRec> recoLcCandSig = nabs(o2::aod::a3_mc_truth::flagMcRec) == static_cast<int>(o2::constants::physics::Pdg::kLambdaCPlus);
293+
Partition<CandsLcRec> recoLcCandSig = nabs(o2::aod::a3_mc_truth::flagMcRec) == static_cast<int>(CharmHadAlice3::Lc);
294294
Partition<CandsLcRec> recoLcCandBkg = nabs(o2::aod::a3_mc_truth::flagMcRec) == 0;
295-
Partition<CandsLcRecWMl> recoLcCandSigWMl = nabs(o2::aod::a3_mc_truth::flagMcRec) == static_cast<int>(o2::constants::physics::Pdg::kLambdaCPlus);
295+
Partition<CandsLcRecWMl> recoLcCandSigWMl = nabs(o2::aod::a3_mc_truth::flagMcRec) == static_cast<int>(CharmHadAlice3::Lc);
296296
Partition<CandsLcRecWMl> recoLcCandBkgWMl = nabs(o2::aod::a3_mc_truth::flagMcRec) == 0;
297297

298298
void init(InitContext const&)

‎ALICE3/Tasks/alice3-task3prong.cxx‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -80,7 +80,7 @@ struct Alice3Task3Prong {
8080

8181
Filter filterSelectCandidates = aod::a3_hf_sel_3prong::isSel == true;
8282

83-
Partition<Cands3PGen> candsGenLcs = nabs(aod::a3_mc_truth::flagMcGen) == static_cast<int>(o2::constants::physics::Pdg::kLambdaCPlus);
83+
Partition<Cands3PGen> candsGenLcs = nabs(aod::a3_mc_truth::flagMcGen) == static_cast<int>(CharmHadAlice3::Lc);
8484

8585
ConfigurableAxis thnConfigAxisPt{"thnConfigAxisPt", {72, 0, 36}, ""};
8686
ConfigurableAxis thnConfigAxisMass{"thnConfigAxisMass", {300, 1.98, 2.58}, ""};
@@ -110,7 +110,7 @@ struct Alice3Task3Prong {
110110
}
111111

112112
if (doprocessLc || doprocessLcWMl) {
113-
selectedPdg = o2::constants::physics::Pdg::kLambdaCPlus;
113+
selectedPdg = CharmHadAlice3::Lc;
114114
}
115115

116116
auto addHistogramsRec = [&](const std::string& histoName, const std::string& xAxisTitle, const std::string& yAxisTitle, const HistogramConfigSpec& configSpec) {
@@ -263,7 +263,7 @@ struct Alice3Task3Prong {
263263

264264
const auto pt = candidate.pt();
265265
const auto originType = candidate.originMcRec();
266-
// const auto ptRecB = candidate.ptBhadMotherPart();
266+
const auto ptRecB = candidate.bHadMotherPtRec();
267267

268268
/// Fill histograms
269269
fillHistogramsRecSig<CharmHad, Signal>(candidate);
@@ -332,7 +332,7 @@ struct Alice3Task3Prong {
332332
fillHistogramsGen<CharmHad, Prompt>(particle);
333333
} else if (particle.originMcGen() == RecoDecay::OriginType::NonPrompt) {
334334
LOG(info) << "B-hadron pt still to be implemented!";
335-
// float ptGenB = isPrompt ? -1. : mcParticles.rawIteratorAt(particle.idxBhadMotherPart()).pt();
335+
ptGenB = particle.bHadMotherPtGen();
336336
fillHistogramsGen<CharmHad, NonPrompt>(particle);
337337
}
338338

‎PWGHF/Core/SelectorCuts.h‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1574,7 +1574,7 @@ constexpr double BinsPt[NBinsPt + 1] = {
15741574
36.};
15751575
const auto vecBinsPt = std::vector<double>{BinsPt, BinsPt + NBinsPt + 1};
15761576

1577-
// default values for the cuts m, ptP, ptK, ptPi, chi2PCA, dL, cosp, dLXY, NdLXY, ImpParXY
1577+
// default values for the cuts m, ptP, ptK, ptPi, chi2PCA, cosp, dL, dLXY, NdLXY, ImpParXY
15781578
constexpr double Cuts[NBinsPt][NCutVars] = {{0.4, 0.4, 0.4, 0.4, 0., 0.005, 0., 0., 0., 1e+10}, /* 0 < pT < 1 */
15791579
{0.4, 0.4, 0.4, 0.4, 0., 0.005, 0., 0., 0., 1e+10}, /* 1 < pT < 2 */
15801580
{0.4, 0.4, 0.4, 0.4, 0., 0.005, 0., 0., 0., 1e+10}, /* 2 < pT < 3 */
@@ -1600,7 +1600,7 @@ static const std::vector<std::string> labelsPt = {
16001600
"pT bin 9"};
16011601

16021602
// column labels
1603-
static const std::vector<std::string> labelsCutVar = {"m", "pT prong 0", "pT prong 1", "pT prong 2", "Chi2PCA", "decay length", "cos pointing angle", "decLengthXY", "normDecLXY", "impParXY"};
1603+
static const std::vector<std::string> labelsCutVar = {"m", "pT prong 0", "pT prong 1", "pT prong 2", "Chi2PCA", "cos pointing angle", "decay length", "decLengthXY", "normDecLXY", "impParXY"};
16041604
} // namespace hf_cuts_3prongs_alice3
16051605

16061606
} // namespace o2::analysis

0 commit comments

Comments
 (0)