IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 8fb7bfe5 authored by ralfulrich's avatar ralfulrich
Browse files

new setting, so that ObservationLevel process sees a few particles

parent da7bee3a
No related branches found
No related tags found
No related merge requests found
...@@ -240,6 +240,9 @@ public: ...@@ -240,6 +240,9 @@ public:
// The example main program for a particle cascade // The example main program for a particle cascade
// //
int main() { int main() {
const LengthType height_atmosphere = 112.8_km;
feenableexcept(FE_INVALID); feenableexcept(FE_INVALID);
// initialize random number sequence(s) // initialize random number sequence(s)
random::RNGManager::GetInstance().RegisterRandomStream("cascade"); random::RNGManager::GetInstance().RegisterRandomStream("cascade");
...@@ -281,7 +284,7 @@ int main() { ...@@ -281,7 +284,7 @@ int main() {
// random::RNGManager::GetInstance().RegisterRandomStream("pythia"); // random::RNGManager::GetInstance().RegisterRandomStream("pythia");
// process::pythia::Decay decay(trackedHadrons); // process::pythia::Decay decay(trackedHadrons);
ProcessCut cut(20_GeV); ProcessCut cut(20_GeV);
ObservationLevel obsLevel(1400_m); ObservationLevel obsLevel(height_atmosphere - 2000_m);//1400_m);
// random::RNGManager::GetInstance().RegisterRandomStream("HadronicElasticModel"); // random::RNGManager::GetInstance().RegisterRandomStream("HadronicElasticModel");
// process::HadronicElasticModel::HadronicElasticInteraction // process::HadronicElasticModel::HadronicElasticInteraction
...@@ -327,7 +330,7 @@ int main() { ...@@ -327,7 +330,7 @@ int main() {
cout << "input angles: theta=" << theta << " phi=" << phi << endl; cout << "input angles: theta=" << theta << " phi=" << phi << endl;
cout << "input momentum: " << plab.GetComponents() / 1_GeV << endl; cout << "input momentum: " << plab.GetComponents() / 1_GeV << endl;
Point pos(rootCS, 0_m, 0_m, Point pos(rootCS, 0_m, 0_m,
112.8_km); // this is the CORSIKA 7 start of atmosphere/universe height_atmosphere); // this is the CORSIKA 7 start of atmosphere/universe
stack.AddParticle(std::tuple<particles::Code, units::si::HEPEnergyType, stack.AddParticle(std::tuple<particles::Code, units::si::HEPEnergyType,
corsika::stack::MomentumVector, geometry::Point, corsika::stack::MomentumVector, geometry::Point,
units::si::TimeType, unsigned short, unsigned short>{ units::si::TimeType, unsigned short, unsigned short>{
......
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