Skip to content

[PWGEM/Dilepton] fix TTCA for global muons #15026

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
dsekihat merged 1 commit into from
Feb 18, 2026
Merged

[PWGEM/Dilepton] fix TTCA for global muons #15026

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
29 changes: 20 additions & 9 deletions PWGEM/Dilepton/TableProducer/skimmerPrimaryMuon.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -595,11 +595,16 @@ struct skimmerPrimaryMuon {

// LOGF(info, "muon_tmp.globalIndex() = %d, muon_tmp.matchMCHTrackId() = %d, muon_tmp.matchMFTTrackId() = %d, muon_tmp.chi2MatchMCHMFT() = %f", muon_tmp.globalIndex(), muon_tmp.matchMCHTrackId(), muon_tmp.matchMFTTrackId(), muon_tmp.chi2MatchMCHMFT());

if (cfgApplyPreselectionInBestMatch && !isSelected(pt, eta, muon_tmp.rAtAbsorberEnd(), pDCA, muon_tmp.chi2() / ndf, muon_tmp.trackType(), dcaXY)) {
continue;
}
if (std::sqrt(std::pow(deta / maxDEta, 2) + std::pow(dphi / maxDPhi, 2)) > 1.f) {
continue;
if (cfgApplyPreselectionInBestMatch) {
if (!isSelected(pt, eta, muon_tmp.rAtAbsorberEnd(), pDCA, muon_tmp.chi2() / ndf, muon_tmp.trackType(), dcaXY)) {
continue;
}
if (std::sqrt(std::pow(deta / maxDEta, 2) + std::pow(dphi / maxDPhi, 2)) > 1.f) {
continue;
}
if (muon_tmp.chi2MatchMCHMFT() > maxMatchingChi2MCHMFT) {
continue;
}
}

if (0.f < muon_tmp.chi2MatchMCHMFT() && muon_tmp.chi2MatchMCHMFT() < min_chi2MatchMCHMFT) {
Expand Down Expand Up @@ -762,11 +767,14 @@ struct skimmerPrimaryMuon {
}

const auto& fwdtrackIdsThisCollision = fwdtrackIndices.sliceBy(fwdtrackIndicesPerCollision, collision.globalIndex());
// LOGF(info, "collision.globalIndex() = %d, fwdtrackIdsThisCollision.size() = %d", collision.globalIndex(), fwdtrackIdsThisCollision.size());
for (const auto& fwdtrackId : fwdtrackIdsThisCollision) {
const auto& fwdtrack = fwdtrackId.template fwdtrack_as<MyFwdTracks>();
// LOGF(info, "fwdtrack.globalIndex() = %d, fwdtrack.matchMCHTrackId() = %d, fwdtrack.matchMFTTrackId() = %d, fwdtrack.trackType() = %d", fwdtrack.globalIndex(), fwdtrack.matchMCHTrackId(), fwdtrack.matchMFTTrackId(), fwdtrack.trackType());
if (fwdtrack.trackType() != o2::aod::fwdtrack::ForwardTrackTypeEnum::GlobalMuonTrack && fwdtrack.trackType() != o2::aod::fwdtrack::ForwardTrackTypeEnum::MuonStandaloneTrack) {
continue;
}

if (fwdtrack.trackType() == o2::aod::fwdtrack::ForwardTrackTypeEnum::GlobalMuonTrack && std::find(vec_min_chi2MatchMCHMFT.begin(), vec_min_chi2MatchMCHMFT.end(), std::make_tuple(fwdtrack.globalIndex(), fwdtrack.matchMCHTrackId(), fwdtrack.matchMFTTrackId())) == vec_min_chi2MatchMCHMFT.end()) {
continue;
}
Expand Down Expand Up @@ -1476,6 +1484,7 @@ struct associateSameMFT {
PresliceUnsorted<aod::EMPrimaryMuons> perMFTTrack = o2::aod::emprimarymuon::mfttrackId;
std::vector<int> self_Ids;

// Multiple MCH-MID tracks can match with the same MFTsa. This function is to reject such global muons.
void process(aod::EMPrimaryMuons const& muons)
{
for (const auto& muon : muons) {
Expand All @@ -1486,17 +1495,19 @@ struct associateSameMFT {
if (global_muon.globalIndex() == muon.globalIndex()) { // don't store myself.
continue;
}
self_Ids.emplace_back(global_muon.globalIndex());
// self_Ids.emplace_back(global_muon.globalIndex());

// if (global_muon.collisionId() == muon.collisionId()) {
// self_Ids.emplace_back(global_muon.globalIndex());
// }
if (global_muon.collisionId() == muon.collisionId()) { // the same global muon is repeatedly stored and associated to different collisions if FTTCA is used.
self_Ids.emplace_back(global_muon.globalIndex());
}
}
em_same_mft_ids(self_Ids);
self_Ids.clear();
self_Ids.shrink_to_fit();
} else {
em_same_mft_ids(std::vector<int>{}); // empty for standalone muons
self_Ids.clear();
self_Ids.shrink_to_fit();
}
} // end of muon loop
}
Expand Down
8 changes: 3 additions & 5 deletions PWGEM/Dilepton/Utils/EMTrackUtilities.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -144,11 +144,9 @@ bool isBestMatch(TTrack const& track, TCut const& cut, TTracks const& tracks)
map_chi2MCHMFT[track.globalIndex()] = track.chi2MatchMCHMFT(); // add myself
for (const auto& glmuonId : track.globalMuonsWithSameMFTIds()) {
const auto& candidate = tracks.rawIteratorAt(glmuonId);
if (track.mchtrackId() != candidate.mchtrackId() && track.mfttrackId() == candidate.mfttrackId()) {
if (candidate.trackType() == o2::aod::fwdtrack::ForwardTrackTypeEnum::GlobalMuonTrack) {
if (cut.template IsSelectedTrack<is_wo_acc>(candidate)) {
map_chi2MCHMFT[candidate.globalIndex()] = candidate.chi2MatchMCHMFT();
}
if (candidate.trackType() == o2::aod::fwdtrack::ForwardTrackTypeEnum::GlobalMuonTrack && candidate.emeventId() == track.emeventId()) {
if (cut.template IsSelectedTrack<is_wo_acc>(candidate)) {
map_chi2MCHMFT[candidate.globalIndex()] = candidate.chi2MatchMCHMFT();
}
}
} // end of glmuonId
Expand Down
Loading