IAP GITLAB

Skip to content
Snippets Groups Projects
Commit ac63e87e authored by Alan Coleman's avatar Alan Coleman
Browse files

Merge branch 'update_synchrotron_examples' into 'master'

Small update synchrotron examples

See merge request !633
parents a3ea32ed d3198e44
No related branches found
No related tags found
1 merge request!633Small update synchrotron examples
Pipeline #13865 passed
......@@ -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");
......
......@@ -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),
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment