Non-perturbative renormalization in coordinate space for Nf=2N_f=2 maximally twisted mass fermions with tree-level Symanzik improved gauge action

We present results of a lattice QCD application of a coordinate space renormalization scheme for the extraction of renormalization constants for flavour non-singlet bilinear quark operators. The method consists in the analysis of the small-distance behaviour of correlation functions in Euclidean space and has several theoretical and practical advantages, in particular: it is gauge invariant, easy to implement and has relatively low computational cost. The values of renormalization constants in the X-space scheme can be converted to the MSbar scheme via 4-loop continuum perturbative formulae. Our results for $N_f=2$ maximally twisted mass fermions with tree-level Symanzik improved gauge action are compared to the ones from the RI-MOM scheme and show full agreement with this method.Comment: 29 pages, 7 figures, 6 table

The use of Ethernet in the DataFlow of the ATLAS Trigger & DAQ

The article analyzes a proposed network topology for the ATLAS DAQ DataFlow, and identifies the Ethernet features required for a proper operation of the network: MAC address table size, switch performance in terms of throughput and latency, the use of Flow Control, Virtual LANs and Quality of Service. We investigate these features on some Ethernet switches, and conclude on their usefulness for the ATLAS DataFlow network.Comment: Talk from the 2003 Computing in High Energy and Nuclear Physics (CHEP03), La Jolla, Ca, March 2003, 10 pages, LaTeX, 10 eps figures. PSN MOGT01

J-PET Framework: Software platform for PET tomography data reconstruction and analysis

J-PET Framework is an open-source software platform for data analysis, written in C++ and based on the ROOT package. It provides a common environment for implementation of reconstruction, calibration and filtering procedures, as well as for user-level analyses of Positron Emission Tomography data. The library contains a set of building blocks that can be combined by users with even little programming experience, into chains of processing tasks through a convenient, simple and well-documented API. The generic input-output interface allows processing the data from various sources: low-level data from the tomography acquisition system or from diagnostic setups such as digital oscilloscopes, as well as high-level tomography structures e.g. sinograms or a list of lines-of-response. Moreover, the environment can be interfaced with Monte Carlo simulation packages such as GEANT and GATE, which are commonly used in the medical scientific community.Comment: 14 pages, 5 figure

Step scaling in X-space : running of the quark mass

We perform a benchmark study of the step scaling procedure for the ratios of renormalization constants extracted from position space correlation functions. We work in the quenched approximation and consider the pseudoscalar, scalar, vector and axial vector bilinears. The pseudoscalar/scalar cases allow us to obtain the non-perturbative running of the quark mass over a wide range of energy scales - from around 17 GeV to below 1.5 GeV - which agrees well with the 4-loop prediction of continuum perturbation theory. We find that step scaling is feasible in X-space and we discuss its advantages and potential problems

Selected Issues on Histograming on GPUs

The contemporary large scale measuring systems in the real-time environment make extensive use of histogramming as a tool for the experimental data quality monitoring. The processing of a large number of data channels requires a suitable computing power where the graphical processors seem to be well suited. Histogramming operations run on the central and graphics processing units are discussed. Results of the performance measurements including various configurations of the allocation of the histograms in various parts of the memory of used devices are presented

Non-perturbative running of renormalization constants from correlators in coordinate space using step scaling

Working in a quenched setup with Wilson twisted mass valence fermions, we explore the possibility to compute non-perturbatively the step scaling function using the coordinate (X-space) renormalization scheme. This scheme has the advantage of being on-shell and gauge invariant. The step scaling method allows us to calculate the running of the renormalization constants of quark bilinear operators. We describe here the details of this calculation. The aim of this exploratory study is to identify the feasibility of the X-space scheme when used in small volume simulations required by the step scaling technique. Eventually, we translate our final results to the continuum MSscheme and compare against four-loop analytic formulae finding satisfactory agreement. (C) 2016 The Authors. Published by Elsevier B. V. This is an open access article under the CC BY license

Simulating NEMA characteristics of the modular total-body J-PET scanner -- an economic total-body PET from plastic scintillators

The purpose of the presented research is estimation of the performance characteristics of the economic Total-Body Jagiellonian-PET system (TB-J-PET) constructed from plastic scintillators. The characteristics are estimated according to the NEMA NU-2-2018 standards utilizing the GATE package. The simulated detector consists of 24 modules, each built out of 32 plastic scintillator strips (each with cross section of 6 mm times 30 mm and length of 140 cm or 200 cm) arranged in two layers in regular 24-sided polygon circumscribing a circle with the diameter of 78.6 cm. For the TB-J-PET with an axial field-of-view (AFOV) of 200 cm, a spatial resolutions of 3.7 mm (transversal) and 4.9 mm (axial) are achieved. The NECR peak of 630 kcps is expected at 30 kBq/cc activity concentration and the sensitivity at the center amounts to 38 cps/kBq. The SF is estimated to 36.2 %. The values of SF and spatial resolution are comparable to those obtained for the state-of-the-art clinical PET scanners and the first total-body tomographs: uExplorer and PennPET. With respect to the standard PET systems with AFOV in the range from 16 cm to 26 cm, the TB-J-PET is characterized by an increase in NECR approximately by factor of 4 and by the increase of the whole-body sensitivity by factor of 12.6 to 38. The TOF resolution for the TB-J-PET is expected to be at the level of CRT=240 ps (FWHM). For the TB-J-PET with an axial field-of-view (AFOV) of 140 cm, an image quality of the reconstructed images of a NEMA IEC phantom was presented with a contrast recovery coefficient (CRC) and a background variability parameters. The increase of the whole-body sensitivity and NECR estimated for the TB-J-PET with respect to current commercial PET systems makes the TB-J-PET a promising cost-effective solution for the broad clinical applications of total-body PET scanners.Comment: 31 pages, 11 figures, 6 tables, submitted to Physics in Medicine and Biology 202

Optimisation of the event-based TOF filtered back-projection for online imaging in total-body J-PET

We perform a parametric study of the newly developed time-of-flight (TOF) image reconstruction algorithm, proposed for the real-time imaging in total-body Jagiellonian PET (J-PET) scanners. The asymmetric 3D filtering kernel is applied at each most likely position of electron-positron annihilation, estimated from the emissions of back-to-back $\gamma$-photons. The optimisation of its parameters is studied using Monte Carlo simulations of a 1-mm spherical source, NEMA IEC and XCAT phantoms inside the ideal J-PET scanner. The combination of high-pass filters which included the TOF filtered back-projection (FBP), resulted in spatial resolution, 1.5 $\times$ higher in the axial direction than for the conventional 3D FBP. For realistic $10$-minute scans of NEMA IEC and XCAT, which require a trade-off between the noise and spatial resolution, the need for Gaussian TOF kernel components, coupled with median post-filtering, is demonstrated. The best sets of 3D filter parameters were obtained by the Nelder-Mead minimisation of the mean squared error between the resulting and reference images. The approach allows training the reconstruction algorithm for custom scans, using the IEC phantom, when the temporal resolution is below 50 ps. The image quality parameters, estimated for the best outcomes, were systematically better than for the non-TOF FBP

Feasibility study of the time reversal symmetry tests in decay of metastable positronium atoms with the J-PET detector

This article reports on the feasibility of testing of the symmetry under reversal in time in a purely leptonic system constituted by positronium atoms using the J-PET detector. The present state of T symmetry tests is discussed with an emphasis on the scarcely explored sector of leptonic systems. Two possible strategies of searching for manifestations of T violation in nonvanishing angular correlations of final state observables in the decay of metastable triplet states of positronium available with J-PET are proposed and discussed. Results of a pilot measurement with J-PET and assessment of its performance in reconstruction of three-photon decays are shown along with an analysis of its impact on the sensitivity of the detector for the determination of T-violation sensitive observables

Realistic Total-Body J-PET Geometry Optimization -- Monte Carlo Study

Total-Body PET is one of the most promising medical diagnostics modalities. The high sensitivity provided by Total-Body technology can be advantageous for novel tomography methods like positronium imaging. Several efforts are ongoing to lower the price of the TB-PET systems. Among the alternatives, the Jagiellonian PET (J-PET) technology, based on plastic scintillator strips, offers a low-cost alternative. The work aimed to compare five Total-Body J-PET geometries as a possible next generation J-PET scanner design. We present comparative studies of performance characteristics of the cost-effective Total-Body PET scanners using J-PET technology. We investigated in silico five Total-Body scanner geometries. Monte Carlo simulations of the XCAT phantom, the 2-meter sensitivity line source and positronium sensitivity phantoms were performed. We compared the sensitivity profiles for 2-gamma and 3-gamma tomography, relative cost of the setups and performed quantitative analysis of the reconstructed images. The analysis of the reconstructed XCAT images reveals the superiority of the seven-ring scanners over the three-ring setups. However, the three-ring scanners would be approximately 2-3 times cheaper. The peak sensitivity values for two-gamma vary from 20 to 34 cps/kBq. The sensitivity curves for the positronium tomography have a similar shape to the two-gamma sensitivity profiles. The peak values are lower compared to the two-gamma cases, from about 20-28 times, with a maximum of 1.66 cps/kBq. The results show the feasibility of multi-organ imaging of all the systems to be considered for the next generation of TB J-PET designs. The relative cost for all the scanners is about 10-4 times lower compared to the cost of the uExplorer. These properties coupled together with J-PET cost-effectiveness, make the J-PET technology an attractive solution for broad application in clinics