diff --git a/Documentation/Examples/CMakeLists.txt b/Documentation/Examples/CMakeLists.txt index 30073cfdc3928921ad74ea735b23b0d2b65bb27c..d1bb8d1d9e1f7e9ef755f4dbe30c5c0566c284f7 100644 --- a/Documentation/Examples/CMakeLists.txt +++ b/Documentation/Examples/CMakeLists.txt @@ -25,7 +25,7 @@ add_executable (cascade_example cascade_example.cc) target_compile_options(cascade_example PRIVATE -g) # do not skip asserts target_link_libraries (cascade_example SuperStupidStack CORSIKAunits CORSIKAlogging CORSIKArandom - #~ ProcessSibyll + ProcessSibyll CORSIKAcascade ProcessStackInspector ProcessTrackWriter diff --git a/Documentation/Examples/cascade_example.cc b/Documentation/Examples/cascade_example.cc index 5fa635dfd92f2eb1dd51fb8f5b42a207b1b5ca2b..f157ba9f8a53d6de7b85cbcf1d08176a6f9c800a 100644 --- a/Documentation/Examples/cascade_example.cc +++ b/Documentation/Examples/cascade_example.cc @@ -10,9 +10,9 @@ #include <corsika/cascade/Cascade.h> #include <corsika/process/ProcessSequence.h> +#include <corsika/process/hadronic_elastic_model/HadronicElasticModel.h> #include <corsika/process/stack_inspector/StackInspector.h> #include <corsika/process/tracking_line/TrackingLine.h> -#include <corsika/process/hadronic_elastic_model/HadronicElasticModel.h> #include <corsika/setup/SetupStack.h> #include <corsika/setup/SetupTrajectory.h> @@ -230,20 +230,19 @@ int main() { // setup processes, decays and interactions tracking_line::TrackingLine<setup::Stack> tracking(env); stack_inspector::StackInspector<setup::Stack> p0(true); - + corsika::random::RNGManager::GetInstance().RegisterRandomStream("s_rndm"); corsika::process::sibyll::Interaction sibyll(env); corsika::process::sibyll::Decay decay; ProcessCut cut(8_GeV); - + corsika::random::RNGManager::GetInstance().RegisterRandomStream("HadronicElasticModel"); corsika::process::HadronicElasticModel::HadronicElasticInteraction hadronicElastic(env); - + corsika::process::TrackWriter::TrackWriter trackWriter("tracks.dat"); // assemble all processes into an ordered process list - //~ auto sequence = p0 << sibyll << decay << cut << trackWriter; - auto sequence = hadronicElastic << trackWriter; + auto sequence = p0 << sibyll << decay << hadronicElastic << cut << trackWriter; // cout << "decltype(sequence)=" << type_id_with_cvr<decltype(sequence)>().pretty_name() // << "\n"; @@ -251,7 +250,7 @@ int main() { // setup particle stack, and add primary particle setup::Stack stack; stack.Clear(); - const HEPEnergyType E0 = 10_GeV; + const HEPEnergyType E0 = 1000_GeV; double theta = 0.; double phi = 0.; {