Skip to content

Skip Chi0c and Chi1c transport by default for EPOS4 #2078

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
jackal1-66 merged 1 commit into from
Aug 2, 2025
Merged

Skip Chi0c and Chi1c transport by default for EPOS4 #2078

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
9 changes: 4 additions & 5 deletions MC/config/examples/epos4/generator_EPOS4.C
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -15,16 +15,15 @@ class GeneratorEPOS4 : public o2::eventgen::GeneratorHepMC
LOG(error) << "Failed to import particles from HepMC event!";
return false;
}
// Remove charmonia chi_0c and chi_1c from the particles list (incompatible with default G4 physics list)
// Skip transport of charmonia chi_0c and chi_1c (incompatible with default G4 physics list)
// These are not decayed by EPOS4 (no daughters)
if (!mEnChi)
{
for (int a = 0; a < mParticles.size(); ++a) {
if (mParticles[a].GetPdgCode() == 10441 || mParticles[a].GetPdgCode() == 20443)
{
LOG(debug) << "Removing charmonium state " << mParticles[a].GetPdgCode() << " from particles list";
mParticles.erase(mParticles.begin() + a);
--a; // Adjust index after erasing
LOG(debug) << "Setting charmonium state " << mParticles[a].GetPdgCode() << " with index " << a << " not to be transported";
mParticles[a].SetBit(ParticleStatus::kToBeDone, false);
}
}
}
Expand All @@ -34,7 +33,7 @@ class GeneratorEPOS4 : public o2::eventgen::GeneratorHepMC
void setChiFlag(bool &flag) {
mEnChi = flag;
if (!flag) {
LOG(info) << "Charmonium states chi_0c and chi_1c will be removed from the particles list";
LOG(info) << "Charmonium states chi_0c and chi_1c will not be transported.";
}
}

Expand Down