Skip to content

Commit 814256a

Browse files
fchinufchinu
authored
[PWGHF] Add more D variables to B0 reduced task and MLResponse (#17532)
1 parent e4940f1 commit 814256a

5 files changed

Lines changed: 108 additions & 11 deletions

File tree

‎PWGHF/Core/HfMlResponseB0ToDPi.h‎

Lines changed: 35 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -30,10 +30,9 @@
3030
// Fill the map of available input features
3131
// the key is the feature's name (std::string)
3232
// the value is the corresponding value in EnumInputFeatures
33-
#define FILL_MAP_B0(FEATURE) \
34-
{ \
35-
#FEATURE, static_cast<uint8_t>(InputFeaturesB0ToDPi::FEATURE) \
36-
}
33+
#define FILL_MAP_B0(FEATURE) \
34+
{ \
35+
#FEATURE, static_cast<uint8_t>(InputFeaturesB0ToDPi::FEATURE)}
3736

3837
// Check if the index of mCachedIndices (index associated to a FEATURE)
3938
// matches the entry in EnumInputFeatures associated to this FEATURE
@@ -90,6 +89,14 @@ enum class InputFeaturesB0ToDPi : uint8_t {
9089
decayLengthXYNormalised,
9190
cpa,
9291
cpaXY,
92+
decayLengthDToPv,
93+
decayLengthDToPvXY,
94+
decayLengthDToB,
95+
decayLengthDToBXY,
96+
cpaDToPv,
97+
cpaDToPvXY,
98+
cpaDToB,
99+
cpaDToBXY,
93100
maxNormalisedDeltaIP,
94101
prong0MlScoreBkg,
95102
prong0MlScorePrompt,
@@ -139,6 +146,14 @@ class HfMlResponseB0ToDPi : public HfMlResponse<TypeOutputScore>
139146
CHECK_AND_FILL_VEC_B0(decayLengthXYNormalised);
140147
CHECK_AND_FILL_VEC_B0(cpa);
141148
CHECK_AND_FILL_VEC_B0(cpaXY);
149+
CHECK_AND_FILL_VEC_B0(decayLengthDToPv);
150+
CHECK_AND_FILL_VEC_B0(decayLengthDToPvXY);
151+
CHECK_AND_FILL_VEC_B0(decayLengthDToB);
152+
CHECK_AND_FILL_VEC_B0(decayLengthDToBXY);
153+
CHECK_AND_FILL_VEC_B0(cpaDToPv);
154+
CHECK_AND_FILL_VEC_B0(cpaDToPvXY);
155+
CHECK_AND_FILL_VEC_B0(cpaDToB);
156+
CHECK_AND_FILL_VEC_B0(cpaDToBXY);
142157
CHECK_AND_FILL_VEC_B0(maxNormalisedDeltaIP);
143158
// TPC PID variable
144159
CHECK_AND_FILL_VEC_B0_FULL(prongBachPi, tpcNSigmaPi1, tpcNSigmaPi);
@@ -202,6 +217,14 @@ class HfMlResponseB0ToDPi : public HfMlResponse<TypeOutputScore>
202217
CHECK_AND_FILL_VEC_B0(decayLengthXYNormalised);
203218
CHECK_AND_FILL_VEC_B0(cpa);
204219
CHECK_AND_FILL_VEC_B0(cpaXY);
220+
CHECK_AND_FILL_VEC_B0(decayLengthDToPv);
221+
CHECK_AND_FILL_VEC_B0(decayLengthDToPvXY);
222+
CHECK_AND_FILL_VEC_B0(decayLengthDToB);
223+
CHECK_AND_FILL_VEC_B0(decayLengthDToBXY);
224+
CHECK_AND_FILL_VEC_B0(cpaDToPv);
225+
CHECK_AND_FILL_VEC_B0(cpaDToPvXY);
226+
CHECK_AND_FILL_VEC_B0(cpaDToB);
227+
CHECK_AND_FILL_VEC_B0(cpaDToBXY);
205228
CHECK_AND_FILL_VEC_B0(maxNormalisedDeltaIP);
206229
// TPC PID variable
207230
CHECK_AND_FILL_VEC_B0_FULL(prongBachPi, tpcNSigmaPiBachPi, tpcNSigmaPi);
@@ -257,6 +280,14 @@ class HfMlResponseB0ToDPi : public HfMlResponse<TypeOutputScore>
257280
FILL_MAP_B0(decayLengthXYNormalised),
258281
FILL_MAP_B0(cpa),
259282
FILL_MAP_B0(cpaXY),
283+
FILL_MAP_B0(decayLengthDToPv),
284+
FILL_MAP_B0(decayLengthDToPvXY),
285+
FILL_MAP_B0(decayLengthDToB),
286+
FILL_MAP_B0(decayLengthDToBXY),
287+
FILL_MAP_B0(cpaDToPv),
288+
FILL_MAP_B0(cpaDToPvXY),
289+
FILL_MAP_B0(cpaDToB),
290+
FILL_MAP_B0(cpaDToBXY),
260291
FILL_MAP_B0(maxNormalisedDeltaIP),
261292
FILL_MAP_B0(prong0MlScoreBkg),
262293
FILL_MAP_B0(prong0MlScorePrompt),

‎PWGHF/D2H/TableProducer/candidateCreatorB0Reduced.cxx‎

Lines changed: 4 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -225,7 +225,8 @@ struct HfCandidateCreatorB0Reduced {
225225
pVecSoftPion[0], pVecSoftPion[1], pVecSoftPion[2],
226226
pVecBachPion[0], pVecBachPion[1], pVecBachPion[2],
227227
dcaD.getY(), dcaSoftPion.getY(), dcaBachPion.getY(),
228-
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaSoftPion.getSigmaY2()), std::sqrt(dcaBachPion.getSigmaY2()));
228+
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaSoftPion.getSigmaY2()), std::sqrt(dcaBachPion.getSigmaY2()),
229+
candD.xSecondaryVertex(), candD.ySecondaryVertex(), candD.zSecondaryVertex());
229230

230231
rowCandidateProngsDStar(candD.globalIndex(), trackBachPion.globalIndex(), trackSoftPion.globalIndex());
231232

@@ -327,7 +328,8 @@ struct HfCandidateCreatorB0Reduced {
327328
pVecD[0], pVecD[1], pVecD[2],
328329
pVecPion[0], pVecPion[1], pVecPion[2],
329330
dcaD.getY(), dcaPion.getY(),
330-
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaPion.getSigmaY2()));
331+
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaPion.getSigmaY2()),
332+
candD.xSecondaryVertex(), candD.ySecondaryVertex(), candD.zSecondaryVertex());
331333

332334
rowCandidateProngs(candD.globalIndex(), trackPion.globalIndex());
333335

‎PWGHF/D2H/Tasks/taskB0Reduced.cxx‎

Lines changed: 28 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -90,11 +90,17 @@ DECLARE_SOA_COLUMN(DecayLengthNormalised, decayLengthNormalised, float);
9090
DECLARE_SOA_COLUMN(DecayLengthXYNormalised, decayLengthXYNormalised, float); //! Normalised transverse decay length of candidate
9191
DECLARE_SOA_COLUMN(DecayLengthD, decayLengthD, float); //! Decay length of D-meson daughter candidate (cm)
9292
DECLARE_SOA_COLUMN(DecayLengthXYD, decayLengthXYD, float); //! Transverse decay length of D-meson daughter candidate (cm)
93+
DECLARE_SOA_COLUMN(DecayLengthDToB, decayLengthDToB, float); //! Decay length of D-meson daughter candidate from B decay vertex (cm)
94+
DECLARE_SOA_COLUMN(DecayLengthXYDToB, decayLengthXYDToB, float); //! Transverse decay length of D-meson daughter candidate from B decay vertex (cm)
9395
DECLARE_SOA_COLUMN(ImpactParameterD, impactParameterD, float); //! Impact parameter product of D-meson daughter candidate
9496
DECLARE_SOA_COLUMN(ImpactParameterBach, impactParameterBach, float); //! Impact parameter product of bachelor pion
9597
DECLARE_SOA_COLUMN(ImpactParameterProduct, impactParameterProduct, float); //! Impact parameter product of daughters
9698
DECLARE_SOA_COLUMN(Cpa, cpa, float); //! Cosine pointing angle of candidate
9799
DECLARE_SOA_COLUMN(CpaXY, cpaXY, float); //! Cosine pointing angle of candidate in transverse plane
100+
DECLARE_SOA_COLUMN(CpaD, cpaD, float); //! Cosine pointing angle of D candidate to PV
101+
DECLARE_SOA_COLUMN(CpaDXY, cpaDXY, float); //! Cosine pointing angle of D candidate to PV in transverse plane
102+
DECLARE_SOA_COLUMN(CpaDToB, cpaDToB, float); //! Cosine pointing angle of D candidate to B decay vertex
103+
DECLARE_SOA_COLUMN(CpaDToBXY, cpaDToBXY, float); //! Cosine pointing angle of D candidate to B decay vertex in transverse plane
98104
DECLARE_SOA_COLUMN(MaxNormalisedDeltaIP, maxNormalisedDeltaIP, float); //! Maximum normalized difference between measured and expected impact parameter of candidate prongs
99105
DECLARE_SOA_COLUMN(MlScoreSig, mlScoreSig, float); //! ML score for signal class
100106
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); //! Flag for association with wrong collision
@@ -123,6 +129,12 @@ DECLARE_SOA_TABLE(HfRedCandB0Lites, "AOD", "HFREDCANDB0LITE", //! Table with som
123129
hf_cand_b0_lite::PtD,
124130
hf_cand_b0_lite::DecayLengthD,
125131
hf_cand_b0_lite::DecayLengthXYD,
132+
hf_cand_b0_lite::DecayLengthDToB,
133+
hf_cand_b0_lite::DecayLengthXYDToB,
134+
hf_cand_b0_lite::CpaD,
135+
hf_cand_b0_lite::CpaDXY,
136+
hf_cand_b0_lite::CpaDToB,
137+
hf_cand_b0_lite::CpaDToBXY,
126138
hf_cand_b0_lite::ImpactParameterD,
127139
hf_cand_b0_lite::PtDmesProngMin,
128140
hf_cand_b0_lite::AbsEtaDmesProngMin,
@@ -621,8 +633,14 @@ struct HfTaskB0Reduced {
621633
// D-meson features
622634
invMassD,
623635
ptD,
624-
decLenD,
625-
decLenXyD,
636+
candidate.decayLengthDToPv(),
637+
candidate.decayLengthDToPvXY(),
638+
candidate.decayLengthDToB(),
639+
candidate.decayLengthDToBXY(),
640+
candidate.cpaDToPv(),
641+
candidate.cpaDToPvXY(),
642+
candidate.cpaDToB(),
643+
candidate.cpaDToBXY(),
626644
candidate.impactParameter0(),
627645
candD.ptProngMin(),
628646
candD.absEtaProngMin(),
@@ -878,8 +896,14 @@ struct HfTaskB0Reduced {
878896
// D-meson features
879897
invMassD,
880898
ptD,
881-
decLenD,
882-
decLenXyD,
899+
candidate.decayLengthDToPv(),
900+
candidate.decayLengthDToPvXY(),
901+
candidate.decayLengthDToB(),
902+
candidate.decayLengthDToBXY(),
903+
candidate.cpaDToPv(),
904+
candidate.cpaDToPvXY(),
905+
candidate.cpaDToB(),
906+
candidate.cpaDToBXY(),
883907
candidate.impactParameter0(),
884908
candD.ptProngMin(),
885909
candD.absEtaProngMin(),

‎PWGHF/DataModel/CandidateReconstructionTables.h‎

Lines changed: 39 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1687,6 +1687,9 @@ DECLARE_SOA_TABLE(HfCandLbMcGen, "AOD", "HFCANDLBMCGEN", //!
16871687
// specific B0 candidate properties
16881688
namespace hf_cand_b0
16891689
{
1690+
DECLARE_SOA_COLUMN(XSecondaryVertexD, xSecondaryVertexD, float);
1691+
DECLARE_SOA_COLUMN(YSecondaryVertexD, ySecondaryVertexD, float);
1692+
DECLARE_SOA_COLUMN(ZSecondaryVertexD, zSecondaryVertexD, float);
16901693
DECLARE_SOA_INDEX_COLUMN_FULL(Prong0, prong0, int, HfCand3Prong, "_0"); // D index
16911694
DECLARE_SOA_DYNAMIC_COLUMN(ImpactParameterProduct, impactParameterProduct, // Impact parameter product for B0 -> J/Psi K*0
16921695
[](float pxJpsiDauPos, float pyJpsiDauPos, float pzJpsiDauPos, float pxJpsiDauNeg, float pyJpsiDauNeg, float pzJpsiDauNeg, float pxLfTrack0, float pyLfTrack0, float pzLfTrack0, float pxLfTrack1, float pyLfTrack1, float pzLfTrack1, float xVtxP, float yVtxP, float zVtxP, float xVtxS, float yVtxS, float zVtxS) -> float {
@@ -1698,6 +1701,22 @@ DECLARE_SOA_DYNAMIC_COLUMN(ImpactParameterProductJpsi, impactParameterProductJps
16981701
[](float dcaDauPos, float dcaDauNeg) -> float { return dcaDauPos * dcaDauNeg; });
16991702
DECLARE_SOA_DYNAMIC_COLUMN(ImpactParameterProductK0Star, impactParameterProductK0Star, // K*0 impact parameter product for B0 -> J/Psi K*0
17001703
[](float dcaLfTrack0, float dcaLfTrack1) -> float { return dcaLfTrack0 * dcaLfTrack1; });
1704+
DECLARE_SOA_DYNAMIC_COLUMN(CpaDToPv, cpaDToPv, //!
1705+
[](float xVtxP, float yVtxP, float zVtxP, float xVtxS, float yVtxS, float zVtxS, float px, float py, float pz) -> float { return RecoDecay::cpa(std::array{xVtxP, yVtxP, zVtxP}, std::array{xVtxS, yVtxS, zVtxS}, std::array{px, py, pz}); });
1706+
DECLARE_SOA_DYNAMIC_COLUMN(CpaDToPvXY, cpaDToPvXY, //!
1707+
[](float xVtxP, float yVtxP, float xVtxS, float yVtxS, float px, float py) -> float { return RecoDecay::cpaXY(std::array{xVtxP, yVtxP}, std::array{xVtxS, yVtxS}, std::array{px, py}); });
1708+
DECLARE_SOA_DYNAMIC_COLUMN(CpaDToB, cpaDToB, //!
1709+
[](float xVtxB, float yVtxB, float zVtxB, float xVtxS, float yVtxS, float zVtxS, float px, float py, float pz) -> float { return RecoDecay::cpa(std::array{xVtxB, yVtxB, zVtxB}, std::array{xVtxS, yVtxS, zVtxS}, std::array{px, py, pz}); });
1710+
DECLARE_SOA_DYNAMIC_COLUMN(CpaDToBXY, cpaDToBXY, //!
1711+
[](float xVtxB, float yVtxB, float xVtxS, float yVtxS, float px, float py) -> float { return RecoDecay::cpaXY(std::array{xVtxB, yVtxB}, std::array{xVtxS, yVtxS}, std::array{px, py}); });
1712+
DECLARE_SOA_DYNAMIC_COLUMN(DecayLengthDToPv, decayLengthDToPv, //!
1713+
[](float xVtxP, float yVtxP, float zVtxP, float xVtxS, float yVtxS, float zVtxS) -> float { return RecoDecay::distance(std::array{xVtxP, yVtxP, zVtxP}, std::array{xVtxS, yVtxS, zVtxS}); });
1714+
DECLARE_SOA_DYNAMIC_COLUMN(DecayLengthDToPvXY, decayLengthDToPvXY, //!
1715+
[](float xVtxP, float yVtxP, float xVtxS, float yVtxS) -> float { return RecoDecay::distanceXY(std::array{xVtxP, yVtxP}, std::array{xVtxS, yVtxS}); });
1716+
DECLARE_SOA_DYNAMIC_COLUMN(DecayLengthDToB, decayLengthDToB, //!
1717+
[](float xVtxB, float yVtxB, float zVtxB, float xVtxS, float yVtxS, float zVtxS) -> float { return RecoDecay::distance(std::array{xVtxB, yVtxB, zVtxB}, std::array{xVtxS, yVtxS, zVtxS}); });
1718+
DECLARE_SOA_DYNAMIC_COLUMN(DecayLengthDToBXY, decayLengthDToBXY, //!
1719+
[](float xVtxB, float yVtxB, float xVtxS, float yVtxS) -> float { return RecoDecay::distanceXY(std::array{xVtxB, yVtxB}, std::array{xVtxS, yVtxS}); });
17011720

17021721
enum DecayTypeMc : uint8_t { B0ToDplusPiToPiKPiPi = 0,
17031722
B0ToDsPiToKKPiPi,
@@ -1718,7 +1737,17 @@ DECLARE_SOA_TABLE(HfCandB0Base, "AOD", "HFCANDB0BASE",
17181737
hf_cand::PxProng1, hf_cand::PyProng1, hf_cand::PzProng1,
17191738
hf_cand::ImpactParameter0, hf_cand::ImpactParameter1,
17201739
hf_cand::ErrorImpactParameter0, hf_cand::ErrorImpactParameter1,
1740+
// D meson decay vertex position
1741+
hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD,
17211742
/* dynamic columns */
1743+
hf_cand_b0::CpaDToPv<collision::PosX, collision::PosY, collision::PosZ, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0, hf_cand::PzProng0>,
1744+
hf_cand_b0::CpaDToPvXY<collision::PosX, collision::PosY, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0>,
1745+
hf_cand_b0::CpaDToB<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand::ZSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0, hf_cand::PzProng0>,
1746+
hf_cand_b0::CpaDToBXY<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0>,
1747+
hf_cand_b0::DecayLengthDToPv<collision::PosX, collision::PosY, collision::PosZ, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD>,
1748+
hf_cand_b0::DecayLengthDToPvXY<collision::PosX, collision::PosY, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD>,
1749+
hf_cand_b0::DecayLengthDToB<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand::ZSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD>,
1750+
hf_cand_b0::DecayLengthDToBXY<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD>,
17221751
hf_cand_2prong::M<hf_cand::PxProng0, hf_cand::PyProng0, hf_cand::PzProng0, hf_cand::PxProng1, hf_cand::PyProng1, hf_cand::PzProng1>,
17231752
hf_cand_2prong::M2<hf_cand::PxProng0, hf_cand::PyProng0, hf_cand::PzProng0, hf_cand::PxProng1, hf_cand::PyProng1, hf_cand::PzProng1>,
17241753
hf_cand_2prong::ImpactParameterProduct<hf_cand::ImpactParameter0, hf_cand::ImpactParameter1>,
@@ -1755,7 +1784,17 @@ DECLARE_SOA_TABLE(HfCandB0DStar, "AOD", "HFCANDB0DSTAR",
17551784
hf_cand::PxProng2, hf_cand::PyProng2, hf_cand::PzProng2,
17561785
hf_cand::ImpactParameter0, hf_cand::ImpactParameter1, hf_cand::ImpactParameter2,
17571786
hf_cand::ErrorImpactParameter0, hf_cand::ErrorImpactParameter1, hf_cand::ErrorImpactParameter2,
1787+
// D meson decay vertex position
1788+
hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD,
17581789
/* dynamic columns */
1790+
hf_cand_b0::CpaDToPv<collision::PosX, collision::PosY, collision::PosZ, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0, hf_cand::PzProng0>,
1791+
hf_cand_b0::CpaDToPvXY<collision::PosX, collision::PosY, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0>,
1792+
hf_cand_b0::CpaDToB<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand::ZSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0, hf_cand::PzProng0>,
1793+
hf_cand_b0::CpaDToBXY<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0>,
1794+
hf_cand_b0::DecayLengthDToPv<collision::PosX, collision::PosY, collision::PosZ, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD>,
1795+
hf_cand_b0::DecayLengthDToPvXY<collision::PosX, collision::PosY, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD>,
1796+
hf_cand_b0::DecayLengthDToB<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand::ZSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD>,
1797+
hf_cand_b0::DecayLengthDToBXY<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD>,
17591798
hf_cand_3prong::M<hf_cand::PxProng0, hf_cand::PyProng0, hf_cand::PzProng0, hf_cand::PxProng1, hf_cand::PyProng1, hf_cand::PzProng1, hf_cand::PxProng2, hf_cand::PyProng2, hf_cand::PzProng2>,
17601799
hf_cand_3prong::M2<hf_cand::PxProng0, hf_cand::PyProng0, hf_cand::PzProng0, hf_cand::PxProng1, hf_cand::PyProng1, hf_cand::PzProng1, hf_cand::PxProng2, hf_cand::PyProng2, hf_cand::PzProng2>,
17611800
hf_cand_2prong::ImpactParameterProduct<hf_cand::ImpactParameter0, hf_cand::ImpactParameter2>,

‎PWGHF/TableProducer/candidateCreatorB0.cxx‎

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -373,7 +373,8 @@ struct HfCandidateCreatorB0 {
373373
pVecD[0], pVecD[1], pVecD[2],
374374
pVecPion[0], pVecPion[1], pVecPion[2],
375375
dcaD.getY(), dcaPion.getY(),
376-
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaPion.getSigmaY2()));
376+
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaPion.getSigmaY2()),
377+
candD.xSecondaryVertex(), candD.ySecondaryVertex(), candD.zSecondaryVertex());
377378

378379
rowCandidateProngs(candD.globalIndex(), trackPion.globalIndex());
379380
} // pi loop

0 commit comments

Comments
 (0)