Skip to content

[PWGHF] Add hadron and electron PID cuts in HFE–hadron correlation #14767

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

Open
rashigupt wants to merge 17 commits into
base: master
Choose a base branch
from
+270 −9
Open

[PWGHF] Add hadron and electron PID cuts in HFE–hadron correlation #14767

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
17 commits
Select commit Hold shift + click to select a range
127ef4c
Add 3D histogram before electron Selection
rashigupt Nov 11, 2025
734c26e
remove typo error
rashigupt Nov 11, 2025
e830592
use Capitalise properly
rashigupt Nov 11, 2025
c03d12c
change TPC nSigma Info to TPC nSigma
rashigupt Nov 13, 2025
683556a
Add LSmass and ULmass seprate column
rashigupt Dec 1, 2025
250f093
Add two separate columns, LSmass and ULmass.
rashigupt Dec 1, 2025
2a75245
We use the mass column correctly by separating it into LSmass and ULm…
rashigupt Dec 1, 2025
e1579ec
Merge branch 'master' into rashi
rashigupt Dec 1, 2025
398771c
remove clang error
rashigupt Dec 1, 2025
0cc3808
remove clang error
rashigupt Dec 1, 2025
8cf1be2
Remove formate error
rashigupt Dec 2, 2025
b826c36
Update CorrelationTables.h
rashigupt Dec 3, 2025
602d04c
Fix formatting issues to pass MegaLinter
rashi234 Dec 4, 2025
d060b20
Merge branch 'AliceO2Group:master' into rashi
rashigupt Dec 4, 2025
fabfa0d
Merge branch 'AliceO2Group:master' into rashi
rashigupt Feb 3, 2026
7e69ec8
Update electronSelectionWithTpcEmcal.cxx
rashigupt Feb 3, 2026
2758ffb
Update correlatorHfeHadrons.cxx
rashigupt Feb 3, 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
82 changes: 74 additions & 8 deletions PWGHF/HFC/TableProducer/correlatorHfeHadrons.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -138,8 +138,21 @@ struct HfCorrelatorHfeHadrons {
registry.add("hMCgenInclusiveEHCorrl", "Sparse for Delta phi and Delta eta for McGen Electron pair with Hadron;p_{T}^{e} (GeV#it{/c});p_{T}^{h} (GeV#it{/c});#Delta#varphi;#Delta#eta;", {HistType::kTHnSparseF, {{axisPt}, {axisPt}, {axisDeltaPhi}, {axisDeltaEta}}});
registry.add("hptElectron", "hptElectron", {HistType::kTH1D, {axisPt}});
registry.add("hptHadron", "hptHadron", {HistType::kTH1D, {axisPt}});
registry.add("hptHadronMcRec", "hptHadronMcRec", {HistType::kTH1D, {axisPt}});
registry.add("hptHadronMcRecPrimary", "hptHadronMcRecPrimary", {HistType::kTH1D, {axisPt}});
registry.add("hMCRecptEleHadron", "hMCRecptEleHadron", {HistType::kTH1D, {axisPt}});
registry.add("hMCRecptMuonHadron", "hMCRecptMuonHadron", {HistType::kTH1D, {axisPt}});
registry.add("hMCRecptPionHadron", "hMCRecptPionHadron", {HistType::kTH1D, {axisPt}});
registry.add("hMCRecptKaonHadron", "hMCRecptKaonHadron", {HistType::kTH1D, {axisPt}});
registry.add("hMCRecptProtonHadron", "hMCRecptProtonHadron", {HistType::kTH1D, {axisPt}});

registry.add("hMCgenptHadron", "hMCgenptHadron", {HistType::kTH1D, {axisPt}});
registry.add("hMCgenptHadronprimary", "hMCgenptHadronprimary", {HistType::kTH1D, {axisPt}});
registry.add("hMCgenptEleHadron", "hMCgenptEleHadron", {HistType::kTH1D, {axisPt}});
registry.add("hMCgenptMuonHadron", "hMCgenptMuonHadron", {HistType::kTH1D, {axisPt}});
registry.add("hMCgenptPionHadron", "hMCgenptPionHadron", {HistType::kTH1D, {axisPt}});
registry.add("hMCgenptKaonHadron", "hMCgenptKaonHadron", {HistType::kTH1D, {axisPt}});
registry.add("hMCgenptProtonHadron", "hMCgenptProtonHadron", {HistType::kTH1D, {axisPt}});

registry.add("hMixEventInclusiveEHCorrl", "Sparse for mix event Delta phi and Delta eta Inclusive Electron with Hadron;p_{T}^{e} (GeV#it{/c});p_{T}^{h} (GeV#it{/c});#Delta#varphi;#Delta#eta;", {HistType::kTHnSparseF, {{axisPt}, {axisPt}, {axisDeltaPhi}, {axisDeltaEta}}});
registry.add("hMixEventLSEHCorrel", "Sparse for mix event Delta phi and Delta eta Like sign Electron pair with Hadron;p_{T}^{e} (GeV#it{/c});p_{T}^{h} (GeV#it{/c});#Delta#varphi;#Delta#eta;", {HistType::kTHnSparseF, {{axisPt}, {axisPt}, {axisDeltaPhi}, {axisDeltaEta}}});
Expand Down Expand Up @@ -173,8 +186,8 @@ struct HfCorrelatorHfeHadrons {
}

// Electron-hadron Correlation
template <typename TracksType, typename ElectronType, typename CollisionType, typename BcType>
void fillCorrelation(CollisionType const& collision, ElectronType const& electrons, TracksType const& tracks, BcType const&)
template <bool isMc, typename TracksType, typename ElectronType, typename CollisionType, typename BcType, typename McParticles>
void fillCorrelation(CollisionType const& collision, ElectronType const& electrons, TracksType const& tracks, BcType const&, McParticles const&)
{
if (!(isRun3 ? collision.sel8() : (collision.sel7() && collision.alias_bit(kINT7)))) {
return;
Expand All @@ -192,15 +205,50 @@ struct HfCorrelatorHfeHadrons {
}

registry.fill(HIST("hNevents"), 1);

// Add hadron Table For Mix Event Electron Hadron correlation
if (!skipEventTableFilling) {
registry.fill(HIST("hZvertex"), collision.posZ());
for (const auto& hTrack : tracks) {
if (!selAssoHadron(hTrack)) {
continue;
}

// Mc rec hadron efficiency
if constexpr (isMc) {
if (hTrack.has_mcParticle()) {
auto mcParticle = hTrack.template mcParticle_as<aod::McParticles>();
if ((std::abs(mcParticle.pdgCode()) != kElectron) && (std::abs(mcParticle.pdgCode()) != kMuonMinus) && (std::abs(mcParticle.pdgCode()) != kPiPlus) && (std::abs(mcParticle.pdgCode()) != kKPlus) && (std::abs(mcParticle.pdgCode()) != kProton)) {
continue;
}

registry.fill(HIST("hptHadronMcRec"), hTrack.pt());
if (mcParticle.isPhysicalPrimary()) {

registry.fill(HIST("hptHadronMcRecPrimary"), hTrack.pt());

if (std::abs(mcParticle.pdgCode()) != kElectron) {
registry.fill(HIST("hMCRecptEleHadron"), hTrack.pt());
}
if (std::abs(mcParticle.pdgCode()) != kMuonMinus) {
registry.fill(HIST("hMCRecptMuonHadron"), hTrack.pt());
}
if (std::abs(mcParticle.pdgCode()) != kPiPlus) {
registry.fill(HIST("hMCRecptPionHadron"), hTrack.pt());
}
if (std::abs(mcParticle.pdgCode()) != kKPlus) {
registry.fill(HIST("hMCRecptKaonHadron"), hTrack.pt());
}
if (std::abs(mcParticle.pdgCode()) != kProton) {
registry.fill(HIST("hMCRecptProtonHadron"), hTrack.pt());
}
}
}
}
registry.fill(HIST("hTracksBin"), poolBin);

registry.fill(HIST("hptHadron"), hTrack.pt());

entryHadron(hTrack.phi(), hTrack.eta(), hTrack.pt(), poolBin, gCollisionId, timeStamp);
}
}
Expand Down Expand Up @@ -381,22 +429,22 @@ struct HfCorrelatorHfeHadrons {
TableTracks const& tracks,
aod::BCsWithTimestamps const& bcs)
{
fillCorrelation(collision, electrons, tracks, bcs);
fillCorrelation<false>(collision, electrons, tracks, bcs, 0);
}

PROCESS_SWITCH(HfCorrelatorHfeHadrons, processData, "Process for Data", true);
PROCESS_SWITCH(HfCorrelatorHfeHadrons, processData, "Process for Data", false);

// ======= Process starts for McRec, Same event ============

void processMcRec(McTableCollision const& mcCollision,
aod::HfCorrSelEl const& mcElectrons,
McTableTracks const& mcTracks,
aod::BCsWithTimestamps const& bcs)
aod::BCsWithTimestamps const& bcs, aod::McParticles const& mcParticle)
{
fillCorrelation(mcCollision, mcElectrons, mcTracks, bcs);
fillCorrelation<true>(mcCollision, mcElectrons, mcTracks, bcs, mcParticle);
}

PROCESS_SWITCH(HfCorrelatorHfeHadrons, processMcRec, "Process MC Reco mode", false);
PROCESS_SWITCH(HfCorrelatorHfeHadrons, processMcRec, "Process MC Reco mode", true);

void processMcGen(McGenTableCollision const& mcCollision, aod::McParticles const& mcParticles, aod::HfMcGenSelEl const& electrons)
{
Expand All @@ -411,11 +459,29 @@ struct HfCorrelatorHfeHadrons {
if (particleMc.pt() < ptTrackMin) {
continue;
}

if ((std::abs(particleMc.pdgCode()) != kElectron) && (std::abs(particleMc.pdgCode()) != kMuonMinus) && (std::abs(particleMc.pdgCode()) != kPiPlus) && (std::abs(particleMc.pdgCode()) != kKPlus) && (std::abs(particleMc.pdgCode()) != kProton)) {
continue;
}
registry.fill(HIST("hMCgenptHadron"), particleMc.pt());
if (particleMc.isPhysicalPrimary()) {

registry.fill(HIST("hMCgenptHadronprimary"), particleMc.pt());

if (std::abs(particleMc.pdgCode()) != kElectron) {
registry.fill(HIST("hMCgenptEleHadron"), particleMc.pt());
}
if (std::abs(particleMc.pdgCode()) != kMuonMinus) {
registry.fill(HIST("hMCgenptMuonHadron"), particleMc.pt());
}
if (std::abs(particleMc.pdgCode()) != kPiPlus) {
registry.fill(HIST("hMCgenptPionHadron"), particleMc.pt());
}
if (std::abs(particleMc.pdgCode()) != kKPlus) {
registry.fill(HIST("hMCgenptKaonHadron"), particleMc.pt());
}
if (std::abs(particleMc.pdgCode()) != kProton) {
registry.fill(HIST("hMCgenptProtonHadron"), particleMc.pt());
}
}
}

Expand Down
Loading
Loading