diff --git a/examples/boundary_example.cpp b/examples/boundary_example.cpp index bca45481a3841f690e80ea5d936063a4ffe91acf..161b875ae3a2bae69e773c33a7b3fdc7160ce9f6 100644 --- a/examples/boundary_example.cpp +++ b/examples/boundary_example.cpp @@ -13,6 +13,7 @@ #include <corsika/framework/random/RNGManager.hpp> #include <corsika/framework/utility/CorsikaFenv.hpp> #include <corsika/framework/core/Logging.hpp> +#include <corsika/output/DummyOutputManager.hpp> #include <corsika/setup/SetupEnvironment.hpp> #include <corsika/setup/SetupStack.hpp> @@ -42,6 +43,7 @@ #include <iostream> #include <limits> #include <typeinfo> +#include <fstream> using namespace corsika; using namespace std; @@ -121,7 +123,7 @@ int main() { ParticleCut cut(50_GeV, true, true); - TrackWriter trackWriter("boundary_tracks.dat"); + TrackWriter trackWriter; MyBoundaryCrossingProcess<true> boundaryCrossing("crossings.dat"); // assemble all processes into an ordered process list @@ -164,7 +166,8 @@ int main() { } // define air shower object, run simulation - Cascade EAS(env, tracking, sequence, stack); + DummyOutputManager output; + Cascade EAS(env, tracking, sequence, output, stack); EAS.run(); diff --git a/examples/cascade_example.cpp b/examples/cascade_example.cpp index eb0938d9fe08e93ac3a3f9d972362d9c58a0ef48..bee5869abfc43e7d938e2c34d94b7d6f508383fb 100644 --- a/examples/cascade_example.cpp +++ b/examples/cascade_example.cpp @@ -14,6 +14,7 @@ #include <corsika/framework/utility/CorsikaFenv.hpp> #include <corsika/framework/core/Logging.hpp> +#include <corsika/output/DummyOutputManager.hpp> #include <corsika/media/Environment.hpp> #include <corsika/media/HomogeneousMedium.hpp> @@ -140,7 +141,7 @@ int main() { // cascade with only HE model ==> HE cut ParticleCut cut(80_GeV, true, true); - TrackWriter trackWriter("tracks.dat"); + TrackWriter trackWriter; BetheBlochPDG eLoss{showerAxis}; // assemble all processes into an ordered process list @@ -148,7 +149,8 @@ int main() { make_sequence(stackInspect, sibyll, sibyllNuc, decay, eLoss, cut, trackWriter); // define air shower object, run simulation - Cascade EAS(env, tracking, sequence, stack); + DummyOutputManager output; + Cascade EAS(env, tracking, sequence, output, stack); EAS.run(); diff --git a/examples/cascade_proton_example.cpp b/examples/cascade_proton_example.cpp index 4b24f3b0148045a2db97546c37c6a2916e3c1d1e..7c17d97eb19a84585a6a10e6019ccb55ebb93a62 100644 --- a/examples/cascade_proton_example.cpp +++ b/examples/cascade_proton_example.cpp @@ -14,6 +14,7 @@ #include <corsika/framework/utility/CorsikaFenv.hpp> #include <corsika/framework/core/Logging.hpp> +#include <corsika/output/DummyOutputManager.hpp> #include <corsika/media/Environment.hpp> #include <corsika/media/HomogeneousMedium.hpp> @@ -22,6 +23,7 @@ #include <corsika/media/MediumPropertyModel.hpp> #include <corsika/media/UniformMagneticField.hpp> + #include <corsika/setup/SetupEnvironment.hpp> #include <corsika/setup/SetupStack.hpp> #include <corsika/setup/SetupTrajectory.hpp> @@ -130,7 +132,7 @@ int main() { // HadronicElasticModel::HadronicElasticInteraction // hadronicElastic(env); - TrackWriter trackWriter("tracks.dat"); + TrackWriter trackWriter; ShowerAxis const showerAxis{injectionPos, Vector{rootCS, 0_m, 0_m, -100_km}, env}; BetheBlochPDG eLoss{showerAxis}; @@ -141,7 +143,8 @@ int main() { auto sequence = make_sequence(pythia, decay, eLoss, cut, trackWriter, stackInspect); // define air shower object, run simulation - Cascade EAS(env, tracking, sequence, stack); + DummyOutputManager output; + Cascade EAS(env, tracking, sequence, output, stack); EAS.run(); cout << "Result: E0=" << E0 / 1_GeV << endl; diff --git a/examples/em_shower.cpp b/examples/em_shower.cpp index e4d8b3cc432a1852f96dbababca6a0c024fc037b..8de0de39d4f249227741f323ddc91f4752c1eb0e 100644 --- a/examples/em_shower.cpp +++ b/examples/em_shower.cpp @@ -17,6 +17,7 @@ #include <corsika/framework/utility/CorsikaFenv.hpp> #include <corsika/framework/process/InteractionCounter.hpp> #include <corsika/framework/core/Logging.hpp> +#include <corsika/output/DummyOutputManager.hpp> #include <corsika/media/Environment.hpp> #include <corsika/media/LayeredSphericalAtmosphereBuilder.hpp> @@ -146,7 +147,7 @@ int main(int argc, char** argv) { corsika::proposal::ContinuousProcess emContinuous(env); InteractionCounter emCascadeCounted(emCascade); - TrackWriter trackWriter("tracks.dat"); + TrackWriter trackWriter; // long. profile; columns for gamma, e+, e- still need to be added LongitudinalProfile longprof{showerAxis}; @@ -159,7 +160,8 @@ int main(int argc, char** argv) { observationLevel, trackWriter); // define air shower object, run simulation setup::Tracking tracking; - Cascade EAS(env, tracking, sequence, stack); + DummyOutputManager output; + Cascade EAS(env, tracking, sequence, output, stack); // to fix the point of first interaction, uncomment the following two lines: // EAS.setNodes(); diff --git a/examples/hybrid_MC.cpp b/examples/hybrid_MC.cpp index 83ca7bd57098a6af751f430077dffc7414fa31e1..0444935039776970b24a3c1e6180001b8c996042 100644 --- a/examples/hybrid_MC.cpp +++ b/examples/hybrid_MC.cpp @@ -23,6 +23,7 @@ #include <corsika/framework/core/PhysicalUnits.hpp> #include <corsika/framework/utility/CorsikaFenv.hpp> #include <corsika/framework/core/Cascade.hpp> +#include <corsika/output/DummyOutputManager.hpp> #include <corsika/framework/geometry/PhysicalGeometry.hpp> #include <corsika/media/Environment.hpp> @@ -243,7 +244,8 @@ int main(int argc, char** argv) { // define air shower object, run simulation setup::Tracking tracking; - Cascade EAS(env, tracking, sequence, stack); + DummyOutputManager output; + Cascade EAS(env, tracking, sequence, output, stack); // to fix the point of first interaction, uncomment the following two lines: // EAS.SetNodes(); diff --git a/examples/vertical_EAS.cpp b/examples/vertical_EAS.cpp index 552bc2bce98d9bd49fbe569be973e5a7e71a1bb6..3b759abfadb2547282e1115010f2b137254fe2be 100644 --- a/examples/vertical_EAS.cpp +++ b/examples/vertical_EAS.cpp @@ -357,3 +357,4 @@ int main(int argc, char** argv) { output.endOfLibrary(); } +}