IAP GITLAB

Skip to content
Snippets Groups Projects

History

Merged Ralf Ulrich requested to merge history into master
Compare and Show latest version
6 files
+ 18
15
Compare changes
  • Side-by-side
  • Inline
Files
6
@@ -37,8 +37,6 @@
#include <corsika/units/PhysicalUnits.h>
#include <corsika/utl/CorsikaFenv.h>
#include <corsika/stack/history/HistoryObservationPlane.hpp>
#include <iomanip>
#include <iostream>
#include <limits>
@@ -204,9 +202,8 @@ int main(int argc, char** argv) {
process::longitudinal_profile::LongitudinalProfile longprof{showerAxis};
Plane const obsPlane(showerCore, Vector<dimensionless_d>(rootCS, {0., 0., 1.}));
//~ process::observation_plane::ObservationPlane observationLevel(obsPlane,
//~ "particles.dat");
corsika::history::HistoryObservationPlane observationLevel{stack, obsPlane};
process::observation_plane::ObservationPlane observationLevel(obsPlane,
"particles.dat");
process::UrQMD::UrQMD urqmd;
process::interaction_counter::InteractionCounter urqmdCounted{urqmd};
@@ -246,8 +243,6 @@ int main(int argc, char** argv) {
cout << "total dEdX energy (GeV): " << eLoss.GetTotal() / 1_GeV << endl
<< "relative difference (%): " << eLoss.GetTotal() / E0 * 100 << endl;
observationLevel.print();
auto const hists = sibyllCounted.GetHistogram() + sibyllNucCounted.GetHistogram() +
urqmdCounted.GetHistogram();
hists.saveLab("inthist_lab.txt");
Loading