Skip to content

[PWGLF] Added ITS-standalone tracks in dNdEta study with phi #11291

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
romainschotter merged 8 commits into from
May 21, 2025
Merged

[PWGLF] Added ITS-standalone tracks in dNdEta study with phi #11291

Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
cf01149
Added mult histos in dNdEta process functions
scannito May 20, 2025
2c3ece9
Add ITS-standalone tracks
scannito May 21, 2025
32dd8c9
Add ITS-standalone tracks 2
scannito May 21, 2025
5beb338
Changed track filter
scannito May 21, 2025
f74cb19
Fix clang format
scannito May 21, 2025
ca23b9d
Changed track filter 2
scannito May 21, 2025
a21ec32
Fix clang format 2
scannito May 21, 2025
733769d
Fix namespace issue
scannito May 21, 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
17 changes: 15 additions & 2 deletions PWGLF/Tasks/Strangeness/phik0shortanalysis.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -58,6 +58,7 @@
using namespace o2;
using namespace o2::framework;
using namespace o2::framework::expressions;
using namespace o2::aod::track;

struct Phik0shortanalysis {
// Histograms are defined with HistogramRegistry
Expand Down Expand Up @@ -186,6 +187,11 @@
double massK0S = o2::constants::physics::MassK0Short;
double massLambda = o2::constants::physics::MassLambda0;

// Defining track flags
static constexpr TrackSelectionFlags::flagtype TrackSelectionITS = TrackSelectionFlags::kITSNCls | TrackSelectionFlags::kITSChi2NDF | TrackSelectionFlags::kITSHits;
static constexpr TrackSelectionFlags::flagtype TrackSelectionTPC = TrackSelectionFlags::kTPCNCls | TrackSelectionFlags::kTPCCrossedRowsOverNCls | TrackSelectionFlags::kTPCChi2NDF;
static constexpr TrackSelectionFlags::flagtype TrackSelectionDCA = TrackSelectionFlags::kDCAz | TrackSelectionFlags::kDCAxy;

// Defining filters for events (event selection)
// Processed events will be already fulfilling the event selection requirements
Filter eventFilter = (o2::aod::evsel::sel8 == true);
Expand All @@ -194,8 +200,12 @@
// Defining filters on V0s (cannot filter on dynamic columns)
Filter preFilterV0 = (nabs(aod::v0data::dcapostopv) > v0Configs.v0SettingDCAPosToPV && nabs(aod::v0data::dcanegtopv) > v0Configs.v0SettingDCANegToPV && aod::v0data::dcaV0daughters < v0Configs.v0SettingDCAV0Dau);

// Defining filters on tracks (cannot filter on dynamic columns)
Filter trackFilter = requireGlobalTrackWoDCAInFilter();
// Defining filters on tracks
Filter trackFilter = ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::ITS) &&
ncheckbit(aod::track::trackCutFlag, TrackSelectionITS) &&
ifnode(ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::TPC), ncheckbit(aod::track::trackCutFlag, TrackSelectionTPC), true) &&
ncheckbit(aod::track::trackCutFlag, TrackSelectionDCA) &&
ncheckbit(aod::track::trackCutFlag, TrackSelectionFlags::kInAcceptanceTracks);

// Defining the type of the collisions for data and MC
using SelCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::PVMults>;
Expand Down Expand Up @@ -1540,7 +1550,7 @@
if (mcTrack.isPhysicalPrimary()) {
mcPionHist.fill(HIST("h3RecMCDCAxyPrimPi"), track.pt(), track.dcaXY());
} else {
if (mcTrack.getProcess() == 4) { // Selection of secondary pions from weak decay

Check failure on line 1553 in PWGLF/Tasks/Strangeness/phik0shortanalysis.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.
mcPionHist.fill(HIST("h3RecMCDCAxySecWeakDecayPi"), track.pt(), track.dcaXY());
} else { // Selection of secondary pions from material interactions
mcPionHist.fill(HIST("h3RecMCDCAxySecMaterialPi"), track.pt(), track.dcaXY());
Expand Down Expand Up @@ -2013,7 +2023,7 @@
if (mcParticle1.pdgCode() != o2::constants::physics::Pdg::kPhi)
continue;
auto kDaughters = mcParticle1.daughters_as<aod::McParticles>();
if (kDaughters.size() != 2)

Check failure on line 2026 in PWGLF/Tasks/Strangeness/phik0shortanalysis.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.
continue;
bool isPosKaon = false, isNegKaon = false;
for (const auto& kDaughter : kDaughters) {
Expand Down Expand Up @@ -2150,7 +2160,7 @@
continue;
if (cfgisGenMCForClosure) {
auto kDaughters = mcParticle2.daughters_as<aod::McParticles>();
if (kDaughters.size() != 2)

Check failure on line 2163 in PWGLF/Tasks/Strangeness/phik0shortanalysis.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.
continue;
bool isPosKaon = false, isNegKaon = false;
for (const auto& kDaughter : kDaughters) {
Expand Down Expand Up @@ -2227,7 +2237,7 @@
continue;
if (cfgisGenMCForClosure) {
auto kDaughters = mcParticle2.daughters_as<aod::McParticles>();
if (kDaughters.size() != 2)

Check failure on line 2240 in PWGLF/Tasks/Strangeness/phik0shortanalysis.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.
continue;
bool isPosKaon = false, isNegKaon = false;
for (const auto& kDaughter : kDaughters) {
Expand Down Expand Up @@ -2278,6 +2288,7 @@
return;

float multiplicity = collision.centFT0M();
dataEventHist.fill(HIST("hMultiplicityPercent"), multiplicity);

for (const auto& track : filteredTracks)
dataEventHist.fill(HIST("h2EtaDistribution"), multiplicity, track.eta());
Expand All @@ -2300,6 +2311,7 @@
continue;

float genmultiplicity = mcCollision.centFT0M();
mcEventHist.fill(HIST("hRecMCGenMultiplicityPercent"), genmultiplicity);

auto mcTracksThisColl = filteredMCTracks.sliceByCached(aod::track::collisionId, collision.globalIndex(), cache);

Expand Down Expand Up @@ -2336,6 +2348,7 @@
}

float genmultiplicity = mcCollision.centFT0M();
mcEventHist.fill(HIST("hGenMCMultiplicityPercent"), genmultiplicity);

for (const auto& mcParticle : mcParticles) {
if (!mcParticle.isPhysicalPrimary() || std::abs(mcParticle.eta()) > trackConfigs.etaMax)
Expand Down
Loading