We analyse the newest diffractive deep inelastic scattering data from the DESY collider HERA with the help of dipole models. We find good agreement with the data on the diffractive structure functions provided the diffractive open charm contribution
is taken into account. However, the region of large diffractive mass (small values of a parameter beta) needs some refinement with the help of an additional gluon radiation.
This work covers methodology of solving QCD evolution equation of the parton distribution using Markovian Monte Carlo (MMC) algorithms in a class of models ranging from DGLAP to CCFM. One of the purposes of the above MMCs is to test the other more so
phisticated Monte Carlo programs, the so-called Constrained Monte Carlo (CMC) programs, which will be used as a building block in the parton shower MC. This is why the mapping of the evolution variables (eikonal variable and evolution time) into four-momenta is also defined and tested. The evolution time is identified with the rapidity variable of the emitted parton. The presented MMCs are tested independently, with ~0.1% precision, against the non-MC program APCheb especially devised for this purpose.
In this contribution we briefly review the current status of the dipole models and parton saturation on the basis of results presented at the HERA-LHC workshops in the years 2006-2008. The problem of foundations of the dipole models is addressed with
in the QCD formalism. Some limitations of the models and open problems are pointed out. Furthermore, we review and compare the currently used dipole models and summarise the applications to describe various sets of HERA data. Finally we outline some of the theoretical approaches to the problem of multiple scattering and saturation.
Diffractive parton distributions of the proton are determined from fits to diffractive data from HERA. In addition to the twist--2 contribution, the twist--4 contribution from longitudinally polarised virtual photons is considered, which is important
in the region of small diffractive masses. A new prediction for the longitudinal diffractive structure function is presented which differs significantly from that obtained in the pure twist--2 analyses.
