Skip to content

[PWGHF] replace eta selection with rapidity selection #15334

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

Draft
wuctlby wants to merge 3 commits into
base: master
Choose a base branch
from
+22 −3
Draft

[PWGHF] replace eta selection with rapidity selection #15334

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
f57fb58
replace eta selection with rapidity selection
wuctlby Mar 10, 2026
497726f
Please consider the following formatting changes
alibuild Mar 10, 2026
0a9f31d
Merge pull request #19 from alibuild/alibot-cleanup-15334
wuctlby Mar 10, 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
25 changes: 22 additions & 3 deletions PWGHF/HFC/TableProducer/derivedDataCreatorCorrelationsReduced.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -88,7 +88,6 @@ struct HfDerivedDataCreatorCorrelationsReduced {
Configurable<float> ptCandMax{"ptCandMax", 24., "max. cand. pT"};
Configurable<int> tpcNClsCrossedRowsMin{"tpcNClsCrossedRowsMin", 70, "min. TPC crossed rows for associated tracks"};
Configurable<float> etaTrkMax{"etaTrkMax", 1., "max. track eta"};
Configurable<float> etaCandMax{"etaCandMax", 1., "max. trigger candidate eta"};
Configurable<float> ptTrkMin{"ptTrkMin", 0.2, "min. track pT"};
Configurable<float> ptTrkMax{"ptTrkMax", 3., "max. track pT"};
Configurable<float> dcaXYTrkMax{"dcaXYTrkMax", 1., "max. track DCA XY"};
Expand Down Expand Up @@ -264,6 +263,26 @@ struct HfDerivedDataCreatorCorrelationsReduced {
return outputMl;
}

/// Cut on rapidity of the candidate
/// \param candidate is the charm hadron candidate
template <CandidateType CandType, typename TCand>
bool cutCandRapidity(const TCand& candidate)
{
double y = 0.0;
if constexpr (CandType == CandidateType::DsToKKPi || CandType == CandidateType::DsToPiKK) {
y = HfHelper::yDs(candidate);
} else if constexpr (CandType == CandidateType::DplusToPiKPi) {
y = HfHelper::yDplus(candidate);
} else if constexpr (CandType == CandidateType::D0ToPiK || CandType == CandidateType::D0ToKPi) {
y = HfHelper::yD0(candidate);
} else if constexpr (CandType == CandidateType::LcToPKPi) {
y = HfHelper::yLc(candidate);
} else {
return true;
}
return std::fabs(y) < yCandMax;
}

/// Check event selections for collision and fill the collision table
/// \param collision is the collision
template <typename Coll>
Expand Down Expand Up @@ -336,7 +355,7 @@ struct HfDerivedDataCreatorCorrelationsReduced {
const float collCentrality)
{
for (const auto& trigCand : trigCands) {
if (std::fabs(trigCand.eta()) >= etaCandMax) {
if (!cutCandRapidity<CandType>(trigCand)) {
continue;
}
double trigCandPt = trigCand.pt();
Expand Down Expand Up @@ -392,7 +411,7 @@ struct HfDerivedDataCreatorCorrelationsReduced {
void fillCharmMixedEvent(TTrigCands const& trigCands)
{
for (const auto& trigCand : trigCands) {
if (std::fabs(trigCand.eta()) >= etaCandMax) {
if (!cutCandRapidity<CandType>(trigCand)) {
continue;
}
registry.fill(HIST("hPhiVsPtTrig"), RecoDecay::constrainAngle(trigCand.phi(), -o2::constants::math::PIHalf), trigCand.pt());
Expand Down
Loading