Skip to content

Commit 991b2b6

Browse files
vkuceravkucera
committed
Fix member constants
1 parent 16347b2 commit 991b2b6

File tree

9 files changed

+52
-93
lines changed

9 files changed

+52
-93
lines changed

‎PWGHF/HFC/TableProducer/correlatorD0Hadrons.cxx‎

Lines changed: 14 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -239,9 +239,6 @@ struct HfCorrelatorD0Hadrons {
239239
Configurable<bool> useCentrality{"useCentrality", false, "Flag for centrality dependent analyses"};
240240

241241
int leadingIndex = 0;
242-
double massD0{0.};
243-
double massPi{0.};
244-
double massK{0.};
245242
double softPiMass = 0.14543; // pion mass + Q-value of the D*->D0pi decay
246243

247244
SliceCache cache;
@@ -274,10 +271,6 @@ struct HfCorrelatorD0Hadrons {
274271

275272
void init(InitContext&)
276273
{
277-
massD0 = MassD0;
278-
massPi = MassPiPlus;
279-
massK = MassKPlus;
280-
281274
AxisSpec axisMassD = {binsMassD, "inv. mass (#pi K) (GeV/#it{c}^{2})"};
282275
AxisSpec const axisEta = {binsEta, "#it{#eta}"};
283276
AxisSpec const axisPhi = {binsPhi, "#it{#varphi}"};
@@ -415,8 +408,8 @@ struct HfCorrelatorD0Hadrons {
415408
}
416409

417410
// ========================== Define parameters for soft pion removal ================================
418-
auto ePiK = RecoDecay::e(candidate.pVectorProng0(), massPi) + RecoDecay::e(candidate.pVectorProng1(), massK);
419-
auto eKPi = RecoDecay::e(candidate.pVectorProng0(), massK) + RecoDecay::e(candidate.pVectorProng1(), massPi);
411+
auto ePiK = RecoDecay::e(candidate.pVectorProng0(), MassPiPlus) + RecoDecay::e(candidate.pVectorProng1(), MassKPlus);
412+
auto eKPi = RecoDecay::e(candidate.pVectorProng0(), MassKPlus) + RecoDecay::e(candidate.pVectorProng1(), MassPiPlus);
420413

421414
// ========================== trigger efficiency ================================
422415
double efficiencyWeight = 1.;
@@ -479,7 +472,7 @@ struct HfCorrelatorD0Hadrons {
479472
// ========== soft pion removal ===================================================
480473
double invMassDstar1 = 0., invMassDstar2 = 0.;
481474
auto pSum2 = RecoDecay::p2(candidate.pVector(), track.pVector());
482-
auto ePion = track.energy(massPi);
475+
auto ePion = track.energy(MassPiPlus);
483476
invMassDstar1 = std::sqrt((ePiK + ePion) * (ePiK + ePion) - pSum2);
484477
invMassDstar2 = std::sqrt((eKPi + ePion) * (eKPi + ePion) - pSum2);
485478

@@ -655,8 +648,8 @@ struct HfCorrelatorD0Hadrons {
655648
entryD0CandGenInfo(isD0Prompt);
656649

657650
// ===================== Define parameters for soft pion removal ========================
658-
auto ePiK = RecoDecay::e(candidate.pVectorProng0(), massPi) + RecoDecay::e(candidate.pVectorProng1(), massK);
659-
auto eKPi = RecoDecay::e(candidate.pVectorProng0(), massK) + RecoDecay::e(candidate.pVectorProng1(), massPi);
651+
auto ePiK = RecoDecay::e(candidate.pVectorProng0(), MassPiPlus) + RecoDecay::e(candidate.pVectorProng1(), MassKPlus);
652+
auto eKPi = RecoDecay::e(candidate.pVectorProng0(), MassKPlus) + RecoDecay::e(candidate.pVectorProng1(), MassPiPlus);
660653

661654
// ============== D-h correlation dedicated section ====================================
662655

@@ -684,7 +677,7 @@ struct HfCorrelatorD0Hadrons {
684677
double invMassDstar1 = 0, invMassDstar2 = 0;
685678
bool isSoftPiD0 = false, isSoftPiD0bar = false;
686679
auto pSum2 = RecoDecay::p2(candidate.pVector(), track.pVector());
687-
auto ePion = track.energy(massPi);
680+
auto ePion = track.energy(MassPiPlus);
688681
invMassDstar1 = std::sqrt((ePiK + ePion) * (ePiK + ePion) - pSum2);
689682
invMassDstar2 = std::sqrt((eKPi + ePion) * (eKPi + ePion) - pSum2);
690683

@@ -856,7 +849,7 @@ struct HfCorrelatorD0Hadrons {
856849
poolBin,
857850
correlationStatus,
858851
cent);
859-
entryD0HadronRecoInfo(massD0, massD0, 0); // dummy info
852+
entryD0HadronRecoInfo(MassD0, MassD0, 0); // dummy info
860853
entryD0HadronGenInfo(isD0Prompt, particleAssoc.isPhysicalPrimary(), trackOrigin);
861854
} // end inner loop (Tracks)
862855
}
@@ -897,12 +890,12 @@ struct HfCorrelatorD0Hadrons {
897890
// soft pion removal, signal status 1,3 for D0 and 2,3 for D0bar (SoftPi removed), signal status 11,13 for D0 and 12,13 for D0bar (only SoftPi)
898891
const auto invMassD0 = HfHelper::invMassD0ToPiK(candidate);
899892
const auto invMassD0bar = HfHelper::invMassD0barToKPi(candidate);
900-
auto ePiK = RecoDecay::e(candidate.pVectorProng0(), massPi) + RecoDecay::e(candidate.pVectorProng1(), massK);
901-
auto eKPi = RecoDecay::e(candidate.pVectorProng0(), massK) + RecoDecay::e(candidate.pVectorProng1(), massPi);
893+
auto ePiK = RecoDecay::e(candidate.pVectorProng0(), MassPiPlus) + RecoDecay::e(candidate.pVectorProng1(), MassKPlus);
894+
auto eKPi = RecoDecay::e(candidate.pVectorProng0(), MassKPlus) + RecoDecay::e(candidate.pVectorProng1(), MassPiPlus);
902895
double invMassDstar1 = 0., invMassDstar2 = 0.;
903896
bool isSoftPiD0 = false, isSoftPiD0bar = false;
904897
auto pSum2 = RecoDecay::p2(candidate.pVector(), particleAssoc.pVector());
905-
auto ePion = particleAssoc.energy(massPi);
898+
auto ePion = particleAssoc.energy(MassPiPlus);
906899
invMassDstar1 = std::sqrt((ePiK + ePion) * (ePiK + ePion) - pSum2);
907900
invMassDstar2 = std::sqrt((eKPi + ePion) * (eKPi + ePion) - pSum2);
908901
std::vector<float> outputMlD0 = {-1., -1., -1.};
@@ -1004,12 +997,12 @@ struct HfCorrelatorD0Hadrons {
1004997
// soft pion removal
1005998
const auto invMassD0 = HfHelper::invMassD0ToPiK(candidate);
1006999
const auto invMassD0bar = HfHelper::invMassD0barToKPi(candidate);
1007-
auto ePiK = RecoDecay::e(candidate.pVectorProng0(), massPi) + RecoDecay::e(candidate.pVectorProng1(), massK);
1008-
auto eKPi = RecoDecay::e(candidate.pVectorProng0(), massK) + RecoDecay::e(candidate.pVectorProng1(), massPi);
1000+
auto ePiK = RecoDecay::e(candidate.pVectorProng0(), MassPiPlus) + RecoDecay::e(candidate.pVectorProng1(), MassKPlus);
1001+
auto eKPi = RecoDecay::e(candidate.pVectorProng0(), MassKPlus) + RecoDecay::e(candidate.pVectorProng1(), MassPiPlus);
10091002
double invMassDstar1 = 0., invMassDstar2 = 0.;
10101003
bool isSoftPiD0 = false, isSoftPiD0bar = false;
10111004
auto pSum2 = RecoDecay::p2(candidate.pVector(), particleAssoc.pVector());
1012-
auto ePion = particleAssoc.energy(massPi);
1005+
auto ePion = particleAssoc.energy(MassPiPlus);
10131006
invMassDstar1 = std::sqrt((ePiK + ePion) * (ePiK + ePion) - pSum2);
10141007
invMassDstar2 = std::sqrt((eKPi + ePion) * (eKPi + ePion) - pSum2);
10151008

@@ -1128,7 +1121,7 @@ struct HfCorrelatorD0Hadrons {
11281121
int trackOrigin = RecoDecay::getCharmHadronOrigin(mcParticles, particleAssoc, true);
11291122
bool isD0Prompt = particleTrigg.originMcGen() == RecoDecay::OriginType::Prompt;
11301123
entryD0HadronPair(getDeltaPhi(particleAssoc.phi(), particleTrigg.phi()), particleAssoc.eta() - particleTrigg.eta(), particleTrigg.pt(), particleAssoc.pt(), poolBin, correlationStatus, cent);
1131-
entryD0HadronRecoInfo(massD0, massD0, 0); // dummy info
1124+
entryD0HadronRecoInfo(MassD0, MassD0, 0); // dummy info
11321125
entryD0HadronGenInfo(isD0Prompt, particleAssoc.isPhysicalPrimary(), trackOrigin);
11331126
}
11341127
}

‎PWGHF/HFC/TableProducer/correlatorDMesonPairs.cxx‎

Lines changed: 1 addition & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -984,9 +984,6 @@ struct HfCorrelatorDMesonPairs {
984984
registry.fill(HIST("hInputCheckD0OrD0barMcGen"), nDorDbarevent);
985985
}
986986

987-
auto massD = MassD0;
988-
auto massDbar = MassD0Bar;
989-
990987
for (const auto& particle1 : mcParticles) {
991988
// check if the particle is D0 or D0bar
992989
if (std::abs(particle1.pdgCode()) != Pdg::kD0) {
@@ -1112,7 +1109,7 @@ struct HfCorrelatorDMesonPairs {
11121109
}
11131110

11141111
// Fill pair Selection Status
1115-
entryD0PairMcGen(particle1.pt(), particle2.pt(), particle1.y(), particle2.y(), particle1.phi(), particle2.phi(), massD, massDbar, massD, massDbar, pairType, particleType1, particleType2);
1112+
entryD0PairMcGen(particle1.pt(), particle2.pt(), particle1.y(), particle2.y(), particle1.phi(), particle2.phi(), MassD0, MassD0Bar, MassD0, MassD0Bar, pairType, particleType1, particleType2);
11161113
entryD0PairMcGenInfo(originGen1, originGen2, matchedGen1, matchedGen2);
11171114

11181115
} // end inner loop

‎PWGHF/TableProducer/candidateCreator2Prong.cxx‎

Lines changed: 12 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -116,12 +116,12 @@ struct HfCandidateCreator2Prong {
116116

117117
HfEventSelection hfEvSel; // event selection and monitoring
118118
o2::vertexing::DCAFitterN<2> df; // 2-prong vertex fitter
119-
Service<o2::ccdb::BasicCCDBManager> ccdb;
119+
Service<o2::ccdb::BasicCCDBManager> ccdb{};
120120

121121
int runNumber{0};
122122
double bz{0.};
123123

124-
const float toMicrometers = 10000.; // from cm to µm
124+
static constexpr float CentiToMicro = 10000.; // from cm to µm
125125

126126
std::shared_ptr<TH1> hCandidates;
127127

@@ -313,13 +313,13 @@ struct HfCandidateCreator2Prong {
313313
o2::dataformats::DCA impactParameter1;
314314
trackParVar0.propagateToDCA(primaryVertex, bz, &impactParameter0);
315315
trackParVar1.propagateToDCA(primaryVertex, bz, &impactParameter1);
316-
registry.fill(HIST("hDcaXYProngs"), track0.pt(), impactParameter0.getY() * toMicrometers);
317-
registry.fill(HIST("hDcaXYProngs"), track1.pt(), impactParameter1.getY() * toMicrometers);
318-
registry.fill(HIST("hDcaZProngs"), track0.pt(), impactParameter0.getZ() * toMicrometers);
319-
registry.fill(HIST("hDcaZProngs"), track1.pt(), impactParameter1.getZ() * toMicrometers);
316+
registry.fill(HIST("hDcaXYProngs"), track0.pt(), impactParameter0.getY() * CentiToMicro);
317+
registry.fill(HIST("hDcaXYProngs"), track1.pt(), impactParameter1.getY() * CentiToMicro);
318+
registry.fill(HIST("hDcaZProngs"), track0.pt(), impactParameter0.getZ() * CentiToMicro);
319+
registry.fill(HIST("hDcaZProngs"), track1.pt(), impactParameter1.getZ() * CentiToMicro);
320320

321321
// get uncertainty of the decay length
322-
double phi, theta;
322+
double phi{}, theta{};
323323
getPointDirection(std::array{primaryVertex.getX(), primaryVertex.getY(), primaryVertex.getZ()}, secondaryVertex, phi, theta);
324324
auto errorDecayLength = std::sqrt(getRotatedCovMatrixXX(covMatrixPV, phi, theta) + getRotatedCovMatrixXX(covMatrixPCA, phi, theta));
325325
auto errorDecayLengthXY = std::sqrt(getRotatedCovMatrixXX(covMatrixPV, phi, 0.) + getRotatedCovMatrixXX(covMatrixPCA, phi, 0.));
@@ -439,15 +439,15 @@ struct HfCandidateCreator2Prong {
439439

440440
float impactParameter0XY = 0., errImpactParameter0XY = 0., impactParameter1XY = 0., errImpactParameter1XY = 0.;
441441
if (!kfPosPion.GetDistanceFromVertexXY(kfpV, impactParameter0XY, errImpactParameter0XY)) {
442-
registry.fill(HIST("hDcaXYProngs"), track0.pt(), impactParameter0XY * toMicrometers);
443-
registry.fill(HIST("hDcaZProngs"), track0.pt(), std::sqrt(kfPosPion.GetDistanceFromVertex(kfpV) * kfPosPion.GetDistanceFromVertex(kfpV) - impactParameter0XY * impactParameter0XY) * toMicrometers);
442+
registry.fill(HIST("hDcaXYProngs"), track0.pt(), impactParameter0XY * CentiToMicro);
443+
registry.fill(HIST("hDcaZProngs"), track0.pt(), std::sqrt(kfPosPion.GetDistanceFromVertex(kfpV) * kfPosPion.GetDistanceFromVertex(kfpV) - impactParameter0XY * impactParameter0XY) * CentiToMicro);
444444
} else {
445445
registry.fill(HIST("hDcaXYProngs"), track0.pt(), -999.f);
446446
registry.fill(HIST("hDcaZProngs"), track0.pt(), -999.f);
447447
}
448448
if (!kfNegPion.GetDistanceFromVertexXY(kfpV, impactParameter1XY, errImpactParameter1XY)) {
449-
registry.fill(HIST("hDcaXYProngs"), track1.pt(), impactParameter1XY * toMicrometers);
450-
registry.fill(HIST("hDcaZProngs"), track1.pt(), std::sqrt(kfNegPion.GetDistanceFromVertex(kfpV) * kfNegPion.GetDistanceFromVertex(kfpV) - impactParameter1XY * impactParameter1XY) * toMicrometers);
449+
registry.fill(HIST("hDcaXYProngs"), track1.pt(), impactParameter1XY * CentiToMicro);
450+
registry.fill(HIST("hDcaZProngs"), track1.pt(), std::sqrt(kfNegPion.GetDistanceFromVertex(kfpV) * kfNegPion.GetDistanceFromVertex(kfpV) - impactParameter1XY * impactParameter1XY) * CentiToMicro);
451451
} else {
452452
registry.fill(HIST("hDcaXYProngs"), track1.pt(), -999.f);
453453
registry.fill(HIST("hDcaZProngs"), track1.pt(), -999.f);
@@ -472,7 +472,7 @@ struct HfCandidateCreator2Prong {
472472

473473
auto* covMatrixSV = kfCandD0.CovarianceMatrix();
474474

475-
double phi, theta;
475+
double phi{}, theta{};
476476
getPointDirection(std::array{kfpV.GetX(), kfpV.GetY(), kfpV.GetZ()}, std::array{kfCandD0.GetX(), kfCandD0.GetY(), kfCandD0.GetZ()}, phi, theta);
477477
auto errorDecayLength = std::sqrt(getRotatedCovMatrixXX(covMatrixPV, phi, theta) + getRotatedCovMatrixXX(covMatrixSV, phi, theta));
478478
auto errorDecayLengthXY = std::sqrt(getRotatedCovMatrixXX(covMatrixPV, phi, 0.) + getRotatedCovMatrixXX(covMatrixSV, phi, 0.));

‎PWGHF/TableProducer/candidateCreator3Prong.cxx‎

Lines changed: 17 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -133,13 +133,13 @@ struct HfCandidateCreator3Prong {
133133

134134
HfEventSelection hfEvSel; // event selection and monitoring
135135
o2::vertexing::DCAFitterN<3> df; // 3-prong vertex fitter
136-
Service<o2::ccdb::BasicCCDBManager> ccdb;
136+
Service<o2::ccdb::BasicCCDBManager> ccdb{};
137137

138138
int runNumber{0};
139139
double bz{0.};
140140

141-
const float toMicrometers = 10000.; // from cm to µm
142-
constexpr static float UndefValueFloat{-999.f};
141+
static constexpr float CentiToMicro = 10000.; // from cm to µm
142+
static constexpr float UndefValueFloat{-999.f};
143143

144144
using FilteredHf3Prongs = soa::Filtered<aod::Hf3Prongs>;
145145
using FilteredPvRefitHf3Prongs = soa::Filtered<soa::Join<aod::Hf3Prongs, aod::HfPvRefit3Prong>>;
@@ -371,15 +371,15 @@ struct HfCandidateCreator3Prong {
371371
trackParVar0.propagateToDCA(primaryVertex, bz, &impactParameter0);
372372
trackParVar1.propagateToDCA(primaryVertex, bz, &impactParameter1);
373373
trackParVar2.propagateToDCA(primaryVertex, bz, &impactParameter2);
374-
registry.fill(HIST("hDcaXYProngs"), track0.pt(), impactParameter0.getY() * toMicrometers);
375-
registry.fill(HIST("hDcaXYProngs"), track1.pt(), impactParameter1.getY() * toMicrometers);
376-
registry.fill(HIST("hDcaXYProngs"), track2.pt(), impactParameter2.getY() * toMicrometers);
377-
registry.fill(HIST("hDcaZProngs"), track0.pt(), impactParameter0.getZ() * toMicrometers);
378-
registry.fill(HIST("hDcaZProngs"), track1.pt(), impactParameter1.getZ() * toMicrometers);
379-
registry.fill(HIST("hDcaZProngs"), track2.pt(), impactParameter2.getZ() * toMicrometers);
374+
registry.fill(HIST("hDcaXYProngs"), track0.pt(), impactParameter0.getY() * CentiToMicro);
375+
registry.fill(HIST("hDcaXYProngs"), track1.pt(), impactParameter1.getY() * CentiToMicro);
376+
registry.fill(HIST("hDcaXYProngs"), track2.pt(), impactParameter2.getY() * CentiToMicro);
377+
registry.fill(HIST("hDcaZProngs"), track0.pt(), impactParameter0.getZ() * CentiToMicro);
378+
registry.fill(HIST("hDcaZProngs"), track1.pt(), impactParameter1.getZ() * CentiToMicro);
379+
registry.fill(HIST("hDcaZProngs"), track2.pt(), impactParameter2.getZ() * CentiToMicro);
380380

381381
// get uncertainty of the decay length
382-
double phi, theta;
382+
double phi{}, theta{};
383383
getPointDirection(std::array{primaryVertex.getX(), primaryVertex.getY(), primaryVertex.getZ()}, secondaryVertex, phi, theta);
384384
auto errorDecayLength = std::sqrt(getRotatedCovMatrixXX(covMatrixPV, phi, theta) + getRotatedCovMatrixXX(covMatrixPCA, phi, theta));
385385
auto errorDecayLengthXY = std::sqrt(getRotatedCovMatrixXX(covMatrixPV, phi, 0.) + getRotatedCovMatrixXX(covMatrixPCA, phi, 0.));
@@ -513,22 +513,22 @@ struct HfCandidateCreator3Prong {
513513

514514
float impactParameter0XY = 0., errImpactParameter0XY = 0., impactParameter1XY = 0., errImpactParameter1XY = 0., impactParameter2XY = 0., errImpactParameter2XY = 0.;
515515
if (!kfFirstProton.GetDistanceFromVertexXY(kfpV, impactParameter0XY, errImpactParameter0XY)) {
516-
registry.fill(HIST("hDcaXYProngs"), track0.pt(), impactParameter0XY * toMicrometers);
517-
registry.fill(HIST("hDcaZProngs"), track0.pt(), std::sqrt(kfFirstProton.GetDistanceFromVertex(kfpV) * kfFirstProton.GetDistanceFromVertex(kfpV) - impactParameter0XY * impactParameter0XY) * toMicrometers);
516+
registry.fill(HIST("hDcaXYProngs"), track0.pt(), impactParameter0XY * CentiToMicro);
517+
registry.fill(HIST("hDcaZProngs"), track0.pt(), std::sqrt(kfFirstProton.GetDistanceFromVertex(kfpV) * kfFirstProton.GetDistanceFromVertex(kfpV) - impactParameter0XY * impactParameter0XY) * CentiToMicro);
518518
} else {
519519
registry.fill(HIST("hDcaXYProngs"), track0.pt(), UndefValueFloat);
520520
registry.fill(HIST("hDcaZProngs"), track0.pt(), UndefValueFloat);
521521
}
522522
if (!kfSecondKaon.GetDistanceFromVertexXY(kfpV, impactParameter1XY, errImpactParameter1XY)) {
523-
registry.fill(HIST("hDcaXYProngs"), track1.pt(), impactParameter1XY * toMicrometers);
524-
registry.fill(HIST("hDcaZProngs"), track1.pt(), std::sqrt(kfSecondKaon.GetDistanceFromVertex(kfpV) * kfSecondKaon.GetDistanceFromVertex(kfpV) - impactParameter1XY * impactParameter1XY) * toMicrometers);
523+
registry.fill(HIST("hDcaXYProngs"), track1.pt(), impactParameter1XY * CentiToMicro);
524+
registry.fill(HIST("hDcaZProngs"), track1.pt(), std::sqrt(kfSecondKaon.GetDistanceFromVertex(kfpV) * kfSecondKaon.GetDistanceFromVertex(kfpV) - impactParameter1XY * impactParameter1XY) * CentiToMicro);
525525
} else {
526526
registry.fill(HIST("hDcaXYProngs"), track1.pt(), UndefValueFloat);
527527
registry.fill(HIST("hDcaZProngs"), track1.pt(), UndefValueFloat);
528528
}
529529
if (!kfThirdProton.GetDistanceFromVertexXY(kfpV, impactParameter2XY, errImpactParameter2XY)) {
530-
registry.fill(HIST("hDcaXYProngs"), track2.pt(), impactParameter2XY * toMicrometers);
531-
registry.fill(HIST("hDcaZProngs"), track2.pt(), std::sqrt(kfThirdProton.GetDistanceFromVertex(kfpV) * kfThirdProton.GetDistanceFromVertex(kfpV) - impactParameter2XY * impactParameter2XY) * toMicrometers);
530+
registry.fill(HIST("hDcaXYProngs"), track2.pt(), impactParameter2XY * CentiToMicro);
531+
registry.fill(HIST("hDcaZProngs"), track2.pt(), std::sqrt(kfThirdProton.GetDistanceFromVertex(kfpV) * kfThirdProton.GetDistanceFromVertex(kfpV) - impactParameter2XY * impactParameter2XY) * CentiToMicro);
532532
} else {
533533
registry.fill(HIST("hDcaXYProngs"), track2.pt(), UndefValueFloat);
534534
registry.fill(HIST("hDcaZProngs"), track2.pt(), UndefValueFloat);
@@ -624,7 +624,7 @@ struct HfCandidateCreator3Prong {
624624
registry.fill(HIST("hCovSVZZ"), kfCandPKPi.Covariance(2, 2));
625625

626626
auto* covMatrixSV = kfCandPKPi.CovarianceMatrix();
627-
double phi, theta;
627+
double phi{}, theta{};
628628
getPointDirection(std::array{kfpV.GetX(), kfpV.GetY(), kfpV.GetZ()}, std::array{kfCandPKPi.GetX(), kfCandPKPi.GetY(), kfCandPKPi.GetZ()}, phi, theta);
629629
auto errorDecayLength = std::sqrt(getRotatedCovMatrixXX(covMatrixPV, phi, theta) + getRotatedCovMatrixXX(covMatrixSV, phi, theta));
630630
auto errorDecayLengthXY = std::sqrt(getRotatedCovMatrixXX(covMatrixPV, phi, 0.) + getRotatedCovMatrixXX(covMatrixSV, phi, 0.));

‎PWGHF/TableProducer/candidateCreatorB0.cxx‎

Lines changed: 4 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -97,14 +97,11 @@ struct HfCandidateCreatorB0 {
9797
Configurable<std::string> ccdbPathGrp{"ccdbPathGrp", "GLO/GRP/GRP", "Path of the grp file (Run 2)"};
9898
Configurable<std::string> ccdbPathGrpMag{"ccdbPathGrpMag", "GLO/Config/GRPMagField", "CCDB path of the GRPMagField object (Run 3)"};
9999

100-
Service<o2::ccdb::BasicCCDBManager> ccdb;
100+
Service<o2::ccdb::BasicCCDBManager> ccdb{};
101101
o2::base::MatLayerCylSet* lut{};
102102
o2::base::Propagator::MatCorrType matCorr = o2::base::Propagator::MatCorrType::USEMatCorrLUT;
103103
int runNumber{};
104104

105-
double massPi{0.};
106-
double massD{0.};
107-
double massB0{0.};
108105
double invMass2DPiMin{0.};
109106
double invMass2DPiMax{0.};
110107
double bz{0.};
@@ -136,11 +133,8 @@ struct HfCandidateCreatorB0 {
136133
void init(InitContext const&)
137134
{
138135
// invariant-mass window cut
139-
massPi = MassPiPlus;
140-
massD = MassDMinus;
141-
massB0 = MassB0;
142-
invMass2DPiMin = (massB0 - invMassWindowB0) * (massB0 - invMassWindowB0);
143-
invMass2DPiMax = (massB0 + invMassWindowB0) * (massB0 + invMassWindowB0);
136+
invMass2DPiMin = (MassB0 - invMassWindowB0) * (MassB0 - invMassWindowB0);
137+
invMass2DPiMax = (MassB0 + invMassWindowB0) * (MassB0 + invMassWindowB0);
144138

145139
// Initialise fitter for B vertex (2-prong vertex filter)
146140
df2.setPropagateToPCA(propagateToPCA);
@@ -352,7 +346,7 @@ struct HfCandidateCreatorB0 {
352346
df2.getTrack(1).getPxPyPzGlo(pVecPion); // momentum of Pi at the B0 vertex
353347

354348
// calculate invariant mass square and apply selection
355-
auto invMass2DPi = RecoDecay::m2(std::array{pVecD, pVecPion}, std::array{massD, massPi});
349+
auto invMass2DPi = RecoDecay::m2(std::array{pVecD, pVecPion}, std::array{MassDMinus, MassPiPlus});
356350
if ((invMass2DPi < invMass2DPiMin) || (invMass2DPi > invMass2DPiMax)) {
357351
continue;
358352
}

0 commit comments

Comments
 (0)