Skip to content

[PWGLF] Try to use PV and SV to correct momentum of short tracks #12024

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

Merged
maciacco merged 4 commits into from
Jul 12, 2025
Merged

[PWGLF] Try to use PV and SV to correct momentum of short tracks #12024

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
8be92e0
Try to use PV and SV to correct track momentum
wang-yuanzhe Jul 11, 2025
bdd798a
Do propagation before introducing new point
wang-yuanzhe Jul 11, 2025
82978de
Fix bug
wang-yuanzhe Jul 11, 2025
d815da5
Fix comparison of different type
wang-yuanzhe Jul 12, 2025
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
5 changes: 5 additions & 0 deletions PWGLF/DataModel/LFHyperhelium4sigmaTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -36,6 +36,9 @@ DECLARE_SOA_COLUMN(ZMoth, zMoth, float); // Z o
DECLARE_SOA_COLUMN(PxMoth, pxMoth, float); //! Px of the mother track at the decay vertex
DECLARE_SOA_COLUMN(PyMoth, pyMoth, float); //! Py of the mother track at the decay vertex
DECLARE_SOA_COLUMN(PzMoth, pzMoth, float); //! Pz of the mother track at the decay vertex
DECLARE_SOA_COLUMN(RefitPxMoth, refitPxMoth, float); //! Refit Px of the mother track at the decay vertex
DECLARE_SOA_COLUMN(RefitPyMoth, refitPyMoth, float); //! Refit Py of the mother track at the decay vertex
DECLARE_SOA_COLUMN(RefitPzMoth, refitPzMoth, float); //! Refit Pz of the mother track at the decay vertex
DECLARE_SOA_COLUMN(PxAlpha, pxAlpha, float); //! Px of the daughter alpha track at the decay vertex
DECLARE_SOA_COLUMN(PyAlpha, pyAlpha, float); //! Py of the daughter alpha track at the decay vertex
DECLARE_SOA_COLUMN(PzAlpha, pzAlpha, float); //! Pz of the daughter alpha track at the decay vertex
Expand Down Expand Up @@ -77,6 +80,7 @@ DECLARE_SOA_TABLE(He4S2BCands, "AOD", "HE4S2BCANDS",
he4scand::IsMatter,
he4scand::XMoth, he4scand::YMoth, he4scand::ZMoth,
he4scand::PxMoth, he4scand::PyMoth, he4scand::PzMoth,
he4scand::RefitPxMoth, he4scand::RefitPyMoth, he4scand::RefitPzMoth,
he4scand::PxAlpha, he4scand::PyAlpha, he4scand::PzAlpha,
he4scand::DcaMothPv, he4scand::DcaAlphaPv, he4scand::DcaKinkTopo,
he4scand::ItsChi2Moth, he4scand::ItsClusterSizesMoth, he4scand::ItsClusterSizesAlpha,
Expand All @@ -89,6 +93,7 @@ DECLARE_SOA_TABLE(MCHe4S2BCands, "AOD", "MCHE4S2BCANDS",
he4scand::IsMatter,
he4scand::XMoth, he4scand::YMoth, he4scand::ZMoth,
he4scand::PxMoth, he4scand::PyMoth, he4scand::PzMoth,
he4scand::RefitPxMoth, he4scand::RefitPyMoth, he4scand::RefitPzMoth,
he4scand::PxAlpha, he4scand::PyAlpha, he4scand::PzAlpha,
he4scand::DcaMothPv, he4scand::DcaAlphaPv, he4scand::DcaKinkTopo,
he4scand::ItsChi2Moth, he4scand::ItsClusterSizesMoth, he4scand::ItsClusterSizesAlpha,
Expand Down
52 changes: 52 additions & 0 deletions PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -80,6 +80,8 @@ constexpr std::array<int, kNDaughterType> kDaughterPDG = {
PDG_t::kNeutron,
PDG_t::kPi0};

const std::array<float, 6> covPosSV{6.4462712107237135f, 0.1309793068144521f, 6.626654155592929f, -0.4510297694023185f, 0.16996629627762413f, 4.109195981415627f};

std::shared_ptr<TH1> hMotherCounter;
std::shared_ptr<TH1> hMother2BCounter;
std::shared_ptr<TH1> hDauCounter[kNChargedDaughterType];
Expand Down Expand Up @@ -224,6 +226,41 @@ float getTPCNSigma(const TTrack& track, const int daughterType)
return nSigma;
}

//--------------------------------------------------------------
// Refit the momentum of the mother track
template <typename TCollision, typename TTrack>
std::array<float, 3> refitMotherTrack(TCollision& collision, TTrack& track, std::array<float, 3> posSV)
{
auto trackPar = getTrackParCov(track);

float trackIUPos[2] = {track.y(), track.z()};
float trackIUCov[3] = {track.cYY(), track.cZY(), track.cZZ()};
o2::dataformats::VertexBase primaryVtx = {{collision.posX(), collision.posY(), collision.posZ()}, {collision.covXX(), collision.covXY(), collision.covYY(), collision.covXZ(), collision.covYZ(), collision.covZZ()}};
o2::dataformats::VertexBase secondaryVtx = {{posSV[0], posSV[1], posSV[2]}, {covPosSV[0], covPosSV[1], covPosSV[2], covPosSV[3], covPosSV[4], covPosSV[5]}};

o2::base::Propagator::Instance()->propagateToDCABxByBz(primaryVtx, trackPar, 2.f, o2::base::Propagator::MatCorrType::USEMatCorrLUT);

trackPar.resetCovariance(999.f);
std::array<float, 3> refitP = {-999.f, -999.f, -999.f};
if (!trackPar.update(primaryVtx, 999.f)) {
return refitP;
}

trackPar.rotate(track.alpha());
o2::base::Propagator::Instance()->PropagateToXBxByBz(trackPar, track.x());
if (!trackPar.update(trackIUPos, trackIUCov)) {
return refitP;
}

o2::base::Propagator::Instance()->propagateToDCABxByBz(secondaryVtx, trackPar, 2.f, o2::base::Propagator::MatCorrType::USEMatCorrLUT);
if (!trackPar.update(secondaryVtx, 999.f)) {
return refitP;
}

trackPar.getPxPyPzGlo(refitP);
return refitP;
}

//--------------------------------------------------------------
struct Hyphe4sCandidate {

Expand Down Expand Up @@ -475,12 +512,19 @@ struct Hyperhelium4sigmaRecoTask {
Hyphe4sCandidate hyphe4sCand;
fillCandidate(hyphe4sCand, collision, kinkCand, motherTrack, dauTrack);

std::array<float, 3> posDecVtx = {kinkCand.xDecVtx() + collision.posX(), kinkCand.yDecVtx() + collision.posY(), kinkCand.zDecVtx() + collision.posZ()};
auto refitP = refitMotherTrack(collision, motherTrack, posDecVtx);
for (size_t i = 0; i < refitP.size(); ++i) {
refitP[i] *= 2.f;
}

outputDataTable(
hyphe4sCand.primVtx[0], hyphe4sCand.primVtx[1], hyphe4sCand.primVtx[2],
hyphe4sCand.decVtx[0], hyphe4sCand.decVtx[1], hyphe4sCand.decVtx[2],
hyphe4sCand.isMatter,
hyphe4sCand.lastPosMoth[0], hyphe4sCand.lastPosMoth[1], hyphe4sCand.lastPosMoth[2],
hyphe4sCand.momMoth[0], hyphe4sCand.momMoth[1], hyphe4sCand.momMoth[2],
refitP[0], refitP[1], refitP[2],
hyphe4sCand.momDaug[0], hyphe4sCand.momDaug[1], hyphe4sCand.momDaug[2],
hyphe4sCand.dcaXYMothPv, hyphe4sCand.dcaXYDauPv, hyphe4sCand.dcaKinkTopo,
hyphe4sCand.chi2ITSMoth, hyphe4sCand.itsClusterSizeMoth, hyphe4sCand.itsClusterSizeDau,
Expand Down Expand Up @@ -554,6 +598,12 @@ struct Hyperhelium4sigmaRecoTask {
Hyphe4sCandidate hyphe4sCand;
fillCandidate(hyphe4sCand, collision, kinkCand, motherTrack, dauTrack);

std::array<float, 3> posDecVtx = {kinkCand.xDecVtx() + collision.posX(), kinkCand.yDecVtx() + collision.posY(), kinkCand.zDecVtx() + collision.posZ()};
auto refitP = refitMotherTrack(collision, motherTrack, posDecVtx);
for (size_t i = 0; i < refitP.size(); ++i) {
refitP[i] *= 2.f;
}

// qa for true signal
if (isTrueSignal) {
auto mcMotherTrack = motherTrack.mcParticle_as<aod::McParticles>();
Expand Down Expand Up @@ -620,6 +670,7 @@ struct Hyperhelium4sigmaRecoTask {
hyphe4sCand.isMatter,
hyphe4sCand.lastPosMoth[0], hyphe4sCand.lastPosMoth[1], hyphe4sCand.lastPosMoth[2],
hyphe4sCand.momMoth[0], hyphe4sCand.momMoth[1], hyphe4sCand.momMoth[2],
refitP[0], refitP[1], refitP[2],
hyphe4sCand.momDaug[0], hyphe4sCand.momDaug[1], hyphe4sCand.momDaug[2],
hyphe4sCand.dcaXYMothPv, hyphe4sCand.dcaXYDauPv, hyphe4sCand.dcaKinkTopo,
hyphe4sCand.chi2ITSMoth, hyphe4sCand.itsClusterSizeMoth, hyphe4sCand.itsClusterSizeDau,
Expand Down Expand Up @@ -661,6 +712,7 @@ struct Hyperhelium4sigmaRecoTask {
-1, -1, -1,
-1, -1, -1,
-1, -1, -1,
-1, -1, -1,
-1, -1,
true, false, isReconstructedMCCollisions[mcparticle.mcCollisionId()], isSelectedMCCollisions[mcparticle.mcCollisionId()],
hyphe4sCand.trueDecVtx[0], hyphe4sCand.trueDecVtx[1], hyphe4sCand.trueDecVtx[2],
Expand Down
Loading