|
Oliver, J. F., & Rafecas, M. (2010). Improving the singles rate method for modeling accidental coincidences in high-resolution PET. Phys. Med. Biol., 55(22), 6951–6971.
Abstract: Random coincidences ('randoms') are one of the main sources of image degradation in PET imaging. In order to correct for this effect, an accurate method to estimate the contribution of random events is necessary. This aspect becomes especially relevant for high-resolution PET scanners where the highest image quality is sought and accurate quantitative analysis is undertaken. One common approach to estimate randoms is the so-called singles rate method (SR) widely used because of its good statistical properties. SR is based on the measurement of the singles rate in each detector element. However, recent studies suggest that SR systematically overestimates the correct random rate. This overestimation can be particularly marked for low energy thresholds, below 250 keV used in some applications and could entail a significant image degradation. In this work, we investigate the performance of SR as a function of the activity, geometry of the source and energy acceptance window used. We also investigate the performance of an alternative method, which we call 'singles trues' (ST) that improves SR by properly modeling the presence of true coincidences in the sample. Nevertheless, in any real data acquisition the knowledge of which singles are members of a true coincidence is lost. Therefore, we propose an iterative method, STi, that provides an estimation based on ST but which only requires the knowledge of measurable quantities: prompts and singles. Due to inter-crystal scatter, for wide energy windows ST only partially corrects SR overestimations. While SR deviations are in the range 86-300% (depending on the source geometry), the ST deviations are systematically smaller and contained in the range 4-60%. STi fails to reproduce the ST results, although for not too high activities the deviation with respect to ST is only a few percent. For conventional energy windows, i.e. those without inter-crystal scatter, the ST method corrects the SR overestimations, and deviations from the true random rate are of the order of 1% or less. In addition, in the case of conventional energy window STi results reproduce ST results and therefore the former can be used to obtain the true random rate.
|
|
|
BABAR Collaboration(del Amo Sanchez, P. et al), Lopez-March, N., Martinez-Vidal, F., Milanes, D. A., & Oyanguren, A. (2010). Search for B+ meson decay to a(1)(+)(1260)K*(0)(892). Phys. Rev. D, 82(9), 091101–8pp.
Abstract: We present a search for the decay B+ -> a(1)(+)(1260)K*(0)(892). The data, collected with the BABAR detector at the SLAC National Accelerator Laboratory, represent 465 X 10(6)B (B) over bar pairs produced in e(+)e(-) annihilation at the energy of the gamma(4S). We find no significant signal and set an upper limit at 90% confidence level on the product of branching fractions B(B+ -> a(1)(+) (1260)K*(0)(892)) X B(a(1)(+)(1260) -> pi(+) pi(-) pi(+)) of 1.8 X 10(-6).
|
|
|
CDF Collaboration(Aaltonen, T. et al), & Cabrera, S. (2010). Measurement of the WW plus WZ Production Cross Section Using the lepton plus jets Final State at CDF II. Phys. Rev. Lett., 104(10), 101801–8pp.
Abstract: We report two complementary measurements of the WW + WZ cross section in the final state consisting of an electron or muon, missing transverse energy, and jets, performed using p (p) over bar collision data at root s = 1.96 TeV collected by the CDF II detector. The first method uses the dijet invariant mass distribution while the second more sensitive method uses matrix-element calculations. The result from the second method has a signal significance of 5.4 sigma and is the first observation of WW + WZ production using this signature. Combining the results gives sigma(WW+WZ) = 16.0 +/- 3.3 pb, in agreement with the standard model prediction.
|
|
|
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.
|
|