Choose timezone
Your profile timezone:
CORSIKA Cosmic Ray Simulation Workshop Karlsruhe
→
Europe/Berlin
Hotel ACHAT PLAZA Karlsruhe
Hotel ACHAT PLAZA Karlsruhe
Mendelssohnplatz, 76131 Karlsruhe
https://achat-hotels.com/hotel/karlsruhe-city
Ralf Ulrich
(KIT)
Description
General and open workshop to discuss technology and science of air shower simulations, and meeting of the CORSIKA 8 Project.
Support
Participants
-
- 9:30 AM → 10:30 AM
-
10:30 AM
→
11:00 AM
Break 30m
- 11:00 AM → 11:10 AM
-
11:10 AM
→
12:25 PM
Status and progress of air shower simulations¶
The current state-of-the-art in air shower simulations, with a focus on the status and progress of CORSIKA 8.
- 11:10 AM
-
11:25 AM
Speaker: Ralf Ulrich (KIT)
-
11:45 AM
Speaker: Maximilian Reininghaus (KIT / IKP)
-
12:05 PM
This talk describes the lepton propagator PROPOSAL, its current capabilities and its possible use in future versions of CORSIKA, and necessary extensions to function as a replacement for the EGS4 electromagnetic shower simulation code.
Speaker: Alexander Sandrock (Technische Universität Dortmund)
-
12:30 PM
→
2:00 PM
Lunch 1h 30m
-
2:00 PM
→
3:00 PM
Applications and requirements of CORSIKA¶
Contributions from experimental collaborations, phenomenologist or theorists about important application of CORSIKA, and future wishes for new functionality.
-
2:00 PM
The KM3NeT Collaboration has successfully deployed its first detection units in the Mediterranean Sea in 2016 (ARCA) and 2017 (ORCA). The sample of data collected between December 2016 and December 2017 has been used to measure the atmospheric muon flux at two different depths under the sea level: 3.5 km with ARCA and 2.5 km with ORCA. The atmospheric muon flux represents an abundant signal for a neutrino telescope and can be used to test the reliability of the Monte Carlo simulation chain. In this case, the measurements are compared to Monte Carlo simulations based on MUPAGE and CORSIKA codes. MUPAGE is fast, parametric simulation code developed for the ANTARES experiment that generates muon events, induced by cosmic rays impinging the Earth atmosphere, at different undersea depths. The atmospheric muon distributions at the sea level obtained with MUPAGE/CORSIKA are propagated to the detector with the MUSIC code, a 3D program that takes into account all muon energy loss processes. Finally, all muons arriving at the detector are transported through the active volume, Cherenkov light is generated and the same triggering algorithms used during data taking are applied. The main features of the simulation and reconstruction chain are discussed and presented.
Speaker: Piotr Kalaczyński (NCBJ Warsaw) -
2:20 PM
The Baikal experiment was proposed to register cosmological neutrinos and map the high-energy neutrino sky in the Southern Hemisphere. It will use a km3-scale high-energy neutrino telescope in Lake Baikal. The selection of events from neutrino induced upward going muons, suggests a fairly reliable estimation of the expected background from atmospheric muons. The procedure for simulating background events from atmospheric muons in an array is performed in several step. The CORSIKA7.64 was used to simulate the flux of atmospheric muons at sea level with appropriate chemical composition of the primary cosmic rays. The muon propagation through water and rock to the array level was then simulated with the MUM code. As the last step of simulation chain simulation of the detector response to the Cherenkov radiation of muons by taking into account the features of array measuring systems was performed. The main features of the CORSIKA Monte Carlo code and following steps of the full simulation chain will be briefly reported.
Speaker: Blahoslav Pastircak (IEP SAS Kosice) -
2:35 PM
New measurements at the LHC are needed to resolve the Muon Puzzle in cosmic-ray induced air showers and to reduce the systematic uncertainties of air shower simulations overall. In particular, the very forward production of light hadrons has to be studied with LHCb and other forward detectors. The LHC has accelerated proton and lead beams so far (and a bit of Xenon). The hadronic interaction models used to simulate air showers are primarily tuned to p-p data, while the most common collision in an air shower is between a high energy hadron and a oxygen or nitrogen nucleus. The predictions for these interactions vary between models by 50 % in multiplicity. This variance needs to be reduced to 10 % and the fraction of energy which ends up in neutral pions need to be clarified. This can be achieved by studying proton-oxygen collisions at the LHC, which have been proposed for 2023. We review the science case for measuring proton-oxygen collisions and report on the status of the proposal.
Speaker: Hans Dembinski (Max Planck Institute for Nuclear Physics, Heidelberg)
-
2:00 PM
-
3:00 PM
→
3:30 PM
Coffee 30m
-
3:30 PM
→
5:05 PM
Status and progress of air shower simulations¶
The current state-of-the-art in air shower simulations, with a focus on the status and progress of CORSIKA 8.
-
3:30 PM
Speaker: Takashi Sako (ICRR, University of Tokyo)
-
3:50 PM
Tutorial part (DIY) and introduction, overview, results
Speaker: Ralf Ulrich (KIT) -
4:45 PM
Template meta-programming (TMP) is a fancy and mysterious sounding topic. CORSIKA 8 is using TMP a lot, so in this talk I will explain when it is useful and demonstrate a few key techniques and best practices, based on my experience as the author of Boost.Histogram, which is another TMP-heavy C++ library.
Speaker: Hans Dembinski (Max Planck Institute for Nuclear Physics, Heidelberg)
-
3:30 PM
-
-
9:00 AM
→
12:45 PM
Celebration: 30years of CORSIKA¶
The first and original version of the program "COsmic Ray SImulation for KAscade" was released in 1989 - 30 years in the past. We want to celebrate the extraordinary development CORSIKA has taken until today.
- 9:00 AM
- 9:10 AM
-
9:30 AM
From Quarks to Air Showers: A Historic Review of the Development of the Simulation Program System ASICO, Results and Impact, and its Evolution to CORSIKA¶ 30m
I set the scene with a brief presentation of the status of our knowledge of high energy cosmic rays and particle physics around 1960. This is followed by an overview of the quark hunt in the early 1960s, when it was believed that free quarks exist, and the search that was undertaken at the Niels Bohr Institute in Copenhagen in 1965, and elsewhere. Instead of quarks, the unexpected copious production of
-
10:05 AM
Speaker: Capdevielle Jean-Noel
-
10:25 AM
Coffee 30m
-
11:00 AM
Speaker: Johannes Knapp
- 11:25 AM
-
11:50 AM
Speaker: Konrad Bernlöhr
- 12:15 PM
-
12:45 PM
→
2:15 PM
Lunch 1h 30m
-
2:15 PM
→
5:40 PM
Applications and requirements of CORSIKA¶
Contributions from experimental collaborations, phenomenologist or theorists about important application of CORSIKA, and future wishes for new functionality.
-
2:15 PM
Speaker: Sergey Ostapchenko
-
2:55 PM
Speaker: Anna Schäfer
-
3:30 PM
Coffee 30m
-
4:00 PM
Zero degree calorimeters at LHC like LHCf or in CMS provide the opportunity to identify pion-exchange reactions at very high energies p+p --> p+pi+n where the neutron escapes at zero degree and the pion and the proton subsequently collide in the detector. This provides a unique opportunity to learn about pion collisions, which are one of the largest uncertainties in air showers.
Speaker: Tanguy Pierog (IKP, KIT) - 4:20 PM
-
4:40 PM
Speaker: Takashi Sako (ICRR, University of Tokyo)
- 5:00 PM
-
5:15 PM
Speaker: Jean-Noel Capdevielle
-
2:15 PM
-
9:00 AM
→
12:45 PM
-
-
9:00 AM
→
10:30 AM
Applications and requirements of CORSIKA¶
Contributions from experimental collaborations, phenomenologist or theorists about important application of CORSIKA, and future wishes for new functionality.
- 9:00 AM
-
9:15 AM
The detector is currently in an R&D phase for an upgrade that will include adding an array of scintillator panels and radio antennas, to the existing array of ice Cherenkov detectors. Understanding the potential physics output of such an array requires a large library of CoREAS simulations. I will present how we are currently using CORSIKA to develop the radio component of the IceTop upgrade as well as how our needs could be better met with improvements to the simulation software.
Speaker: Alan Coleman (University of Delaware) -
9:35 AM
There are many reasons for producing millions of equitype CORSIKA simulations, e.g. to cover necessary parameter space for the model or to fit experimental data. Moreover, sometimes it is necessary to perform fast pre-simulations, e.g. with CONEX or other computation-reduced options. We present a SIMulation Manager (SIMM) framework for steering, running and evaluating status of CORSIKA simulations. The main features of SIMM is the interface for CORSIKA configuration, scheduling the simulations for running on distributed resources, integration with Auger Offline software and flexible plugin-based framework for custom simulation sequences (i.e. CONEX-CORSIKA-CoREAS with evaluation at each step). The software is implemented on Python with SQLAlchemy framework, which allows one to use any of SQL engines for metadata storage. We are successfully using SIMM more than five years for the handling simulations for Tunka-Rex and TAIGA experiments. The ideas developed in the frame of SIMM can be further used in future Astroparticle Data Life Cycle frameworks.
Speaker: Dmitriy Kostunin (DESY) -
9:55 AM
CORELib is a collection of simulated showers induced by primary cosmic rays and a common framework to produce and manage such large production. Indeed, CORELib provides high statics background simulation for various astroparticles and astronomy experiments. The simulations are obtained employing widely used CORSIKA simulation framework and the GRID infrastructure as computational resources; in the current productions the observation level is placed at sea level, while the Standard European Atmospheric Model is used. The parameters of simulations are available and easily accessible to the users. Two different primary energy spectra are considered: one with -2 spectral index and one with flat spectrum. Four different high energy models are evaluated: QGSJET01 with CHARM option, QGSJET01 with TAULEP option, QGSJETII with TAULEP option and EPOS with TAULEP option. The common low energy model chosen is GHEISHA. Protons and Heavy Nuclei (He, C, N, O, Fe) induced showers are simulated. Proton induced showers two production are available: with and without Cherenkov radiation. All the information about every single shower is stored in SQLite database queryable by the user. In order to simplify the access at the library, a dedicate tool named CORANT is used to translate from binary to ASCII the CORSIKA output files. The information about showers is placed in separated files, one for each particles class (EM, Hadrons, Muons, Neutrinos ). The whole production of CORELib will be stored at the CNAF, the Information Technology national center of INFN (Italian Institute for Nuclear Physics) . It can be downloaded through gridFTP with X509 certificate.
Speaker: Simona Maria Stellacci
-
10:30 AM
→
11:00 AM
Coffee 30m
-
11:00 AM
→
12:35 PM
Status and progress of air shower simulations¶
The current state-of-the-art in air shower simulations, with a focus on the status and progress of CORSIKA 8.
-
11:00 AM
This talk intends to give technical insights into the propagation process of PROPOSAL, a tool to propagate leptons through media. Furthermore, we are going to discuss the requirements and challenges regarding possible use cases for PROPOSAL in CORSIKA in the future.
Speaker: Jean-Marco Alameddine -
11:20 AM
Emca is a new python simulation package, based on MCEq, for the calculation of electromagnetic particle fluxes in various materials. The calculations are based on a cascade equation approach, allowing for an iterative calculation of showers. This allows for fast and efficient calculations of cascades, allowing easier testing and prototyping compared to a Monte Carlo approach. The methodology of EmCa and comparisons between it and CONEX will be shown.
Speaker: Stephan Meighen-Berger -
11:40 AM
Tutorial part for DIY, introduction and concepts
Speaker: felix riehn (LIP, Lisbon)
-
11:00 AM
-
12:35 PM
→
2:00 PM
Lunch 1h 25m
-
2:00 PM
→
5:00 PM
Status and progress of air shower simulations¶
The current state-of-the-art in air shower simulations, with a focus on the status and progress of CORSIKA 8.
-
2:00 PM
In high-energy physics (HEP), C++ is still the dominant language, but Python is continuously growing and should overtake C++ in the near future. In the industry, Python is already the dominant language for data science and driving the leading frameworks for machine learning. In this talk, I will show these and other interesting facts and explain how a "slow" interpreted language like Python is able to beat a "fast" language like C++, in a field where code execution speed actually matters. I will present the Scikit-HEP project that aims to provide key functionality in Python for HEP analyses that's currently missing, like suitable histograms. Finally, I will argue why the upcoming ROOT 7, the first backward incompatible change in ROOT's history and a bold step forward, will not win the users back.
Speaker: Hans Dembinski (Max Planck Institute for Nuclear Physics, Heidelberg) - 2:30 PM
- 2:50 PM
-
3:10 PM
Speaker: Tim Huege (KIT)
-
3:25 PM
Coffee 30m
-
2:00 PM
-
6:00 PM
→
9:00 PM
Workshop Dinner 3h Schützenhaus am Turmberg
Schützenhaus am Turmberg
Schützenhaus auf dem Turmberg Jean-Ritzert-Strasse 8 | 76227 KarlsruheTelefon (0721) 94318122
info@schuetzenhaus-turmberg.de
https://www.schuetzenhaus-turmberg.de/see also: https://indico.kit.edu/event/529/page/475-circular-5th-june
-
9:00 AM
→
10:30 AM
-
-
9:00 AM
→
9:40 AM
Applications and requirements of CORSIKA¶
Contributions from experimental collaborations, phenomenologist or theorists about important application of CORSIKA, and future wishes for new functionality.
-
9:00 AM
The Cherenkov Telescope Array (CTA) is the next generation ground-based array of IACTs, the world’s largest and most sensitive high-energy gamma-ray observatory. The full setup will consist of more than 100 telescopes located in the northern (CTA-N) and southern hemispheres (CTA-S). The Observatorio del Roque de Los Muchachos (ORM) on the LaPalma island, Spain is selected as a site for the northern array (CTA-N).
In order to determine the influence of clouds on the atmospheric transparency and performance of CTA-N, transmittance simulations and simulations of instrument response in such an atmosphere are being performed.
Speaker: Mario Pecimotika (University of Rijeka, Department of Physics, Croatia) - 9:20 AM
-
9:00 AM
-
9:40 AM
→
10:20 AM
Status and progress of air shower simulations¶
The current state-of-the-art in air shower simulations, with a focus on the status and progress of CORSIKA 8.
-
9:40 AM
Speaker: Maximilian Reininghaus (KIT / IKP)
-
9:40 AM
-
10:20 AM
→
10:50 AM
Coffee 30m
-
10:50 AM
→
12:20 PM
Status and progress of air shower simulations¶
The current state-of-the-art in air shower simulations, with a focus on the status and progress of CORSIKA 8.
-
10:50 AM
Overview of different features and possibilites of an extensive Stack/Memory managment implementation.
Speaker: Dominik Baack (TU Dortmund) -
11:10 AM
As part of this, I would also discuss "Workspace environment"
https://gitlab.ikp.kit.edu/AirShowerPhysics/corsika/issues/198Speaker: Dr Lukas Nellen (I de CIencias Nucleares, Universidad Nacional Autonoma de Mexico) -
11:55 AM
The time complexity of extensive air shower simulations rises approximately linearly with the incident particle energy for the CORSIKA 7 framework. The range of cosmic ray energies observed on earth covers several orders of magnitude. In order to simulate the highest energies in the cosmic ray spectrum, one has to introduce some sort of heuristic (e.g. thinning) which reduces runtime and preserves the shower properties to leading order. The physical content on higher order effects, like shower-to-shower fluctuations, is usually reduced. In this talk I am going to present my ideas on how to supplement current heuristics by training neural networks on CORSIKA simulations.
Speaker: Marcel Köpke (Karlsruhe Institute of Technology)
-
10:50 AM
- 12:25 PM → 12:45 PM
-
12:45 PM
→
2:10 PM
Lunch 1h 25m
-
2:10 PM
→
4:05 PM
This is a post-workshop event. Interested people should contact ralf.ulrich@kit.edu so that I can organize a room large enough for us. This session extends to Friday, and will also continue the following complete week (Max is visiting from Argentina).
Post-workshop CORSIKA discussions, 20th June 2019
Stephan, Alexander, Lukas, Felix, Max, Hans, Ralf
Stack
- First public demonstration of Stack via dynstack
~1 week by Doiminik
- Suggest to add "Particle" class, evtl. only as facade
- rename IsEmpty() -> empty(), see also below
- Particles delete:
Only set flag, physically purge particles when needed,
see https://plflib.org/colony.htm
General
- keep gcc7 as minimum now, gcc8 crystalizes as standard in a while
- look at "ranges", but don't overdo interface work
- Use std compatible naming: NOT CamelCase ! →decided
- Variable names: small →decided
- Member variables with postfix “_” → decided
- Getters: GetX() -> x() → decided
- Setters: SetX(const&) -> set_x(const&) → decided
- Method/function parameter names: no special signature ?? not clear yet
- Typenames in templates: keep the "T" prefix, not clear yet.
Optimization
Vectorisierung, is important, but requires later discussion
Profiling needed
-> time scale 6...9month to start
Configuration
Nothing about configuration has been finalized or decided yet. This is still all under discussion. Example would be extremely useful.
- Need global baseline definition of Particle cuts and thresholds (property of particle-properties??)
- Configuration: singleton class, map of map (several layers), python access,
- Init functions should disappear -> constructors are better
- Constructor -> single-config search in config maps -> fill once locally, during construction
- Config object with main functionality in c++. Not in python. Python is just the interface.
- Config object must be able to stream data in/out. Validation when? With Schema probably not; directly in c++?
- Problem with unit testing?
- Config object may count how often entries are read. If zero: might indicate typo or mistake
Preliminary Conclusions:
-
Need test implementation.
-
Python needs to handle units.
-
Need for validation procedures, need for debugging
-
Check impact on unit testing
-
Global state (singleton) is bad by default… investigate different options
-> Marcel will send configuration code example
Coordinates
Coordinate system pointers change in reference systems -> template types combined with dynamic configuration
-> Lukas will look at this
-
https://gitlab.ikp.kit.edu/AirShowerPhysics/corsika/issues?label_name%5B%5D=workshop+discussion
- First public demonstration of Stack via dynstack
-
9:00 AM
→
9:40 AM
