secondaries produced in Pythia interaction are configured in the wrong coordinate system
When converting particles from Pythia::Particle to Corsika particle, it seems that we have used a wrong coordinate system. In this line:
MomentumVector const pyPlab(labCS, {p8p.px() * 1_GeV, p8p.py() * 1_GeV, p8p.pz() * 1_GeV});
I think the right coordinate should be the rotatedCS
which can be obtained from boost.getRotatedCS()
.