Ions in final_state()

[michael-pitt]

Hello,

I'm looking for a ion remnants from OO collisions, but I cannot find a stable particles with ion pid.

import chromo
from chromo.kinematics import CenterOfMass
from chromo.constants import GeV, TeV
import chromo.models as im

ekin = CenterOfMass(2*8*6.8*TeV, "O", "O")

models = [
    im.EposLHC(ekin),
    im.QGSJetII04(ekin),
]

for model in models:
    print(model)
    for event in model(10):
        print('number of stable particles :',len(event.final_state().pid))
        print('list of ions',[pid for pid in event.final_state().pid if pid > 99999])

How can I find the ion remnants?

Thank you in advance, Michael

[afedynitch]

Hi Michael,

most models won’t produce realistic remnant since these high energy models only do Glauber and ignore other nuclear effects. So the remnants won’t be realistic.

We could consistently provide whatever the models return though but it’s not implemented at this moment. Please leave this issue and we will follow this up asap.

For now maybe Pythia8/Argamtyr could return remnants (but also not realistic).

[michael-pitt]

Hi @afedynitch,

Thanks for the quick feedback, I was interested in DPMJet-III output. I see that some particles have pid of 99999, maybe those are the lost ions?

For pythia, ions are not allowed as input particles

Pythia8.projectiles
{<PDGID: 2112>, <PDGID: 321>, <PDGID: 130>, <PDGID: 211>, <PDGID: 2212>, <PDGID: 310>, <PDGID: 22>}

I guess this will require integrating the Angantyr settings.

regards, Michael