|
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.
|
|
|
ANTARES Collaboration(Aguilar, J. A. et al), Bigongiari, C., Dornic, D., Emanuele, U., Gomez-Gonzalez, J. P., Hernandez-Rey, J. J., et al. (2010). Performance of the front-end electronics of the ANTARES neutrino telescope. Nucl. Instrum. Methods Phys. Res. A, 622(1), 59–73.
Abstract: ANTARES is a high-energy neutrino telescope installed in the Mediterranean Sea at a depth of 2475 m. It consists of a three-dimensional array of optical modules, each containing a large photomultiplier tube. A total of 2700 front-end ASICs named analogue ring samplers (ARS) process the phototube signals, measure their arrival time, amplitude and shape as well as perform monitoring and calibration tasks. The ARS chip processes the analogue signals from the optical modules and converts information into digital data. All the information is transmitted to shore through further multiplexing electronics and an optical link. This paper describes the performance of the ARS chip: results from the functionality and characterization tests in the laboratory are summarized and the long-term performance in the apparatus is illustrated.
|
|
|
Pierre Auger Collaboration(Abreu, P. et al), & Pastor, S. (2010). Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter. Astropart Phys., 34(5), 314–326.
Abstract: Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, 6 x 10(19) eV. The anisotropy was measured by the fraction of arrival directions that are less than 3.1 degrees from the position of an active galactic nucleus within 75 Mpc (using the Veron-Cetty and Veron 12th catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is (38(-6)(+7))%, compared with 21% expected for isotropic cosmic rays. This is down from the early estimate of (69-(+11)(13))%. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation.
|
|
|
Aguilar, A. C., Binosi, D., & Papavassiliou, J. (2010). QCD effective charges from lattice data. J. High Energy Phys., 07(7), 002–24pp.
Abstract: We use recent lattice data on the gluon and ghost propagators, as well as the Kugo-Ojima function, in order to extract the non-perturbative behavior of two particular definitions of the QCD effective charge, one based on the pinch technique construction, and one obtained from the standard ghost-gluon vertex. The construction relies crucially on the definition of two dimensionful quantities, which are invariant under the renormalization group, and are built out of very particular combinations of the aforementioned Green's functions. The main non-perturbative feature of both effective charges, encoded in the infrared finiteness of the gluon propagator and ghost dressing function used in their definition, is the freezing at a common finite (non-vanishing) value, in agreement with a plethora of theoretical and phenomenological expectations. We discuss the sizable discrepancy between the freezing values obtained from the present lattice analysis and the corresponding estimates derived from several phenomenological studies, and attribute its origin to the difference in the gauges employed. A particular toy calculation suggests that the modifications induced to the non-perturbative gluon propagator by the gauge choice may indeed account for the observed deviation of the freezing values.
|
|
|
Cervera-Villanueva, A., Laing, A., Martin-Albo, J., & Soler, F. J. P. (2010). Performance of the MIND detector at a Neutrino Factory using realistic muon reconstruction. Nucl. Instrum. Methods Phys. Res. A, 624(3), 601–614.
Abstract: A Neutrino Factory producing an intense beam composed of v(e)((v) over bar (e)) and (v) over bar (mu)(v(mu)) from muon decays has been shown to have the greatest sensitivity to the two currently unmeasured neutrino mixing parameters theta(13) and delta(CP) Using the wrong-sign muon signal to measure v(e)-> v(mu)((v) over bar (e) ->(v) over bar (mu)) oscillations in a 50kt Magnetised Iron Neutrino Detector (MIND) sensitivity to delta(CP) could be maintained down to small values of theta(13) However the detector efficiencies used in these previous studies were calculated assuming perfect pattern recognition In this paper MIND is reassessed taking into account for the first time a realistic pattern recognition for the muon candidate Reoptimisation of the analysis utilises a combination of methods including a multivariate analysis similar to the one used in MINOS to maintain high efficiency while suppressing backgrounds ensuring that the signal selection efficiency and the background levels are comparable or better than the ones in previous analyses As a result MIND remains the most sensitive future facility for the discovery of CP violation from neutrino oscillations.
|
|
|
Carbone, C., Mena, O., & Verde, L. (2010). Cosmological parameters degeneracies and non-Gaussian halo bias. J. Cosmol. Astropart. Phys., 07(7), 020–17pp.
Abstract: We study the impact of the cosmological parameters uncertainties on the measurements of primordial non-Gaussianity through the large-scale non-Gaussian halo bias effect. While this is not expected to be an issue for the standard Lambda CDM model, it may not be the case for more general models that modify the large-scale shape of the power spectrum. We consider the so-called local non-Gaussianity model, parametrized by the f(NL) non-Gaussianity parameter which is zero for a Gaussian case, and make forecasts on f(NL) from planned surveys, alone and combined with a Planck CMB prior. In particular, we consider EUCLID- and LSST-like surveys and forecast the correlations among f(NL) and the running of the spectral index alpha(s), the dark energy equation of state w, the effective sound speed of dark energy perturbations c(s)(2), the total mass of massive neutrinos M-nu = Sigma m(nu), and the number of extra relativistic degrees of freedom N-nu(rel). Neglecting CMB information on f(NL) and scales k > 0.03h/Mpc, we find that, if N-nu(rel) is assumed to be known, the uncertainty on cosmological parameters increases the error on f(NL) by 10 to 30% depending on the survey. Thus the f(NL) constraint is remarkable robust to cosmological model uncertainties. On the other hand, if N-nu(rel) is simultaneously constrained from the data, the f(NL) error increases by similar to 80%. Finally, future surveys which provide a large sample of galaxies or galaxy clusters over a volume comparable to the Hubble volume can measure primordial non-Gaussianity of the local form with a marginalized 1-sigma error of the order Delta f(NL) similar to 2 – 5, after combination with CMB priors for the remaining cosmological parameters. These results are competitive with CMB bispectrum constraints achievable with an ideal CMB experiment.
|
|
|
Belver, D., Cabanelas, P., Castro, E., Garzon, J. A., Gil, A., Gonzalez-Diaz, D., et al. (2010). Performance of the Low-Jitter High-Gain/Bandwidth Front-End Electronics of the HADES tRPC Wall. IEEE Trans. Nucl. Sci., 57(5), 2848–2856.
Abstract: A front-end electronics (FEE) chain for accurate time measurements has been developed for the new Resistive Plate Chamber (RPC)-based Time-of-Flight (TOF) wall of the High Acceptance Di-Electron Spectrometer (HADES). The wall covers an area of around 8 m(2) divided in 6 sectors. In total, 1122 4-gap timing RPC cells are read-out by 2244 time and charge sensitive channels. The FEE chain consists of 2 custom-made boards: a 4-channel Daughter BOard(DBO) and a 32-channel MotherBOard (MBO). The DBO uses a fast 2 GHz amplifier feeding a dual high-speed discriminator. The time and charge information are encoded, respectively, in the leading edge and the width of an LVDS signal. Each MBO houses up to 8 DBOs providing them regulated voltage supply, threshold values via DACs, test signals and, additionally, routing out a signal proportional to the channel multiplicity needed for a 1st level trigger decision. The MBO delivers LVDS signals to a multi-purpose Trigger Readout Board (TRB) for data acquisition. The FEE allows achieving a system resolution around 75 ps fulfilling comfortably the requirements of the HADES upgrade [1]. The commissioning of the whole RPC wall is finished and the 6 sectors are already mounted in their final position in the HADES spectrometer and ready to take data during the beam-times foreseen for 2010.
|
|
|
Horvat, S., Magas, V. K., Strottman, D. D., & Csernai, L. P. (2010). Entropy development in ideal relativistic fluid dynamics with the Bag Model equation of state. Phys. Lett. B, 692(4), 277–280.
Abstract: We consider an idealized situation where the Quark-Gluon Plasma (QGP) is described by a perfect, (3 + 1)-dimensional fluid dynamic model starting from an initial state and expanding until a final state where freeze-out and/or hadronization takes place. We study the entropy production with attention to effects of (i) numerical viscosity, (ii) late stages of flow where the Bag Constant and the partonic pressure are becoming similar, (iii) and the consequences of final freeze-out and constituent quark matter formation.
|
|
|
Aguiar, P., Rafecas, M., Ortuño, J. E., Kontaxakis, G., Santos, A., Pavia, J., et al. (2010). Geometrical and Monte Carlo projectors in 3D PET reconstruction. Med. Phys., 37(11), 5691–5702.
Abstract: Purpose: In the present work, the authors compare geometrical and Monte Carlo projectors in detail. The geometrical projectors considered were the conventional geometrical Siddon ray-tracer (S-RT) and the orthogonal distance-based ray-tracer (OD-RT), based on computing the orthogonal distance from the center of image voxel to the line-of-response. A comparison of these geometrical projectors was performed using different point spread function (PSF) models. The Monte Carlo-based method under consideration involves an extensive model of the system response matrix based on Monte Carlo simulations and is computed off-line and stored on disk. Methods: Comparisons were performed using simulated and experimental data of the commercial small animal PET scanner rPET. Results: The results demonstrate that the orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions yield better images in terms of contrast and spatial resolution than those obtained after using the conventional method and the multiray-based S-RT. Furthermore, the Monte Carlo-based method yields slight improvements in terms of contrast and spatial resolution with respect to these geometrical projectors. Conclusions: The orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions represent satisfactory alternatives to factorizing the system matrix or to the conventional on-the-fly ray-tracing methods for list-mode reconstruction, where an extensive modeling based on Monte Carlo simulations is unfeasible.
|
|
|
Wagner, C., Verde, L., & Boubekeur, L. (2010). N-body simulations with generic non-Gaussian initial conditions I: power spectrum and halo mass function. J. Cosmol. Astropart. Phys., 10(10), 022–24pp.
Abstract: We address the issue of setting up generic non-Gaussian initial conditions for N-body simulations. We consider inflationary-motivated primordial non-Gaussianity where the perturbations in the Bardeen potential are given by a dominant Gaussian part plus a non-Gaussian part specified by its bispectrum. The approach we explore here is suitable for any bispectrum, i.e. it does not have to be of the so-called separable or factorizable form. The procedure of generating a non-Gaussian field with a given bispectrum (and a given power spectrum for the Gaussian component) is not univocal, and care must be taken so that higher-order corrections do not leave a too large signature on the power spectrum. This is so far a limiting factor of our approach. We then run N-body simulations for the most popular inflationary-motivated non-Gaussian shapes. The halo mass function and the non-linear power spectrum agree with theoretical analytical approximations proposed in the literature, even if they were so far developed and tested only for a particular shape (the local one). We plan to make the simulations outputs available to the community via the non-Gaussian simulations comparison project web site http://icc.ub.edu/similar to liciaverde/NGSCP.html.
|
|