From 5ddaccd3dfb3dee0f32a5c6bb0e3d132bd328dd5 Mon Sep 17 00:00:00 2001
From: Maximilian Sackel <maximilian.sackel@tu-dortmund.de>
Date: Tue, 23 Jun 2020 10:00:50 +0000
Subject: [PATCH] try to modify cascade example so that proposal is used.
---
Documentation/Examples/cascade_example.cc | 12 +++++++-----
1 file changed, 7 insertions(+), 5 deletions(-)
diff --git a/Documentation/Examples/cascade_example.cc b/Documentation/Examples/cascade_example.cc
index 8ff257cc5..7aa0561c0 100644
--- a/Documentation/Examples/cascade_example.cc
+++ b/Documentation/Examples/cascade_example.cc
@@ -137,6 +137,7 @@ int main() {
random::RNGManager::GetInstance().RegisterRandomStream("sibyll");
random::RNGManager::GetInstance().RegisterRandomStream("pythia");
+ random::RNGManager::GetInstance().RegisterRandomStream("proposal");
process::sibyll::Interaction sibyll;
process::sibyll::NuclearInteraction sibyllNuc(sibyll, env);
process::sibyll::Decay decay;
@@ -148,21 +149,22 @@ int main() {
process::energy_loss::EnergyLoss eLoss{showerAxis};
// assemble all processes into an ordered process list
- auto sequence = stackInspect << sibyll << sibyllNuc << decay << eLoss << cut
- << trackWriter;
+ auto sequence = stackInspect << sibyll << sibyllNuc << proposal << decay
+ /* << eLoss */
+ << cut << trackWriter;
// define air shower object, run simulation
cascade::Cascade EAS(env, tracking, sequence, stack);
EAS.Run();
- eLoss.PrintProfile(); // print longitudinal profile
+ /* eLoss.PrintProfile(); // print longitudinal profile */
cut.ShowResults();
const HEPEnergyType Efinal =
cut.GetCutEnergy() + cut.GetInvEnergy() + cut.GetEmEnergy();
cout << "total cut energy (GeV): " << Efinal / 1_GeV << endl
<< "relative difference (%): " << (Efinal / E0 - 1) * 100 << endl;
- cout << "total dEdX energy (GeV): " << eLoss.GetTotal() / 1_GeV << endl
- << "relative difference (%): " << eLoss.GetTotal() / E0 * 100 << endl;
+ /* cout << "total dEdX energy (GeV): " << eLoss.GetTotal() / 1_GeV << endl */
+ /* << "relative difference (%): " << eLoss.GetTotal() / E0 * 100 << endl; */
}
--
GitLab