Skip to content
Merged
Show file tree
Hide file tree
Changes from 5 commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 2 additions & 4 deletions PWGHF/D2H/TableProducer/candidateCreatorB0Reduced.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -192,8 +192,7 @@
registry.fill(HIST("hCovSVXX"), covMatrixPCA[0]);
registry.fill(HIST("hCovPVXX"), covMatrixPV[0]);

// propagate D and Pi to the B0 vertex
df3.propagateTracksToVertex();
// get D and Pi tracks (propagated to the B0 vertex if propagateToPCA==true)
// track.getPxPyPzGlo(pVec) modifies pVec of track
df3.getTrack(0).getPxPyPzGlo(pVecD0); // momentum of D at the B0 vertex
df3.getTrack(1).getPxPyPzGlo(pVecSoftPion); // momentum of SoftPi at the B0 vertex
Expand Down Expand Up @@ -299,8 +298,7 @@
registry.fill(HIST("hCovSVXX"), covMatrixPCA[0]);
registry.fill(HIST("hCovPVXX"), covMatrixPV[0]);

// propagate D and Pi to the B0 vertex
df2.propagateTracksToVertex();
// get D and Pi tracks (propagated to the B0 vertex if propagateToPCA==true)
// track.getPxPyPzGlo(pVec) modifies pVec of track
df2.getTrack(0).getPxPyPzGlo(pVecD); // momentum of D at the B0 vertex
df2.getTrack(1).getPxPyPzGlo(pVecPion); // momentum of Pi at the B0 vertex
Expand Down Expand Up @@ -365,7 +363,7 @@
auto candsDThisColl = candsD.sliceBy(candsDplusPerCollision, thisCollId);
auto tracksPionThisCollision = tracksPion.sliceBy(tracksPionPerCollision, thisCollId);
runCandidateCreation<false>(collision, candsDThisColl, tracksPionThisCollision, invMass2DPiMin, invMass2DPiMax);
if (ncol % 10000 == 0) {

Check failure on line 366 in PWGHF/D2H/TableProducer/candidateCreatorB0Reduced.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
LOGP(debug, "collisions parsed {}", ncol);
}
ncol++;
Expand Down Expand Up @@ -399,7 +397,7 @@
auto candsDThisColl = candsD.sliceBy(candsDplusPerCollision, thisCollId);
auto tracksPionThisCollision = tracksPion.sliceBy(tracksPionPerCollision, thisCollId);
runCandidateCreation<true>(collision, candsDThisColl, tracksPionThisCollision, invMass2DPiMin, invMass2DPiMax);
if (ncol % 10000 == 0) {

Check failure on line 400 in PWGHF/D2H/TableProducer/candidateCreatorB0Reduced.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
LOGP(debug, "collisions parsed {}", ncol);
}
ncol++;
Expand Down Expand Up @@ -434,7 +432,7 @@
auto candsDThisColl = candsD.sliceBy(candsDstarPerCollision, thisCollId);
auto tracksPionThisCollision = tracksPion.sliceBy(tracksPionPerCollision, thisCollId);
runCandidateCreationDStar<false>(collision, candsDThisColl, softPions, tracksPionThisCollision, invMass2DPiMin, invMass2DPiMax);
if (ncol % 10000 == 0) {

Check failure on line 435 in PWGHF/D2H/TableProducer/candidateCreatorB0Reduced.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
LOGP(debug, "collisions parsed {}", ncol);
}
ncol++;
Expand Down Expand Up @@ -469,7 +467,7 @@
auto candsDThisColl = candsD.sliceBy(candsDstarWithMlPerCollision, thisCollId);
auto tracksPionThisCollision = tracksPion.sliceBy(tracksPionPerCollision, thisCollId);
runCandidateCreationDStar<true>(collision, candsDThisColl, softPions, tracksPionThisCollision, invMass2DPiMin, invMass2DPiMax);
if (ncol % 10000 == 0) {

Check failure on line 470 in PWGHF/D2H/TableProducer/candidateCreatorB0Reduced.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
LOGP(debug, "collisions parsed {}", ncol);
}
ncol++;
Expand Down
6 changes: 2 additions & 4 deletions PWGHF/D2H/TableProducer/candidateCreatorBToJpsiReduced.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -238,8 +238,7 @@ struct HfCandidateCreatorBToJpsiReduced {
registry.fill(HIST("hCovSVXX"), covMatrixPCA[0]);
registry.fill(HIST("hCovPVXX"), covMatrixPV[0]);

// propagate Jpsi daugthers and K to the B+ vertex
df3.propagateTracksToVertex();
// get JPsi daughters and K tracks (propagated to the B+ vertex if propagateToPCA==true)
// track.getPxPyPzGlo(pVec) modifies pVec of track
df3.getTrack(0).getPxPyPzGlo(pVecDauPos); // momentum of positive Jpsi daughter at the B+ vertex
df3.getTrack(1).getPxPyPzGlo(pVecDauNeg); // momentum of negative Jpsi daughter at the B+ vertex
Expand Down Expand Up @@ -309,8 +308,7 @@ struct HfCandidateCreatorBToJpsiReduced {
registry.fill(HIST("hCovSVXX"), covMatrixPCA[0]);
registry.fill(HIST("hCovPVXX"), covMatrixPV[0]);

// propagate Jpsi and phi to the Bs vertex
df4.propagateTracksToVertex();
// get JPsi daughters and K tracks (propagated to the Bs vertex if propagateToPCA==true)
// track.getPxPyPzGlo(pVec) modifies pVec of track
df4.getTrack(0).getPxPyPzGlo(pVecDauPos); // momentum of Jpsi at the B+ vertex
df4.getTrack(1).getPxPyPzGlo(pVecDauNeg); // momentum of Jpsi at the B+ vertex
Expand Down
3 changes: 1 addition & 2 deletions PWGHF/D2H/TableProducer/candidateCreatorBplusReduced.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -169,8 +169,7 @@
registry.fill(HIST("hCovSVXX"), covMatrixPCA[0]);
registry.fill(HIST("hCovPVXX"), covMatrixPV[0]);

// propagate D0 and Pi to the B+ vertex
df2.propagateTracksToVertex();
// get D and Pi tracks (propagated to the B+ vertex if propagateToPCA==true)
// track.getPxPyPzGlo(pVec) modifies pVec of track
df2.getTrack(0).getPxPyPzGlo(pVecD0); // momentum of D0 at the B+ vertex
df2.getTrack(1).getPxPyPzGlo(pVecPion); // momentum of Pi at the B+ vertex
Expand Down Expand Up @@ -240,7 +239,7 @@
auto candsDThisColl = candsD.sliceBy(candsDPerCollision, thisCollId);
auto tracksPionThisCollision = tracksPion.sliceBy(tracksPionPerCollision, thisCollId);
runCandidateCreation<false>(collision, candsDThisColl, tracksPionThisCollision, invMass2D0PiMin, invMass2D0PiMax);
if (ncol % 10000 == 0) {

Check failure on line 242 in PWGHF/D2H/TableProducer/candidateCreatorBplusReduced.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
LOG(debug) << ncol << " collisions parsed";
}
ncol++;
Expand Down Expand Up @@ -274,7 +273,7 @@
auto candsDThisColl = candsD.sliceBy(candsDPerCollision, thisCollId);
auto tracksPionThisCollision = tracksPion.sliceBy(tracksPionPerCollision, thisCollId);
runCandidateCreation<true>(collision, candsDThisColl, tracksPionThisCollision, invMass2D0PiMin, invMass2D0PiMax);
if (ncol % 10000 == 0) {

Check failure on line 276 in PWGHF/D2H/TableProducer/candidateCreatorBplusReduced.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
LOGP(debug, "collisions parsed {}", ncol);
}
ncol++;
Expand Down
3 changes: 1 addition & 2 deletions PWGHF/D2H/TableProducer/candidateCreatorBsReduced.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -169,8 +169,7 @@
registry.fill(HIST("hCovSVXX"), covMatrixPCA[0]);
registry.fill(HIST("hCovPVXX"), covMatrixPV[0]);

// propagate Ds and Pi to the Bs vertex
df2.propagateTracksToVertex();
// get Ds and Pi tracks (propagated to the B0 vertex if propagateToPCA==true)
// track.getPxPyPzGlo(pVec) modifies pVec of track
df2.getTrack(0).getPxPyPzGlo(pVecD); // momentum of Ds at the Bs vertex
df2.getTrack(1).getPxPyPzGlo(pVecPion); // momentum of Pi at the Bs vertex
Expand Down Expand Up @@ -240,7 +239,7 @@
auto candsDThisColl = candsD.sliceBy(candsDPerCollision, thisCollId);
auto tracksPionThisCollision = tracksPion.sliceBy(tracksPionPerCollision, thisCollId);
runCandidateCreation<false>(collision, candsDThisColl, tracksPionThisCollision, invMass2DPiMin, invMass2DPiMax);
if (ncol % 10000 == 0) {

Check failure on line 242 in PWGHF/D2H/TableProducer/candidateCreatorBsReduced.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
LOGP(debug, "collisions parsed {}", ncol);
}
ncol++;
Expand Down Expand Up @@ -274,7 +273,7 @@
auto candsDThisColl = candsD.sliceBy(candsDPerCollision, thisCollId);
auto tracksPionThisCollision = tracksPion.sliceBy(tracksPionPerCollision, thisCollId);
runCandidateCreation<true>(collision, candsDThisColl, tracksPionThisCollision, invMass2DPiMin, invMass2DPiMax);
if (ncol % 10000 == 0) {

Check failure on line 276 in PWGHF/D2H/TableProducer/candidateCreatorBsReduced.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
LOGP(debug, "collisions parsed {}", ncol);
}
ncol++;
Expand Down
3 changes: 1 addition & 2 deletions PWGHF/D2H/TableProducer/candidateCreatorLbReduced.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -189,8 +189,7 @@ struct HfCandidateCreatorLbReduced {
registry.fill(HIST("hCovSVXX"), covMatrixPCA[0]);
registry.fill(HIST("hCovPVXX"), covMatrixPV[0]);

// propagate Lc and Pi to the Lb vertex
df2.propagateTracksToVertex();
// get Lc and Pi tracks (propagated to the Lb vertex if propagateToPCA==true)
// track.getPxPyPzGlo(pVec) modifies pVec of track
df2.getTrack(0).getPxPyPzGlo(pVecLc); // momentum of Lc at the Lb vertex
df2.getTrack(1).getPxPyPzGlo(pVecPion); // momentum of Pi at the Lb vertex
Expand Down
2 changes: 0 additions & 2 deletions PWGHF/D2H/TableProducer/dataCreatorJpsiHadReduced.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -881,7 +881,6 @@ struct HfDataCreatorJpsiHadReduced {
std::array<float, 3> pVecBPlus{}, pVec0{}, pVec1{}, pVec2{};

auto secondaryVertexBPlus = df3.getPCACandidate();
df3.propagateTracksToVertex();
df3.getTrack(0).getPxPyPzGlo(pVec0);
df3.getTrack(1).getPxPyPzGlo(pVec1);
df3.getTrack(2).getPxPyPzGlo(pVec2);
Expand Down Expand Up @@ -969,7 +968,6 @@ struct HfDataCreatorJpsiHadReduced {
std::array<float, 3> pVecBS{}, pVec0{}, pVec1{}, pVecPhi{};

auto secondaryVertexBS = df4.getPCACandidate();
df4.propagateTracksToVertex();
df4.getTrack(0).getPxPyPzGlo(pVec0);
df4.getTrack(1).getPxPyPzGlo(pVec1);
df4.getTrack(2).getPxPyPzGlo(pVec2);
Expand Down
3 changes: 1 addition & 2 deletions PWGHF/TableProducer/candidateCreatorB0.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -207,7 +207,7 @@
auto primaryVertex = getPrimaryVertex(collision);
auto covMatrixPV = primaryVertex.getCov();

if (ncol % 10000 == 0) {

Check failure on line 210 in PWGHF/TableProducer/candidateCreatorB0.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
LOG(debug) << ncol << " collisions parsed";
}
ncol++;
Expand Down Expand Up @@ -345,8 +345,7 @@
hCovSVXX->Fill(covMatrixPCA[0]);
hCovPVXX->Fill(covMatrixPV[0]);

// propagate D and Pi to the B0 vertex
df2.propagateTracksToVertex();
// get D and Pi tracks (propagated to the B0 vertex if propagateToPCA==true)
// track.getPxPyPzGlo(pVec) modifies pVec of track
df2.getTrack(0).getPxPyPzGlo(pVecD); // momentum of D at the B0 vertex
df2.getTrack(1).getPxPyPzGlo(pVecPion); // momentum of Pi at the B0 vertex
Expand Down
2 changes: 1 addition & 1 deletion PWGHF/TableProducer/candidateCreatorBplus.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -326,7 +326,7 @@ struct HfCandidateCreatorBplus {
}
hCandidatesB->Fill(SVFitting::FitOk);

dfB.propagateTracksToVertex(); // propagate the bachelor and D0 to the B+ vertex
// get D and Pi tracks (propagated to the B+ vertex if propagateToPCA==true)
trackD0.getPxPyPzGlo(pVecD0); // momentum of D0 at the B+ vertex
trackParCovPi.getPxPyPzGlo(pVecBach); // momentum of pi+ at the B+ vertex

Expand Down
3 changes: 1 addition & 2 deletions PWGHF/TableProducer/candidateCreatorBs.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -324,8 +324,7 @@
auto chi2PCA = df2.getChi2AtPCACandidate();
auto covMatrixPCA = df2.calcPCACovMatrixFlat();

// propagate Ds and Pi to the Bs vertex
df2.propagateTracksToVertex();
// get Ds and Pi tracks (propagated to the Bs vertex if propagateToPCA==true)
// track.getPxPyPzGlo(pVec) modifies pVec of track
df2.getTrack(0).getPxPyPzGlo(pVecDs); // momentum of Ds at the Bs vertex
df2.getTrack(1).getPxPyPzGlo(pVecPion); // momentum of Pi at the Bs vertex
Expand Down Expand Up @@ -418,7 +417,7 @@
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughtersDs, Pdg::kDSBar, std::array{-kKPlus, +kKPlus, -kPiPlus}, true, &sign, 2);
if (indexRec > -1) {
RecoDecay::getDaughters(mcParticles.rawIteratorAt(indexRec), &arrDaughDsIndex, std::array{0}, 1);
if (arrDaughDsIndex.size() == 2) {

Check failure on line 420 in PWGHF/TableProducer/candidateCreatorBs.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
for (auto iProng = 0u; iProng < arrDaughDsIndex.size(); ++iProng) {
auto daughI = mcParticles.rawIteratorAt(arrDaughDsIndex[iProng]);
arrPDGDaughDs[iProng] = std::abs(daughI.pdgCode());
Expand Down
1 change: 0 additions & 1 deletion PWGHF/TableProducer/candidateCreatorCascade.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -277,7 +277,6 @@ struct HfCandidateCreatorCascade {
auto chi2PCA = df.getChi2AtPCACandidate();
auto covMatrixPCA = df.calcPCACovMatrixFlat();
registry.fill(HIST("hCovSVXX"), covMatrixPCA[0]); // FIXME: Calculation of errorDecayLength(XY) gives wrong values without this line.
// do I have to call "df.propagateTracksToVertex();"?
auto trackParVarV0 = df.getTrack(0);
auto trackParVarBach = df.getTrack(1);

Expand Down
2 changes: 0 additions & 2 deletions PWGHF/TableProducer/candidateCreatorDstar.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -286,8 +286,6 @@ struct HfCandidateCreatorDstar {
registry.fill(HIST("Refit/hCovSVXZ"), covMatrixPCA[3]);
registry.fill(HIST("Refit/hCovSVZZ"), covMatrixPCA[5]);

// Doubt:................Below, track object are at secondary vertex!
// < track param propagated to V0 candidate (no check for the candidate validity). propagateTracksToVertex must be called in advance
auto trackD0ProngParVar0 = df.getTrack(0);
auto trackD0ProngParVar1 = df.getTrack(1);

Expand Down
2 changes: 1 addition & 1 deletion PWGHF/TableProducer/candidateCreatorLb.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -215,7 +215,7 @@ struct HfCandidateCreatorLb {
auto chi2PCA = df2.getChi2AtPCACandidate();
auto covMatrixPCA = df2.calcPCACovMatrixFlat();

df2.propagateTracksToVertex();
// get Lc and Pi tracks (propagated to the Lb vertex if propagateToPCA==true)
df2.getTrack(0).getPxPyPzGlo(pvecLc);
df2.getTrack(1).getPxPyPzGlo(pvecPion);

Expand Down
7 changes: 1 addition & 6 deletions PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -105,7 +105,7 @@ struct HfCandidateCreatorXic0Omegac0 {
Produces<aod::HfCandToXiPiKfQa> rowKfXic0Qa;
Produces<aod::HfCandToOmegaKaKf> kfCandidateOmegaKaData;

Configurable<bool> propagateToPCA{"propagateToPCA", false, "create tracks version propagated to PCA"};
Configurable<bool> propagateToPCA{"propagateToPCA", true, "create tracks version propagated to PCA"};
Configurable<bool> useAbsDCA{"useAbsDCA", true, "Minimise abs. distance rather than chi2"};
Configurable<bool> useWeightedFinalPCA{"useWeightedFinalPCA", true, "Recalculate vertex position using track covariances, effective only if useAbsDCA is true"};
Configurable<double> maxR{"maxR", 200., "reject PCA's above this radius"};
Expand All @@ -114,7 +114,6 @@ struct HfCandidateCreatorXic0Omegac0 {
Configurable<double> minParamChange{"minParamChange", 1.e-3, "stop iterations if largest change of any X is smaller than this"};
Configurable<double> minRelChi2Change{"minRelChi2Change", 0.9, "stop iterations is chi2/chi2old > this"};
Configurable<double> maxChi2{"maxChi2", 100., "discard vertices with chi2/Nprongs > this (or sum{DCAi^2}/Nprongs for abs. distance minimization)"};
Configurable<bool> refitWithMatCorr{"refitWithMatCorr", true, "when doing propagateTracksToVertex, propagate tracks to vtx with material corrections and rerun minimization"};
Configurable<bool> rejDiffCollTrack{"rejDiffCollTrack", true, "Reject tracks coming from different collisions"};
Configurable<bool> fillAllHist{"fillAllHist", true, "Fill additional KF histograms to check selector cuts"};
Configurable<bool> doCascadePreselection{"doCascadePreselection", true, "Use invariant mass and dcaXY cuts to preselect cascade candidates"};
Expand Down Expand Up @@ -545,10 +544,6 @@ struct HfCandidateCreatorXic0Omegac0 {
auto vertexCharmBaryonFromFitter = df.getPCACandidate();
std::array<float, 3> pVecCascAsD{};
std::array<float, 3> pVecCharmBachelorAsD{};
df.propagateTracksToVertex();
if (!df.isPropagateTracksToVertexDone()) {
continue;
}
df.getTrack(0).getPxPyPzGlo(pVecCascAsD);
df.getTrack(1).getPxPyPzGlo(pVecCharmBachelorAsD);
std::array<float, 3> pVecCharmBaryon = {pVecCascAsD[0] + pVecCharmBachelorAsD[0], pVecCascAsD[1] + pVecCharmBachelorAsD[1], pVecCascAsD[2] + pVecCharmBachelorAsD[2]};
Expand Down
1 change: 0 additions & 1 deletion PWGHF/TableProducer/candidateCreatorXicc.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -173,7 +173,6 @@ struct HfCandidateCreatorXicc {
auto chi2PCA = df2.getChi2AtPCACandidate();
auto covMatrixPCA = df2.calcPCACovMatrixFlat();

df2.propagateTracksToVertex();
df2.getTrack(0).getPxPyPzGlo(pvecxic);
df2.getTrack(1).getPxPyPzGlo(pvecpion);

Expand Down
3 changes: 0 additions & 3 deletions PWGHF/TableProducer/trackIndexSkimCreator.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -3265,7 +3265,6 @@ struct HfTrackIndexSkimCreatorCascades {
if (nCand2 == 0) {
continue;
}
df2.propagateTracksToVertex(); // propagate the bach and V0 to the Lc vertex
df2.getTrack(0).getPxPyPzGlo(pVecV0); // take the momentum at the Lc vertex
df2.getTrack(1).getPxPyPzGlo(pVecBach);
}
Expand Down Expand Up @@ -3333,7 +3332,6 @@ struct HfTrackIndexSkimCreatorLfCascades {
Configurable<double> minParamChange{"minParamChange", 1.e-3, "stop iterations if largest change of any X is smaller than this"};
Configurable<double> minRelChi2Change{"minRelChi2Change", 0.9, "stop iterations is chi2/chi2old > this"};
Configurable<double> maxChi2{"maxChi2", 100., "discard vertices with chi2/Nprongs > this (or sum{DCAi^2}/Nprongs for abs. distance minimization)"};
Configurable<bool> refitWithMatCorr{"refitWithMatCorr", true, "when doing propagateTracksToVertex, propagate tracks to vtx with material corrections and rerun minimization"};

// Selection criteria
// selections have been set to run2 lambda dedicated cuts
Expand Down Expand Up @@ -3604,7 +3602,6 @@ struct HfTrackIndexSkimCreatorLfCascades {
const auto magneticField = o2::base::Propagator::Instance()->getNominalBz(); // z component

df2.setBz(magneticField);
df2.setRefitWithMatCorr(config.refitWithMatCorr);

// cascade loop
const auto thisCollId = collision.globalIndex();
Expand Down
1 change: 0 additions & 1 deletion PWGHF/TableProducer/treeCreatorTccToD0D0Pi.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -599,7 +599,6 @@ struct HfTreeCreatorTccToD0D0Pi {
}
hCandidatesTcc->Fill(SVFitting::FitOk);

dfTcc.propagateTracksToVertex(); // propagate the softpi and D0 pair to the Tcc vertex
trackD1.getPxPyPzGlo(pVecD1); // momentum of D1 at the Tcc vertex
trackD2.getPxPyPzGlo(pVecD2); // momentum of D2 at the Tcc vertex
trackParCovPi.getPxPyPzGlo(pVecSoftPi); // momentum of pi at the Tcc vertex
Expand Down
Loading