IAP GITLAB

Skip to content
Snippets Groups Projects
Commit b44993dd authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan: Committed by ralfulrich
Browse files

first try real-world history with vertical_EAS

parent dcc9a333
No related branches found
No related tags found
No related merge requests found
...@@ -115,6 +115,7 @@ if (Pythia8_FOUND) ...@@ -115,6 +115,7 @@ if (Pythia8_FOUND)
CORSIKAgeometry CORSIKAgeometry
CORSIKAenvironment CORSIKAenvironment
CORSIKAprocesssequence CORSIKAprocesssequence
CORSIKAhistory # for HistoryObservationPlane
) )
endif() endif()
......
...@@ -37,6 +37,9 @@ ...@@ -37,6 +37,9 @@
#include <corsika/units/PhysicalUnits.h> #include <corsika/units/PhysicalUnits.h>
#include <corsika/utl/CorsikaFenv.h> #include <corsika/utl/CorsikaFenv.h>
#include <corsika/history/HistoryObservationPlane.hpp>
#include <corsika/history/HistorySecondaryView.hpp>
#include <iomanip> #include <iomanip>
#include <iostream> #include <iostream>
#include <limits> #include <limits>
...@@ -63,9 +66,9 @@ void registerRandomStreams(const int seed) { ...@@ -63,9 +66,9 @@ void registerRandomStreams(const int seed) {
random::RNGManager::GetInstance().RegisterRandomStream("urqmd"); random::RNGManager::GetInstance().RegisterRandomStream("urqmd");
random::RNGManager::GetInstance().RegisterRandomStream("proposal"); random::RNGManager::GetInstance().RegisterRandomStream("proposal");
if (seed==0) if (seed == 0)
random::RNGManager::GetInstance().SeedAll(); random::RNGManager::GetInstance().SeedAll();
else else
random::RNGManager::GetInstance().SeedAll(seed); random::RNGManager::GetInstance().SeedAll(seed);
} }
...@@ -81,8 +84,7 @@ int main(int argc, char** argv) { ...@@ -81,8 +84,7 @@ int main(int argc, char** argv) {
feenableexcept(FE_INVALID); feenableexcept(FE_INVALID);
int seed = 0; int seed = 0;
if (argc>4) if (argc > 4) seed = std::stoi(std::string(argv[4]));
seed = std::stoi(std::string(argv[4]));
// initialize random number sequence(s) // initialize random number sequence(s)
registerRandomStreams(seed); registerRandomStreams(seed);
...@@ -205,8 +207,9 @@ int main(int argc, char** argv) { ...@@ -205,8 +207,9 @@ int main(int argc, char** argv) {
process::longitudinal_profile::LongitudinalProfile longprof{showerAxis}; process::longitudinal_profile::LongitudinalProfile longprof{showerAxis};
Plane const obsPlane(showerCore, Vector<dimensionless_d>(rootCS, {0., 0., 1.})); Plane const obsPlane(showerCore, Vector<dimensionless_d>(rootCS, {0., 0., 1.}));
process::observation_plane::ObservationPlane observationLevel(obsPlane, //~ process::observation_plane::ObservationPlane observationLevel(obsPlane,
"particles.dat"); //~ "particles.dat");
corsika::history::HistoryObservationPlane observationLevel{stack, obsPlane};
process::UrQMD::UrQMD urqmd; process::UrQMD::UrQMD urqmd;
process::interaction_counter::InteractionCounter urqmdCounted{urqmd}; process::interaction_counter::InteractionCounter urqmdCounted{urqmd};
...@@ -223,7 +226,9 @@ int main(int argc, char** argv) { ...@@ -223,7 +226,9 @@ int main(int argc, char** argv) {
// define air shower object, run simulation // define air shower object, run simulation
tracking_line::TrackingLine tracking; tracking_line::TrackingLine tracking;
cascade::Cascade EAS(env, tracking, sequence, stack); cascade::Cascade<decltype(tracking), decltype(sequence), decltype(stack),
corsika::history::HistorySecondaryView<corsika::setup::StackView> >
EAS(env, tracking, sequence, stack);
// to fix the point of first interaction, uncomment the following two lines: // to fix the point of first interaction, uncomment the following two lines:
// EAS.SetNodes(); // EAS.SetNodes();
......
...@@ -20,7 +20,7 @@ ...@@ -20,7 +20,7 @@
namespace corsika::history { namespace corsika::history {
namespace detail { namespace detail {
auto hist_factory() { inline auto hist_factory() {
/*auto h = boost::histogram::make_histogram( /*auto h = boost::histogram::make_histogram(
boost::histogram::axis::regular<double, boost::histogram::axis::transform::log>{ boost::histogram::axis::regular<double, boost::histogram::axis::transform::log>{
130, 1e8, 1e21, "muon energy/eV"}, 130, 1e8, 1e21, "muon energy/eV"},
......
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