Skip to content

[PWGHF] Add centrality information to OmegaC tree creator #13423

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
fcatalan92 merged 4 commits into from
Oct 20, 2025
Merged

[PWGHF] Add centrality information to OmegaC tree creator #13423

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
0a18212
Add centrality column to trees + cleanup keepOnlyMcSignal logic
fcatalan92 Oct 15, 2025
ebabcdb
Add process functions with centrality (FT0M) information
fcatalan92 Oct 15, 2025
b9e2e14
Address Vit's comments
fcatalan92 Oct 17, 2025
b4cd776
Remove not needed function from task
fcatalan92 Oct 17, 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
15 changes: 3 additions & 12 deletions PWGHF/D2H/Tasks/taskOmegac0ToOmegapi.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -168,15 +168,6 @@ struct HfTaskOmegac0ToOmegapi {
registry.get<THnSparse>(HIST("hReco"))->Sumw2();
}

/// Evaluate centrality/multiplicity percentile (centrality estimator is automatically selected based on the used table)
/// \param candidate is candidate
/// \return centrality/multiplicity percentile of the collision
template <typename Coll>
float evaluateCentralityColl(const Coll& collision)
{
return o2::hf_centrality::getCentralityColl<Coll>(collision);
}

template <bool ApplyMl, typename CandType>
void processData(const CandType& candidates)
{
Expand Down Expand Up @@ -214,7 +205,7 @@ struct HfTaskOmegac0ToOmegapi {
continue;
}

float const cent = evaluateCentralityColl(collision);
float const cent = o2::hf_centrality::getCentralityColl(collision);

if constexpr (ApplyMl) {
registry.fill(HIST("hReco"), candidate.invMassCharmBaryon(), candidate.ptCharmBaryon(), candidate.kfRapOmegac(),
Expand Down Expand Up @@ -283,7 +274,7 @@ struct HfTaskOmegac0ToOmegapi {

auto collision = candidate.template collision_as<CollisionsWithMcLabels>();
uint16_t const numPvContributors = collision.numContrib();
float const mcCent = evaluateCentralityColl(collision.template mcCollision_as<McCollisionWithCents>());
float const mcCent = o2::hf_centrality::getCentralityColl(collision.template mcCollision_as<McCollisionWithCents>());

if constexpr (ApplyMl) {
registry.fill(HIST("hReco"), candidate.invMassCharmBaryon(), candidate.ptCharmBaryon(), candidate.kfRapOmegac(), mcCent, numPvContributors, candidate.ptBhadMotherPart(), candidate.originMcRec(), candidate.flagMcMatchRec(), candidate.mlProbOmegac()[0]);
Expand All @@ -309,7 +300,7 @@ struct HfTaskOmegac0ToOmegapi {
maxNumContrib = recCol.numContrib() > maxNumContrib ? recCol.numContrib() : maxNumContrib;
}

float const mcCent = evaluateCentralityColl(mcCollision);
float const mcCent = o2::hf_centrality::getCentralityColl(mcCollision);

if (particle.originMcGen() == RecoDecay::OriginType::Prompt) {
registry.fill(HIST("hMcGen"), ptGen, -1., yGen, RecoDecay::OriginType::Prompt, mcCent, maxNumContrib);
Expand Down
Loading
Loading