A descent variant minus the number of inversions.
This statistic is defined for general finite crystallographic root system $\Phi$ with Weyl group $W$ as follows: Let $2\rho = \sum_{\beta \in \Phi^+} \beta = \sum_{\alpha\in\Delta}b_\alpha \alpha$ be the sum of the positive roots expressed in the simple roots.
For $w \in W$ this statistic is then
$$\operatorname{stat}(w) = \sum_{\alpha\in\Delta\,:\,w(\alpha) \in \Phi^-}b_\alpha - \ell(w)\,,$$
where the sum ranges over all descents of $w$ and $\ell(w)$ is the Coxeter length.
It was shown in [1], that for irreducible groups, it holds that
$$\sum_{w\in W} q^{\operatorname{stat}(w)} = f\prod_{\alpha \in \Delta} \frac{1-q^{b_\alpha}}{1-q^{e_\alpha}}\,,$$
where $\{ e_\alpha \mid \alpha \in \Delta\}$ are the exponents of the group and $f$ is its index of connection, i.e., the index of the root lattice inside the weight lattice.
For a permutation $\sigma \in S_n$, this becomes
$$\operatorname{stat}(\sigma) = \sum_{i \in \operatorname{Des}(\sigma)}i\cdot(n-i) - \operatorname{inv}(\sigma)\,.$$

This map sends a permutation $\pi$ to $\pi^{-1} \star \pi$ where $\star$ denotes the Demazure product on permutations.
This map is a surjection onto the set of involutions, i.e., the set of permutations $\pi$ for which $\pi = \pi^{-1}$.

Return the permutation whose Lehmer code equals the major code of the preimage.
This map sends the major index to the number of inversions.