Skip to content

[PWGLF] add isPhysiscalPrimary option #11432

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 1 commit into from
Jun 2, 2025
+3 −0
Merged

[PWGLF] add isPhysiscalPrimary option #11432

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
3 changes: 3 additions & 0 deletions PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -80,6 +80,7 @@
Configurable<int> nPtBins{"nPtBins", 200, "number of pT bins"};
Configurable<float> rapidityCut{"rapidityCut", 0.5, "max (absolute) rapidity of generated cascade"};
Configurable<int> nRapidityBins{"nRapidityBins", 200, "number of pT bins"};
Configurable<bool> requirePhysicalPrimary{"requirePhysicalPrimary", false, "require the generated cascade to be a physical primary"};

void init(InitContext const&)
{
Expand All @@ -106,8 +107,10 @@
// Count all generated MC particles
// WARNING: event-level losses have to be understood too
for (auto& particle : mcparts) {
if (TMath::Abs(particle.y()) > rapidityCut)

Check failure on line 110 in PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity] 

Replace ROOT entities with equivalents from standard C++ or from O2.
continue;
if (requirePhysicalPrimary && !particle.isPhysicalPrimary())
continue;
if (particle.pdgCode() == 3312) {
registry.fill(HIST("hPtXiMinus"), particle.pt());
registry.fill(HIST("h2DXiMinus"), particle.pt(), particle.y());
Expand Down Expand Up @@ -157,12 +160,12 @@
void process(aod::Collision const& collision, LabeledCascades const& Cascades, aod::McParticles const&)
{
for (auto& casc : Cascades) {
Int_t lPDG = 0;

Check failure on line 163 in PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity] 

Replace ROOT entities with equivalents from standard C++ or from O2.
if (assocMC) {
if (!casc.has_mcParticle())
continue;
auto cascmc = casc.mcParticle();
if (TMath::Abs(cascmc.pdgCode()) == 3312 || TMath::Abs(cascmc.pdgCode()) == 3334)

Check failure on line 168 in PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity] 

Replace ROOT entities with equivalents from standard C++ or from O2.
lPDG = cascmc.pdgCode();
}

Expand Down Expand Up @@ -269,21 +272,21 @@
auto negTrack = lCascade.template negTrack_as<TCascTracksTo>();

// Bachelor: depends on type
if (TMath::Abs(bachTrack.tpcNSigmaPi()) > tpcNsigmaBachelor && tpcNsigmaBachelor < 9.99)

Check failure on line 275 in PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity] 

Replace ROOT entities with equivalents from standard C++ or from O2.
lConsistentWithXi = false;
if (TMath::Abs(bachTrack.tpcNSigmaKa()) > tpcNsigmaBachelor && tpcNsigmaBachelor < 9.99)

Check failure on line 277 in PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity] 

Replace ROOT entities with equivalents from standard C++ or from O2.
lConsistentWithOm = false;

// Proton check: depends on cascade sign
if (lCascade.sign() < 0 && TMath::Abs(posTrack.tpcNSigmaPr()) > tpcNsigmaProton && tpcNsigmaProton < 9.99)

Check failure on line 281 in PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity] 

Replace ROOT entities with equivalents from standard C++ or from O2.
lConsistentWithLambda = false;
if (lCascade.sign() > 0 && TMath::Abs(negTrack.tpcNSigmaPr()) > tpcNsigmaProton && tpcNsigmaProton < 9.99)

Check failure on line 283 in PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity] 

Replace ROOT entities with equivalents from standard C++ or from O2.
lConsistentWithLambda = false;

// Pion check: depends on cascade sign
if (lCascade.sign() < 0 && TMath::Abs(negTrack.tpcNSigmaPi()) > tpcNsigmaPion && tpcNsigmaPion < 9.99)

Check failure on line 287 in PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity] 

Replace ROOT entities with equivalents from standard C++ or from O2.
lConsistentWithLambda = false;
if (lCascade.sign() > 0 && TMath::Abs(posTrack.tpcNSigmaPi()) > tpcNsigmaPion && tpcNsigmaPion < 9.99)

Check failure on line 289 in PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity] 

Replace ROOT entities with equivalents from standard C++ or from O2.
lConsistentWithLambda = false;

// bit-packing (first bit -> consistent with Xi, second bit -> consistent with Omega)
Expand All @@ -297,7 +300,7 @@
registry.fill(HIST("hCandidateCounter"), 0.5); // all candidates

// check mc association if requested
Int_t lPDG = 0;

Check failure on line 303 in PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity] 

Replace ROOT entities with equivalents from standard C++ or from O2.
if (assocMC) {
if (!casc.has_mcParticle())
return;
Expand Down
Loading