Skip to content

Check for the tertiary tracks #7028

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
njacazio merged 4 commits into from
Jul 30, 2024
Merged

Check for the tertiary tracks #7028

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
b655eb8
Update qaEfficiency.cxx
upasanasharma31 Jul 29, 2024
3a1ca1e
Update qaEfficiency.cxx
upasanasharma31 Jul 29, 2024
b262b6d
Update qaEfficiency.cxx
upasanasharma31 Jul 29, 2024
5ee38fa
Update qaEfficiency.cxx
upasanasharma31 Jul 29, 2024
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
23 changes: 5 additions & 18 deletions DPG/Tasks/AOTTrack/qaEfficiency.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -925,24 +925,11 @@ struct QaEfficiency {
}
bool isFinal(const o2::aod::McParticles::iterator& mcParticle)
{
// Example conditions to determine if a particle is final (tertiary)
// Here, we assume that final state particles are those not originating from primary vertex
// and not further decaying into other particles
// Check if the particle has no daughters
if (!mcParticle.has_daughters()) {

// Check if the particle is not a primary particle
if (!mcParticle.isPhysicalPrimary()) {
// Check if the particle is produced in a secondary decay
if (mcParticle.getProcess() == 4) {
// Get the mother particle's index and the mother particle itself
auto mothers = mcParticle.mothers_as<o2::aod::McParticles>();
for (const auto& mother : mothers) {
// Check if the mother particle is not primary and produced in a weak decay
if (!mother.isPhysicalPrimary() && mother.getProcess() == 4) {
return true; // Consider it as a tertiary particle
}
}
if (!mcParticle.has_daughters() && !mcParticle.isPhysicalPrimary() && mcParticle.getProcess() == 4) {
auto mothers = mcParticle.mothers_as<o2::aod::McParticles>();
for (const auto& mother : mothers) {
if (!mother.isPhysicalPrimary() && mother.getProcess() == 4) {
return true;
}
}
}
Expand Down