CORSIKA Cosmic Ray Simulation Workshop Karlsruhe
Hotel ACHAT PLAZA Karlsruhe
General and open workshop to discuss technology and science of air shower simulations, and meeting of the CORSIKA 8 Project.
Abdul Rehman
Alan Coleman
Alexander Sandrock
Anna Schäfer
Bernardino Spisso
Blahoslav Pastirčák
Darko Veberic
Dieter Heck
Dmitriy Kostunin
Dominik Baack
Doris Wochele
Eva Santos
felix riehn
Frank Schröder
G. Schatz
Gabor Veres
Hans Dembinski
Ioana Alexandra Caracas
Jan Ebr
Jean-Marco Alameddine
Jean-Noel Capdevielle
Johannes Blümer
Johannes Knapp
Juergen Oehlschlaeger
Konrad Bernlöhr
Luisa Arrabito
Lukas Nellen
Marcel Koepke
Mario Pecimotika
Markus Roth
Maximilian Reininghaus
Michael Unger
Oliver Suranyi
Peter Grieder
Piotr Kalaczyński
Ralf Ulrich
Ralph Engel
Roxanne Turcotte
Sabine Bucher
Sergey Ostapchenko
Simona Maria Stellacci
Stephan Meighen-Berger
Suren Chilingaryan
Takashi Sako
Tanguy Pierog
Tim Huege
Yating Chai
-
-
1
Arrival, Registration, CoffeeSpeaker: Organizers
-
10:30 AM
Break
- 2
-
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
IntroductionSpeaker: Ralph Engel
-
4
The CORSIKA 8 ProjectSpeaker: Ralf Ulrich (KIT)
-
5
Status of capabilities of CORSIKA 8, first results and comparisonsSpeaker: Maximilian Reininghaus (KIT / IKP)
-
6
The lepton propagator PROPOSAL and the air shower simulation code CORSIKA
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)
-
3
-
12:30 PM
Lunch
-
Applications and requirements of CORSIKA
Contributions from experimental collaborations, phenomenologist or theorists about important application of CORSIKA, and future wishes for new functionality.
-
7
Measured and simulated muon flux for the first KM3NeT Detection Units
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) -
8
Corsika simulations of background in Baikal experiment
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) -
9
Proton-Oxygen collisions at the LHC for air shower research
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)
-
7
-
3:00 PM
Coffee
-
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
The COSMOS air shower simulation programSpeaker: Takashi Sako (ICRR, University of Tokyo)
-
11
Hands-On: Simulate your first CORSIKA 8 air shower
Tutorial part (DIY) and introduction, overview, results
Speaker: Ralf Ulrich (KIT) -
12
Introduction to template meta-programming in C++
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)
-
10
-
1
-
-
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.
-
13
IntroductionSpeaker: Johannes Blümer (KIT)
-
14
Original motivation for CORSIKASpeaker: Gerd Schatz
-
15
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
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 $N\overline{N}$ in air showers was discovered, which prompted me to develop what later on became the air shower simulation program system ASICO. I discuss briefly the contemporary shower and high energy hadronic interaction models, the sources of information to construct such models, the architecture of the initial cascade simulation process, and the first results. This is followed by an outline of the unique Echo Lake experiment, a landmark in cosmic ray physics, its results and the problem it caused. Subsequently, a series of significant discoveries and new theoretical developments of the epoch are summarized, which culminated in the coming into operation of the CERN ISR. This collider with its approximately 1800\,GeV laboratory equivalent energy yielded an enormous wealth of new information and insight into a hitherto "uncharted" energy domain. Its results confirmed several cosmic ray discoveries and solved the Echo Lake problem. The presentation closes with a potpourri of results of ASICO after over twenty years of development, and its transition to CORSIKA.
Speaker: Prof. Peter Grieder (University of Bern, Switzerland) - 16
-
10:25 AM
Coffee
-
17
The link to particle physicsSpeaker: Johannes Knapp
- 18
-
19
CORSIKA and Imaging Air-Cherenkov TelescopesSpeaker: Konrad Bernlöhr
- 20
-
13
-
12:45 PM
Lunch
-
Applications and requirements of CORSIKA
Contributions from experimental collaborations, phenomenologist or theorists about important application of CORSIKA, and future wishes for new functionality.
-
21
Modeling of high energy cosmic ray interactions: selected topicsSpeaker: Sergey Ostapchenko
-
22
SMASH - a novel transport model to simulate low-energy hadronic interactionsSpeaker: Anna Schäfer
-
3:30 PM
Coffee
-
23
Opportunities with Calorimeters at LHC: How to solve the muon puzzle?
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) - 24
-
25
Results from the LHCf calorimeterSpeaker: Takashi Sako (ICRR, University of Tokyo)
- 26
-
27
GHOST - a new generator for nuclear collisionsSpeaker: Jean-Noel Capdevielle
-
21
-
-
-
Applications and requirements of CORSIKA
Contributions from experimental collaborations, phenomenologist or theorists about important application of CORSIKA, and future wishes for new functionality.
- 28
-
29
Use of CORSIKA for the IceTop Radio Enhancement
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) -
30
Handling the metadata for large amount of CORSIKA simulations
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) -
31
CORELib: COsmic Ray Event Library
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
Coffee
-
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.
-
32
Particle propagation for CORSIKA in PROPOSAL
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 -
33
EmCa - Electromagnetic Cascades Simulation Package
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 -
34
Hands-On: Physics processes
Tutorial part for DIY, introduction and concepts
Speaker: felix riehn (LIP, Lisbon)
-
32
-
12:35 PM
Lunch
-
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.
-
35
Python in High-Energy Physics
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) - 36
- 37
-
38
Simulation of radio emission from particle showers in CORSIKA 7 and CORSIKA 8Speaker: Tim Huege (KIT)
-
3:25 PM
Coffee
-
35
-
6:00 PM
Workshop Dinner 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.scc.kit.edu/event/529/page/475-circular-5th-june
-
-
-
Applications and requirements of CORSIKA
Contributions from experimental collaborations, phenomenologist or theorists about important application of CORSIKA, and future wishes for new functionality.
-
39
Simulations for the atmosphere with clouds
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) - 40
-
39
-
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.
-
41
Hands-On: Geometry and environment setupSpeaker: Maximilian Reininghaus (KIT / IKP)
-
41
-
10:20 AM
Coffee
-
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.
-
42
DynStack / Memory managment - implementation and features
Overview of different features and possibilites of an extensive Stack/Memory managment implementation.
Speaker: Dominik Baack (TU Dortmund) -
43
Run control and steering
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) -
44
Speedup of extensive air shower simulations with neural networks
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)
-
42
-
45
Conclusion and end of workshopSpeaker: Ralf Ulrich (KIT)
-
12:45 PM
Lunch
-
46
Working session on CORSIKA 8 (optional, post-workshop)
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
-
a) Open issues
https://gitlab.ikp.kit.edu/AirShowerPhysics/corsika/issues?label_name%5B%5D=workshop+discussion
-
b) Further planning and tasks
- First public demonstration of Stack via dynstack
-