  The theoretical description of observables at collider experiments relies on
factorization theorems separating perturbative dynamics from universal
non-perturbative matrix elements. Despite significant recent progress in
extending these factorization theorems to increasingly differential jet
substructure observables, the focus has been primarily on infrared safe
observables sensitive only to correlations in the energy of final state
hadrons. However, significant information about the dynamics of the underlying
collision is encoded in how energy is correlated between hadrons of different
quantum numbers. In this paper we extend the class of calculable correlations
by deriving factorization theorems for a broad class of correlations, $\langle
\mathcal{E}_{R_1}(n_1) \cdots \mathcal{E}_{R_k}(n_k) \rangle$, between the
energy flux carried by hadrons specified by quantum numbers, $R_1, \cdots,
R_k$. We show that these factorization theorems involve moments of a new class
of universal non-perturbative functions, the "joint track functions", which
extend the track function formalism to describe the fraction of energy carried
by hadrons of multiple quantum numbers arising from the fragmentation of quarks
or gluons. We study the general properties of these functions, and then apply
this to the specific case of joint track functions for positive and negative
electromagnetic charges. We extract these from parton shower Monte Carlo
programs and use them to calculate correlations in electromagnetically charged
energy flux. We additionally propose and study a C-odd
$\mathcal{E}_{\mathcal{Q}}$ detector, which results in a qualitatively distinct
scaling behavior compared to the standard energy correlators. Our formalism
significantly extends the class of observables that can be computed at hadron
colliders, with a wide range of applications from particle to nuclear physics.
