Source author record

Martin Erdmann

Martin Erdmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Fast Columnar Physics Analyses of Terabyte-Scale LHC Data on a Cache-Aware Dask Cluster

The development of an LHC physics analysis involves numerous investigations that require the repeated processing of terabytes of data. Thus, a rapid completion of each of these analysis cycles is central to mastering the science project. We present a solution to efficiently handle and accelerate physics analyses on small-size institute clusters. Our solution is based on three key concepts: Vectorized processing of collision events, the "MapReduce" paradigm for scaling out on computing clusters, and efficiently utilized SSD caching to reduce latencies in IO operations. Using simulations from a Higgs pair production physics analysis as an example, we achieve an improvement factor of $6.3$ in runtime after one cycle and even an overall speedup of a factor of $14.9$ after $10$ cycles.

preprint2021arXiv

Inference of cosmic-ray source properties by conditional invertible neural networks

The inference of physical parameters from measured distributions constitutes a core task in physics data analyses. Among recent deep learning methods, so-called conditional invertible neural networks provide an elegant approach owing to their probability-preserving bijective mapping properties. They enable training the parameter-observation correspondence in one mapping direction and evaluating the parameter posterior distributions in the reverse direction. Here, we study the inference of cosmic-ray source properties from cosmic-ray observations on Earth using extensive astrophysical simulations. We compare the performance of conditional invertible neural networks (cINNs) with the frequently used Markov Chain Monte Carlo (MCMC) method. While cINNs are trained to directly predict the parameters' posterior distributions, the MCMC method extracts the posterior distributions through a likelihood function that matches simulations with observations. Overall, we find good agreement between the physics parameters derived by the two different methods. As a result of its computational efficiency, the cINN method allows for a swift assessment of inference quality.

preprint2021arXiv

Towards extracting cosmic magnetic field structures from cosmic-ray arrival directions

We present a novel method to search for structures of coherently aligned patterns in ultra-high energy cosmic-ray arrival directions simultaneously across the entire sky. This method can be used to obtain information on the Galactic magnetic field, in particular the integrated component perpendicular to the line of sight, from cosmic-ray data only. Using a likelihood-ratio approach, neighboring cosmic rays are related by rotatable, elliptically shaped density distributions and the significance of their alignment with respect to circular distributions is evaluated. In this way, a vector field tangential to the celestial sphere is fitted which approximates the local deflections in cosmic magnetic fields if significant deflection structures are detected. The sensitivity of the method is evaluated on the basis of astrophysical simulations of the ultra-high energy cosmic-ray sky, where a discriminative power between isotropic and signal-induced scenarios is found.

preprint2020arXiv

Precise simulation of electromagnetic calorimeter showers using a Wasserstein Generative Adversarial Network

Simulations of particle showers in calorimeters are computationally time-consuming, as they have to reproduce both energy depositions and their considerable fluctuations. A new approach to ultra-fast simulations are generative models where all calorimeter energy depositions are generated simultaneously. We use GEANT4 simulations of an electron beam impinging on a multi-layer electromagnetic calorimeter for adversarial training of a generator network and a critic network guided by the Wasserstein distance. The generator is constraint during the training such that the generated showers show the expected dependency on the initial energy and the impact position. It produces realistic calorimeter energy depositions, fluctuations and correlations which we demonstrate in distributions of typical calorimeter observables. In most aspects, we observe that generated calorimeter showers reach the level of showers as simulated with the GEANT4 program.

preprint2016arXiv

CRPropa 3 - a Public Astrophysical Simulation Framework for Propagating Extraterrestrial Ultra-High Energy Particles

We present the simulation framework CRPropa version 3 designed for efficient development of astrophysical predictions for ultra-high energy particles. Users can assemble modules of the most relevant propagation effects in galactic and extragalactic space, include their own physics modules with new features, and receive on output primary and secondary cosmic messengers including nuclei, neutrinos and photons. In extension to the propagation physics contained in a previous CRPropa version, the new version facilitates high-performance computing and comprises new physical features such as an interface for galactic propagation using lensing techniques, an improved photonuclear interaction calculation, and propagation in time dependent environments to take into account cosmic evolution effects in anisotropy studies and variable sources. First applications using highlighted features are presented as well.

preprint2016arXiv

Energy Estimation of Cosmic Rays with the Engineering Radio Array of the Pierre Auger Observatory

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

preprint2016arXiv

Measurement of the Radiation Energy in the Radio Signal of Extensive Air Showers as a Universal Estimator of Cosmic-Ray Energy

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

preprint2016arXiv

Simulation of Radiation Energy Release in Air Showers

A simulation study of the energy released by extensive air showers in the form of MHz radiation is performed using the CoREAS simulation code. We develop an efficient method to extract this radiation energy from air-shower simulations. We determine the longitudinal profile of the radiation energy release and compare it to the longitudinal profile of the energy deposit by the electromagnetic component of the air shower. We find that the radiation energy corrected for the geometric dependence of the geomagnetic emission scales quadratically with the energy in the electromagnetic component of the air shower with a second-order dependence on the atmospheric density at the position of the maximum shower development $X_\mathrm{max}$. In a measurement where $X_\mathrm{max}$ is not accessible, this second order dependence can be approximated using the zenith angle of the incoming direction of the air shower with only a minor loss in accuracy. Our method results in an intrinsic uncertainty of 4% in the determination of the energy in the electromagnetic air-shower component, which is well below current experimental uncertainties.

preprint2016arXiv

The Nuclear Window to the Extragalactic Universe

We investigate two recent parameterizations of the galactic magnetic field with respect to their impact on cosmic nuclei traversing the field. We present a comprehensive study of the size of angular deflections, dispersion in the arrival probability distributions, multiplicity in the images of arrival on Earth, variance in field transparency, and influence of the turbulent field components. To remain restricted to ballistic deflections, a cosmic nucleus with energy E and charge Z should have a rigidity above E/Z=6 EV. In view of the differences resulting from the two field parameterizations as a measure of current knowledge in the galactic field, this rigidity threshold may have to be increased. For a point source search with E/Z>60 EV, field uncertainties increase the required signal events for discovery moderately for sources in the northern and southern regions, but substantially for sources near the galactic disk.

preprint2015arXiv

An improved limit to the diffuse flux of ultra-high energy neutrinos from the Pierre Auger Observatory

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_ν^2 dN_ν/dE_ν< 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.

preprint2015arXiv

Investigation of the Galactic Magnetic Field with Ultra-High Energy Cosmic Rays

We present a method to correct for deflections of ultra-high energy cosmic rays in the galactic magnetic field. We perform these corrections by simulating the expected arrival directions of protons using a parameterization of the field derived from Faraday rotation and synchrotron emission measurements. To evaluate the method we introduce a simulated astrophysical scenario and two observables designed for testing cosmic ray deflections. We show that protons can be identified by taking advantage of the galactic magnetic field pattern. Consequently, cosmic ray deflection in the galactic field can be verified experimentally. The method also enables searches for directional correlations of cosmic rays with source candidates.

preprint2015arXiv

Measurement of the cosmic ray spectrum above $4{\times}10^{18}$ eV using inclined events detected with the Pierre Auger Observatory

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-γ}$ with index $γ=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.

preprint2014arXiv

A field study of data analysis exercises in a bachelor physics course using the internet platform VISPA

Bachelor physics lectures on particle physics and astrophysics were complemented by exercises related to data analysis and data interpretation at the RWTH Aachen University recently. The students performed these exercises using the internet platform VISPA, which provides a development environment for physics data analyses. We describe the platform and its application within the physics course, and present the results of a student survey. The students acceptance of the learning project was positive. The level of acceptance was related to their individual preference for learning with a computer. Furthermore, students with good programming skills favor working individually, while students who attribute themselves having low programming abilities favor working in teams. The students appreciated approaching actual research through the data analysis tasks.

preprint2014arXiv

Cosmic ray propagation with CRPropa 3

Solving the question of the origin of ultra-high energy cosmic rays (UHECRs) requires the development of detailed simulation tools in order to interpret the experimental data and draw conclusions on the UHECR universe. CRPropa is a public Monte Carlo code for the galactic and extragalactic propagation of cosmic ray nuclei above $\sim 10^{17}$ eV, as well as their photon and neutrino secondaries. In this contribution the new algorithms and features of CRPropa 3, the next major release, are presented. CRPropa 3 introduces time-dependent scenarios to include cosmic evolution in the presence of cosmic ray deflections in magnetic fields. The usage of high resolution magnetic fields is facilitated by shared memory parallelism, modulated fields and fields with heterogeneous resolution. Galactic propagation is enabled through the implementation of galactic magnetic field models, as well as an efficient forward propagation technique through transformation matrices. To make use of the large Python ecosystem in astrophysics CRPropa 3 can be steered and extended in Python.

preprint2013arXiv

CRPropa 3.0 - a Public Framework for Propagating UHE Cosmic Rays through Galactic and Extragalactic Space

The interpretation of experimental data of ultra-high energy cosmic rays (UHECRs) above 10^17 eV is still under controversial debate. The development and improvement of numerical tools to propagate UHECRs in galactic and extragalactic space is a crucial ingredient to interpret data and to draw conclusions on astrophysical parameters. In this contribution the next major release of the publicly available code CRPropa (3.0) is presented. It reflects a complete redesign of the code structure to facilitate high performance computing and comprises new physical features such as an interface for galactic propagation using lensing techniques and inclusion of cosmological effects in a three-dimensional environment. The performance is benchmarked and first applications are presented.

preprint2013arXiv

Detecting Local Deflection Patterns of Ultra-high Energy Cosmic Rays using the Principal Axes of the Directional Energy Distribution

From deflections in galactic and extragalactic magnetic fields energy dependent structures in the arrival directions of ultra-high energy cosmic rays (UHECR) are expected. We propose to characterize these structures by the strength of collimation of energy along the principal axes in selected regions in the sky. While the strength of collimation are indicators of anisotropy in the arrival distribution of UHECR, the orientation of the principal system holds information about the direction of the deflections of UHECR. We discuss the method and present expected limits on the strength of deflection and density of sources using simulated scenarios of UHECR proton propagation.

preprint2013arXiv

PARSEC: A Parametrized Simulation Engine for Ultra-High Energy Cosmic Ray Protons

We present a new simulation engine for fast generation of ultra-high energy cosmic ray data based on parametrizations of common assumptions of UHECR origin and propagation. Implemented are deflections in unstructured turbulent extragalactic fields, energy losses for protons due to photo-pion production and electron-pair production, as well as effects from the expansion of the universe. Additionally, a simple model to estimate propagation effects from iron nuclei is included. Deflections in galactic magnetic fields are included using a matrix approach with precalculated lenses generated from backtracked cosmic rays. The PARSEC program is based on object oriented programming paradigms enabling users to extend the implemented models and is steerable with a graphical user interface.

preprint2012arXiv

A new Monte Carlo Generator for Ultra-High Energy Cosmic Rays from the Local and Distant Universe

For the understanding of the origin and propagation of ultra-high energy cosmic rays (UHECR) we developed a new approach to simulating UHECRs from an arbitrary number of sources based on Monte Carlo technique. The method consists of a combination of three steps. For distant sources we apply commonly accepted parameterizations to calculate the contribution to the observed cosmic ray flux. For sources of the local universe we use forward tracking through realistic matter distributions and magnetic fields resulting from explicit simulations of large-scale structure formation. From the calculations and the forward tracking we generate maps of the probability to observe a particle with a given energy from a discrete direction. To account for deflections in the galactic field, these probability maps are transformed by matrices calculated from backtracking of antiparticles through field parameterizations. Based on the combined probability maps, Monte Carlo production of individual UHECR data is performed which are then used in comparisons with UHECR measurements. The simulated UHECR data serves for optimizing the analysis techniques used in UHECR measurements as well as for constraining the parameter space of the underlying source and magnetic field models.

preprint2010arXiv

A Method of Measuring Cosmic Magnetic Fields with Ultra High Energy Cosmic Ray Data

We present a method to measure cosmic magnetic fields with ultra high energy cosmic rays (UHECRs). We apply an advanced autocorrelation method to simulated UHECRs which includes their directional as well as energy information. Without explicit knowledge of the UHECR sources, such measurements are sensitive to the number of sources and to the magnetic field strength subjected to the UHECRs. Using a UHECR Monte Carlo model including sources, random walk propagation and a coherent deflection, we explain the procedure of reconstructing the allowed phase space of the model parameters from a simulated autocorrelation measurement.