diff --git a/Processes/Sibyll/Interaction.cc b/Processes/Sibyll/Interaction.cc index e508482d5039c808b5f19c9443684dd729fd9c02..00ac5976813da70061aa9feca2f42d7df84651b8 100644 --- a/Processes/Sibyll/Interaction.cc +++ b/Processes/Sibyll/Interaction.cc @@ -312,13 +312,15 @@ namespace corsika::process::sibyll { auto const pCoM = Vector<hepmomentum_d>(csPrime, tmp); HEPEnergyType const eCoM = psib.GetEnergy(); auto const Plab = boost.fromCoM(FourVector(eCoM, pCoM)); + auto const p3lab = Plab.GetSpaceLikeComponents(); + assert(p3lab.GetCoordinateSystem() == originalCS); // just to be sure! // add to corsika stack auto pnew = vP.AddSecondary( tuple<particles::Code, units::si::HEPEnergyType, stack::MomentumVector, geometry::Point, units::si::TimeType>{ - pid, Plab.GetTimeLikeComponent(), Plab.GetSpaceLikeComponents(), pOrig, - tOrig}); + process::sibyll::ConvertFromSibyll(psib.GetPID()), + Plab.GetTimeLikeComponent(), p3lab, pOrig, tOrig}); Plab_final += pnew.GetMomentum(); Elab_final += pnew.GetEnergy();