In Warren Siegel's Fields page 131 , He says that in order to get rid of the negative norm zero component of the vector wavefunction , we impose the constraint SabPa+ωPb=0 in analogy to the on-shell constraint P2+m2=0. I don't see how this can help us get rid of negative probability wavefunctions.