An extension of MacMahon's Equidistribution Theorem to ordered set partitions

We prove a conjecture of Haglund which can be seen as an extension of the equidistribution of the inversion number and the major index over permutations to ordered set partitions. Haglund's conjecture implicitly defines two statistics on ordered set partitions and states that they are equidistributed. The implied inversion statistic is equivalent to a statistic on ordered set partitions studied by Steingr\'imsson, Ishikawa, Kasraoui, and Zeng, and is known to have a nice distribution in terms of $q$-Stirling numbers. The resulting major index exhibits a combinatorial relationship between $q$-Stirling numbers and the Euler-Mahonian distribution on the symmetric group, solving a problem posed by Steingr\'imsson.