code: def chi_a(mass1, mass2, spin1z, spin2z):
      
      return (spin2z * mass2 - spin1z * mass1) / (mass2 + mass1)
docs: Returns the aligned mass-weighted spin difference from mass1, mass2,
      spin1z, and spin2z.