Skip to content

Commit e37977e

Browse files
vkuceravkucera
committed
Use namespace where useful
1 parent 44d55eb commit e37977e

24 files changed

+255
-231
lines changed

‎PWGHF/ALICE3/TableProducer/candidateCreatorChic.cxx‎

Lines changed: 10 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -31,6 +31,7 @@
3131
using namespace o2;
3232
using namespace o2::analysis;
3333
using namespace o2::aod;
34+
using namespace o2::constants::physics;
3435
using namespace o2::framework;
3536
using namespace o2::framework::expressions;
3637

@@ -80,7 +81,7 @@ struct HfCandidateCreatorChic {
8081

8182
void init(InitContext const&)
8283
{
83-
massJpsi = o2::constants::physics::MassJPsi;
84+
massJpsi = MassJPsi;
8485
}
8586

8687
void process(aod::Collision const& collision,
@@ -244,19 +245,19 @@ struct HfCandidateCreatorChicMc {
244245
auto arrayJpsiDaughters = std::array{daughterPosJpsi, daughterNegJpsi};
245246

246247
// chi_c → J/ψ gamma
247-
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayJpsiDaughters, o2::constants::physics::Pdg::kJPsi, std::array{+kMuonPlus, -kMuonPlus}, true);
248+
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayJpsiDaughters, Pdg::kJPsi, std::array{+kMuonPlus, -kMuonPlus}, true);
248249
if (indexRec > -1) {
249250
hMassJpsiToMuMuMatched->Fill(hfHelper.invMassJpsiToMuMu(candidate.prong0()));
250251

251-
int indexMother = RecoDecay::getMother(mcParticles, mcParticles.rawIteratorAt(indexRec), o2::constants::physics::Pdg::kChiC1);
252-
int indexMotherGamma = RecoDecay::getMother(mcParticles, mcParticles.rawIteratorAt(candidate.prong1().mcparticleId()), o2::constants::physics::Pdg::kChiC1);
252+
int indexMother = RecoDecay::getMother(mcParticles, mcParticles.rawIteratorAt(indexRec), Pdg::kChiC1);
253+
int indexMotherGamma = RecoDecay::getMother(mcParticles, mcParticles.rawIteratorAt(candidate.prong1().mcparticleId()), Pdg::kChiC1);
253254
if (indexMother > -1 && indexMotherGamma == indexMother && candidate.prong1().mcparticle().pdgCode() == kGamma) {
254255
auto particleMother = mcParticles.rawIteratorAt(indexMother);
255256
hEphotonMatched->Fill(candidate.prong1().e());
256257
hMassEMatched->Fill(sqrt(candidate.prong1().px() * candidate.prong1().px() + candidate.prong1().py() * candidate.prong1().py() + candidate.prong1().pz() * candidate.prong1().pz()));
257258
if (particleMother.has_daughters()) {
258259
std::vector<int> arrAllDaughtersIndex;
259-
RecoDecay::getDaughters(particleMother, &arrAllDaughtersIndex, std::array{static_cast<int>(kGamma), static_cast<int>(o2::constants::physics::Pdg::kJPsi)}, 1);
260+
RecoDecay::getDaughters(particleMother, &arrAllDaughtersIndex, std::array{static_cast<int>(kGamma), static_cast<int>(Pdg::kJPsi)}, 1);
260261
if (arrAllDaughtersIndex.size() == 2) {
261262
flag = 1 << hf_cand_chic::DecayType::ChicToJpsiToMuMuGamma;
262263
hMassChicToJpsiToMuMuGammaMatched->Fill(hfHelper.invMassChicToJpsiGamma(candidate));
@@ -278,17 +279,17 @@ struct HfCandidateCreatorChicMc {
278279
channel = 0;
279280

280281
// chi_c → J/ψ gamma
281-
if (RecoDecay::isMatchedMCGen(mcParticles, particle, o2::constants::physics::Pdg::kChiC1, std::array{static_cast<int>(o2::constants::physics::Pdg::kJPsi), static_cast<int>(kGamma)}, true)) {
282+
if (RecoDecay::isMatchedMCGen(mcParticles, particle, Pdg::kChiC1, std::array{static_cast<int>(Pdg::kJPsi), static_cast<int>(kGamma)}, true)) {
282283
// Match J/psi --> e+e-
283284
std::vector<int> arrDaughter;
284-
RecoDecay::getDaughters(particle, &arrDaughter, std::array{static_cast<int>(o2::constants::physics::Pdg::kJPsi)}, 1);
285+
RecoDecay::getDaughters(particle, &arrDaughter, std::array{static_cast<int>(Pdg::kJPsi)}, 1);
285286
auto jpsiCandMC = mcParticles.rawIteratorAt(arrDaughter[0]);
286-
if (RecoDecay::isMatchedMCGen(mcParticles, jpsiCandMC, o2::constants::physics::Pdg::kJPsi, std::array{+kElectron, -kElectron}, true)) {
287+
if (RecoDecay::isMatchedMCGen(mcParticles, jpsiCandMC, Pdg::kJPsi, std::array{+kElectron, -kElectron}, true)) {
287288
flag = 1 << hf_cand_chic::DecayType::ChicToJpsiToEEGamma;
288289
}
289290

290291
if (flag == 0) {
291-
if (RecoDecay::isMatchedMCGen(mcParticles, jpsiCandMC, o2::constants::physics::Pdg::kJPsi, std::array{+kMuonPlus, -kMuonPlus}, true)) {
292+
if (RecoDecay::isMatchedMCGen(mcParticles, jpsiCandMC, Pdg::kJPsi, std::array{+kMuonPlus, -kMuonPlus}, true)) {
292293
flag = 1 << hf_cand_chic::DecayType::ChicToJpsiToMuMuGamma;
293294
}
294295
}

‎PWGHF/ALICE3/TableProducer/candidateCreatorX.cxx‎

Lines changed: 7 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -31,6 +31,7 @@
3131
using namespace o2;
3232
using namespace o2::analysis;
3333
using namespace o2::aod;
34+
using namespace o2::constants::physics;
3435
using namespace o2::framework;
3536
using namespace o2::framework::expressions;
3637

@@ -80,8 +81,8 @@ struct HfCandidateCreatorX {
8081

8182
void init(InitContext const&)
8283
{
83-
massPi = o2::constants::physics::MassPiPlus;
84-
massJpsi = o2::constants::physics::MassJPsi;
84+
massPi = MassPiPlus;
85+
massJpsi = MassJPsi;
8586
}
8687

8788
void process(aod::Collision const& collision,
@@ -276,8 +277,8 @@ struct HfCandidateCreatorXMc {
276277
aod::McParticles const& mcParticles)
277278
{
278279
int indexRec = -1;
279-
int pdgCodeX = o2::constants::physics::Pdg::kX3872;
280-
int pdgCodeJpsi = o2::constants::physics::Pdg::kJPsi;
280+
int pdgCodeX = Pdg::kX3872;
281+
int pdgCodeJpsi = Pdg::kJPsi;
281282
int8_t sign = 0;
282283
int8_t flag = 0;
283284
int8_t origin = 0;
@@ -300,7 +301,7 @@ struct HfCandidateCreatorXMc {
300301
// X → J/ψ π+ π−
301302

302303
// J/ψ → e+ e−
303-
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayJpsiDaughters, o2::constants::physics::Pdg::kJPsi, std::array{+kElectron, -kElectron}, true);
304+
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayJpsiDaughters, Pdg::kJPsi, std::array{+kElectron, -kElectron}, true);
304305
// X → π+ π− e+ e−
305306
if (indexRec > -1) {
306307
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughters, pdgCodeX, std::array{+kPiPlus, -kPiPlus, +kElectron, -kElectron}, true, &sign, 2);
@@ -311,7 +312,7 @@ struct HfCandidateCreatorXMc {
311312

312313
// J/ψ → μ+ μ−
313314
if (flag == 0) {
314-
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayJpsiDaughters, o2::constants::physics::Pdg::kJPsi, std::array{+kMuonPlus, -kMuonPlus}, true);
315+
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayJpsiDaughters, Pdg::kJPsi, std::array{+kMuonPlus, -kMuonPlus}, true);
315316
// X → π+ π− μ+ μ−
316317
if (indexRec > -1) {
317318
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughters, pdgCodeX, std::array{+kPiPlus, -kPiPlus, +kMuonPlus, -kMuonPlus}, true, &sign, 2);

‎PWGHF/D2H/TableProducer/dataCreatorD0PiReduced.cxx‎

Lines changed: 9 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -37,6 +37,7 @@
3737
using namespace o2;
3838
using namespace o2::analysis;
3939
using namespace o2::aod;
40+
using namespace o2::constants::physics;
4041
using namespace o2::framework;
4142
using namespace o2::framework::expressions;
4243

@@ -162,9 +163,9 @@ struct HfDataCreatorD0PiReduced {
162163
runNumber = 0;
163164

164165
// invariant-mass window cut
165-
massPi = o2::constants::physics::MassPiPlus;
166-
massD0 = o2::constants::physics::MassD0;
167-
massBplus = o2::constants::physics::MassBPlus;
166+
massPi = MassPiPlus;
167+
massD0 = MassD0;
168+
massBplus = MassBPlus;
168169
invMass2D0PiMin = (massBplus - invMassWindowD0Pi) * (massBplus - invMassWindowD0Pi);
169170
invMass2D0PiMax = (massBplus + invMassWindowD0Pi) * (massBplus + invMassWindowD0Pi);
170171
}
@@ -348,18 +349,18 @@ struct HfDataCreatorD0PiReduced {
348349
int8_t flag{0};
349350
// B+ → D0(bar) π+ → (K+ π-) π+
350351
// Printf("Checking B+ → D0bar π+");
351-
auto indexRec = RecoDecay::getMatchedMCRec(particlesMc, arrayDaughtersBplus, o2::constants::physics::Pdg::kBPlus, std::array{+kPiPlus, +kKPlus, -kPiPlus}, true, &sign, 2);
352+
auto indexRec = RecoDecay::getMatchedMCRec(particlesMc, arrayDaughtersBplus, Pdg::kBPlus, std::array{+kPiPlus, +kKPlus, -kPiPlus}, true, &sign, 2);
352353
if (indexRec > -1) {
353354
// D0bar → K+ π-
354355
// Printf("Checking D0bar → K+ π-");
355-
indexRec = RecoDecay::getMatchedMCRec(particlesMc, arrayDaughtersD0, o2::constants::physics::Pdg::kD0, std::array{+kPiPlus, -kKPlus}, true, &sign, 1);
356+
indexRec = RecoDecay::getMatchedMCRec(particlesMc, arrayDaughtersD0, Pdg::kD0, std::array{+kPiPlus, -kKPlus}, true, &sign, 1);
356357
if (indexRec > -1) {
357358
flag = sign * BIT(hf_cand_bplus::DecayType::BplusToD0Pi);
358359
} else {
359360
LOGF(info, "WARNING: B+ decays in the expected final state but the condition on the intermediate state is not fulfilled");
360361
}
361362
}
362-
auto indexMother = RecoDecay::getMother(particlesMc, trackPion.template mcParticle_as<P>(), o2::constants::physics::Pdg::kBPlus, true);
363+
auto indexMother = RecoDecay::getMother(particlesMc, trackPion.template mcParticle_as<P>(), Pdg::kBPlus, true);
363364
auto particleMother = particlesMc.rawIteratorAt(indexMother);
364365

365366
rowHfD0PiMcRecReduced(indexHfCand2Prong, selectedTracksPion[trackPion.globalIndex()], flag, particleMother.pt());
@@ -405,11 +406,11 @@ struct HfDataCreatorD0PiReduced {
405406
int8_t sign{0};
406407
int8_t flag{0};
407408
// B+ → D0bar π+
408-
if (RecoDecay::isMatchedMCGen(particlesMc, particle, o2::constants::physics::Pdg::kBPlus, std::array{static_cast<int>(o2::constants::physics::Pdg::kD0), +kPiPlus}, true)) {
409+
if (RecoDecay::isMatchedMCGen(particlesMc, particle, Pdg::kBPlus, std::array{static_cast<int>(Pdg::kD0), +kPiPlus}, true)) {
409410
// Match D0bar -> π- K+
410411
auto candD0MC = particlesMc.rawIteratorAt(particle.daughtersIds().front());
411412
// Printf("Checking D0bar -> π- K+");
412-
if (RecoDecay::isMatchedMCGen(particlesMc, candD0MC, static_cast<int>(o2::constants::physics::Pdg::kD0), std::array{+kPiPlus, -kKPlus}, true, &sign)) {
413+
if (RecoDecay::isMatchedMCGen(particlesMc, candD0MC, static_cast<int>(Pdg::kD0), std::array{+kPiPlus, -kKPlus}, true, &sign)) {
413414
flag = sign * BIT(hf_cand_bplus::DecayType::BplusToD0Pi);
414415
}
415416
}

‎PWGHF/D2H/TableProducer/dataCreatorDplusPiReduced.cxx‎

Lines changed: 9 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -37,6 +37,7 @@
3737
using namespace o2;
3838
using namespace o2::analysis;
3939
using namespace o2::aod;
40+
using namespace o2::constants::physics;
4041
using namespace o2::framework;
4142
using namespace o2::framework::expressions;
4243

@@ -161,9 +162,9 @@ struct HfDataCreatorDplusPiReduced {
161162
runNumber = 0;
162163

163164
// invariant-mass window cut
164-
massPi = o2::constants::physics::MassPiPlus;
165-
massD = o2::constants::physics::MassDMinus;
166-
massB0 = o2::constants::physics::MassB0;
165+
massPi = MassPiPlus;
166+
massD = MassDMinus;
167+
massB0 = MassB0;
167168
invMass2DPiMin = (massB0 - invMassWindowDPi) * (massB0 - invMassWindowDPi);
168169
invMass2DPiMax = (massB0 + invMassWindowDPi) * (massB0 + invMassWindowDPi);
169170
}
@@ -351,19 +352,19 @@ struct HfDataCreatorDplusPiReduced {
351352
int8_t flag{0};
352353
int8_t debug{0};
353354
// B0 → D- π+ → (π- K+ π-) π+
354-
auto indexRec = RecoDecay::getMatchedMCRec(particlesMc, arrayDaughtersB0, o2::constants::physics::Pdg::kB0, std::array{-kPiPlus, +kKPlus, -kPiPlus, +kPiPlus}, true, &sign, 3);
355+
auto indexRec = RecoDecay::getMatchedMCRec(particlesMc, arrayDaughtersB0, Pdg::kB0, std::array{-kPiPlus, +kKPlus, -kPiPlus, +kPiPlus}, true, &sign, 3);
355356
if (indexRec > -1) {
356357
// D- → π- K+ π-
357358
// Printf("Checking D- → π- K+ π-");
358-
indexRec = RecoDecay::getMatchedMCRec(particlesMc, arrayDaughtersD, o2::constants::physics::Pdg::kDMinus, std::array{-kPiPlus, +kKPlus, -kPiPlus}, true, &sign, 2);
359+
indexRec = RecoDecay::getMatchedMCRec(particlesMc, arrayDaughtersD, Pdg::kDMinus, std::array{-kPiPlus, +kKPlus, -kPiPlus}, true, &sign, 2);
359360
if (indexRec > -1) {
360361
flag = sign * BIT(hf_cand_b0::DecayType::B0ToDPi);
361362
} else {
362363
debug = 1;
363364
LOGF(debug, "B0 decays in the expected final state but the condition on the intermediate state is not fulfilled");
364365
}
365366
}
366-
auto indexMother = RecoDecay::getMother(particlesMc, trackPion.template mcParticle_as<P>(), o2::constants::physics::Pdg::kB0, true);
367+
auto indexMother = RecoDecay::getMother(particlesMc, trackPion.template mcParticle_as<P>(), Pdg::kB0, true);
367368
auto particleMother = particlesMc.rawIteratorAt(indexMother);
368369
rowHfDPiMcRecReduced(indexHfCand3Prong, selectedTracksPion[trackPion.globalIndex()], flag, debug, particleMother.pt());
369370
}
@@ -409,11 +410,11 @@ struct HfDataCreatorDplusPiReduced {
409410
int8_t sign{0};
410411
int8_t flag{0};
411412
// B0 → D- π+
412-
if (RecoDecay::isMatchedMCGen(particlesMc, particle, o2::constants::physics::Pdg::kB0, std::array{-static_cast<int>(o2::constants::physics::Pdg::kDPlus), +kPiPlus}, true)) {
413+
if (RecoDecay::isMatchedMCGen(particlesMc, particle, Pdg::kB0, std::array{-static_cast<int>(Pdg::kDPlus), +kPiPlus}, true)) {
413414
// Match D- -> π- K+ π-
414415
auto candDMC = particlesMc.rawIteratorAt(particle.daughtersIds().front());
415416
// Printf("Checking D- -> π- K+ π-");
416-
if (RecoDecay::isMatchedMCGen(particlesMc, candDMC, -static_cast<int>(o2::constants::physics::Pdg::kDPlus), std::array{-kPiPlus, +kKPlus, -kPiPlus}, true, &sign)) {
417+
if (RecoDecay::isMatchedMCGen(particlesMc, candDMC, -static_cast<int>(Pdg::kDPlus), std::array{-kPiPlus, +kKPlus, -kPiPlus}, true, &sign)) {
417418
flag = sign * BIT(hf_cand_b0::DecayType::B0ToDPi);
418419
}
419420
}

‎PWGHF/HFC/TableProducer/correlatorD0D0bar.cxx‎

Lines changed: 8 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -29,6 +29,7 @@
2929

3030
using namespace o2;
3131
using namespace o2::analysis;
32+
using namespace o2::constants::physics;
3233
using namespace o2::framework;
3334
using namespace o2::framework::expressions;
3435

@@ -397,10 +398,10 @@ struct HfCorrelatorD0D0bar {
397398
// MC gen level
398399
for (const auto& particle1 : mcParticles) {
399400
// check if the particle is D0 or D0bar (for general plot filling and selection, so both cases are fine) - NOTE: decay channel is not probed!
400-
if (std::abs(particle1.pdgCode()) != o2::constants::physics::Pdg::kD0) {
401+
if (std::abs(particle1.pdgCode()) != Pdg::kD0) {
401402
continue;
402403
}
403-
double yD = RecoDecay::y(std::array{particle1.px(), particle1.py(), particle1.pz()}, o2::constants::physics::MassD0);
404+
double yD = RecoDecay::y(std::array{particle1.px(), particle1.py(), particle1.pz()}, MassD0);
404405
if (yCandMax >= 0. && std::abs(yD) > yCandMax) {
405406
continue;
406407
}
@@ -415,15 +416,15 @@ struct HfCorrelatorD0D0bar {
415416

416417
// D-Dbar correlation dedicated section
417418
// if it's a D0 particle, search for D0bar and evaluate correlations
418-
if (particle1.pdgCode() != o2::constants::physics::Pdg::kD0) { // just checking the particle PDG, not the decay channel (differently from Reco: you have a BR factor btw such levels!)
419+
if (particle1.pdgCode() != Pdg::kD0) { // just checking the particle PDG, not the decay channel (differently from Reco: you have a BR factor btw such levels!)
419420
continue;
420421
}
421422
registry.fill(HIST("hCountD0triggersMCGen"), 0, particle1.pt()); // to count trigger D0 (for normalisation)
422423
for (const auto& particle2 : mcParticles) {
423-
if (particle2.pdgCode() != o2::constants::physics::Pdg::kD0Bar) { // check that inner particle is D0bar
424+
if (particle2.pdgCode() != Pdg::kD0Bar) { // check that inner particle is D0bar
424425
continue;
425426
}
426-
if (yCandMax >= 0. && std::abs(RecoDecay::y(std::array{particle2.px(), particle2.py(), particle2.pz()}, o2::constants::physics::MassD0Bar)) > yCandMax) {
427+
if (yCandMax >= 0. && std::abs(RecoDecay::y(std::array{particle2.px(), particle2.py(), particle2.pz()}, MassD0Bar)) > yCandMax) {
427428
continue;
428429
}
429430
if (ptCandMin >= 0. && particle2.pt() < ptCandMin) {
@@ -498,7 +499,7 @@ struct HfCorrelatorD0D0bar {
498499
continue;
499500
}
500501
counterCCbarBeforeEtasel++; // count c or cbar (before kinematic selection)
501-
double yC = RecoDecay::y(std::array{particle1.px(), particle1.py(), particle1.pz()}, o2::constants::physics::MassCharm);
502+
double yC = RecoDecay::y(std::array{particle1.px(), particle1.py(), particle1.pz()}, MassCharm);
502503
if (yCandMax >= 0. && std::abs(yC) > yCandMax) {
503504
continue;
504505
}
@@ -522,7 +523,7 @@ struct HfCorrelatorD0D0bar {
522523
if (particle2.pdgCode() != PDG_t::kCharmBar) { // check that inner particle is a cbar
523524
continue;
524525
}
525-
if (yCandMax >= 0. && std::abs(RecoDecay::y(std::array{particle2.px(), particle2.py(), particle2.pz()}, o2::constants::physics::MassCharmBar)) > yCandMax) {
526+
if (yCandMax >= 0. && std::abs(RecoDecay::y(std::array{particle2.px(), particle2.py(), particle2.pz()}, MassCharmBar)) > yCandMax) {
526527
continue;
527528
}
528529
if (ptCandMin >= 0. && particle2.pt() < ptCandMin) {

‎PWGHF/HFC/TableProducer/correlatorD0D0barBarrelFullPid.cxx‎

Lines changed: 8 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -29,6 +29,7 @@
2929

3030
using namespace o2;
3131
using namespace o2::analysis;
32+
using namespace o2::constants::physics;
3233
using namespace o2::framework;
3334
using namespace o2::framework::expressions;
3435

@@ -398,10 +399,10 @@ struct HfCorrelatorD0D0barBarrelFullPid {
398399
// MC gen level
399400
for (const auto& particle1 : mcParticles) {
400401
// check if the particle is D0 or D0bar (for general plot filling and selection, so both cases are fine) - NOTE: decay channel is not probed!
401-
if (std::abs(particle1.pdgCode()) != o2::constants::physics::Pdg::kD0) {
402+
if (std::abs(particle1.pdgCode()) != Pdg::kD0) {
402403
continue;
403404
}
404-
double yD = RecoDecay::y(std::array{particle1.px(), particle1.py(), particle1.pz()}, o2::constants::physics::MassD0);
405+
double yD = RecoDecay::y(std::array{particle1.px(), particle1.py(), particle1.pz()}, MassD0);
405406
if (yCandMax >= 0. && std::abs(yD) > yCandMax) {
406407
continue;
407408
}
@@ -416,15 +417,15 @@ struct HfCorrelatorD0D0barBarrelFullPid {
416417

417418
// D-Dbar correlation dedicated section
418419
// if it's a D0 particle, search for D0bar and evaluate correlations
419-
if (particle1.pdgCode() != o2::constants::physics::Pdg::kD0) { // just checking the particle PDG, not the decay channel (differently from Reco: you have a BR factor btw such levels!)
420+
if (particle1.pdgCode() != Pdg::kD0) { // just checking the particle PDG, not the decay channel (differently from Reco: you have a BR factor btw such levels!)
420421
continue;
421422
}
422423
registry.fill(HIST("hCountD0triggersMCGen"), 0, particle1.pt()); // to count trigger D0 (for normalisation)
423424
for (const auto& particle2 : mcParticles) {
424-
if (particle2.pdgCode() != o2::constants::physics::Pdg::kD0Bar) { // check that inner particle is D0bar
425+
if (particle2.pdgCode() != Pdg::kD0Bar) { // check that inner particle is D0bar
425426
continue;
426427
}
427-
if (yCandMax >= 0. && std::abs(RecoDecay::y(std::array{particle2.px(), particle2.py(), particle2.pz()}, o2::constants::physics::MassD0Bar)) > yCandMax) {
428+
if (yCandMax >= 0. && std::abs(RecoDecay::y(std::array{particle2.px(), particle2.py(), particle2.pz()}, MassD0Bar)) > yCandMax) {
428429
continue;
429430
}
430431
if (ptCandMin >= 0. && particle2.pt() < ptCandMin) {
@@ -499,7 +500,7 @@ struct HfCorrelatorD0D0barBarrelFullPid {
499500
continue;
500501
}
501502
counterCCbarBeforeEtasel++; // count c or cbar (before kinematic selection)
502-
double yC = RecoDecay::y(std::array{particle1.px(), particle1.py(), particle1.pz()}, o2::constants::physics::MassCharm);
503+
double yC = RecoDecay::y(std::array{particle1.px(), particle1.py(), particle1.pz()}, MassCharm);
503504
if (yCandMax >= 0. && std::abs(yC) > yCandMax) {
504505
continue;
505506
}
@@ -523,7 +524,7 @@ struct HfCorrelatorD0D0barBarrelFullPid {
523524
if (particle2.pdgCode() != PDG_t::kCharmBar) { // check that inner particle is a cbar
524525
continue;
525526
}
526-
if (yCandMax >= 0. && std::abs(RecoDecay::y(std::array{particle2.px(), particle2.py(), particle2.pz()}, o2::constants::physics::MassCharmBar)) > yCandMax) {
527+
if (yCandMax >= 0. && std::abs(RecoDecay::y(std::array{particle2.px(), particle2.py(), particle2.pz()}, MassCharmBar)) > yCandMax) {
527528
continue;
528529
}
529530
if (ptCandMin >= 0. && particle2.pt() < ptCandMin) {

0 commit comments

Comments
 (0)