|
|
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.
|
|
|
|
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.
|
|
|
|
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.
|
|
|
|
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.
|
|
|
|
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.
|
|
