The Eulerian distribution on the involutions of the hyperoctahedral group is unimodal

The Eulerian distribution on the involutions of the symmetric group is unimodal, as shown by Guo and Zeng. In this paper we prove that the Eulerian distribution on the involutions of the hyperoctahedral group, when viewed as a colored permutation group, is unimodal in a similar way and we compute its generating function, using signed quasisymmetric functions.