The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10
8
GeV. Its goal is to solve the ...long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10000 radio antennas deployed over 10000 km
2
. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND.
TURTLE is a C library providing utilities allowing to navigate through a topography described by a Digital Elevation Model (DEM). The library has been primarily designed for the Monte Carlo transport ...of particles scattering over medium to long ranges, e.g. atmospheric muons. But, it can also efficiently handle ray tracing problems with very large DEMs (109 nodes or more), e.g. for neutrino simulations. The TURTLE library was built on an optimistic ray tracing algorithm, detailed in the present paper. This algorithm proceeds by trials and errors, approximating the topography within the modelling uncertainties of the DEM data. This allows to traverse a topography in constant time, i.e. independently of the number of grid nodes, and with no added memory. Detailed performance studies are provided by comparison to other ray tracing algorithms and as an application to muon transport in a Monte Carlo simulation.
Program Title: The TURTLE library
Program Files doi:http://dx.doi.org/10.17632/pbnsctnbk8.1
Licensing provisions: LGPL-3.0
Programming language: C99
Nature of problem: Transport long range Monte Carlo particles through a topography described by a high resolution Digital Elevation Model (DEM).
Solution method: Mapping the ground surface on the fly using a heuristic requiring only the original DEM data.
The PUMAS library is a transport engine for muon and tau leptons in matter. It can operate with a configurable level of details, from a fast deterministic CSDA mode to a detailed Monte Carlo ...simulation. A peculiarity of PUMAS is that it is revertible, i.e. it can run in forward or in backward mode. Thus, the PUMAS library is particularly well suited for muography applications. In the present document, we provide a detailed description of PUMAS, of its physics and of its implementation.
Program Title: The PUMAS library
CPC Library link to program files:https://doi.org/10.17632/dnry8jwkd7.1
Developer's repository link:https://niess.github.io/pumas-pages
Code Ocean capsule:https://codeocean.com/capsule/4593807
Licensing provisions: LGPL-3.0
Programming language: C99
Nature of problem: Transport of high energy muon or tau leptons in matter.
Solution method: Transport engine with a configurable level of details, from a fast deterministic CSDA mode to a detailed Monte Carlo simulation. The transport engine can operate in both forward and backward modes.
We present an algorithm for simulating reverse Monte Carlo decays given an existing forward Monte Carlo decay engine. This algorithm is implemented in the Alouette library, a TAUOLA thin wrapper for ...simulating decays of τ-leptons. We provide a detailed description of Alouette, as well as validation results.
Program title: Alouette
CPC Library link to program files:https://doi.org/10.17632/w2z5w546c5.1
Developer's repository link:https://github.com/niess/alouette
Code Ocean capsule:https://codeocean.com/capsule/9167021
Licensing provisions: LGPL-3.0
Programming language: C, Fortran and Python
Nature of problem: Perform reverse Monte Carlo decays.
Solution method: Invert an existing forward Monte Carlo engine using the Jacobian backward method. Apply the algorithm to τ decays generated by TAUOLA.
SUMMARY
Gravimetry is a technique widely used to image the structure of the Earth. However, inversions are ill-posed and the imaging power of the technique rapidly decreases with depth. To overcome ...this limitation, muography, a new imaging technique relying on high energy atmospheric muons, has recently been developed. Because muography only provides integrated densities above the detector from a limited number of observation points, inversions are also ill-posed. Previous studies have shown that joint muographic and gravimetric inversions better reconstruct the 3-D density structure of volcanic edifices than independent density inversions. These studies address the ill-posedness of the joint problem by regularizing the solution with respect to a prior density model. However, the obtained solutions depend on some hyperparameters, which are either determined relative to a single test case or rely on ad-hoc parameters. This can lead to inaccurate retrieved models, sometimes associated with artefacts linked to the muon data acquisition. In this study, we use a synthetic example based on the Puy de Dôme volcano to determine a robust method to obtain the resulting model closest to the synthetic model and devoid of acquisition artefacts. We choose a Bayesian approach to include an a priori density model and a smoothing by a Gaussian spatial correlation function relying on two hyperparameters: an a priori density standard deviation and an isotropic spatial correlation length. This approach has the advantage to provide a posteriori standard deviations on the resulting densities. Using our synthetic volcano, we investigate the most reliable criterion to determine the hyperparameters. Our results suggest that k-fold Cross-Validation Sum of Squares and the Leave One Out methods are more robust criteria than the classically used L-curves. The determined hyperparameters allow to overcome the artefacts linked to the data acquisition geometry, even when only a limited number of muon telescopes is available. We also illustrate the behaviour of the inversion in case of offsets in the a priori density or in the data and show that they lead to recognizable structures that help identify them.
The determination of the fundamental parameters of the Standard Model (and its extensions) is often limited by the presence of statistical and theoretical uncertainties. We present several models for ...the latter uncertainties (random, nuisance, external) in the frequentist framework, and we derive the corresponding
p
values. In the case of the nuisance approach where theoretical uncertainties are modeled as biases, we highlight the important, but arbitrary, issue of the range of variation chosen for the bias parameters. We introduce the concept of adaptive
p
value, which is obtained by adjusting the range of variation for the bias according to the significance considered, and which allows us to tackle metrology and exclusion tests with a single and well-defined unified tool, which exhibits interesting frequentist properties. We discuss how the determination of fundamental parameters is impacted by the model chosen for theoretical uncertainties, illustrating several issues with examples from quark flavor physics.
The Giant Radio Array for Neutrino Detection Martineau-Huynh, Olivier; Kotera, Kumiko; Bustamente, Mauricio ...
EPJ Web of Conferences,
01/2016, Letnik:
116
Journal Article, Conference Proceeding
Recenzirano
Odprti dostop
High-energy neutrino astronomy will probe the working of the most violent phenomena in the Universe. The Giant Radio Array for Neutrino Detection (GRAND) project consists of an array of ∼ 105 radio ...antennas deployed over ∼ 200 000 km2 in a mountainous site. It aims at detecting high-energy neutrinos via the measurement of air showers induced by the decay in the atmosphere of τ leptons produced by the interaction of cosmic neutrinos under the Earth surface. Our objective with GRAND is to reach a neutrino sensitivity of 5 × 10−11E−2 GeV−1 cm−2 s−1 sr−1 above 3 × 1016 eV. This sensitivity ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and up to 100 events per year are expected for the standard models. GRAND would also probe the neutrino signals produced at the potential sources of UHECRs.
The Giant Radio Array for Neutrino Detection Martineau-Huynh, Olivier; Bustamante, Mauricio; Carvalho, Washington ...
EPJ Web of Conferences,
01/2017, Letnik:
135
Journal Article, Conference Proceeding
Recenzirano
Odprti dostop
The Giant Radio Array for Neutrino Detection (GRAND) is a planned array of ~ 2·105 radio antennas deployed over ~ 200 000 km2 in a mountainous site. It aims primarly at detecting high-energy ...neutrinos via the observation of extensive air showers induced by the decay in the atmosphere of taus produced by the interaction of cosmic neutrinos under the Earth surface. GRAND aims at reaching a neutrino sensitivity of 5 · 10−11 E−2 GeV−1 cm−2 s−1 sr−1 above 3 · 1016 eV. This ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and ~50 events per year are expected for the standard models. The instrument will also detect UHECRs and possibly FRBs. Here we show how our preliminary design should enable us to reach our sensitivity goals, and discuss the steps to be taken to achieve GRAND, while the compelling science case for GRAND is discussed in more details in 1.
Over the last decades, radio detection of air showers has been established as a promising detection technique for ultrahigh-energy cosmic rays and neutrinos. Very large or dense antenna arrays are ...necessary to be proficient at collecting information about these particles and understanding their properties accurately. The exploitation of such arrays requires to run massive air-shower simulations to evaluate the radio signal at each antenna position, taking into account features such as the ground topology. In order to reduce computational costs, we have developed a fast computation of the emitted radio signal on the basis of generic shower simulations, called Radio Morphing. The method consists in the calculation of the radio signal of any air-shower by i) a scaling of the electric-field amplitude of a reference air shower to the target shower, ii) an isometry on the simulated positions and iii) an interpolation of the radio pulse at the desired position. This technique enables one to gain many orders of magnitude in CPU time compared to a standard computation. In this contribution, we present this novel tool and explain its methodology. In particular, Radio Morphing will be a key element for the simulation chain of the Giant Radio Array for Neutrino Detection (GRAND) project, that aims at detecting ultra-high-energy neutrinos with an array of 200 000 radio antennas in mountainous regions.
Isospin analysis of charmless B-meson decays Charles, J.; Deschamps, O.; Descotes-Genon, S. ...
The European physical journal. C, Particles and fields,
08/2017, Letnik:
77, Številka:
8
Journal Article
Recenzirano
Odprti dostop
We discuss the determination of the CKM angle
α
using the non-leptonic two-body decays
B
→
π
π
,
B
→
ρ
ρ
and
B
→
ρ
π
using the latest data available. We illustrate the methods used in each case and ...extract the corresponding value of
α
. Combining all these elements, we obtain the determination
α
dir
=
(
86.2
-
4.0
+
4.4
∪
178.4
-
5.1
+
3.9
)
∘
. We assess the uncertainties associated to the breakdown of the isospin hypothesis and the choice of the statistical framework in detail. We also determine the hadronic amplitudes (tree and penguin) describing the QCD dynamics involved in these decays, briefly comparing our results with theoretical expectations. For each observable of interest in the
B
→
π
π
,
B
→
ρ
ρ
and
B
→
ρ
π
systems, we perform an indirect determination based on the constraints from all the other observables available and we discuss the compatibility between indirect and direct determinations. Finally, we review the impact of future improved measurements on the determination of
α
.