diff --git a/examples/cascade_examples/radio_em_shower.cpp b/examples/cascade_examples/radio_em_shower.cpp
index 7b1740e1f9413eb4b7546dea26baeb36db3cebe4..9631ce895caeb73f2a24b798760ca190d31481a1 100644
--- a/examples/cascade_examples/radio_em_shower.cpp
+++ b/examples/cascade_examples/radio_em_shower.cpp
@@ -119,7 +119,7 @@ int main(int argc, char** argv) {
   Point const center{rootCS, 0_m, 0_m, 0_m};
 
   double const refractive_index = 1.000327;
-  MagneticFieldVector bField{rootCS, 50_uT, 0_T, 0_T};
+  MagneticFieldVector const bField{rootCS, 50_uT, 0_T, 0_T};
   create_5layer_atmosphere<EnvironmentInterface, MyExtraEnv>(
       env, AtmosphereId::LinsleyUSStd, center, refractive_index, Medium::AirDry1Atm,
       bField);
@@ -132,19 +132,19 @@ int main(int argc, char** argv) {
   for (auto [pcode, energy] : energy_resolution)
     set_energy_production_threshold(pcode, energy);
 
-  const Code beamCode = Code::Electron;
+  Code const beamCode = Code::Electron;
   auto const mass = get_mass(beamCode);
-  const HEPEnergyType E0 = 1_GeV * std::stof(std::string(argv[1]));
-  double theta = 0.;
+  HEPEnergyType const E0 = 1_GeV * std::stof(std::string(argv[1]));
+  double const theta = 0.;
   auto const thetaRad = theta / 180. * M_PI;
 
-  HEPMomentumType P0 = calculate_momentum(E0, mass);
+  HEPMomentumType const P0 = calculate_momentum(E0, mass);
   auto momentumComponents = [](double thetaRad, HEPMomentumType ptot) {
     return std::make_tuple(ptot * sin(thetaRad), 0_eV, -ptot * cos(thetaRad));
   };
 
   auto const [px, py, pz] = momentumComponents(thetaRad, P0);
-  auto plab = MomentumVector(rootCS, {px, py, pz});
+  auto const plab = MomentumVector(rootCS, {px, py, pz});
 
   auto const observationHeight = 1.4_km + constants::EarthRadius::Mean;
   auto const injectionHeight = 112.75_km + constants::EarthRadius::Mean;
@@ -163,8 +163,8 @@ int main(int argc, char** argv) {
 
   // Radio antennas and relevant information
   // the antenna time variables
-  const TimeType duration{1e-6_s};
-  const InverseTimeType sampleRate{1e+9_Hz};
+  TimeType const duration{1e-6_s};
+  InverseTimeType const sampleRate{1e+9_Hz};
 
   // the detector (aka antenna collection) for CoREAS and ZHS
   AntennaCollection<TimeDomainAntenna> detectorCoREAS;