Skip to content

[PWGLF] adding only primary in mc for pT dependent cut in DCA #17142

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
maciacco merged 10 commits into from
Jul 21, 2026
+19 −17
Merged

[PWGLF] adding only primary in mc for pT dependent cut in DCA #17142

Changes from all commits
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
e3dd8bd
Modify track settings and DCA parameters for proton
jtanwar2212 Jul 16, 2026
9e9fd82
Set default values for member variables
jtanwar2212 Jul 16, 2026
898812d
Format code for consistency in nucleitpcpbpb.cxx
jtanwar2212 Jul 16, 2026
5315c1b
Add mass²/z² plots for MC
jtanwar2212 Jul 16, 2026
536b78a
Please consider the following formatting changes
alibuild Jul 16, 2026
9923be5
Merge pull request #8 from alibuild/alibot-cleanup-17084
jtanwar2212 Jul 16, 2026
add5c6a
Merge branch 'AliceO2Group:master' into master
jtanwar2212 Jul 20, 2026
83aa59a
Merge branch 'AliceO2Group:master' into master
jtanwar2212 Jul 21, 2026
6f6ce7a
Add Delta Pt histograms for particles and anti-particles
jtanwar2212 Jul 21, 2026
3e8a180
Add decayType check for track histogram filling
jtanwar2212 Jul 21, 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
36 changes: 19 additions & 17 deletions PWGLF/Tasks/Nuspex/nucleitpcpbpb.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.

Check failure on line 1 in PWGLF/Tasks/Nuspex/nucleitpcpbpb.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/workflow-file] 

Name of a workflow file must match the name of the main struct in it (without the PWG prefix). (Class implementation files should be in "Core" directories.)
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
Expand Down Expand Up @@ -918,9 +918,24 @@
histomc.fill(HIST("hSpectramc"), particleAnti, collision.centFT0C(), ptReco, ptTOF);

// Basic track histograms
histos.fill(HIST("dcaXY"), ptReco, track.dcaXY());
histos.fill(HIST("dcaZ"), ptReco, track.dcaZ());
histos.fill(HIST("Tpcsignal"), getRigidity(track) * track.sign(), track.tpcSignal());
if (decayType == 0) {
histos.fill(HIST("dcaXY"), ptReco, track.dcaXY());
histos.fill(HIST("dcaZ"), ptReco, track.dcaZ());
histos.fill(HIST("Tpcsignal"), getRigidity(track) * track.sign(), track.tpcSignal());
}
// Delta Pt histograms
float ptGen = matchedMCParticle.pt();
float deltaPt = ptReco - ptGen;

if (pdg == -particlePdgCodes.at(i) && decayType == 0) { // Anti-particle
histomc.fill(HIST("histDeltaPtVsPtGenanti"), ptReco, deltaPt);
histomc.fill(HIST("histPIDtrackanti"), ptReco, track.pidForTracking());
}

if (pdg == particlePdgCodes.at(i) && decayType == 0) { // Particle
histomc.fill(HIST("histDeltaPtVsPtGen"), ptReco, deltaPt);
histomc.fill(HIST("histPIDtrack"), ptReco, track.pidForTracking());
}

// Fill mass²/z² for MC - separate for particles and anti-particles
if (track.hasTOF()) {
Expand All @@ -940,7 +955,7 @@
}
}

if (!skipHe4) {
if (!skipHe4 && decayType == 0) {
if (pdg > 0) {
histomc.fill(HIST("hMassVsPtMC"), ptReco, massSquareOverChargeSquare, collision.centFT0C());
} else {
Expand All @@ -949,19 +964,6 @@
}
}
}

// Delta Pt histograms
float ptGen = matchedMCParticle.pt();
float deltaPt = ptReco - ptGen;

if (pdg == -particlePdgCodes.at(i)) { // Anti-particle
histomc.fill(HIST("histDeltaPtVsPtGenanti"), ptReco, deltaPt);
histomc.fill(HIST("histPIDtrackanti"), ptReco, track.pidForTracking());
}
if (pdg == particlePdgCodes.at(i)) { // Particle
histomc.fill(HIST("histDeltaPtVsPtGen"), ptReco, deltaPt);
histomc.fill(HIST("histPIDtrack"), ptReco, track.pidForTracking());
}
}
}
break; // Found the matching collision, break out of collision loop
Expand Down Expand Up @@ -1131,7 +1133,7 @@
// ---- Primary particles ----
decayType = 0;
if (matchedMCParticle.has_mothers()) {
for (auto& motherparticle : matchedMCParticle.mothers_as<aod::McParticles>()) {

Check failure on line 1136 in PWGLF/Tasks/Nuspex/nucleitpcpbpb.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop] 

Use constant references for non-modified iterators in range-based for loops.
if (std::find(hfMothCodes.begin(), hfMothCodes.end(), std::abs(motherparticle.pdgCode())) != hfMothCodes.end()) {
decayType = 1;
break;
Expand All @@ -1146,7 +1148,7 @@
decayType = 1;

// Check if it's from HF decay
for (auto& motherparticle : matchedMCParticle.mothers_as<aod::McParticles>()) {

Check failure on line 1151 in PWGLF/Tasks/Nuspex/nucleitpcpbpb.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop] 

Use constant references for non-modified iterators in range-based for loops.
if (std::find(hfMothCodes.begin(), hfMothCodes.end(), std::abs(motherparticle.pdgCode())) != hfMothCodes.end()) {
break;
}
Expand All @@ -1169,7 +1171,7 @@
if (matchedMCParticle.isPhysicalPrimary()) {
histomc.fill(HIST("DCAxy_vs_pT_primary"), ptDCA, track.dcaXY(), collision.centFT0C());
}
if (decayType == 2) { // Transport/Material

Check failure on line 1174 in PWGLF/Tasks/Nuspex/nucleitpcpbpb.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
histomc.fill(HIST("DCAxy_vs_pT_transport"), ptDCA, track.dcaXY(), collision.centFT0C());
} else if (decayType == 1) { // Weak decay (including HF)
histomc.fill(HIST("DCAxy_vs_pT_weakdecay"), ptDCA, track.dcaXY(), collision.centFT0C());
Expand Down
Loading