Skip to content

[PWGUD] UPC Cand Tabel producer GlobalMuon changed propagate method #16549

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
mochimochino wants to merge 1 commit into
base: master
Choose a base branch
from
+48 −3
Open

[PWGUD] UPC Cand Tabel producer GlobalMuon changed propagate method #16549

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
51 changes: 48 additions & 3 deletions PWGUD/TableProducer/upcCandProducerGlobalMuon.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -90,7 +90,7 @@ struct UpcCandProducerGlobalMuon {
Configurable<int> fBcWindowMCH{"fBcWindowMCH", 20, "Time window for MCH-MID to MCH-only matching for Muon candidates"};
Configurable<float> fMaxFV0Amp{"fMaxFV0Amp", 100.f, "Max FV0 amplitude in the same BC"};

// NEW: MFT/Global track support configurables
// MFT/Global track support configurables
Configurable<bool> fEnableMFT{"fEnableMFT", true, "Enable MFT/global track processing"};
Configurable<bool> fSaveMFTClusters{"fSaveMFTClusters", true, "Save MFT cluster information"};

Expand Down Expand Up @@ -369,6 +369,7 @@ struct UpcCandProducerGlobalMuon {
scrollBackForth(inGbc, maxDbc, mapBcs, fillAmps);
}

// Propagate MCH track to z = 0
auto propagateToZero(ForwardTracks::iterator const& muon)
{
using SMatrix55 = ROOT::Math::SMatrix<double, 5, 5, ROOT::Math::MatRepSym<double, 5>>;
Expand All @@ -392,6 +393,45 @@ struct UpcCandProducerGlobalMuon {
return propmuon;
}

// Propagate MCH track to target z point (set MFT z)
auto propagateMCHtoZ(ForwardTracks::iterator const& muon, double targetZ)
{
using SMatrix55 = ROOT::Math::SMatrix<double, 5, 5, ROOT::Math::MatRepSym<double, 5>>;
using SMatrix5 = ROOT::Math::SVector<double, 5>;
SMatrix5 tpars(muon.x(), muon.y(), muon.phi(), muon.tgl(), muon.signed1Pt());
std::vector<double> v1{muon.cXX(), muon.cXY(), muon.cYY(), muon.cPhiX(), muon.cPhiY(),
muon.cPhiPhi(), muon.cTglX(), muon.cTglY(), muon.cTglPhi(), muon.cTglTgl(),
muon.c1PtX(), muon.c1PtY(), muon.c1PtPhi(), muon.c1PtTgl(), muon.c1Pt21Pt2()};
SMatrix55 tcovs(v1.begin(), v1.end());
o2::dataformats::GlobalFwdTrack propmuon;
o2::dataformats::GlobalFwdTrack track;
track.setParameters(tpars);
track.setZ(muon.z());
track.setCovariances(tcovs);
auto mchTrack = fMatching.FwdtoMCH(track);
// Here is different point with propagateToZero
o2::mch::TrackExtrap::extrapToVertex(mchTrack, 0., 0., targetZ, 0., 0.);
auto proptrack = fMatching.MCHtoFwd(mchTrack);
o2::dataformats::GlobalFwdTrack propMuonAtZ;
propMuonAtZ.setParameters(proptrack.getParameters());
propMuonAtZ.setZ(proptrack.getZ());
propMuonAtZ.setCovariances(proptrack.getCovariances());
return propMuonAtZ;
}

// Propagate Global track to z = 0
auto propagateGlobalFwdToZero(const o2::dataformats::GlobalFwdTrack& globalTrack)
{
auto mchTrack = fMatching.FwdtoMCH(globalTrack);
o2::mch::TrackExtrap::extrapToZCov(mchTrack, 0.);
auto proptrack = fMatching.MCHtoFwd(mchTrack);
o2::dataformats::GlobalFwdTrack propMuonAtZero;
propMuonAtZero.setParameters(proptrack.getParameters());
propMuonAtZero.setZ(proptrack.getZ());
propMuonAtZero.setCovariances(proptrack.getCovariances());
return propMuonAtZero;
}

bool addToFwdTable(int64_t candId, int64_t trackId, uint64_t gbc, float trackTime,
ForwardTracks const& fwdTracks, o2::aod::MFTTracks const& /*mftTracks*/,
const o2::aod::McFwdTrackLabels* mcFwdTrackLabels)
Expand All @@ -414,14 +454,19 @@ struct UpcCandProducerGlobalMuon {
if (isGlobal && fEnableMFT) {
// Refit global muon: propagate MCH component to vertex, combine with MFT spatial info
auto mchTrack = track.matchMCHTrack_as<ForwardTracks>();
auto propMuon = propagateToZero(mchTrack);
// Propgate to Z of MFT -> z =0
auto mfttrack = track.matchMFTTrack_as<o2::aod::MFTTracks>();
double zMFT = mfttrack.z();

auto propMuonAtMFT = propagateMCHtoZ(mchTrack, zMFT);
using SMatrix5 = ROOT::Math::SVector<double, 5>;
using SMatrix55 = ROOT::Math::SMatrix<double, 5, 5, ROOT::Math::MatRepSym<double, 5>>;
SMatrix5 tpars(mfttrack.x(), mfttrack.y(), mfttrack.phi(), mfttrack.tgl(), mfttrack.signed1Pt());
SMatrix55 tcovs{};
o2::track::TrackParCovFwd mft{mfttrack.z(), tpars, tcovs, mfttrack.chi2()};
pft = o2::aod::fwdtrackutils::refitGlobalMuonCov(propMuon, mft);
auto refitAtMFT = o2::aod::fwdtrackutils::refitGlobalMuonCov(propMuonAtMFT, mft);
// To z = 0
pft = propagateGlobalFwdToZero(refitAtMFT);
} else {
pft = propagateToZero(track);
}
Expand Down
Loading