From 9bcb3367122fc223e909dd4d8f56730d478a5350 Mon Sep 17 00:00:00 2001
From: ralfulrich <ralf.ulrich@kit.edu>
Date: Wed, 13 Feb 2019 15:58:38 +0100
Subject: [PATCH] fixed SIBYLL: cannot Delete particle before calling
AddSecondaries....
---
Processes/Sibyll/Decay.cc | 4 ++--
Processes/Sibyll/Interaction.cc | 5 ++---
2 files changed, 4 insertions(+), 5 deletions(-)
diff --git a/Processes/Sibyll/Decay.cc b/Processes/Sibyll/Decay.cc
index 1e0d1ca2..eec1f27f 100644
--- a/Processes/Sibyll/Decay.cc
+++ b/Processes/Sibyll/Decay.cc
@@ -157,8 +157,6 @@ namespace corsika::process::sibyll {
// remember position
Point const decayPoint = p.GetPosition();
TimeType const t0 = p.GetTime();
- // remove original particle from corsika stack
- p.Delete();
// set all particles/hadrons unstable
// setHadronsUnstable();
setUnstable(pCode);
@@ -184,6 +182,8 @@ namespace corsika::process::sibyll {
}
// empty sibyll stack
ss.Clear();
+ // remove original particle from corsika stack
+ p.Delete();
}
} // namespace corsika::process::sibyll
diff --git a/Processes/Sibyll/Interaction.cc b/Processes/Sibyll/Interaction.cc
index 5ccc1b74..4c5f6044 100644
--- a/Processes/Sibyll/Interaction.cc
+++ b/Processes/Sibyll/Interaction.cc
@@ -320,9 +320,6 @@ namespace corsika::process::sibyll {
sib_list_(print_unit);
fNucCount += get_nwounded() - 1;
- // delete current particle
- p.Delete();
-
// add particles from sibyll to stack
// link to sibyll stack
SibStack ss;
@@ -356,6 +353,8 @@ namespace corsika::process::sibyll {
<< ", Plab_final=" << (Plab_final / 1_GeV).GetComponents() << endl;
}
}
+ // delete current particle
+ p.Delete();
return process::EProcessReturn::eOk;
}
--
GitLab