IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 8d91d2bf authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan: Committed by Felix Riehn
Browse files

fixed coordinate systems in Sibyll::Interaction

parent 69aa3dd5
No related branches found
No related tags found
2 merge requests!234WIP: Initial example of python as script language from C++,!204Resolve "boost & coordinate system in process::sibyll::Interaction"
...@@ -315,12 +315,13 @@ namespace corsika::process::sibyll { ...@@ -315,12 +315,13 @@ namespace corsika::process::sibyll {
auto const p3lab = Plab.GetSpaceLikeComponents(); auto const p3lab = Plab.GetSpaceLikeComponents();
assert(p3lab.GetCoordinateSystem() == originalCS); // just to be sure! assert(p3lab.GetCoordinateSystem() == originalCS); // just to be sure!
auto const pid = process::sibyll::ConvertFromSibyll(psib.GetPID());
// add to corsika stack // add to corsika stack
auto pnew = vP.AddSecondary( auto pnew = vP.AddSecondary(
tuple<particles::Code, units::si::HEPEnergyType, stack::MomentumVector, tuple<particles::Code, units::si::HEPEnergyType, stack::MomentumVector,
geometry::Point, units::si::TimeType>{ geometry::Point, units::si::TimeType>{
process::sibyll::ConvertFromSibyll(psib.GetPID()), pid, Plab.GetTimeLikeComponent(), p3lab, pOrig, tOrig});
Plab.GetTimeLikeComponent(), p3lab, pOrig, tOrig});
Plab_final += pnew.GetMomentum(); Plab_final += pnew.GetMomentum();
Elab_final += pnew.GetEnergy(); Elab_final += pnew.GetEnergy();
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment