We calculate three jet cross sections in photoproduction using exact matrix
elements for the direct and resolved contributions. Numerical distribution
s are presented in a generic, irreducible set of variables that allows to d
isentangle the dynamics of partonic QCD subprocesses from each other and fr
om pure phase space distributions. The results are compared to preliminary
data from the ZEUS collaboration at HERA. It is found that the largest cont
ribution comes from photon-gluon fusion in the mass range 36 GeV < M3-jet <
80 GeV. The measured leading jet scattering angle distribution is consiste
nt with the t-channel exchange of a massless fermion in 2 --> 2 scattering,
where the third parton is assumed to arise from soft bremsstrahlung. The d
ata are inconsistent with pure phase space and Rutherford scattering distri
butions.