Include parsing of spin_1z and spin_2z from posterior samples files

Currently the source_classification_pefunction only allows for reading total spin and tilt. We need to make sure that if the samples provide spin_1z and spin_1z instead then the code is able to use those too instead of defaulting the spins to zero.