In this paper, we develop sample path large deviations for multivariate Hawkes processes with heavy-tailed mutual excitation rates. Our results address a broad class of rare events in Hawkes processes at the sample path level and, via the cluster representation of Hawkes processes and a recent result on the tail asymptotics of the cluster sizes, unravel the most likely configurations of (multiple) large clusters that could trigger the target events. Our proof hinges on establishing the asymptotic equivalence, in terms of M-convergence, between a suitably scaled multivariate Hawkes process and a coupled Lévy process with multivariate hidden regular variation. Hence, along the way, we derive a sample path large deviations principle for a class of Lévy processes with multivariate hidden regular variation, which not only plays an auxiliary role in our analysis but is also of independent interest.