We present a new approach to computing event shape distributions or, more precisely, charge flow correlations in a generic conformal field theory (CFT). These infrared finite observables are familiar from collider physics studies and describe the angular distribution of global charges in outgoing radiation created from the vacuum by some source. The charge flow correlations can be expressed in terms of Wightman correlation functions in a certain limit. We explain how to compute these quantities starting from their Euclidean analogues by means of a nontrivial analytic continuation which, in the framework of CFT, can be performed elegantly in Mellin space. The relation between the charge flow correlations and Euclidean correlation functions can be reformulated directly in configuration space, bypassing the Mellin representation, as a certain Lorentzian double discontinuity of the correlation function integrated along the cuts. We illustrate the general formalism in N = 4 SYM, making use of the well-known results on the four-point correlation function of half-BPS scalar operators. We compute the double scalar flow correlation in N = 4 SYM, at weak and strong coupling and show that it agrees with known results obtained by different techniques. One of the remarkable features of the N = 4 theory is that the scalar and energy flow correlations are proportional to each other. Imposing natural physical conditions on the energy flow correlations (finiteness, positivity and regularity), we formulate additional constraints on the four-point correlation functions in N = 4SYM that should be valid at any coupling and away from the planar limit.

We study event shapes in N = 4SYM describing the angular distribution of energy and R-charge in the final states created by the simplest half-BPS scalar operator. Applying the approach developed in the companion paper arXiv:1309.0769, we compute these observables using the correlation functions of certain components of the N = 4 stress-tensor supermultiplet: the half-BPS operator itself, the R-symmetry current and the stress tensor. We present master formulas for the all-order event shapes as convolutions of the Mellin amplitude defining the correlation function of the half-BPS operators, with a coupling-independent kernel determined by the choice of the observable. We find remarkably simple relations between various event shapes following from N = 4 superconformal symmetry. We perform thorough checks at leading order in the weak coupling expansion and show perfect agreement with the conventional calculations based on amplitude techniques. We extend our results to strong coupling using the correlation function of half-BPS operators obtained from the AdS/CFT correspondence.