We present a detailed analysis of QED corrections to B¯¯¯¯→K¯¯¯¯¯ℓ+ℓ− decays at the double-differential level. Cancellations of soft and collinear divergences are demonstrated analytically using the phase space slicing method. Whereas soft divergences are found to cancel at the differential level, the cancellation of the hard-collinear logs ln mℓ require, besides photon-inclusiveness, a specific choice of kinematic variables. In particular, hard-collinear logs in the lepton-pair invariant mass distribution (q2), are sizeable and need to be treated with care when comparing with experiment. Virtual and real amplitudes are evaluated using an effective mesonic Lagrangian. Crucially, we show that going beyond this approximation does not introduce any further infrared sensitive terms. All analytic computations are performed for generic charges and are therefore adaptable to semileptonic decays such as B¯¯¯¯→Dℓν¯¯¯.