IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 5f34a359 authored by ralfulrich's avatar ralfulrich
Browse files

new setting, so that ObservationLevel process sees a few particles

parent 979a7a7c
No related branches found
No related tags found
No related merge requests found
......@@ -236,6 +236,9 @@ public:
// The example main program for a particle cascade
//
int main() {
const LengthType height_atmosphere = 112.8_km;
feenableexcept(FE_INVALID);
// initialize random number sequence(s)
random::RNGManager::GetInstance().RegisterRandomStream("cascade");
......@@ -283,7 +286,7 @@ int main() {
process::sibyll::NuclearInteraction sibyllNuc(env, sibyll);
process::sibyll::Decay decay(trackedHadrons);
ProcessCut cut(20_GeV);
ObservationLevel obsLevel(1400_m);
ObservationLevel obsLevel(height_atmosphere - 2000_m);//1400_m);
process::TrackWriter::TrackWriter trackWriter("tracks.dat");
process::EnergyLoss::EnergyLoss eLoss;
......@@ -318,7 +321,7 @@ int main() {
cout << "input angles: theta=" << theta << " phi=" << phi << endl;
cout << "input momentum: " << plab.GetComponents() / 1_GeV << endl;
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,
corsika::stack::MomentumVector, geometry::Point,
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