7 km

Documentation/Examples/vertical_EAS.cc

@@ -137,7 +137,7 @@ int main() {
   Line const line(injectionPos, plab.normalized() * 1_m * 1_Hz);
   auto const velocity = line.GetV0().norm();
 
-  auto const observationHeight = 1.452_km;
+  auto const observationHeight = 7_km;
 
   setup::Trajectory const showerAxis(line, (112.8_km - observationHeight) / velocity);
 
@@ -161,7 +161,7 @@ int main() {
   process::pythia::Decay decay(trackedHadrons);
   process::particle_cut::ParticleCut cut(60_GeV);
 
-  process::track_writer::TrackWriter trackWriter("tracks.dat");
+//  process::track_writer::TrackWriter trackWriter("tracks.dat");
   process::energy_loss::EnergyLoss eLoss(showerAxis);
 
   Plane const obsPlane(Point(rootCS, 0_m, 0_m, observationHeight),
@@ -175,8 +175,8 @@ int main() {
 
   auto sibyllSequence = sibyll << sibyllNuc;
   process::switch_process::SwitchProcess switchProcess(urqmd, sibyllSequence, 55_GeV);
-  auto sequence = switchProcess << decay << eLoss << cut << observationLevel
-                               /* << trackWriter */;
+  auto sequence = switchProcess << decay << eLoss << cut << observationLevel;
+//                                << trackWriter;
 
   // define air shower object, run simulation
   tracking_line::TrackingLine tracking;