Wounded nucleon information missing in multiple models #225
Description
Dear expert,
Recently I found out that there are some models that would not generate wounded nucleons successfully in proton-Oxygen collision.
| model name | wounded |
|---|---|
| DPMJET-III-19.1 | Exist |
| DPMJET-III-19.3 | Exist |
| DPMJET-III-3.0-7 | Exist |
| EPOS-LHC | Exist |
| EPOS-LHC-R | Exist |
| QGSJet-01d | Not Exist |
| QGSJet-II-03 | Exist |
| QGSJet-II-04 | Exist |
| QGSJet-III | Exist |
| SIBYLL-2.1 | Not Exist |
| SIBYLL-2.3c | Not Exist |
| SIBYLL-2.3d | Not Exist |
| SIBYLL-2.3dStar-mix | Not Exist |
| SIBYLL-2.3e | Not Exist |
| SIBYLL-2.3eStar-mix | Not Exist |
Let me also attach the script:
from chromo.models import (DpmjetIII191, DpmjetIII193, DpmjetIII307, EposLHC, EposLHCR,
QGSJet01d, QGSJetII03, QGSJetII04, QGSJetIII, Sibyll21,
Sibyll23c, Sibyll23d, Sibyll23dStarMixed, Sibyll23e, Sibyll23eStarMixed,
)
models = [DpmjetIII191, DpmjetIII193, DpmjetIII307, EposLHC, EposLHCR,
QGSJet01d, QGSJetII03, QGSJetII04, QGSJetIII, Sibyll21,
Sibyll23c, Sibyll23d, Sibyll23dStarMixed, Sibyll23e, Sibyll23eStarMixed]
logger_info = []
kin = CenterOfMass(5 * TeV, "p", (16, 8))
for model in models:
try:
gen = model(kin)
gen.set_stable(111)
for event in gen(1):
f = event.final_state()
wounded = f.n_wounded
if wounded[1] != 0:
logger_info.append(f'{gen.label} : Exist')
else:
logger_info.append(f'{gen.label} : Not Exist')
except Exception as e:
logger_info.append(f'Error: {e}')
for info in logger_info:
print(info)Running in chromo: 0.9.0
Kindly confirmed that, thanks!