Analytical calculation of volumes-of-intersection for iterative, fully 3-D PET reconstruction

Use of iterative algorithms to reconstruct three-dimensional (3-D) positron emission tomography (PET) data requires the computation of the system probability matrix. The pure geometrical contribution can easily be approximated by the length-of-intersection (LOI) between lines-of-response (LOR) and individual voxels. However, more accurate geometrical projectors are desirable. Therefore, we have developed a fast method for the analytical calculation of the 3-D shape and volume of volumes-of-intersection (VOI). This method provides an alternative robust projector with a uniformly continuous sampling of the image space. The enhanced calculation effort is facilitated by using several speedup techniques. Exploiting intrinsic symmetry relations and the sparseness of the system matrix allows to create an efficiently compressed matrix which can be precomputed and completely stored in memory. In addition, a new voxel addressing scheme has been implemented. This scheme avoids time-consuming symmetry transformations of voxel addresses by using an octant-wise symmetrically ordered field of voxels. The above methods have been applied for a fully 3-D, iterative reconstruction of 3-D sinograms recorded with a Siemens/CTI ECAT HR+ PET scanner. A comparison of the performance of the reconstruction using LOI weighting and VOI weighting is presented

Massively parallelizable list-mode reconstruction using a Monte Carlo-based elliptical Gaussian model

Purpose: A fully three-dimensional (3D) massively parallelizable list-mode ordered-subsets expectation-maximization (LM-OSEM) reconstruction algorithm has been developed for high-resolution PET cameras. System response probabilities are calculated online from a set of parameters derived from Monte Carlo simulations. The shape of a system response for a given line of response (LOR) has been shown to be asymmetrical around the LOR. This work has been focused on the development of efficient region-search techniques to sample the system response probabilities, which are suitable for asymmetric kernel models, including elliptical Gaussian models that allow for high accuracy and high parallelization efficiency. The novel region-search scheme using variable kernel models is applied in the proposed PET reconstruction algorithm...This work was partially supported by Spain’s Ministry of Science and Innovation through CDTI’s CENIT program (AMIT project) and INNPACTO (PRECISION project), Instituto de Salud Carlos III (PI09/91058 and PI09/91065), and Project Nos. TEC2010-21619-C04-03 and TEC2011-28972- C02-02, Comunidad de Madrid (ARTEMIS S2009/DPI- 1802), and the European Regional Development Funds (FEDER). CIBER-BBN is an initiative funded by the VI National R&D&i Plan 2008–2011, Iniciativa Ingenio 2010, Consolider Program, CIBER Actions, and finance by the Instituto de Salud Carlos III with assistance from the European Regional Development Fund.Publicad

Searches at HERA for Squarks in R-Parity Violating Supersymmetry

A search for squarks in R-parity violating supersymmetry is performed in e^+p collisions at HERA at a centre of mass energy of 300 GeV, using H1 data corresponding to an integrated luminosity of 37 pb^(-1). The direct production of single squarks of any generation in positron-quark fusion via a Yukawa coupling lambda' is considered, taking into account R-parity violating and conserving decays of the squarks. No significant deviation from the Standard Model expectation is found. The results are interpreted in terms of constraints within the Minimal Supersymmetric Standard Model (MSSM), the constrained MSSM and the minimal Supergravity model, and their sensitivity to the model parameters is studied in detail. For a Yukawa coupling of electromagnetic strength, squark masses below 260 GeV are excluded at 95% confidence level in a large part of the parameter space. For a 100 times smaller coupling strength masses up to 182 GeV are excluded.Comment: 32 pages, 14 figures, 3 table

Differential (2+1) Jet Event Rates and Determination of alpha_s in Deep Inelastic Scattering at HERA

Events with a (2+1) jet topology in deep-inelastic scattering at HERA are studied in the kinematic range 200 < Q^2< 10,000 GeV^2. The rate of (2+1) jet events has been determined with the modified JADE jet algorithm as a function of the jet resolution parameter and is compared with the predictions of Monte Carlo models. In addition, the event rate is corrected for both hadronization and detector effects and is compared with next-to-leading order QCD calculations. A value of the strong coupling constant of alpha_s(M_Z^2)= 0.118+- 0.002 (stat.)^(+0.007)_(-0.008) (syst.)^(+0.007)_(-0.006) (theory) is extracted. The systematic error includes uncertainties in the calorimeter energy calibration, in the description of the data by current Monte Carlo models, and in the knowledge of the parton densities. The theoretical error is dominated by the renormalization scale ambiguity.Comment: 25 pages, 6 figures, 3 tables, submitted to Eur. Phys.

Multi-Jet Event Rates in Deep Inelastic Scattering and Determination of the Strong Coupling Constant

Jet event rates in deep inelastic ep scattering at HERA are investigated applying the modified JADE jet algorithm. The analysis uses data taken with the H1 detector in 1994 and 1995. The data are corrected for detector and hadronization effects and then compared with perturbative QCD predictions using next-to-leading order calculations. The strong coupling constant alpha_S(M_Z^2) is determined evaluating the jet event rates. Values of alpha_S(Q^2) are extracted in four different bins of the negative squared momentum transfer~\qq in the range from 40 GeV2 to 4000 GeV2. A combined fit of the renormalization group equation to these several alpha_S(Q^2) values results in alpha_S(M_Z^2) = 0.117+-0.003(stat)+0.009-0.013(syst)+0.006(jet algorithm).Comment: 17 pages, 4 figures, 3 tables, this version to appear in Eur. Phys. J.; it replaces first posted hep-ex/9807019 which had incorrect figure 4

Multiplicity Structure of the Hadronic Final State in Diffractive Deep-Inelastic Scattering at HERA

The multiplicity structure of the hadronic system X produced in deep-inelastic processes at HERA of the type ep -> eXY, where Y is a hadronic system with mass M_Y< 1.6 GeV and where the squared momentum transfer at the pY vertex, t, is limited to |t|<1 GeV^2, is studied as a function of the invariant mass M_X of the system X. Results are presented on multiplicity distributions and multiplicity moments, rapidity spectra and forward-backward correlations in the centre-of-mass system of X. The data are compared to results in e+e- annihilation, fixed-target lepton-nucleon collisions, hadro-produced diffractive final states and to non-diffractive hadron-hadron collisions. The comparison suggests a production mechanism of virtual photon dissociation which involves a mixture of partonic states and a significant gluon content. The data are well described by a model, based on a QCD-Regge analysis of the diffractive structure function, which assumes a large hard gluonic component of the colourless exchange at low Q^2. A model with soft colour interactions is also successful.Comment: 22 pages, 4 figures, submitted to Eur. Phys. J., error in first submission - omitted bibliograph

Measurements of Transverse Energy Flow in Deep-Inelastic Scattering at HERA

Measurements of transverse energy flow are presented for neutral current deep-inelastic scattering events produced in positron-proton collisions at HERA. The kinematic range covers squared momentum transfers Q^2 from 3.2 to 2,200 GeV^2, the Bjorken scaling variable x from 8.10^{-5} to 0.11 and the hadronic mass W from 66 to 233 GeV. The transverse energy flow is measured in the hadronic centre of mass frame and is studied as a function of Q^2, x, W and pseudorapidity. A comparison is made with QCD based models. The behaviour of the mean transverse energy in the central pseudorapidity region and an interval corresponding to the photon fragmentation region are analysed as a function of Q^2 and W.Comment: 26 pages, 8 figures, submitted to Eur. Phys.

Forward pi^0 Production and Associated Transverse Energy Flow in Deep-Inelastic Scattering at HERA

Deep-inelastic positron-proton interactions at low values of Bjorken-x down to x \approx 4.10^-5 which give rise to high transverse momentum pi^0 mesons are studied with the H1 experiment at HERA. The inclusive cross section for pi^0 mesons produced at small angles with respect to the proton remnant (the forward region) is presented as a function of the transverse momentum and energy of the pi^0 and of the four-momentum transfer Q^2 and Bjorken-x. Measurements are also presented of the transverse energy flow in events containing a forward pi^0 meson. Hadronic final state calculations based on QCD models implementing different parton evolution schemes are confronted with the data.Comment: 27 pages, 8 figures and 3 table

Low Q^2 Jet Production at HERA and Virtual Photon Structure

The transition between photoproduction and deep-inelastic scattering is investigated in jet production at the HERA ep collider, using data collected by the H1 experiment. Measurements of the differential inclusive jet cross-sections dsigep/dEt* and dsigmep/deta*, where Et* and eta* are the transverse energy and the pseudorapidity of the jets in the virtual photon-proton centre of mass frame, are presented for 0 < Q2 < 49 GeV2 and 0.3 < y < 0.6. The interpretation of the results in terms of the structure of the virtual photon is discussed. The data are best described by QCD calculations which include a partonic structure of the virtual photon that evolves with Q2.Comment: 20 pages, 5 Figure

Hadron Production in Diffractive Deep-Inelastic Scattering

Characteristics of hadron production in diffractive deep-inelastic positron-proton scattering are studied using data collected in 1994 by the H1 experiment at HERA. The following distributions are measured in the centre-of-mass frame of the photon dissociation system: the hadronic energy flow, the Feynman-x (x_F) variable for charged particles, the squared transverse momentum of charged particles (p_T^{*2}), and the mean p_T^{*2} as a function of x_F. These distributions are compared with results in the gamma^* p centre-of-mass frame from inclusive deep-inelastic scattering in the fixed-target experiment EMC, and also with the predictions of several Monte Carlo calculations. The data are consistent with a picture in which the partonic structure of the diffractive exchange is dominated at low Q^2 by hard gluons.Comment: 16 pages, 6 figures, submitted to Phys. Lett.