Skip to content

Add ITS fake clusters information to the mcMask #14666

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
shahor02 merged 2 commits into from
Sep 25, 2025
Merged

Add ITS fake clusters information to the mcMask #14666

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
be70ea7
Add ITS fake clusters information to the mcMask
Sep 10, 2025
4ae732a
remove redundant check
Sep 10, 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
17 changes: 16 additions & 1 deletion Detectors/AOD/src/AODProducerWorkflowSpec.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1126,7 +1126,7 @@ void AODProducerWorkflowDPL::fillMCTrackLabelsTable(MCTrackLabelCursorType& mcTr
if (!needToStore(mGIDToTableID)) {
continue;
}
if (mcTruth.isValid()) { // if not set, -1 will be stored
if (mcTruth.isValid()) { // if not set, -1 will be stored
labelHolder.labelID = (mToStore[mcTruth.getSourceID()][mcTruth.getEventID()])[mcTruth.getTrackID()]; // defined by TPC if it contributes, otherwise: by ITS
if (mcTruth.isFake()) {
labelHolder.labelMask |= (0x1 << 15);
Expand All @@ -1139,6 +1139,21 @@ void AODProducerWorkflowDPL::fillMCTrackLabelsTable(MCTrackLabelCursorType& mcTr
}
}
}
if (trackIndex.includesDet(DetID::ITS)) {
auto itsGID = data.getITSContributorGID(trackIndex);
auto itsSource = itsGID.getSource();
if (itsSource == GIndex::ITS) {
auto& itsTrack = data.getITSTrack(itsGID);
for (unsigned int iL = 0; iL < 7; ++iL) {
if (itsTrack.isFakeOnLayer(iL)) {
labelHolder.labelMask |= (0x1 << iL);
}
}
} else if (itsSource == GIndex::ITSAB) {
labelHolder.labelMask |= (data.getTrackMCLabel(itsGID).isFake() << 12);
}
}

} else if (mcTruth.isNoise()) {
labelHolder.labelMask |= (0x1 << 14);
}
Expand Down
2 changes: 1 addition & 1 deletion Framework/Core/include/Framework/AnalysisDataModel.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2017,7 +2017,7 @@ namespace aod
namespace mctracklabel
{
DECLARE_SOA_INDEX_COLUMN(McParticle, mcParticle); //! MC particle
DECLARE_SOA_COLUMN(McMask, mcMask, uint16_t); //! Bit mask to indicate detector mismatches (bit ON means mismatch). Bit 0-6: mismatch at ITS layer. Bit 7-9: # of TPC mismatches in the ranges 0, 1, 2-3, 4-7, 8-15, 16-31, 32-63, >64. Bit 10: TRD, bit 11: TOF, bit 15: indicates negative label
DECLARE_SOA_COLUMN(McMask, mcMask, uint16_t); //! Bit mask to indicate detector mismatches (bit ON means mismatch). Bit 0-6: mismatch at ITS layer. Bit 12: ITSAB tracklet mismatch. Bit 13: ITS-TPC mismatch. Bit 14: isNoise == True (global track), Bit 15: isFake == True (global track)
} // namespace mctracklabel

DECLARE_SOA_TABLE(McTrackLabels, "AOD", "MCTRACKLABEL", //! Table joined to the track table containing the MC index
Expand Down