Cone Beam Computerized Tomography (CBCT) and Positron Emission Tomography (PET) Scans are medical imaging devices that require solving ill-posed inverse problems. The models considered come directly from the physics of the acquisition devices, and take into account the specificity of the Poisson noise. We propose various fast numerical schemes to compute the solution. In particular, we show that a new algorithm recently introduced by A. Chambolle and T. Pock is well suited in the PET case when considering non differentiable regularizations such as total variation or wavelet l1-regularization. Numerical experiments indicate that the proposed algorithms compare favorably with respect to well-established methods in tomography. Moreover, the new generation of hybrid pixel detectors working in a photon counting mode makes possible to select photons according to their energy, thus offering the opportunity to go beyond the classical reconstruction techniques that provides an absorption reconstruction volume. We will study the open questions of reconstructing colored X-ray volumes.