Skip to content

[PWGDQ] fix MFT momentum rescaling - #17461

Draft
aferrero2707 wants to merge 2 commits into
AliceO2Group:masterfrom
aferrero2707:PWGDQ-fix-mft-momentum-rescaling
Draft

[PWGDQ] fix MFT momentum rescaling#17461
aferrero2707 wants to merge 2 commits into
AliceO2Group:masterfrom
aferrero2707:PWGDQ-fix-mft-momentum-rescaling

Conversation

@aferrero2707

@aferrero2707 aferrero2707 commented Aug 14, 2026

Copy link
Copy Markdown
Contributor

The MCH tracks are extrapolated to the first measured MFT point using TrackExtrap::extrapToVertex() instead of TrackExtrap::extrapToVertexWithoutBranson(), which cures large inaccuracies in the MCH momentum estimation at the MFT front for some tracks.
Those inaccuracies were contributing to high-mass tails in the di-muon invariant mass forward tracks with rescaled momentum.

Below is a comparison of the J/psi invariant mass before and after the fix, estimated from the same AO2Ds of LHC25i4 (OO MC simulation with injected J/psi and psi(2S)).

Before:
image

After:
image

The MCH tracks are extrapolated to the first measured MFT point using
TrackExtrap::extrapToVertex() instead of TrackExtrap::extrapToVertexWithoutBranson(),
which cures large inaccuracies in the MCH momentum estimation at the MFT
front for some tracks. Those inaccuracies were contributing to high-mass
tails in the di-muon invariant mass distribution reconstructed from global
forward tracks with rescaled momentum.
@github-actions

Copy link
Copy Markdown

O2 linter results: ❌ 0 errors, ⚠️ 61 warnings, 🔕 12 disabled

@mcoquet642

Copy link
Copy Markdown
Collaborator

why not use the common tool for rescaling introduced in fwdtrackUtilities.h ?

@aferrero2707

Copy link
Copy Markdown
Contributor Author

why not use the common tool for rescaling introduced in fwdtrackUtilities.h ?

Mostly to have the freedom of experimenting without affecting the common tools used in analysis.
For this specific case, if I understand correctly one could use refitGlobalMuonCov() for the momentum rescaling, right? However, the issue is in the determination of the MCH momentum to be used for the rescaling, so the common function would not help.

Once the fix is confirmed we can see how to properly propagate it to the common code. What do you think?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

Development

Successfully merging this pull request may close these issues.

2 participants