Rodriguez-Alvarez, M. J., Sanchez, F., Soriano, A., & Iborra, A. (2010). Sparse Givens resolution of large system of linear equations: Applications to image reconstruction. Math. Comput. Model., 52(7-8), 1258–1264.
Abstract: In medicine, computed tomographic images are reconstructed from a large number of measurements of X-ray transmission through the patient (projection data). The mathematical model used to describe a computed tomography device is a large system of linear equations of the form AX = B. In this paper we propose the QR decomposition as a direct method to solve the linear system. QR decomposition can be a large computational procedure. However, once it has been calculated for a specific system, matrices Q and R are stored and used for any acquired projection on that system. Implementation of the QR decomposition in order to take more advantage of the sparsity of the system matrix is discussed.
|
Rodriguez-Alvarez, M. J., Sanchez, F., Soriano, A., Iborra, A., & Mora, C. (2011). Exploiting symmetries for weight matrix design in CT imaging. Math. Comput. Model., 54(7-8), 1655–1664.
Abstract: In this paper we propose several methods of constructing the system matrix (SM) of a Computed Tomography (CT) scanner with two objectives: (1) to construct SMs in the shortest possible time and store them in an ordinary PC without losing quality, (2) to analyze the possible applications of the proposed method to 3D, taking into account SMs' sizes, computing time and reconstructed image quality. In order to build the SM, we propose two new field of view (FOV) pixellation schemes, based on a polar coordinate system (polar grid) by taking advantage of the polar rotation symmetries of CT devices. Comparisons between the SMs proposed are performed using two phantom and a real CT-simulator images. Global error, contrast, noise and homogeneity of the reconstructed images are discussed.
|
Freitas, E. D. C., Monteiro, C. M. B., Ball, M., Gomez-Cadenas, J. J., Lopes, J. A. M., Lux, T., et al. (2010). Secondary scintillation yield in high-pressure xenon gas for neutrinoless double beta decay (0 nu beta beta) search. Phys. Lett. B, 684(4-5), 205–210.
Abstract: The search for neutrinoless double beta decay (0 nu beta beta) is an important topic in contemporary physics with many active experiments. New projects are planning to use high-pressure xenon gas as both source and detection medium. The secondary scintillation processes available in noble gases permit large amplification with negligible statistical fluctuations, offering the prospect of energy resolution approaching the Fano factor limit. This Letter reports results for xenon secondary scintillation yield, at room temperature, as a function of electric field in the gas scintillation gap for pressures ranging from 2 to 10 bar. A Large Area Avalanche Photodiode (LAAPD) collected the VUV secondary scintillation produced in the gas. X-rays directly absorbed in the LAAPD are used as a reference for determining the number of charge carriers produced by the scintillation pulse and, hence, the number of photons impinging the LAAPD. The number of photons produced per drifting electron and per kilovolt, the so-called scintillation amplification parameter, displays a small increase with pressure, ranging from 141 +/- 6 at 2 bar to 170 +/- 10 at 8 bar. In our setup, this Parameter does not increase above 8 bar due to nonnegligible electron attachment. The results are in good agreement with those presented in the literature in the 1 to 3 bar range. The increase of the scintillation amplification parameter with pressure for high gas densities has been also observed in former work at cryogenic temperatures.
|
Bourguille, B., Nieves, J., & Sanchez, F. (2021). Inclusive and exclusive neutrino-nucleus cross sections and the reconstruction of the interaction kinematics. J. High Energy Phys., 04(4), 004–42pp.
Abstract: We present a full kinematic analysis of neutrino-nucleus charged current quasielastic interactions based on the Local Fermi Gas model and the Random Phase Approximation. The model was implemented in the NEUT Monte Carlo framework, which allows us to investigate potentially measurable observables, including hadron distributions. We compare the predictions simultaneously to the most recent T2K and MINERvA charged current (CC) inclusive, CC0 pi and transverse kinematic-imbalance variable results. We pursuit a microscopic interpretation of the relevant reaction mechanisms, with the aim to achieving in neutrino oscillation experiments a correct reconstruction of the incoming neutrino kinematics, free of conceptual biasses. Such study is of the utmost importance for the ambitious experimental program which is underway to precisely determine neutrino properties, test the three-generation paradigm, establish the order of mass eigenstates and investigate leptonic CP violation.
|
Sobczyk, J. E., Nieves, J., & Sanchez, F. (2020). Exclusive-final-state hadron observables from neutrino-nucleus multinucleon knockout. Phys. Rev. C, 102(2), 024601–16pp.
Abstract: We present results of an updated calculation of the two particle two hole (2p2h) contribution to the neutrino-induced charge-current cross section. We provide also some exclusive observables, interesting from the point of view of experimental studies, e.g., distributions of momenta of the outgoing nucleons and of available energy, which we compare with the results obtained within the NEUT generator. We also compute, and separate from the total, the contributions of 3p3h mechanisms. Finally, we discuss the differences between the present results and previous implementations of the model in MC event generators, done at the level of inclusive cross sections, which might significantly influence the experimental analyses, particularly in the cases where the hadronic observables are considered.
|