Unfolding spectral analysis of the Fermi-LAT data

The Large Area Telescope (LAT) onboard the Fermi satellite is observing the gamma-ray sky in the high energy region, above 20 MeV. We have developed a method to reconstruct the energy spectra of the gamma-rays detected by the Fermi LAT instrument based on a Bayesian unfolding approach, that takes into account the energy dispersion introduced by the instrument response. The method has been successfully applied to reconstruct the energy spectra of both steady and pulsating point sources. The analysis technique will be illustrated and the results obtained in some significant test cases will be discussed.Comment: 6 pages, 9 figure

Hadronic interactions of primary cosmic rays with the FLUKA code

The measured fluxes of secondary particles produced by the interactions of cosmic rays with the astronomical environment represent a powerful tool to infer some properties of primary cosmic rays. In this work we investigate the production of secondary particles in inelastic hadronic interactions between several cosmic rays species of projectiles and different target nuclei of the interstellar medium. The yields of secondary particles have been calculated with the FLUKA simulation package, that provides with very good accuracy the energy distributions of secondary products in a large energy range. An application to the propagation and production of secondaries in the Galaxy is presented.Comment: 8 pages, 4 figures; Contribution to the 34th International Cosmic Ray Conference, July 30 to August 6, The Hague, Netherlands; fixing a typo in the y-axis label of Fig.

A Bayesian approach to evaluate confidence intervals in counting experiments with background

In this paper we propose a procedure to evaluate Bayesian confidence intervals in counting experiments where both signal and background fluctuations are described by the Poisson statistics. The results obtained when the method is applied to the calculation of upper limits will also be illustrated.Comment: 19 pages, 6 figure

Markov chain Monte Carlo analyses of the flux ratios of B, Be and Li with the DRAGON2 code

Recent cosmic-ray measurements are challenging our models of propagation in the Galaxy. A good characterization of the secondary cosmic rays (B, Be, Li and sub-iron species) is crucial to constrain these models and exploit the precision of modern CR experiments. In this work, a Markov chain Monte Carlo analysis has been implemented to fit the experimental flux ratios between B, Be and Li and their flux ratios to the primary nuclei C and O. We have fitted the data using two different parametrizations for the spallation cross sections. The uncertainties in the evaluation of the spectra of these secondary cosmic rays, due to spallation cross sections, have been considered by introducing scale factors as nuisance parameters. We have also tested two different formulations for the diffusion coefficient, which differ in the origin of the high energy hardening of cosmic rays. Additionally, two different approaches are used to scale the cross sections, one based on a combined analysis of all the species ("combined" analysis) and the other reproducing the high energy spectra of the secondary-to-secondary flux ratios of Be/B, Li/B, Li/Be ("scaled" analysis). This allows us to make a better comparison between the propagation parameters inferred from the cross sections parametrizations tested in this work. This novel analysis has been successfully implemented using the numerical code DRAGON2 to reproduce the cosmic-ray nuclei data up to $Z=14$ from the AMS-02 experiment. It is found that the ratios of Li favor a harder spectral index of the diffusion coefficient, but compatible with the other ratios inside the observed $2\sigma$ uncertainties. In addition, it is shown that, including these scale factors, the secondary-to-primary flux ratios can be simultaneously reproduced.Comment: 27 pages, 10 figures, 9 table

FLUKA cross sections for cosmic-ray interactions with the DRAGON2 code

Secondary particles produced in spallation reactions of cosmic rays with the interstellar gas provide valuable information that allow us to investigate the injection and transport of charged particles in the Galaxy. A good understanding of the cross sections of production of these particles is crucial to correctly interpret our models, although the existing experimental data is very scarce and uncertain. We have developed a new set of cross sections, both inelastic and inclusive, computed with the FLUKA Monte Carlo nuclear code and tested its compatibility with CR data. Inelastic and inclusive cross sections have been compared to the most up-to-date data and parameterisations finding a general good agreement. Then, these cross sections have been implemented in the DRAGON2 code to characterize the spectra of CR nuclei up to Z = 26 and the secondary-to-primary ratios of B, Be and Li. Interestingly, we find that the FLUKA cross sections allow us to predict an energy-dependence of the B, Be and Li flux ratios which is compatible with AMS-02 data and to reproduce simultaneously these flux ratios with a scaling lower than 20%. Finally, we implement the cross sections of production of gamma rays, calculated with FLUKA, in the Gammasky code and compute diffuse gamma-ray sky maps and the local HI emissivity spectrum, finding a very good agreement with Fermi Large Area Telescope data