s2max = min(spin2_bounds[1], (mtot*chis - mass1*spin1_bounds[0])/mass2)