fn permutation_polarization<const N: usize>( pair: [Jet<N>; 2], pol: u8, helicity: bool, ) -> [Jet<N>; 2]