AbdusSalam, S. S. et al, & Eberhardt, O. (2022). Simple and statistically sound recommendations for analysing physical theories. Rep. Prog. Phys., 85(5), 052201–11pp.
Abstract: Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at Zenodo.
|
NEXT Collaboration(Haefner, J. et al), Benlloch-Rodriguez, J. M., Carcel, S., Carrion, J. V., Martin-Albo, J., Martinez-Vara, M., et al. (2023). Reflectance and fluorescence characteristics of PTFE coated with TPB at visible, UV, and VUV as a function of thickness. J. Instrum., 18(3), P03016–21pp.
Abstract: Polytetrafluoroethylene (PTFE) is an excellent diffuse reflector widely used in light collection systems for particle physics experiments. In noble element systems, it is often coated with tetraphenyl butadiene (TPB) to allow detection of vacuum ultraviolet scintillation light. In this work this dependence is investigated for PTFE coated with TPB in air for light of wavelengths of 200 nm, 260 nm, and 450 nm. The results show that TPB-coated PTFE has a reflectance of approximately 92% for thicknesses ranging from 5 mm to 10 mm at 450 nm, with negligible variation as a function of thickness within this range. A cross-check of these results using an argon chamber supports the conclusion that the change in thickness from 5 mm to 10 mm does not affect significantly the light response at 128 nm. Our results indicate that pieces of TPB-coated PTFE thinner than the typical 10 mm can be used in particle physics detectors without compromising the light signal.
|
Hirn, J., Garcia, J. E., Montesinos-Navarro, A., Sanchez-Martin, R., Sanz, V., & Verdu, M. (2022). A deep Generative Artificial Intelligence system to predict species coexistence patterns. Methods Ecol. Evol., 13, 1052–1061.
Abstract: Predicting coexistence patterns is a current challenge to understand diversity maintenance, especially in rich communities where these patterns' complexity is magnified through indirect interactions that prevent their approximation with classical experimental approaches. We explore cutting-edge Machine Learning techniques called Generative Artificial Intelligence (GenAI) to predict species coexistence patterns in vegetation patches, training generative adversarial networks (GAN) and variational AutoEncoders (VAE) that are then used to unravel some of the mechanisms behind community assemblage. The GAN accurately reproduces real patches' species composition and plant species' affinity to different soil types, and the VAE also reaches a high level of accuracy, above 99%. Using the artificially generated patches, we found that high-order interactions tend to suppress the positive effects of low-order interactions. Finally, by reconstructing successional trajectories, we could identify the pioneer species with larger potential to generate a high diversity of distinct patches in terms of species composition. Understanding the complexity of species coexistence patterns in diverse ecological communities requires new approaches beyond heuristic rules. Generative Artificial Intelligence can be a powerful tool to this end as it allows to overcome the inherent dimensionality of this challenge.
|
LHCb Collaboration(Aaij, R. et al), Martinez-Vidal, F., Oyanguren, A., Ruiz Valls, P., & Sanchez Mayordomo, C. (2014). Study of production and cold nuclear matter effects in pPb collisions at=5 TeV. J. High Energy Phys., 07(7), 094–19pp.
Abstract: Production of mesons in proton-lead collisions at a nucleon-nucleon centre-of-mass energy = 5 TeV is studied with the LHCb detector. The analysis is based on a data sample corresponding to an integrated luminosity of 1.6 nb(-1). The mesons of transverse momenta up to 15 GeV/c are reconstructed in the dimuon decay mode. The rapidity coverage in the centre-of-mass system is 1.5 < y < 4.0 (forward region) and -5.0 < y < -2.5 (backward region). The forward-backward production ratio and the nuclear modification factor for (1S) mesons are determined. The data are compatible with the predictions for a suppression of (1S) production with respect to proton-proton collisions in the forward region, and an enhancement in the backward region. The suppression is found to be smaller than in the case of prompt J/psi mesons.
|
LHCb Collaboration(Aaij, R. et al), Garcia Martin, L. M., Henry, L., Martinez-Vidal, F., Oyanguren, A., Remon Alepuz, C., et al. (2017). Study of prompt D-0 meson production in pPb collisions at root(NN)-N-s=5 TeV. J. High Energy Phys., 10(10), 090–28pp.
Abstract: Production of prompt D-0 mesons is studied in proton-lead and lead-proton collisions recorded at the LHCb detector at the LHC. The data sample corresponds to an integrated luminosity of 1.58 +/- 0.02 nb(-1) recorded at a nucleon-nucleon centre-of-mass energy of root(NN)-N-s = 5 TeV. Measurements of the differential cross-section, the forward-backward production ratio and the nuclear modification factor are reported using D-0 candidates with transverse momenta less than 10 GeV/c and rapidities in the ranges 1.5 < y* < 4.0 and -5.0 < y* < -2.5 in the nucleon-nucleon centre-of-mass system.
|
Briz, J. A., Nerio, A. N., Ballesteros, C., Borge, M. J. G., Martinez, P., Perea, A., et al. (2022). Proton Radiographs Using Position-Sensitive Silicon Detectors and High-Resolution Scintillators. IEEE Trans. Nucl. Sci., 69(4), 696–702.
Abstract: Proton therapy is a cancer treatment technique currently in growth since it offers advantages with respect to conventional X-ray and gamma-ray radiotherapy. In particular, better control of the dose deposition allowing to reach higher conformity in the treatments causing less secondary effects. However, in order to take full advantage of its potential, improvements in treatment planning and dose verification are required. A new prototype of proton computed tomography scanner is proposed to design more accurate and precise treatment plans for proton therapy. Our prototype is formed by double-sided silicon strip detectors and scintillators of LaBr3(Ce) with high energy resolution and fast response. Here, the results obtained from an experiment performed using a 100-MeV proton beam are presented. Proton radiographs of polymethyl methacrylate (PMMA) samples of 50-mm thickness with spatial patterns in aluminum were taken. Their properties were studied, including reproduction of the dimensions, spatial resolution, and sensitivity to different materials. Structures of up to 2 mm are well resolved and the sensitivity of the system was enough to distinguish the thicknesses of 10 mm of aluminum or PMMA. The spatial resolution of the images was 0.3 line pairs per mm (MTF-10%). This constitutes the first step to validate the device as a proton radiography scanner.
|
Assam, I., Vijande, J., Ballester, F., Perez-Calatayud, J., Poppe, B., & Siebert, F. A. (2022). Evaluation of dosimetric effects of metallic artifact reduction and tissue assignment on Monte Carlo dose calculations for I-125 prostate implants. Med. Phys., 49, 6195–6208.
Abstract: Purpose Monte Carlo (MC) simulation studies, aimed at evaluating the magnitude of tissue heterogeneity in I-125 prostate permanent seed implant brachytherapy (BT), customarily use clinical post-implant CT images to generate a virtual representation of a realistic patient model (virtual patient model). Metallic artifact reduction (MAR) techniques and tissue assignment schemes (TAS) are implemented on the post-implant CT images to mollify metallic artifacts due to BT seeds and to assign tissue types to the voxels corresponding to the bright seed spots and streaking artifacts, respectively. The objective of this study is to assess the combined influence of MAR and TAS on MC absorbed dose calculations in post-implant CT-based phantoms. The virtual patient models used for I-125 prostate implant MC absorbed dose calculations in this study are derived from the CT images of an external radiotherapy prostate patient without BT seeds and prostatic calcifications, thus averting the need to implement MAR and TAS. Methods The geometry of the IsoSeed I25.S17plus source is validated by comparing the MC calculated results of the TG-43 parameters for the line source approximation with the TG-43U1S2 consensus data. Four MC absorbed dose calculations are performed in two virtual patient models using the egs_brachy MC code: (1) TG-43-based D-w,w-TG(43), (2) D-w,D-w-MBDC that accounts for interseed scattering and attenuation (ISA), (3) D-m,D-m that examines ISA and tissue heterogeneity by scoring absorbed dose in tissue, and (4) D-w,D-m that unlike D-m,D-m scores absorbed dose in water. The MC absorbed doses (1) and (2) are simulated in a TG-43 patient phantom derived by assigning the densities of every voxel to 1.00 g cm(-3) (water), whereas MC absorbed doses (3) and (4) are scored in the TG-186 patient phantom generated by mapping the mass density of each voxel to tissue according to a CT calibration curve. The MC absorbed doses calculated in this study are compared with VariSeed v8.0 calculated absorbed doses. To evaluate the dosimetric effect of MAR and TAS, the MC absorbed doses of this work (independent of MAR and TAS) are compared to the MC absorbed doses of different I-125 source models from previous studies that were calculated with different MC codes using post-implant CT-based phantoms generated by implementing MAR and TAS on post-implant CT images. Results The very good agreement of TG-43 parameters of this study and the published consensus data within 3% validates the geometry of the IsoSeed I25.S17plus source. For the clinical studies, the TG-43-based calculations show a D-90 overestimation of more than 4% compared to the more realistic MC methods due to ISA and tissue composition. The results of this work generally show few discrepancies with the post-implant CT-based dosimetry studies with respect to the D-90 absorbed dose metric parameter. These discrepancies are mainly Type B uncertainties due to the different I-125 source models and MC codes. Conclusions The implementation of MAR and TAS on post-implant CT images have no dosimetric effect on the I-125 prostate MC absorbed dose calculation in post-implant CT-based phantoms.
|
Trotta, R., Johannesson, G., Moskalenko, I. V., Porter, T. A., Ruiz de Austri, R., & Strong, A. W. (2011). Constraints on Cosmic-Ray Propagation Models from a Global Bayesian Analysis. Astrophys. J., 729(2), 106–16pp.
Abstract: Research in many areas of modern physics such as, e. g., indirect searches for dark matter and particle acceleration in supernova remnant shocks rely heavily on studies of cosmic rays (CRs) and associated diffuse emissions (radio, microwave, X-rays, gamma-rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions. The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, which uses astrophysical information, and nuclear and particle data as inputs to self-consistently predict CRs, gamma-rays, synchrotron, and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.
|
ANTARES Collaboration(Adrian-Martinez, S. et al), Aguilar, J. A., Bigongiari, C., Dornic, D., Emanuele, U., Gomez-Gonzalez, J. P., et al. (2011). First Search For Point Sources Of High-Energy Cosmic Neutrinos With The Antares Neutrino Telescope. Astrophys. J. Lett., 743(1), L14–6pp.
Abstract: Results are presented of a search for cosmic sources of high-energy neutrinos with the ANTARES neutrino telescope. The data were collected during 2007 and 2008 using detector configurations containing between 5 and 12 detection lines. The integrated live time of the analyzed data is 304 days. Muon tracks are reconstructed using a likelihood-based algorithm. Studies of the detector timing indicate a median angular resolution of 0.5 +/- 0.1 deg. The neutrino flux sensitivity is 7.5 x 10(-8)(E(v)/GeV)(-2) GeV(-1) s(-1) cm(-2) for the part of the sky that is always visible (delta < -48 deg), which is better than limits obtained by previous experiments. No cosmic neutrino sources have been observed.
|
ATLAS Collaboration(Aaboud, M. et al), Alvarez Piqueras, D., Bailey, A. J., Barranco Navarro, L., Cabrera Urban, S., Castillo Gimenez, V., et al. (2018). Comparison between simulated and observed LHC beam backgrounds in the ATLAS experiment at E-beam=4 TeV. J. Instrum., 13, P12006–41pp.
Abstract: Results of dedicated Monte Carlo simulations of beam-induced background (BIB) in the ATLAS experiment at the Large Hadron Collider (LHC) are presented and compared with data recorded in 2012. During normal physics operation this background arises mainly from scattering of the 4 TeV protons on residual gas in the beam pipe. Methods of reconstructing the BIB signals in the ATLAS detector, developed and implemented in the simulation chain based on the FLUKA Monte Carlo simulation package, are described. The interaction rates are determined from the residual gas pressure distribution in the LHC ring in order to set an absolute scale on the predicted rates of BIB so that they can be compared quantitatively with data. Through these comparisons the origins of the BIB leading to different observables in the ATLAS detectors are analysed. The level of agreement between simulation results and BIB measurements by ATLAS in 2012 demonstrates that a good understanding of the origin of BIB has been reached.
|