Improve Pauli algebra for matrices with Lorentz indices

This doesn't work
```
PauliTrace[SI[a, b]] // PauliSimplify
```
This sort of works
```
PauliTrace[SI[a, b]] // LorentzToCartesian // PauliSimplify
```