diff --git a/Documentation/Examples/vertical_EAS.cc b/Documentation/Examples/vertical_EAS.cc index 76d782832c78825c266c4f8183be9eb8ac2773cf..8c8ba1a89a08ee32e4c1db3c260b85a7093dc285 100644 --- a/Documentation/Examples/vertical_EAS.cc +++ b/Documentation/Examples/vertical_EAS.cc @@ -31,8 +31,6 @@ #include <corsika/process/sibyll/Interaction.h> #include <corsika/process/sibyll/NuclearInteraction.h> -//#include <corsika/process/pythia/Decay.h> - #include <corsika/process/urqmd/UrQMD.h> #include <corsika/process/particle_cut/ParticleCut.h> @@ -64,7 +62,6 @@ using namespace corsika::units::si; void registerRandomStreams() { random::RNGManager::GetInstance().RegisterRandomStream("cascade"); random::RNGManager::GetInstance().RegisterRandomStream("s_rndm"); - // random::RNGManager::GetInstance().RegisterRandomStream("pythia"); random::RNGManager::GetInstance().RegisterRandomStream("UrQMD"); random::RNGManager::GetInstance().SeedAll(); @@ -103,9 +100,9 @@ int main() { // setup particle stack, and add primary particle setup::Stack stack; stack.Clear(); - const Code beamCode = Code::Proton; + const Code beamCode = Code::RhoMinus; //Code::KStarMinus; auto const mass = particles::GetMass(beamCode); - const HEPEnergyType E0 = 0.1_PeV; + const HEPEnergyType E0 = 0.1_TeV; double theta = 0.; double phi = 0.; @@ -155,7 +152,7 @@ int main() { particles::Code::KMinus, particles::Code::K0Long, particles::Code::K0Short}; process::sibyll::Interaction sibyll; - process::sibyll::NuclearInteraction sibyllNuc(sibyll, env); + //process::sibyll::NuclearInteraction sibyllNuc(sibyll, env); process::sibyll::Decay decay; //process::pythia::Decay decay(trackedHadrons); @@ -173,7 +170,7 @@ int main() { process::UrQMD::UrQMD urqmd; - auto sibyllSequence = sibyll << sibyllNuc; + auto sibyllSequence = sibyll; // << sibyllNuc; process::switch_process::SwitchProcess switchProcess(urqmd, sibyllSequence, 55_GeV); auto sequence = switchProcess << decay << eLoss << cut << observationLevel << trackWriter;