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!
Dear expert,
Recently I found out that there are some models that would not generate wounded nucleons successfully in proton-Oxygen collision.
Let me also attach the script:
Running in chromo: 0.9.0
Kindly confirmed that, thanks!