From d3198e447ade305be3ae5d08bf3faf40a3726bcb Mon Sep 17 00:00:00 2001
From: Nikos Karastathis <nckcoop@gmail.com>
Date: Sun, 15 Sep 2024 17:15:10 +0000
Subject: [PATCH] Small update synchrotron examples
---
examples/physics_examples/synchrotron_test_C8tracking.cpp | 8 +++++++-
.../physics_examples/synchrotron_test_manual_tracking.cpp | 4 ++--
2 files changed, 9 insertions(+), 3 deletions(-)
diff --git a/examples/physics_examples/synchrotron_test_C8tracking.cpp b/examples/physics_examples/synchrotron_test_C8tracking.cpp
index 81d2e56e5..d03d82d5a 100644
--- a/examples/physics_examples/synchrotron_test_C8tracking.cpp
+++ b/examples/physics_examples/synchrotron_test_C8tracking.cpp
@@ -54,7 +54,13 @@ int main() {
logging::set_level(logging::level::warn);
corsika_logger->set_pattern("[%n:%^%-8l%$] custom pattern: %v");
- CORSIKA_LOG_INFO("Synchrotron radiation");
+ CORSIKA_LOG_INFO("Synchrotron radiation -- CORSIKA 8 tracking");
+ CORSIKA_LOG_INFO(" ");
+ CORSIKA_LOG_WARN(
+ "Adjust the resolution of the circular trajectory via the "
+ "maxMagneticDeflectionAngle parameter which can be found in "
+ "~/corsika/corsika/detail/modules/tracking/TrackingLeapFrogCurved.inl and "
+ "recompile example. Recommended value is 0.00001 rad.");
feenableexcept(FE_INVALID);
RNGManager<>::getInstance().registerRandomStream("cascade");
diff --git a/examples/physics_examples/synchrotron_test_manual_tracking.cpp b/examples/physics_examples/synchrotron_test_manual_tracking.cpp
index d3b62ac25..210442899 100644
--- a/examples/physics_examples/synchrotron_test_manual_tracking.cpp
+++ b/examples/physics_examples/synchrotron_test_manual_tracking.cpp
@@ -137,8 +137,8 @@ int main() {
LengthType const radius{100_m};
TimeType timeCounter{0._s};
- // loop over all the tracks twice (this produces 2 pulses)
- for (size_t i = 0; i <= (n_points)*2; i++) {
+ // loop over the circular trajectory (this produces 1 pulse)
+ for (size_t i = 0; i <= (n_points); i++) {
Point const point_1(rootCS, {radius * cos(M_PI * 2 * i / n_points),
radius * sin(M_PI * 2 * i / n_points), 0_m});
Point const point_2(rootCS, {radius * cos(M_PI * 2 * (i + 1) / n_points),
--
GitLab