ATLAS Collaboration(Abat, E. et al), Castillo Gimenez, V., Ferrer, A., Gonzalez, V., Higon-Rodriguez, E., Mitsou, V. A., et al. (2010). Study of energy response and resolution of the ATLAS barrel calorimeter to hadrons of energies from 20 to 350 GeV. Nucl. Instrum. Methods Phys. Res. A, 621(1-3), 134–150.
Abstract: A fully instrumented slice of the ATLAS detector was exposed to test beams from the SPS (Super Proton Synchrotron) at CERN in 2004. In this paper, the results of the measurements of the response of the barrel calorimeter to hadrons with energies in the range 20-350 GeV and beam impact points and angles corresponding to pseudo-rapidity values in the range 0.2-0.65 are reported. The results are compared to the predictions of a simulation program using the Geant 4 toolkit.
|
Agramunt, J. et al, Tain, J. L., Albiol, F., Algora, A., Domingo-Pardo, C., Jordan, M. D., et al. (2016). Characterization of a neutron-beta counting system with beta-delayed neutron emitters. Nucl. Instrum. Methods Phys. Res. A, 807, 69–78.
Abstract: A new detection system for the measurement of beta-delayed neutron emission probabilities has been characterized using fission products with well known beta-delayed neutron emission properties. The setup consists of BELEN-20, a 4 pi-neutron counter with twenty He-3 proportional tubes arranged inside a large polyethylene neutron moderator, a thin Si detector for beta counting and a self-triggering digital data acquisition system. The use of delayed-neutron precursors with different neutron emission windows allowed the study of the effect of energy dependency on neutron, beta and beta-neutron rates. The observed effect is well reproduced by Monte Carlo simulations. The impact of this dependency on the accuracy of neutron emission probabilities is discussed. A new accurate value of the neutron emission probability for the important delayed-neutron precursor I-137 was obtained, P-n = 7.76(14)%.
|
ATLAS Collaboration(Abat, E. et al), Bernabeu Verdu, J., Castillo Gimenez, V., Costa, M. J., Escobar, C., Ferrer, A., et al. (2011). A layer correlation technique for pion energy calibration at the 2004 ATLAS Combined Beam Test. J. Instrum., 6, P06001–35pp.
Abstract: A new method for calibrating the hadron response of a segmented calorimeter is developed and successfully applied to beam test data. It is based on a principal component analysis of energy deposits in the calorimeter layers, exploiting longitudinal shower development information to improve the measured energy resolution. Corrections for invisible hadronic energy and energy lost in dead material in front of and between the calorimeters of the ATLAS experiment were calculated with simulated Geant4 Monte Carlo events and used to reconstruct the energy of pions impinging on the calorimeters during the 2004 Barrel Combined Beam Test at the CERN H8 area. For pion beams with energies between 20 GeV and 180 GeV, the particle energy is reconstructed within 3% and the energy resolution is improved by between 11% and 25% compared to the resolution at the electromagnetic scale.
|
Jay, G., Arnault, P., & Debbasch, F. (2021). Dirac quantum walks with conserved angular momentum. Quantum Stud. Math. Found., 8, 419–430.
Abstract: A quantum walk (QW) simulating the flat (1+2)D Dirac equation on a spatial polar grid is constructed. Because fermions are represented by spinors, which do not constitute a representation of the rotation group SO(3), but rather of its double cover SU(2), the QW can only be defined globally on an extended spacetime where the polar angle extends from 0 to 4 pi. The coupling of the QW with arbitrary electromagnetic fields is also presented. Finally, the cylindrical relativistic Landau levels of the Dirac equation are computed explicitly and simulated by the QW.
|
Jaworski, G., Palacz, M., Nyberg, J., de Angelis, G., de France, G., Di Nitto, A., et al. (2012). Monte Carlo simulation of a single detector unit for the neutron detector array NEDA. Nucl. Instrum. Methods Phys. Res. A, 673, 64–72.
Abstract: A study of the dimensions and performance of a single detector of the future neutron detector array NEDA was performed by means of Monte Carlo simulations, using GEANT4. Two different liquid scintillators were evaluated: the hydrogen based BC501A and the deuterated BC537. The efficiency and the probability that one neutron will trigger a signal in more than one detector were investigated as a function of the detector size. The simulations were validated comparing the results to experimental measurements performed with two existing neutron detectors, with different geometries, based on the liquid scintillator BC501.
|
Asai, M., Cortes-Giraldo, M. A., Gimenez-Alventosa, V., Gimenez, V., & Salvat, F. (2021). The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4. Front. Physics, 9, 738735–20pp.
Abstract: A translation of the penelope physics subroutines to C++, designed as an extension of the Geant4 toolkit, is presented. The Fortran code system penelope performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. Penelope implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of Geant4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original penelope programs. The combined code, named PenG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of Geant4.
|
de Putter, R., Wagner, C., Mena, O., Verde, L., & Percival, W. J. (2012). Thinking outside the box: effects of modes larger than the survey on matter power spectrum covariance. J. Cosmol. Astropart. Phys., 04(4), 019–31pp.
Abstract: Accurate power spectrum (or correlation function) covariance matrices are a crucial requirement for cosmological parameter estimation from large scale structure surveys. In order to minimize reliance on computationally expensive mock catalogs, it is important to have a solid analytic understanding of the different components that make up a covariance matrix. Considering the matter power spectrum covariance matrix, it has recently been found that there is a potentially dominant effect on mildly non-linear scales due to power in modes of size equal to and larger than the survey volume. This beat coupling effect has been derived analytically in perturbation theory and while it has been tested with simulations, some questions remain unanswered. Moreover, there is an additional effect of these large modes, which has so far not been included in analytic studies, namely the effect on the estimated average density which enters the power spectrum estimate. In this article, we work out analytic, perturbation theory based expressions including both the beat coupling and this local average effect and we show that while, when isolated, beat coupling indeed causes large excess covariance in agreement with the literature, in a realistic scenario this is compensated almost entirely by the local average effect, leaving only similar to 10% of the excess. We test our analytic expressions by comparison to a suite of large N-body simulations, using both full simulation boxes and subboxes thereof to study cases without beat coupling, with beat coupling and with both beat coupling and the local average effect. For the variances, we find excellent agreement with the analytic expressions for k < 0.2 hMpc(-1) at z = 0.5, while the correlation coefficients agree to beyond k = 0.4 hMpc(-1). As expected, the range of agreement increases towards higher redshift and decreases slightly towards z = 0. We finish by including the large-mode effects in a full covariance matrix description for arbitrary survey geometry and confirming its validity using simulations. This may be useful as a stepping stone towards building an actual galaxy (or other tracer's) power spectrum covariance matrix.
|
Bernal, N., Forero-Romero, J. E., Garani, R., & Palomares-Ruiz, S. (2014). Systematic uncertainties from halo asphericity in dark matter searches. J. Cosmol. Astropart. Phys., 09(9), 004–30pp.
Abstract: Although commonly assumed to be spherical, dark matter halos are predicted to be non-spherical by N-body simulations and their asphericity has a potential impact on the systematic uncertainties in dark matter searches. The evaluation of these uncertainties is the main aim of this work, where we study the impact of aspherical dark matter density distributions in Milky-Way-like halos on direct and indirect searches. Using data from the large N-body cosmological simulation Bolshoi, we perform a statistical analysis and quantify the systematic uncertainties on the determination of local dark matter density and the so-called J factors for dark matter annihilations and decays from the galactic center. We find that, due to our ignorance about the extent of the non-sphericity of the Milky Way dark matter halo, systematic uncertainties can be as large as 35%, within the 95% most probable region, for a spherically averaged value for the local density of 0.3-0.4 GeV/cm(3). Similarly, systematic uncertainties on the J factors evaluated around the galactic center can be as large as 10% and 15%, within the 95% most probable region, for dark matter annihilations and decays, respectively.
|
BRIKEN Collaboration(Tarifeño-Saldivia, A. et al), Tain, J. L., Domingo-Pardo, C., Agramunt, J., Algora, A., Morales, A. I., et al. (2017). Conceptual design of a hybrid neutron-gamma detector for study of beta-delayed neutrons at the RIB facility of RIKEN. J. Instrum., 12, P04006–22pp.
Abstract: BRIKEN is a complex detection system to be installed at the RIB-facility of the RIKEN Nishina Center. It is aimed at the detection of heavy-ion implants, β-particles, γ-rays and β-delayed neutrons. The whole detection setup involves the Advanced Implantation Detection Array (AIDA), two HPGe Clover detectors and a large set of 166 counters of 3He embedded in a high-density polyethylene matrix. This article reports on a novel methodology developed for the conceptual design and optimisation of the 3He-tubes array, aiming at the best possible performance in terms of neutron detection. The algorithm is based on a geometric representation of two selected parameters of merit, namely, average neutron detection efficiency and efficiency flatness, as a function of a reduced number of geometric variables. The response of the detection system itself, for each configuration, is obtained from a systematic MC-simulation implemented realistically in Geant4. This approach has been found to be particularly useful. On the one hand, due to the different types and large number of 3He-tubes involved and, on the other hand, due to the additional constraints introduced by the ancillary detectors for charged particles and gamma-rays. Empowered by the robustness of the algorithm, we have been able to design a versatile detection system, which can be easily re-arranged into a compact mode in order to maximize the neutron detection performance, at the cost of the gamma-ray sensitivity. In summary, we have designed a system which shows, for neutron energies up to 1(5) MeV, a rather flat and high average efficiency of 68.6%(64%) and 75.7%(71%) for the hybrid and compact modes, respectively. The performance of the BRIKEN system has been also quantified realistically by means of MC-simulations made with different neutron energy distributions.
|
Roser, J., Muñoz, E., Barrientos, L., Barrio, J., Bernabeu, J., Borja-Lloret, M., et al. (2020). Image reconstruction for a multi-layer Compton telescope: an analytical model for three interaction events. Phys. Med. Biol., 65(14), 145005–17pp.
Abstract: Compton Cameras are electronically collimated photon imagers suitable for sub-MeV to few MeV gamma-ray detection. Such features are desirable to enablein vivorange verification in hadron therapy, through the detection of secondary Prompt Gammas. A major concern with this technique is the poor image quality obtained when the incoming gamma-ray energy is unknown. Compton Cameras with more than two detector planes (multi-layer Compton Cameras) have been proposed as a solution, given that these devices incorporate more signal sequences of interactions to the conventional two interaction events. In particular, three interaction events convey more spectral information as they allow inferring directly the incident gamma-ray energy. A three-layer Compton Telescope based on continuous Lanthanum (III) Bromide crystals coupled to Silicon Photomultipliers is being developed at the IRIS group of IFIC-Valencia. In a previous work we proposed a spectral reconstruction algorithm for two interaction events based on an analytical model for the formation of the signal. To fully exploit the capabilities of our prototype, we present here an extension of the model for three interaction events. Analytical expressions of the sensitivity and the System Matrix are derived and validated against Monte Carlo simulations. Implemented in a List Mode Maximum Likelihood Expectation Maximization algorithm, the proposed model allows us to obtain four-dimensional (energy and position) images by using exclusively three interaction events. We are able to recover the correct spectrum and spatial distribution of gamma-ray sources when ideal data are employed. However, the uncertainties associated to experimental measurements result in a degradation when real data from complex structures are employed. Incorrect estimation of the incident gamma-ray interaction positions, and missing deposited energy associated with escaping secondaries, have been identified as the causes of such degradation by means of a detailed Monte Carlo study. As expected, our current experimental resolution and efficiency to three interaction events prevents us from correctly recovering complex structures of radioactive sources. However, given the better spectral information conveyed by three interaction events, we expect an improvement of the image quality of conventional Compton imaging when including such events. In this regard, future development includes the incorporation of the model assessed in this work to the two interaction events model in order to allow using simultaneously two and three interaction events in the image reconstruction.
|