Skip to content
Closed
Changes from all commits
Commits
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
9 changes: 4 additions & 5 deletions PWGHF/HFL/Tasks/taskSingleMuonMult.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -72,7 +72,7 @@ struct HfTaskSingleMuonMult {
Filter globalTrackFilter = ((o2::aod::track::isGlobalTrack == true) && (nabs(o2::aod::track::eta) < etaTrackMax) && ((o2::aod::track::pt) > ptTrackMin));

HistogramRegistry registry{"registry"};

void init(InitContext&)
{
AxisSpec axisCent = {101, -0.5, 100.5, "centrality"};
Expand Down Expand Up @@ -116,7 +116,7 @@ struct HfTaskSingleMuonMult {

registry.add("hTHnTrk", "Muon information with multiplicity", {HistType::kTHnSparseF, {axisCent, axisNCh, axisPt, axisEta, axisSign}, 5});
registry.add("h3MultNchNmu", "Number of muons and multiplicity", {HistType::kTH3F, {axisCent, axisNCh, axisNMu}});
registry.add("hMultNchNmuTrackType", "Number of muons with different types and multiplicity", {HistType::kTHnSparseF, {axisCent, axisNCh, axisNMu, axisTrackType},4});
registry.add("hMultNchNmuTrackType", "Number of muons with different types and multiplicity", {HistType::kTHnSparseF, {axisCent, axisNCh, axisNMu, axisTrackType}, 4});

auto hEvstat = registry.get<TH1>(HIST("hEvent"));
auto* xEv = hEvstat->GetXaxis();
Expand All @@ -131,10 +131,9 @@ struct HfTaskSingleMuonMult {
xMu->SetBinLabel(3, "RAbsorbCut");
xMu->SetBinLabel(4, "pDcaCut");
xMu->SetBinLabel(5, "chi2Cut");

// Number the types of muon tracks
constexpr uint8_t nTrackTypes{static_cast<uint8_t>(ForwardTrackTypeEnum::MCHStandaloneTrack)};

}

void process(MyCollisions::iterator const& collision,
Expand Down Expand Up @@ -172,7 +171,7 @@ struct HfTaskSingleMuonMult {

// muons per event
int nMu{0};
int nMuType[nTrackTypes+1] = {0};
int nMuType[nTrackTypes + 1] = {0};

for (const auto& muon : muons) {
const auto pt{muon.pt()}, eta{muon.eta()}, theta{90.0f - ((std::atan(muon.tgl())) * constants::math::Rad2Deg)}, pDca{muon.pDca()}, rAbsorb{muon.rAtAbsorberEnd()}, chi2{muon.chi2MatchMCHMFT()};
Expand Down
Loading