Skip to content

[PWGHF] Add more D variables to B0 reduced task and MLResponse #17532

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
fchinu wants to merge 1 commit into
base: master
Choose a base branch
from
Open

[PWGHF] Add more D variables to B0 reduced task and MLResponse #17532

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
39 changes: 35 additions & 4 deletions PWGHF/Core/HfMlResponseB0ToDPi.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -30,10 +30,9 @@
// Fill the map of available input features
// the key is the feature's name (std::string)
// the value is the corresponding value in EnumInputFeatures
#define FILL_MAP_B0(FEATURE) \
{ \
#FEATURE, static_cast<uint8_t>(InputFeaturesB0ToDPi::FEATURE) \
}
#define FILL_MAP_B0(FEATURE) \
{ \
#FEATURE, static_cast<uint8_t>(InputFeaturesB0ToDPi::FEATURE)}

// Check if the index of mCachedIndices (index associated to a FEATURE)
// matches the entry in EnumInputFeatures associated to this FEATURE
Expand Down Expand Up @@ -90,6 +89,14 @@ enum class InputFeaturesB0ToDPi : uint8_t {
decayLengthXYNormalised,
cpa,
cpaXY,
decayLengthDToPv,
decayLengthDToPvXY,
decayLengthDToB,
decayLengthDToBXY,
cpaDToPv,
cpaDToPvXY,
cpaDToB,
cpaDToBXY,
maxNormalisedDeltaIP,
prong0MlScoreBkg,
prong0MlScorePrompt,
Expand Down Expand Up @@ -139,6 +146,14 @@ class HfMlResponseB0ToDPi : public HfMlResponse<TypeOutputScore>
CHECK_AND_FILL_VEC_B0(decayLengthXYNormalised);
CHECK_AND_FILL_VEC_B0(cpa);
CHECK_AND_FILL_VEC_B0(cpaXY);
CHECK_AND_FILL_VEC_B0(decayLengthDToPv);
CHECK_AND_FILL_VEC_B0(decayLengthDToPvXY);
CHECK_AND_FILL_VEC_B0(decayLengthDToB);
CHECK_AND_FILL_VEC_B0(decayLengthDToBXY);
CHECK_AND_FILL_VEC_B0(cpaDToPv);
CHECK_AND_FILL_VEC_B0(cpaDToPvXY);
CHECK_AND_FILL_VEC_B0(cpaDToB);
CHECK_AND_FILL_VEC_B0(cpaDToBXY);
CHECK_AND_FILL_VEC_B0(maxNormalisedDeltaIP);
// TPC PID variable
CHECK_AND_FILL_VEC_B0_FULL(prongBachPi, tpcNSigmaPi1, tpcNSigmaPi);
Expand Down Expand Up @@ -202,6 +217,14 @@ class HfMlResponseB0ToDPi : public HfMlResponse<TypeOutputScore>
CHECK_AND_FILL_VEC_B0(decayLengthXYNormalised);
CHECK_AND_FILL_VEC_B0(cpa);
CHECK_AND_FILL_VEC_B0(cpaXY);
CHECK_AND_FILL_VEC_B0(decayLengthDToPv);
CHECK_AND_FILL_VEC_B0(decayLengthDToPvXY);
CHECK_AND_FILL_VEC_B0(decayLengthDToB);
CHECK_AND_FILL_VEC_B0(decayLengthDToBXY);
CHECK_AND_FILL_VEC_B0(cpaDToPv);
CHECK_AND_FILL_VEC_B0(cpaDToPvXY);
CHECK_AND_FILL_VEC_B0(cpaDToB);
CHECK_AND_FILL_VEC_B0(cpaDToBXY);
CHECK_AND_FILL_VEC_B0(maxNormalisedDeltaIP);
// TPC PID variable
CHECK_AND_FILL_VEC_B0_FULL(prongBachPi, tpcNSigmaPiBachPi, tpcNSigmaPi);
Expand Down Expand Up @@ -257,6 +280,14 @@ class HfMlResponseB0ToDPi : public HfMlResponse<TypeOutputScore>
FILL_MAP_B0(decayLengthXYNormalised),
FILL_MAP_B0(cpa),
FILL_MAP_B0(cpaXY),
FILL_MAP_B0(decayLengthDToPv),
FILL_MAP_B0(decayLengthDToPvXY),
FILL_MAP_B0(decayLengthDToB),
FILL_MAP_B0(decayLengthDToBXY),
FILL_MAP_B0(cpaDToPv),
FILL_MAP_B0(cpaDToPvXY),
FILL_MAP_B0(cpaDToB),
FILL_MAP_B0(cpaDToBXY),
FILL_MAP_B0(maxNormalisedDeltaIP),
FILL_MAP_B0(prong0MlScoreBkg),
FILL_MAP_B0(prong0MlScorePrompt),
Expand Down
6 changes: 4 additions & 2 deletions PWGHF/D2H/TableProducer/candidateCreatorB0Reduced.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -225,7 +225,8 @@
pVecSoftPion[0], pVecSoftPion[1], pVecSoftPion[2],
pVecBachPion[0], pVecBachPion[1], pVecBachPion[2],
dcaD.getY(), dcaSoftPion.getY(), dcaBachPion.getY(),
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaSoftPion.getSigmaY2()), std::sqrt(dcaBachPion.getSigmaY2()));
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaSoftPion.getSigmaY2()), std::sqrt(dcaBachPion.getSigmaY2()),
candD.xSecondaryVertex(), candD.ySecondaryVertex(), candD.zSecondaryVertex());

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

Expand Down Expand Up @@ -327,7 +328,8 @@
pVecD[0], pVecD[1], pVecD[2],
pVecPion[0], pVecPion[1], pVecPion[2],
dcaD.getY(), dcaPion.getY(),
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaPion.getSigmaY2()));
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaPion.getSigmaY2()),
candD.xSecondaryVertex(), candD.ySecondaryVertex(), candD.zSecondaryVertex());

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

Expand Down Expand Up @@ -363,7 +365,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 368 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 @@ -397,7 +399,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 402 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 @@ -432,7 +434,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 437 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 @@ -467,7 +469,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 472 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
32 changes: 28 additions & 4 deletions PWGHF/D2H/Tasks/taskB0Reduced.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -90,11 +90,17 @@ DECLARE_SOA_COLUMN(DecayLengthNormalised, decayLengthNormalised, float);
DECLARE_SOA_COLUMN(DecayLengthXYNormalised, decayLengthXYNormalised, float); //! Normalised transverse decay length of candidate
DECLARE_SOA_COLUMN(DecayLengthD, decayLengthD, float); //! Decay length of D-meson daughter candidate (cm)
DECLARE_SOA_COLUMN(DecayLengthXYD, decayLengthXYD, float); //! Transverse decay length of D-meson daughter candidate (cm)
DECLARE_SOA_COLUMN(DecayLengthDToB, decayLengthDToB, float); //! Decay length of D-meson daughter candidate from B decay vertex (cm)
DECLARE_SOA_COLUMN(DecayLengthXYDToB, decayLengthXYDToB, float); //! Transverse decay length of D-meson daughter candidate from B decay vertex (cm)
DECLARE_SOA_COLUMN(ImpactParameterD, impactParameterD, float); //! Impact parameter product of D-meson daughter candidate
DECLARE_SOA_COLUMN(ImpactParameterBach, impactParameterBach, float); //! Impact parameter product of bachelor pion
DECLARE_SOA_COLUMN(ImpactParameterProduct, impactParameterProduct, float); //! Impact parameter product of daughters
DECLARE_SOA_COLUMN(Cpa, cpa, float); //! Cosine pointing angle of candidate
DECLARE_SOA_COLUMN(CpaXY, cpaXY, float); //! Cosine pointing angle of candidate in transverse plane
DECLARE_SOA_COLUMN(CpaD, cpaD, float); //! Cosine pointing angle of D candidate to PV
DECLARE_SOA_COLUMN(CpaDXY, cpaDXY, float); //! Cosine pointing angle of D candidate to PV in transverse plane
DECLARE_SOA_COLUMN(CpaDToB, cpaDToB, float); //! Cosine pointing angle of D candidate to B decay vertex
DECLARE_SOA_COLUMN(CpaDToBXY, cpaDToBXY, float); //! Cosine pointing angle of D candidate to B decay vertex in transverse plane
DECLARE_SOA_COLUMN(MaxNormalisedDeltaIP, maxNormalisedDeltaIP, float); //! Maximum normalized difference between measured and expected impact parameter of candidate prongs
DECLARE_SOA_COLUMN(MlScoreSig, mlScoreSig, float); //! ML score for signal class
DECLARE_SOA_COLUMN(FlagWrongCollision, flagWrongCollision, int8_t); //! Flag for association with wrong collision
Expand Down Expand Up @@ -123,6 +129,12 @@ DECLARE_SOA_TABLE(HfRedCandB0Lites, "AOD", "HFREDCANDB0LITE", //! Table with som
hf_cand_b0_lite::PtD,
hf_cand_b0_lite::DecayLengthD,
hf_cand_b0_lite::DecayLengthXYD,
hf_cand_b0_lite::DecayLengthDToB,
hf_cand_b0_lite::DecayLengthXYDToB,
hf_cand_b0_lite::CpaD,
hf_cand_b0_lite::CpaDXY,
hf_cand_b0_lite::CpaDToB,
hf_cand_b0_lite::CpaDToBXY,
hf_cand_b0_lite::ImpactParameterD,
hf_cand_b0_lite::PtDmesProngMin,
hf_cand_b0_lite::AbsEtaDmesProngMin,
Expand Down Expand Up @@ -621,8 +633,14 @@ struct HfTaskB0Reduced {
// D-meson features
invMassD,
ptD,
decLenD,
decLenXyD,
candidate.decayLengthDToPv(),
candidate.decayLengthDToPvXY(),
candidate.decayLengthDToB(),
candidate.decayLengthDToBXY(),
candidate.cpaDToPv(),
candidate.cpaDToPvXY(),
candidate.cpaDToB(),
candidate.cpaDToBXY(),
candidate.impactParameter0(),
candD.ptProngMin(),
candD.absEtaProngMin(),
Expand Down Expand Up @@ -878,8 +896,14 @@ struct HfTaskB0Reduced {
// D-meson features
invMassD,
ptD,
decLenD,
decLenXyD,
candidate.decayLengthDToPv(),
candidate.decayLengthDToPvXY(),
candidate.decayLengthDToB(),
candidate.decayLengthDToBXY(),
candidate.cpaDToPv(),
candidate.cpaDToPvXY(),
candidate.cpaDToB(),
candidate.cpaDToBXY(),
candidate.impactParameter0(),
candD.ptProngMin(),
candD.absEtaProngMin(),
Expand Down
39 changes: 39 additions & 0 deletions PWGHF/DataModel/CandidateReconstructionTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1687,6 +1687,9 @@ DECLARE_SOA_TABLE(HfCandLbMcGen, "AOD", "HFCANDLBMCGEN", //!
// specific B0 candidate properties
namespace hf_cand_b0
{
DECLARE_SOA_COLUMN(XSecondaryVertexD, xSecondaryVertexD, float);
DECLARE_SOA_COLUMN(YSecondaryVertexD, ySecondaryVertexD, float);
DECLARE_SOA_COLUMN(ZSecondaryVertexD, zSecondaryVertexD, float);
DECLARE_SOA_INDEX_COLUMN_FULL(Prong0, prong0, int, HfCand3Prong, "_0"); // D index
DECLARE_SOA_DYNAMIC_COLUMN(ImpactParameterProduct, impactParameterProduct, // Impact parameter product for B0 -> J/Psi K*0
[](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 {
Expand All @@ -1698,6 +1701,22 @@ DECLARE_SOA_DYNAMIC_COLUMN(ImpactParameterProductJpsi, impactParameterProductJps
[](float dcaDauPos, float dcaDauNeg) -> float { return dcaDauPos * dcaDauNeg; });
DECLARE_SOA_DYNAMIC_COLUMN(ImpactParameterProductK0Star, impactParameterProductK0Star, // K*0 impact parameter product for B0 -> J/Psi K*0
[](float dcaLfTrack0, float dcaLfTrack1) -> float { return dcaLfTrack0 * dcaLfTrack1; });
DECLARE_SOA_DYNAMIC_COLUMN(CpaDToPv, cpaDToPv, //!
[](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}); });
DECLARE_SOA_DYNAMIC_COLUMN(CpaDToPvXY, cpaDToPvXY, //!
[](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}); });
DECLARE_SOA_DYNAMIC_COLUMN(CpaDToB, cpaDToB, //!
[](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}); });
DECLARE_SOA_DYNAMIC_COLUMN(CpaDToBXY, cpaDToBXY, //!
[](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}); });
DECLARE_SOA_DYNAMIC_COLUMN(DecayLengthDToPv, decayLengthDToPv, //!
[](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}); });
DECLARE_SOA_DYNAMIC_COLUMN(DecayLengthDToPvXY, decayLengthDToPvXY, //!
[](float xVtxP, float yVtxP, float xVtxS, float yVtxS) -> float { return RecoDecay::distanceXY(std::array{xVtxP, yVtxP}, std::array{xVtxS, yVtxS}); });
DECLARE_SOA_DYNAMIC_COLUMN(DecayLengthDToB, decayLengthDToB, //!
[](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}); });
DECLARE_SOA_DYNAMIC_COLUMN(DecayLengthDToBXY, decayLengthDToBXY, //!
[](float xVtxB, float yVtxB, float xVtxS, float yVtxS) -> float { return RecoDecay::distanceXY(std::array{xVtxB, yVtxB}, std::array{xVtxS, yVtxS}); });

enum DecayTypeMc : uint8_t { B0ToDplusPiToPiKPiPi = 0,
B0ToDsPiToKKPiPi,
Expand All @@ -1718,7 +1737,17 @@ DECLARE_SOA_TABLE(HfCandB0Base, "AOD", "HFCANDB0BASE",
hf_cand::PxProng1, hf_cand::PyProng1, hf_cand::PzProng1,
hf_cand::ImpactParameter0, hf_cand::ImpactParameter1,
hf_cand::ErrorImpactParameter0, hf_cand::ErrorImpactParameter1,
// D meson decay vertex position
hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD,
/* dynamic columns */
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>,
hf_cand_b0::CpaDToPvXY<collision::PosX, collision::PosY, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0>,
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>,
hf_cand_b0::CpaDToBXY<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0>,
hf_cand_b0::DecayLengthDToPv<collision::PosX, collision::PosY, collision::PosZ, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD>,
hf_cand_b0::DecayLengthDToPvXY<collision::PosX, collision::PosY, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD>,
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>,
hf_cand_b0::DecayLengthDToBXY<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD>,
hf_cand_2prong::M<hf_cand::PxProng0, hf_cand::PyProng0, hf_cand::PzProng0, hf_cand::PxProng1, hf_cand::PyProng1, hf_cand::PzProng1>,
hf_cand_2prong::M2<hf_cand::PxProng0, hf_cand::PyProng0, hf_cand::PzProng0, hf_cand::PxProng1, hf_cand::PyProng1, hf_cand::PzProng1>,
hf_cand_2prong::ImpactParameterProduct<hf_cand::ImpactParameter0, hf_cand::ImpactParameter1>,
Expand Down Expand Up @@ -1755,7 +1784,17 @@ DECLARE_SOA_TABLE(HfCandB0DStar, "AOD", "HFCANDB0DSTAR",
hf_cand::PxProng2, hf_cand::PyProng2, hf_cand::PzProng2,
hf_cand::ImpactParameter0, hf_cand::ImpactParameter1, hf_cand::ImpactParameter2,
hf_cand::ErrorImpactParameter0, hf_cand::ErrorImpactParameter1, hf_cand::ErrorImpactParameter2,
// D meson decay vertex position
hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD,
/* dynamic columns */
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>,
hf_cand_b0::CpaDToPvXY<collision::PosX, collision::PosY, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0>,
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>,
hf_cand_b0::CpaDToBXY<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand::PxProng0, hf_cand::PyProng0>,
hf_cand_b0::DecayLengthDToPv<collision::PosX, collision::PosY, collision::PosZ, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD, hf_cand_b0::ZSecondaryVertexD>,
hf_cand_b0::DecayLengthDToPvXY<collision::PosX, collision::PosY, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD>,
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>,
hf_cand_b0::DecayLengthDToBXY<hf_cand::XSecondaryVertex, hf_cand::YSecondaryVertex, hf_cand_b0::XSecondaryVertexD, hf_cand_b0::YSecondaryVertexD>,
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>,
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>,
hf_cand_2prong::ImpactParameterProduct<hf_cand::ImpactParameter0, hf_cand::ImpactParameter2>,
Expand Down
3 changes: 2 additions & 1 deletion PWGHF/TableProducer/candidateCreatorB0.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -201,7 +201,7 @@
auto primaryVertex = getPrimaryVertex(collision);
auto covMatrixPV = primaryVertex.getCov();

if (ncol % 10000 == 0) {

Check failure on line 204 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 @@ -373,7 +373,8 @@
pVecD[0], pVecD[1], pVecD[2],
pVecPion[0], pVecPion[1], pVecPion[2],
dcaD.getY(), dcaPion.getY(),
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaPion.getSigmaY2()));
std::sqrt(dcaD.getSigmaY2()), std::sqrt(dcaPion.getSigmaY2()),
candD.xSecondaryVertex(), candD.ySecondaryVertex(), candD.zSecondaryVertex());

rowCandidateProngs(candD.globalIndex(), trackPion.globalIndex());
} // pi loop
Expand Down
Loading