Skip to content

[PWGHF] taskSingleMuonMultMc: Minor fixes. Add in CMakeLists. #15147

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
vkucera merged 14 commits into from
Feb 25, 2026
Merged

[PWGHF] taskSingleMuonMultMc: Minor fixes. Add in CMakeLists. #15147

Show file tree
Hide file tree
Changes from 10 commits
Commits
Show all changes
14 commits
Select commit Hold shift + click to select a range
7e57a6e
Histograms updated
mislam17 Feb 10, 2026
5afacdf
task for single muon identification in MC simulation and Axe calculation
mislam17 Feb 10, 2026
c9ebb4a
Delete PWGHF/HFL/Tasks/taskSingleMuonMultMc.cxx
mislam17 Feb 10, 2026
5ec9dd1
Merge branch 'AliceO2Group:master' into singleMuonMult
mislam17 Feb 10, 2026
67cbd74
Merge branch 'AliceO2Group:master' into singleMuonMult
mislam17 Feb 10, 2026
b9e344a
Merge branch 'AliceO2Group:master' into singleMuonMult
mislam17 Feb 23, 2026
a1b958f
Task to identify various sources of single muons and the calculation …
mislam17 Feb 23, 2026
7a0f3c8
Add task-single-muon-mult-mc workflow to CMakeLists
mislam17 Feb 24, 2026
84f6625
Refactor mask parameter type from uint16_t& to uint16_t
mislam17 Feb 24, 2026
b93e06a
Merge branch 'master' into singleMuonMult
mislam17 Feb 24, 2026
d41d1ad
Add task-single-muon-mult-mc workflow
mislam17 Feb 25, 2026
56cda86
Header files are updated
mislam17 Feb 25, 2026
defc916
Fix formatting in CMakeLists.txt
mislam17 Feb 25, 2026
984e87f
Fixed include format
mislam17 Feb 25, 2026
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 PWGHF/HFL/Tasks/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -48,3 +48,8 @@ o2physics_add_dpl_workflow(task-single-muon-source
SOURCES taskSingleMuonSource.cxx
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(task-single-muon-mult-mc
Comment thread
vkucera marked this conversation as resolved.
Outdated
SOURCES taskSingleMuonMultMc.cxx
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)
49 changes: 29 additions & 20 deletions PWGHF/HFL/Tasks/taskSingleMuonMultMc.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -25,6 +25,7 @@
#include "Framework/O2DatabasePDGPlugin.h"
#include "Framework/runDataProcessing.h"
#include "ReconstructionDataFormats/TrackFwd.h"
#include <Framework/ASoA.h>
Comment thread
vkucera marked this conversation as resolved.
#include <Framework/Configurable.h>

#include <TString.h>
Expand Down Expand Up @@ -257,90 +258,90 @@ struct HfTaskSingleMuonMultMc {
}

// particle has an associated MC particle
bool isIdentified(const uint16_t& mask)
bool isIdentified(const uint16_t mask)
{
return (TESTBIT(mask, IsIdentified));
}
// this particle is muon
bool isMuon(const uint16_t& mask)
bool isMuon(const uint16_t mask)
{
return (TESTBIT(mask, IsIdentified) && TESTBIT(mask, IsMuon));
}

// this muon comes from transport
bool isSecondaryMu(const uint16_t& mask)
bool isSecondaryMu(const uint16_t mask)
{
return (isMuon(mask) && TESTBIT(mask, IsSecondary));
}

// this muon comes from light flavor quark decay
bool isLightDecayMu(const uint16_t& mask)
bool isLightDecayMu(const uint16_t mask)
{
return (isMuon(mask) && TESTBIT(mask, HasLightParent) && (!TESTBIT(mask, IsSecondary)));
}

// this muon comes from tau decays
bool isTauDecayMu(const uint16_t& mask)
bool isTauDecayMu(const uint16_t mask)
{
return (isMuon(mask) && TESTBIT(mask, HasTauParent) && (!TESTBIT(mask, HasWParent)) && (!TESTBIT(mask, HasZParent)) && (!TESTBIT(mask, HasBeautyParent)) && (!TESTBIT(mask, HasCharmParent)));
}

// this muon comes from W+- decay
bool isWBosonDecayMu(const uint16_t& mask)
bool isWBosonDecayMu(const uint16_t mask)
{
return (isMuon(mask) && TESTBIT(mask, HasWParent) && (!TESTBIT(mask, HasZParent)) && (!TESTBIT(mask, HasTauParent)) && (!TESTBIT(mask, HasLightParent)) && (!TESTBIT(mask, IsSecondary)));
}

// this muon comes from Z decay
bool isZBosonDecayMu(const uint16_t& mask)
bool isZBosonDecayMu(const uint16_t mask)
{
return (isMuon(mask) && TESTBIT(mask, HasZParent) && (!TESTBIT(mask, HasWParent)) && (!TESTBIT(mask, HasTauParent)) && (!TESTBIT(mask, HasLightParent)) && (!TESTBIT(mask, IsSecondary)));
}

// this muon comes from quarkonium decay
bool isQuarkoniumDecayMu(const uint16_t& mask)
bool isQuarkoniumDecayMu(const uint16_t mask)
{
return (isMuon(mask) && TESTBIT(mask, HasQuarkoniumParent) && (!TESTBIT(mask, HasBeautyParent)) && (!TESTBIT(mask, HasCharmParent)) && (!TESTBIT(mask, HasLightParent)));
}

// this muon comes from beauty decay and does not have light flavor parent
bool isBeautyMu(const uint16_t& mask)
bool isBeautyMu(const uint16_t mask)
{
return (isMuon(mask) && TESTBIT(mask, HasBeautyParent) && (!TESTBIT(mask, HasQuarkoniumParent)) && (!TESTBIT(mask, HasLightParent)) && (!TESTBIT(mask, IsSecondary)));
}

// this muon comes directly from beauty decay
bool isBeautyDecayMu(const uint16_t& mask)
bool isBeautyDecayMu(const uint16_t mask)
{
return (isBeautyMu(mask) && (!TESTBIT(mask, HasCharmParent) && (!TESTBIT(mask, HasQuarkoniumParent))) && (!TESTBIT(mask, HasLightParent)) && (!TESTBIT(mask, IsSecondary)));
}

// this muon comes from non-prompt charm decay and does not have light flavor parent
bool isNonpromptCharmMu(const uint16_t& mask)
bool isNonpromptCharmMu(const uint16_t mask)
{
return (isBeautyMu(mask) && TESTBIT(mask, HasCharmParent) && (!TESTBIT(mask, HasQuarkoniumParent)) && (!TESTBIT(mask, HasLightParent)) && (!TESTBIT(mask, IsSecondary)));
}

// this muon comes from prompt charm decay and does not have light flavor parent
bool isPromptCharmMu(const uint16_t& mask)
bool isPromptCharmMu(const uint16_t mask)
{
return (isMuon(mask) && TESTBIT(mask, HasCharmParent) && (!TESTBIT(mask, HasBeautyParent)) && (!TESTBIT(mask, HasQuarkoniumParent)) && (!TESTBIT(mask, HasLightParent)) && (!TESTBIT(mask, IsSecondary)));
}

// this muon comes from other sources which have not classified above.
bool isOtherMu(const uint16_t& mask)
bool isOtherMu(const uint16_t mask)
{
return (isMuon(mask) && (!isSecondaryMu(mask)) && (!isLightDecayMu(mask)) && (!isTauDecayMu(mask)) && (!isWBosonDecayMu(mask)) && (!isZBosonDecayMu(mask)) && (!isQuarkoniumDecayMu(mask)) && (!isBeautyMu(mask)) && (!isPromptCharmMu(mask)));
}

// this is a hadron
bool isHadron(const uint16_t& mask)
bool isHadron(const uint16_t mask)
{
return (TESTBIT(mask, IsIdentified) && (!TESTBIT(mask, IsMuon)));
}

// this particle is unidentified
bool isUnidentified(const uint16_t& mask)
bool isUnidentified(const uint16_t mask)
{
return ((!TESTBIT(mask, IsIdentified)));
}
Expand Down Expand Up @@ -398,7 +399,10 @@ struct HfTaskSingleMuonMultMc {
}
}

void process(McGenCollisions::iterator const& mccollision, McMuons const& muons, aod::McParticles const&, McRecCollisions const& collisions)
void process(McGenCollisions::iterator const& mccollision,
McMuons const& muons,
aod::McParticles const&,
McRecCollisions const& collisions)
{

// event selections
Expand Down Expand Up @@ -474,7 +478,10 @@ struct HfTaskSingleMuonMultMc {
}
}

void processResTrack(McGenCollisions::iterator const& mccollision, McRecCollisions const& collisions, aod::McParticles const& particles, MyTracks const& tracks)
void processResTrack(McGenCollisions::iterator const& mccollision,
McRecCollisions const& collisions,
aod::McParticles const& particles,
MyTracks const& tracks)
{
// event selections
if (std::abs(mccollision.posZ()) > zVtxMax) {
Expand Down Expand Up @@ -514,17 +521,19 @@ struct HfTaskSingleMuonMultMc {
auto nTrk = 0;
auto tracksample = tracks.sliceBy(perCol, collision.globalIndex());
for (const auto& track : tracksample) {
if (!track.isGlobalTrack())
if (!track.isGlobalTrack()) {
continue;
}
registry.fill(HIST("hParticleRec"), track.pt(), track.eta());
++nTrk;
}
if (nTrk < 1)
if (nTrk < 1) {
continue;
}
registry.fill(HIST("hTrackResponse"), nP, nTrk);
}
}
PROCESS_SWITCH(HfTaskSingleMuonMultMc, processResTrack, "Process Track Reconstruction/Generation", true);
PROCESS_SWITCH(HfTaskSingleMuonMultMc, processResTrack, "Process Track Reconstruction/Generation", false);
};

WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
Expand Down
Loading