IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 04972957 authored by Nikos Karastathis's avatar Nikos Karastathis :ocean: Committed by Maximilian Reininghaus
Browse files

uncomment accordingly to use fluka

parent 1fc28998
No related branches found
No related tags found
1 merge request!534default configuration for ICRC 2023
...@@ -58,6 +58,7 @@ ...@@ -58,6 +58,7 @@
#include <corsika/modules/StackInspector.hpp> #include <corsika/modules/StackInspector.hpp>
#include <corsika/modules/UrQMD.hpp> #include <corsika/modules/UrQMD.hpp>
#include <corsika/modules/thinning/EMThinning.hpp> #include <corsika/modules/thinning/EMThinning.hpp>
//#include <corsika/modules/FLUKA.hpp>
#include <corsika/setup/SetupStack.hpp> #include <corsika/setup/SetupStack.hpp>
#include <corsika/setup/SetupTrajectory.hpp> #include <corsika/setup/SetupTrajectory.hpp>
...@@ -97,6 +98,7 @@ long registerRandomStreams(long seed) { ...@@ -97,6 +98,7 @@ long registerRandomStreams(long seed) {
RNGManager<>::getInstance().registerRandomStream("epos"); RNGManager<>::getInstance().registerRandomStream("epos");
RNGManager<>::getInstance().registerRandomStream("pythia"); RNGManager<>::getInstance().registerRandomStream("pythia");
RNGManager<>::getInstance().registerRandomStream("urqmd"); RNGManager<>::getInstance().registerRandomStream("urqmd");
// RNGManager<>::getInstance().registerRandomStream("fluka");
RNGManager<>::getInstance().registerRandomStream("proposal"); RNGManager<>::getInstance().registerRandomStream("proposal");
RNGManager<>::getInstance().registerRandomStream("thinning"); RNGManager<>::getInstance().registerRandomStream("thinning");
if (seed == 0) { if (seed == 0) {
...@@ -430,6 +432,9 @@ int main(int argc, char** argv) { ...@@ -430,6 +432,9 @@ int main(int argc, char** argv) {
corsika::urqmd::UrQMD urqmd; corsika::urqmd::UrQMD urqmd;
InteractionCounter urqmdCounted(urqmd); InteractionCounter urqmdCounted(urqmd);
// // until the CI containers have fluka, we keep urqmd. Switch to fluka by uncommenting accordingly.
// corsika::fluka::Interaction leInt{env};
// InteractionCounter leIntCounted{leInt};
StackInspector<setup::Stack<EnvType>> stackInspect(10000, false, E0); StackInspector<setup::Stack<EnvType>> stackInspect(10000, false, E0);
// assemble all processes into an ordered process list // assemble all processes into an ordered process list
...@@ -441,6 +446,9 @@ int main(int argc, char** argv) { ...@@ -441,6 +446,9 @@ int main(int argc, char** argv) {
}; };
auto hadronSequence = auto hadronSequence =
make_select(EnergySwitch(heHadronModelThreshold), urqmdCounted, heCounted); make_select(EnergySwitch(heHadronModelThreshold), urqmdCounted, heCounted);
// // uncomment below and comment the above hadron sequence to use fluka
// auto hadronSequence =
// make_select(EnergySwitch(heHadronModelThreshold), leIntCounted, heCounted);
auto decaySequence = make_sequence(decayPythia, decaySibyll); auto decaySequence = make_sequence(decayPythia, decaySibyll);
// observation plane // observation plane
...@@ -512,6 +520,8 @@ int main(int argc, char** argv) { ...@@ -512,6 +520,8 @@ int main(int argc, char** argv) {
observationLevel.getEnergyGround() / 1_GeV, (Efinal / E0 - 1) * 100); observationLevel.getEnergyGround() / 1_GeV, (Efinal / E0 - 1) * 100);
auto const hists = heCounted.getHistogram() + urqmdCounted.getHistogram(); auto const hists = heCounted.getHistogram() + urqmdCounted.getHistogram();
// uncomment to use fluka
// auto const hists = heCounted.getHistogram() + leIntCounted.getHistogram();
save_hist(hists.labHist(), labHist_file, true); save_hist(hists.labHist(), labHist_file, true);
save_hist(hists.CMSHist(), cMSHist_file, true); save_hist(hists.CMSHist(), cMSHist_file, true);
......
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