From 56db68b6561de2692ecb2d394725b4fd49764243 Mon Sep 17 00:00:00 2001
From: Nikos Karastathis <n.karastathis@kit.edu>
Date: Tue, 23 Mar 2021 16:04:02 +0100
Subject: [PATCH] small updates to radio shower 2 example
---
examples/radio_shower2.cpp | 19 +++----------------
1 file changed, 3 insertions(+), 16 deletions(-)
diff --git a/examples/radio_shower2.cpp b/examples/radio_shower2.cpp
index e7c3e0398..8859c9651 100644
--- a/examples/radio_shower2.cpp
+++ b/examples/radio_shower2.cpp
@@ -62,19 +62,20 @@ using namespace corsika;
using namespace std;
//
-// The example main program for a particle cascade
+// A simple shower to get the electric field trace of an electron
//
int main() {
logging::set_level(logging::level::info);
corsika_logger->set_pattern("[%n:%^%-8l%$] custom pattern: %v");
- std::cout << "cascade_proton_example" << std::endl;
+ std::cout << "radio_shower2" << std::endl;
feenableexcept(FE_INVALID);
// initialize random number sequence(s)
RNGManager::getInstance().registerRandomStream("cascade");
+ // This environment needs a hardcoded refractive index in the propagator at the moment
// setup environment, geometry
using EnvType = setup::Environment;
EnvType env;
@@ -145,36 +146,22 @@ int main() {
setup::Tracking tracking;
StackInspector<setup::Stack> stackInspect(1, true, E0);
-
-// ParticleCut cut(60_GeV, true, true);
// put radio process here
RadioProcess<decltype(detector), CoREAS<decltype(detector),
decltype(StraightPropagator(env))>, decltype(StraightPropagator(env))>
coreas(detector, env);
-
TrackWriter trackWriter("tracks.dat");
ShowerAxis const showerAxis{injectionPos, Vector{rootCS, 0_m, 0_m, -100_km}, env};
-// BetheBlochPDG eLoss{showerAxis};
// assemble all processes into an ordered process list
- // auto sequence = make_sequence(sibyll, sibyllNuc, decay, eLoss, cut, trackWriter,
- // stackInspect); auto sequence = make_sequence(sibyll, decay, eLoss, cut, trackWriter,
- // stackInspect);
auto sequence = make_sequence(coreas, trackWriter, stackInspect);
// define air shower object, run simulation
Cascade EAS(env, tracking, sequence, stack);
EAS.run();
- cout << "Result: E0=" << E0 / 1_GeV << endl;
- cut.showResults();
- const HEPEnergyType Efinal =
- cut.getCutEnergy() + cut.getInvEnergy() + cut.getEmEnergy();
- cout << "total energy (GeV): " << Efinal / 1_GeV << endl
- << "relative difference (%): " << (Efinal / E0 - 1.) * 100 << endl;
-
// get radio output
coreas.writeOutput();
}
--
GitLab